matlab中for循环为什么只运行了一次

Question

a=0.89;b=0.6;alpha=4.1*pi/180;m=0.4;n=2.05;
h0=0.01:0.01:1.2;k=size(h0);
for i=1:k
h=h0(i);
if (2*b-m*tan(alpha))<h & h<2*b
x=0:0.00001:((2*b-h)/tan(alpha)+n);
y=h-(x-(2*b-h)./tan(alpha)).*tan(alpha);
s=(a/b).*(pi*b^2.*(1-1/pi.*acos((y-b)./b))+(y-b).*(b^2-(y-b).^2).^(1/2));
f=pi*a*b*(m-(2*b-h)/tan(alpha))+trapz(x,s);
elseif n*tan(alpha)<h & h<(2*b-m*tan(alpha))
x=0:0.00001:(m+n);
y=h-(x-m).*tan(alpha); s=(a/b).*(pi*b^2.*(1-1/pi.*acos((y-b)./b))+(y-b).*(b^2-(y-b).^2).^(1/2));
f=trapz(x,s);
else
x=0:0.00001:(m+h/tan(alpha));
y=h-(x-m).*tan(alpha);
s=(a/b).*(pi*b^2.*(1-1/pi.*acos((y-b)./b))+(y-b).*(b^2-(y-b).^2).^(1/2));
f=trapz(x,s);
end
v(i)=f;
end
v

v =
0.0035

感觉应该显出来的v是个1*120的向量呀，为啥结果就显示第一次循环的值

Answer 1

把第二行的k=size(h0);换成k=length(h0);
因为size返回数组[1,120]，在后面的for里变成了for i=1:1.

追问

v =
Columns 1 through 4
0.0035 + 0.0000i 0.0063 + 0.0000i 0.0100 + 0.0000i 0.0148 + 0.000i
为什么运行的结果会是虚数啊，我单独运行当h0.001时结果显示为0.0035却不是虚数，理论上结果应该是实属的，这是为啥，求大神讲解

追答

虚部为零，你就用real去了吧，应该是计算误差的原因。

Answer 2

size()返回的是两个值，[行数，列数]=size(a)，所以你代码里k其实等于行数1。把size改成length就可以。

追问

v =
Columns 1 through 4
0.0035 + 0.0000i 0.0063 + 0.0000i 0.0100 + 0.0000i 0.0148 + 0.000i
为什么运行的结果会是虚数啊，我单独运行当h0.001时结果显示为0.0035却不是虚数，理论上结果应该是实属的，这是为啥，求大神讲解