| 查看: 409 | 回复: 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 为什么？是初始值设错了？还是...... 有谁帮忙，急急急，谢谢 |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
filecode=后面第一个是大写字母
已经有4人回复
-
时间戳又变了8-15
已经有20人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有58人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
