| 查看: 377 | 回复: 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 为什么？是初始值设错了？还是...... 有谁帮忙，急急急，谢谢 |
回复此楼» 猜你喜欢
材料学硕297已过四六级求调剂推荐
已经有8人回复
-
本人考085602 化学工程 专硕
已经有21人回复
-
265求调剂
已经有12人回复
-
求调剂
已经有5人回复
-
一志愿中南化学337求调剂
已经有6人回复
-
350求调剂
已经有3人回复
-
295复试调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有22人回复
-
287求调剂
已经有7人回复
-
329求调剂
已经有7人回复
