有20个v,那应该有20个结果吧?
而且你的那个函数的定积分是没有解析解的.所以只能求数值解.
建议你利用trapz函数,具体的看help.
其实也可以用quad但是它需要函数的定义,所以有点麻烦,我不推荐.
下面有代码.
>> x=0:0.01:20;v=1:20;
>> for i=1:20
f(:,i)=((1+16*(x-1).^2)./(v(i)^2-40.*x.*(x-2))).^(1/2);
y(i)=trapz(x,f(:,i));
end
>> y
有20个v,那应该有20个结果吧?
而且你的那个函数的定积分是没有解析解的.所以只能求数值解.
建议你利用trapz函数,具体的看help.
其实也可以用quad但是它需要函数的定义,所以有点麻烦,我不推荐.
下面有代码.
>> x=0:0.01:20;v=1:20;
>> for i=1:20
f(:,i)=((1+16*(x-1).^2)./(v(i)^2-40.*x.*(x-2))).^(1/2);
y(i)=trapz(x,f(:,i));
end
>> y