马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本人新手需要仿真混沌蔡氏电路 编写m文件
function dx=chuas(t,x,beta,rho,sigma,a,c)
if x(1)>1
g=-c*x(1)-a+c;
else if absx(1)<=1
g=-a*x(1);
else
g=-c*x(1)+a-c;
end
end
dx=[beta*(x(2)-x(1)+g);x(1)-x(2)+x(3);-rho*x(2)-sigma*x(3)];
然后
>> t_final=100;x0=[0;0;0];b1=10.0000;r1=15.0000;s1=0.0385;a1=-1.2700;c1=-0.6800;
>> [t,x]=ode45(chuas,[0,t_final],x0,[],b1,r1,s1,a1,c1);plot(t,x),figure;plot(x(:,1),x(:,2));axis([-4 4 -1 1]);
为什么没有仿真结果哪?
真诚恳求大虾帮忙指导。
谢谢啦 |