为什么我的程序做不出文章中给出的非自治系统的图像？
C:\Users\lenovo\Desktop\model and parameters.png最近看了一篇文章，里面分析了一个模型，我编了一个程序想做一下文章中给出的图形，但是总是得不到正确的图像。模型和相应的参数以及图像我作为附件上传了。我写的程序如下。请大家帮我看看是怎么回事，谢谢了
rhd=@(t,x,a0,a1,b,alpher,beta,s,gama)); ...
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=;%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);
=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
看图像吧。谢谢各位了 clear;clc
clf
x0=;%initial value
t_final=200;
options = odeset('RelTol',1e-10,'AbsTol',1e-10);
=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);]; 虽能得出图形，跟你的还是不一样，你看看还有什么地方不对，再交流
页:
[1]