| 查看: 379 | 回复: 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 为什么？是初始值设错了？还是...... 有谁帮忙，急急急，谢谢 |
回复此楼» 猜你喜欢
求调剂，一志愿材料科学与工程985，365分，
已经有5人回复
-
271求调剂
已经有10人回复
-
0854调剂
已经有7人回复
-
272分材料子求调剂
已经有44人回复
-
085410 273求调剂
已经有5人回复
-
材料考研调剂
已经有25人回复
-
材料工程281还有调剂机会吗
已经有11人回复
-
303求调剂
已经有12人回复
-
290调剂生物0860
已经有14人回复
-
农学0904 312求调剂
已经有3人回复
