《绘制Duffing振子的分叉图的程序》由会员分享，可在线阅读，更多相关《绘制Duffing振子的分叉图的程序（7页珍藏版）》请在金锄头文库上搜索。
1、这些程序思想有些也许不对的,有问题，自己改善,我不再负责对这些程序解释。由于我都不懂得道理在哪里。但是盼望您能在程序的提示下，进一步的做改善或者改正，以期获得更为精确的成果。别照搬和迷恋别人的程序！% %绘制ufig振子的庞加莱截面图的程序% % buchang:已知鼓励下步长数值的大小,%tend程序仿真达到1个鼓励周期的总时间,% clear；l% gobamckk3 F0 m% m=1；c0.1；k1=0;k31；mga=1;=12% x=3；% tstart=0;Tbusu=60；bucang=（2*pi/meg)/hu;tnd=(2*pi/mega)*150；% sn=tstr：ca
2、ng:tnd;%t,y=de45(dain3,tpan,x)；cont=nd（t(*piomea*40）); % 去掉前40个周期的鼓励时间以消除瞬态响应的影响% y(count,：)；% TDta=Y(:Tbushu,)-(Tbus+1):Tbushu2,1);% maxvalu,nies=ax(ab(TData)）%oitnum=rond(（tnd-2*pi/omga*40)/buchang/Tbuu)-1；%s=zes（oinnuer,1)； velzrs(poinnumer,1);%fori=1:oine% i（i,1)=Y(Tbush*(-1)+indic,）;% velo（i,1）
3、Y(Tuhu*(i-1)indics,);% end% figur,lt（i,elo，b.,mksz,5);% % 绘制Dufng振子的分叉图的程序% lear;cl% lbl m c k1 3F0 omega；%=1；k1=;k3=1;oa1；F0=2; rane.01：0.01:1; YY=；k=0；% for c=range% kk；% 0=3,4;% tpan0:0.01:2;% t,de45（dafin3,tsan，y0)；% ount=find(100)；% Y=(ut,:);% % 画x的分岔图。% =;% =eng(Y(：,)）; for 2:n-1% f Y(i1，1）+es
4、（i，）& Y（i,1)Y(i1，1)eps %简朴的取出局部最大值。% YY（k,j)=Y(i,)； %使最大值计数个数自动增长% =j1;% end% end% if j1% lot（c,Y（,1:j-1)，b,markersize,）； end% od;% index(k)=-1;% d% xabl(c)；% ylae（ mx)；% ite（afinifuratn iram)；% % % 绘制分岔图的程序% cr,clc% obal c k kFomega% =1;=.；k1=0;k3=1;omega=1；0=12;canu=0.0:0.01:1;% for =1:100 =ccashu
5、(k)% x0=3;4;% spn=0:0.01*pi：500;% t,=ode5(dafin3,a,x0);% disros(50，1);% lo=eos(5，1); fr i=:5% ds(i,1）（00*(i+20),);% vel(i,)=y（0*(+),2);% d% smatrix(k,:)dis;% nd% fure,pt（cans，Datri，b.,arkesize,3）;% %线性参数k1的变化产生的分岔图 a；clc% lobal c k k3 F0 omega m1；c=0.1;k3=;oeg=1;F=12; kcanshu001:0.01:;% =1:200 k=can
6、shu（)% x03;4；% san=0:0.01*pi:500;% t，y=d45(daf,tpn,x0); dis=zro（50,1);% veloos（0,1）;% for i=1：0% s(,1）=y(10*(i+20),1); velo(,)y(100*(i+20),2); e Dimatri（，:）dis;% end plot(kcnsh，isaix,b.，masi,5);%ie(参数变化下的分岔图)% xlabel（线性刚度参数1的变化)% ylabl(X值)非线性参数k3的变化产生的分岔图% car；l% loalm c k k F0 oega% m;c=0.;k0;oega1
7、;F0=2;% kcanshu=0.01:0.:2；% f k=1:2% k3=kcshu(k)% x=3；4；% tspn0:01*2pi:500;% t,=o5(di3,tspan，x）;% dis=zeros（50,);% velo=zeros(50,1）;% for=1:50% dis(i,1)=y(0(i+),1)；% vo（i,）=(00（i+20),2）;% end Dsmatri(k,：)=ds； nd plo(kcash,Dismatri,.，marerize,5);% titl(参数变化下的分岔图)% xlabel(非线性参数k3的变化)% ylabel(X值) %鼓励参数
8、0变化产生的分岔图% cler；clc% glal c1 k3 0 oegam1;c=.;k1=0;k31;mega=1;% F0canhu0.：.：0;% fo k=1:20% F0cnh(k）% x0=3;% tsan=0:0.012*pi：500;% t,=od45(dafin3,tsan,0）;% dis=zeros(50,1);% velo=zeos（5，1)； fori=：50% s（i,1)y(100（i+2),1);% velo(i,）（1*(+0),)；% end% Disatrix(k,:)di;% end% plot（F0cansu,Diarix,b.，markrsize
9、，5);% ite(参数变化下的分岔图)% lbe(鼓励参数F0的变化)%ylael(X值)% % %鼓励频率oega变化产生的分岔图% cear;lc glal m k1 3 F0 omea% m=1;=;k1=0；k=1;F0=2； mgacnshu=0.:0.1:0; ork=1：100% omeg=omacnshu(k) x03;4;% tpan=：.1*2*i/omega:50;% t,y=od5(dfi3,tspan，x）;% is=zer(50,)；% ve=eos（50，1);% o =1:50% dis（,1)=(nd(10*omea*(i+20),1)； velo(i，1)
10、(ound（00*mea*(i2)）,）；% end% Dsmatrix（k，:)=dis;% end% (omeacshu，Dimtrix,.,markersz，)； titl(参数变化下的分岔图)%xlbel(鼓励频率oega的变化）% labl(X值)cer；cl% loal c k1 kF0omega% n=3,rh_ext_fcn=dfin_ext2,fcnintegaore45，tstat=0,stpt=0.5,tend200,% ystar= 4 ,iutp=1，% m1;c=.1;10；0=2;k3=1；%omeaans=01:01:10; fok=1:10% oma =omgacansu(,k）,punovzhisu(,：)lapuovfu(n，rh_extn,f_ntgtor,tstr，tep，ted，ystar,ioutp)%en% fgure,po（omgcshu,lyaunozishu）,%tte（Lpunov动力学指数);%xabel(鼓励频率mega变化）; yabel(Lpuno 指数）;% 绘制分岔图的程序% lar;cl% globa m c k1 F0oega%m=1;c=0.1;k=0；k3=1;omeg1；F0=2; ccnshu=0.:.01:;% or :100 c
