What's wrong with my code?

1 次查看（过去 30 天）

显示更早的评论

Ani Asoyan 2020-6-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/543569-what-s-wrong-with-my-code

评论： Star Strider 2020-6-6

采纳的回答： Star Strider

在 MATLAB Online 中打开

here's my code..

syms a
b=2; e=0.9; q=0.4; g=0.6; k=0; d=0.5;
a=1:1:5;
N(a) = (q*b*e/a-d*e)/((q^2)*(e^2)/a+2*e*g)
pi(a)=(b-q*N*e)/a
m=N(a)
n=pi(a)
plot (a,N(a),'b')
hold on 
plot (a, pi(a),'red')

When I run it, it gives me these errors.

Error using  / 
Matrix dimensions must agree.
Error in sev (line 31)
N(a) = (q*b*e/a-d*e)/((q^2)*(e^2)/a+2*e*g)

When I use "./", it plots only pi, but N is on x axis fixed at 0, what can I do to make it work?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Star Strider 2020-6-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/543569-what-s-wrong-with-my-code#answer_447374

在 MATLAB Online 中打开

Try this:

syms a
b=2; e=0.9; q=0.4; g=0.6; k=0; d=0.5;
% a=1:1:5;
N(a) = (q*b*e/a-d*e)/((q^2)*(e^2)/a+2*e*g)
fpi(a)=(b-q*N*e)/a
m=N(a)
n=fpi(a)
figure
fplot (N(a),[1 5],'b')
hold on 
fplot (fpi(a), [1 5],'red')
hold off