一个小程序

gina-yanke · 发表于 2005-11-4 10:59

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

a=1;b=3;c=1;d=5;s=2;r=0.001; z=[0];y=[0];x=[-1.6]; I2=[];I1=[]; t=1; while t<=3 I2(t)=4*0.08*(exp(-(t-1)/3)-exp(-(t-1)/2)); I1(t)=4*(-1)*(exp(-(t-1)/3)-exp(-(t-1)/2)); z(t+1)=r*(s*(s-x(1))-z(t)); y(t+1)=c-d*x(t)^2-y(t); x(t+1)=y(t)+a*x(t)^3-b*x(t)^2-z(t)+I1(t)+I2(t); t=t+1; end 运行之后 >> x x = 1.0e+009 * -0.0000 -0.0000 -0.0000 -8.7706 但是答案不应该是这么多。于是输入 >> x(2) ans = -11.7760 但是之前的x显示第二个成员的值为0啊。 这是为什么啊？？？？

linqus · 发表于 2005-11-4 11:25

你要使用format long格式， 俺的结果如下： To get started, select MATLAB Help or Demos from the Help menu.
>> a=1;b=3;c=1;d=5;s=2;r=0.001; z=[0];y=[0];x=[-1.6]; I2=[];I1=[]; t=1; while t<=3 I2(t)=4*0.08*(exp(-(t-1)/3)-exp(-(t-1)/2)); I1(t)=4*(-1)*(exp(-(t-1)/3)-exp(-(t-1)/2)); z(t+1)=r*(s*(s-x(1))-z(t)); y(t+1)=c-d*x(t)^2-y(t); x(t+1)=y(t)+a*x(t)^3-b*x(t)^2-z(t)+I1(t)+I2(t); t=t+1; end >> >> x x = 1.0e+009 * -0.00000000160000 -0.00000001177600 -0.00000206126163 -8.77063443897300 >> x(2) ans = -11.77600000000000 >> format short >> x(2) ans = -11.7760 >> x x = 1.0e+009 * -0.0000 -0.0000 -0.0000 -8.7706 >> 呵呵

gina-yanke · 发表于 2005-11-4 18:16

太感谢了

账号		自动登录	找回密码
密码			我要加入