请问这样的方程该如何写?在1stOpt也可以
for i=120:0.5:220yy(2*i-120*2+1)=-0.000639168636347769*exp(-0.5*((i-150.661324113566)/9.03565929752912)^2)+4.74344306737858E-5*exp(-0.5*((i-159.693770225366)/4.10900575550029)^2)+0.00560657896944599*exp(-0.5*((i-109.692079888009)/24.6544137370836)^2);
end
xx=1.44*yy-0.28
xxx=diff(xx)
142957.1Y''+506.35Y'+37.08Y=26.68xx+246.35xxx-430
求Y。 本帖最后由 scott198510 于 2011-7-19 19:02 编辑
1# myfisher
楼组问题有点诡异,求出的xx,xxx都是数组,后面等式右边是一个数组,不是一个值,不知道要什么结果,接前面部分作一下变换,循环时候角标只能是整数,最好不用i作为变量,与虚数i易混淆;
a=120:0.5:220;
aa=(2*a-120*2+1); %求出aa=1:201
for t=min(aa):max(aa);
x=(t-1+120*2)/2;
yy(x)=-0.000639168636347769*exp(-0.5*((x-150.661324113566)/9.03565929752912)^2)+4.74344306737858E-5*exp(-0.5*((x-159.693770225366)/4.10900575550029)^2)+0.00560657896944599*exp(-0.5*((x-109.692079888009)/24.6544137370836)^2);
end
xx=1.44*yy-0.28; %201个
xxx=diff(xx) ;% 200个
%xx与xxx维数不匹配,无法进行后续计算
页:
[1]