matlab进行计算是总是出现：在赋值 A(I)=B 中，B 和 I 中的元素数目必须相同

Question

clf
clear all
clc
syms x y k m z w
y=solve('20*((40/(8.8*y))^(-0.25))*(10*x+12.41*(y-2))+340*x^2*(1+0.2*(y/x))-0.3*8.8*(1.5*x^2-12.25)*y-26950=0','y');
k=[8.25:0.05:8.9];
m=[];
z=[];
for i=1:14
w=subs(y,k(i));
m(i)=w(abs(imag(w))<=eps);
z(i)=8.8*(m(i)-2)*1.5*(k(i))^2+400*k(i)+2*pi*1.975*(m(i)-2)*40+5*1.5*(k(i))^2*40*1.1333*(1+0.2*m(i)/k(i));%这是z关于x，y的值
end
x=k;
y=m;
subplot(1,2,1);
plot(x,y,'+-')
title('20*((40/(8.8*y))^(-0.25))*(10*x+12.41*(y-2))+340*x^2*(1+0.2*(y/x))-0.3*8.8*(1.5*x^2-12.25)*y-26950=0')
xlabel('x')
ylabel('y')
grid
subplot(1,2,2);
plot(x,z,'*-')
title('8.8*(y-2)*1.5*x^2+400*x+2*pi*1.975*(y-2)*40+5*1.5*x^2*40*1.1333*(1+0.2*y/x)')
xlabel('x')
ylabel('z')
grid
这是运行的程序 还请指教

Answer 1

循环语句中的

z(i)是z的第i个元素，matlab默认z(i)为一个数字，长度为1。后面赋值是一个行向量，得到的结果长度和x、y相同（长度201）
这造成等号左右的元素数目不同。
考虑z是一个矩阵，矩阵第一行存储i=1的结果。
你把赋值语句改成
for i=1:2
z(i,:)=c.*(d-x.^2................................
end
即可

Answer 2

你就不能把错误的那一行放上来么。。。MATLAB应该提示是哪一句出了问题吧追问

追答

tt = abs(imag(w))<=eps;
disp(tt);
m(i)=w(abs(imag(w))<=eps);

加两句，看看你的tt是多少
还有之前感觉你没有必要把m，w声明为syms类型变量额

追问

还是出现同样的报错。。。程序是别人帮忙写的，之前算过一次，只是x的取值和函数关系式有一点点差别 但是之前的那个就可以出结果，这个就一直报错

追答

当然报错了，但是我想知道有没有输出tt的数值啊。。。

追问

没有

追答

是这样的，你的m(i)显然只能是一个数，但是我不知道你w(abs(imag(w))<=eps);的输出是什么，如果不只是一个数的话，当然左边就不能被赋值了。

tt = (abs(imag(w))<=eps)
disp(tt);
m(i)=w(abs(imag(w))<=eps);

你看这个有输出么？
我需要知道你的右边到底输出了什么

追问

不知你是否方便留下QQ

追答

1365405275
晚上7点之后在，可以探讨一下

本回答被网友采纳