请问为什么我的答案是个数字零而不是一个向量?。

2 次查看(过去 30 天)
下面是我的代码。我想画个c和函数的解的二维图。但是函数的解为零。这是空集的意思吗?我可以证明函数有唯一解呀。
clear all;
close all;
c=0.001:0.01:1;
k=0.001:0.01:1;
q=c.*(1-k)./(c.^2.*(1-k)-k);
sigma_x=0.01:0.1:10;
sigma_delta=0.001:0.01:1;
alpha_z=sigma_x./(sigma_x+sigma_delta);
sigma_u=0.05:0.5:50;
sigma_a=0.01:0.1:10;
sigma_m=0.025:0.25:25;
alpha_m=sigma_u./(sigma_u+sigma_m);
alpha_a=sigma_u./(sigma_u+sigma_m+sigma_a);
syms b;
solve(b==(-q.*(1-c.*q).*b.^2.*alpha_z.*sigma_x./(alpha_a.*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)).*sigma_u./(c.^2.*((alpha_a-alpha_m).*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)+b.^2.*sigma_x+(q.*(1-c.*q).*b.^2.*alpha_z.*sigma_x./(alpha_a.*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)).^2.*(alpha_a.*sigma_u+q.*b.^2.*alpha_z.*sigma_x)))

采纳的回答

开元显示注单数据延迟风控审核失败不给出款怎么处理?【微81923899】
参考《 函数求解(超链接,自己点进去看) 》,写个for循环 然后每一步针对组数值进行求解
至于解方程,请去阅读论坛《 MATLAB最基础教程(超链接,自己点进去看) 》《 论坛常见问题归纳(超链接,自己点进去看) 》等帖子,或者入门书籍如《 邢树军译版的matlab教程(超链接,自己点进去看) 》,里面全有介绍和讲解
比如对第一个方程的数值求解,要写成
c=0.001;
k=0.001;
q=c.*(1-k)./(c.^2.*(1-k)-k);
sigma_x=0.01;
sigma_delta=0.001;
alpha_z=sigma_x./(sigma_x+sigma_delta);
sigma_u=0.05;
sigma_a=0.01;
sigma_m=0.025;
alpha_m=sigma_u./(sigma_u+sigma_m);
alpha_a=sigma_u./(sigma_u+sigma_m+sigma_a);
f = @( b ) -b + (-q.*(1-c.*q).*b.^2.*alpha_z.*sigma_x./(alpha_a.*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)).*sigma_u./(c.^2.*((alpha_a-alpha_m).*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)+b.^2.*sigma_x+(q.*(1-c.*q).*b.^2.*alpha_z.*sigma_x./(alpha_a.*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)).^2.*(alpha_a.*sigma_u+q.*b.^2.*alpha_z.*sigma_x));
X = fzero( f, 1 )
f( X )

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!