| 查看: 368 | 回复: 0 | ||
mayaya26木虫 (小有名气)
|
[求助]
混沌仿真
|
|
子程序 function dy=Bao2(t,y) dy=zeros(3,1); global a b c a=1.812;c=0.5; dy(1)=(b/(a*(1-c)))*(-exp(-y(2))+1-y(3)); dy(2)=b/(a*c)*y(3); dy(3)=-(a*c*(1-c)/b)*(y(1)+y(2))-y(3)/a; 主程序 clear format long global a b c b1 dt=0.001; tspan=(0:dt:50); Y0=[5 5 5]; for i=1:1:1001 b=0.015*(i-1); sol=ode45(@Bao2,[tspan(1) tspan(end)],Y0); ss=deval(sol,tspan)'; s=ss(20000:end,  ;aa=1; bb=0; cc=0; p=-1; func=((aa*s(:,1))+(bb*s(:,2))+(cc*s(:,3))+p)>0; func_n=(func(1:end-1)>0)&(xor(func(1:end-1),func(2:end))); m=find(func_n==1); t_intsec = -(aa*s(m,1)+bb*s(m,2)+cc*s(m,3)+p)./(aa*(s(m+1,1)-s(m,1))+bb*(s(m+1,2)-s(m,2))+cc*(s(m+1,3)-s(m,3))); z_map=s(m,3)+(s(m+1,3)-s(m,3)).*t_intsec; z1_map=log10(z_map); b1=log10(b*ones(size(s(m,3)))) H=plot(b1,z1_map); set(H,'linestyle','none','marker','.','markersize',1) hold on drawnow end 但是出现错误 b1 = Empty matrix: 0-by-1 为什么？是初始值设错了？还是...... 有谁帮忙，急急急，谢谢 |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有11人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
疑惑？
已经有5人回复