yxy 发表于 2006-10-7 08:27

关于差分方程的一个问题,求助

第一次编matlab程序,出现很多问题,还是让高手来看看吧!
matlab程序如下:

function x=yxy21(x0,n,a1,a2,b,c)
x=x0;x1=a1*b*c*x0;x2=a1*b*c*x1+a2*b*(1-a1)*b*c*x0;
for k=1:n
x(k+2)=a1*b*c*x(k+1)+a2*b*(1-a1)*b*c*x(k);
end
k=(0:20)';
y1=yxy21(100,20,0.5,0.25,0.18,10);y2=yxy21(100,20,0.5,0.25,0.19,10);
y3=yxy21(100,20,0.5,0.25,0.20,10);round();
plot(k,y1,k,y2,':',k,y3,'--');

不知为何,那个function部分总是出现各种错误,改了很久也不明所以。

jimin 发表于 2006-10-7 09:39

function y=yxy21(a,b)
x0=1;
y=[];
for k=1:10
y=;
x0=y(end);
end
>> y=yxy21(2,3)

y =

         5
          13
          29
          61
         125
         253
         509
      1021
      2045
      4093
我随便写了个函数,和你的差不多,你可以参考一下,然后写出你自己需要的

yxy 发表于 2006-10-7 10:09

谢谢了
页: [1]
查看完整版本: 关于差分方程的一个问题,求助