为什么我的程序做不出文章中给出的非自治系统的图像？

蓝云天

C:\Users\lenovo\Desktop\model and parameters.png
最近看了一篇文章，里面分析了一个模型，我编了一个程序想做一下文章中给出的图形，但是总是得不到正确的图像。模型和相应的参数以及图像我作为附件上传了。我写的程序如下。请大家帮我看看是怎么回事，谢谢了

[email=rhd=@(t,x,a0,a1,b,alpher,beta,s,gama)[a0*(1+a1*sin(2*pi*t))*x(4)-(b+s*x(4))*x(1)-beta*x(1)*x(2]rhd=@(t,x,a0,a1,b,alpher,beta,s,gama)[a0*(1+a1*sin(2*pi*t))*x(4)-(b+s*x(4))*x(1)-beta*x(1)*x(2[/email]); ...  
    beta*x(1)*x(2)-(alpher+b+gama+s*x(4))*x(2); ...
    gama*x(2)-(b+s*x(4))*x(3); ((a0*(1+a1*sin(2*pi*t))-b)-s*x(4))*x(4)-alpher*x(2)];
t_final=200;  
x0=[0.6;0.1;0.1;0.8];  %initial value
% parameters
a0=0.02 ;
a1= 0.5;
b= 0.01;
alpher=0.475 ;
beta= 0.936;
s= 0.0001;
gama= 0.025;

options = odeset('RelTol',1e-10,'AbsTol',1e-10);
[t x]=ode45(rhd,0:2*pi/(2*pi)/100:t_final,x0,options,a0,a1,b,alpher,beta,s,gama);
figure
plot(t,x(:,2))

蓝云天

为什么图像显示不出来呢？麻烦大家挪步到
http://www.ilovematlab.cn/viewthread.php?tid=156740&page=1&extra=page%3D1
看图像吧。谢谢各位了

WF1987

clear;clc
clf
x0=[0.6;0.1;0.1;0.8];  %initial value
t_final=200;
options = odeset('RelTol',1e-10,'AbsTol',1e-10);
[t,x]=ode45(@rhd,0:0.01:t_final,x0,options);  
figure
plot(t,x(:,2))
%……………………………………………………………………………………
function dx=rhd(t,x)
a0=0.02 ;
a1=0.5;
b=0.01;
alpher=0.475 ;
beta=0.936;
s=0.0001;
gama=0.025;
pi=3.14;
dx=[a0*(1+a1*sin(2*pi*t))*x(4)-(b-s*x(4))*x(1)+beta*x(1)*x(2);
    beta*x(1)*x(2)-(alpher+b+gama+s*x(4))*x(2);
    gama*x(2)-(b+s*x(4))*x(3);
   ((a0*(1+a1*sin(2*pi*t))-b)-s*x(4))*x(4)-alpher*x(2);];

WF1987

虽能得出图形，跟你的还是不一样，你看看还有什么地方不对，再交流