| 查看: 248 | 回复: 0 | |||
[交流]
龙格库塔 四元求解 求大神
|
|
设置 x y z 的范围 设与不设 范围不受控制 这个模拟的是混沌电路 图形 完全不太像 function huatu x0=0.05*randn(4,1); [x0]=a(x0); [xx]=ch(x0,1,20000); plot3(xx(1,  ,xx(2,,xx(3,,'k');%axis([-4,3,-2.5,2,-4,1.5]); xlabel('x'); ylabel('y'); %zlabel('z'); %wlabei('w'); function[x]=map(xx,t) if abs(xx(4))>=1 hx=-40; else hx=-0.8; end A=[hx 4 0 0 ;-1 0 1 0;0 -1 0.65 0;1 0 0 0]; x=A*xx; function[x0]=a(x) h=0.002; a=h/2; aa=h/6; x=[0.198,0.000012,0.018,0.0001]'; for k=1:100 k0=map(x,k); x1=x+k0*a; k1=map(x1,k); x1=x+k1*a; k2=map(x1,k); x1=x+k2*h; k3=map(x1,k); x=x+aa*(k0+2*(k1+k2)+k3); end x0=x; function[xx]=ch(x,t,symbol) h=0.002; a=h/2; aa=h/6; xx=[]; for j=1:symbol k0=map(x,t); x1=x+k0*a; k1=map(x1,t); x1=x+k1*a; k2=map(x1,t); x1=x+k2*h; k3=map(x1,t); x=x+aa*(k0+2*(k1+k2)+k3); xx=[xx x]; end  O5O5CODIQGZ(IF[V1M1IKQJ.png |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有12人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有19人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有24人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有11人回复
20