matlab画图程序.docx

,建立一个 M文件,用来存贮函数,本例题以达芬方程( Duffing )为例,其中 force 为参数 function df=dafen(t,x,flag,force) df=[x(2);force*cos(*t)-x(1)^3+x(1)-*x(2)]; 第二步,建立一个画图的 M文件 clear force=; options=odeset('RelTol',1e-7);% 定义误差精度的,系统默认 1e-3 ,如果改为 1e-3 , X将等于 0 tt=2*pi/ %定义步长的[ t,x]=ode45(***@dafen,[0:tt/100:80*tt],[0,0],options,[],forc e); figure plot(x(2000:end,1),x(2000:end,2),'-') %X=x-xxx %检验 options 的%pojialai hold on i=2000:100:3000 plot(x(i,1),x(i,2),'*') 2. global x; global y; global k; y=1; x=1; p=plot(x,y,'.','EraseMode','none','MarkerSize',3); axis([0 2 -2 2]) hold on for x=1:200 for k=1:500 y=1-x*y*y/100; set(p,'Xdata',x/100,'Ydata',y); drawnow set(p,'Xdata',x/100,'Ydata',y); drawnow end end 由于使用的是动画画图,所以只能使用屏幕截图保存,若此时点击窗口,已经画出的点会全部消失。当然,将命令改用打点画图的话就可以点击窗口保存。说明: Xn+1 = 1-k*Xn*Xn 图形表示的是:对于 0到2之间的一个 k值,任给 X一个初值,经过上式迭代循环 500 次之后,得到的 X的值与 k的关系图。可以看到: k较小的时候(比如 k= ) ,不论最初的 X取什么,最后总是得到一个稳定值。 增大到某个值时,不论最初 X 取什么,最后得到的是在两个 X 值之间跳跃的结果,即图像开始分裂了,有了二周期。 ,出现四周期和八周期,最后混沌。 ,其中上分裂叉有一条隐边界贯穿到下面。 ...... 3. clear force=; options=odeset('RelTol',1e-7); %定义误差精度的,系统默认 1e-3 ,如果改为 1e-3 ,X 将等于 0 tt=2*pi/ %定义步长的[t,x]=ode45(***@dafen,[0:tt/100:80*tt],[0,0],options,[],force); figure plot(x(2000:end,1),x(2000:end,2),'-') %X=x-xxx %检验 options 的%pojialai hold on i=2000:100:3000 plot(x(i,1),x(i,2),'*') M-FILE function dx=duffing(t,x) mu=; F=; w=1; dx=[x(2);F*sin(w*t)+mu*x(1)-x(1).^3-*((x(1)^2)-1)/((x(1)^2)+1)*x(2)+*cos(t)] 程序、[t,x]=ode45(***@duffing,[0,2800],[0,]); x1=x(:,1); x2=x(:,2); x1=mod(x1,2*pi); x1(x1>pi)=x1(x1>pi)-2*pi; plot(t(1:50:end),x1(1:50:end))% 频闪采样图形 figure h=plot(x1,x2) 同宿轨 ezplot('-sech(*t)','-sech(*t)*tanh(*t)',[-4*pi,4*pi]) Lyapunov 指数图 M文件: function dX = Rossler_ly(t,X) k=; B=; x=X(1); y=X(2); z=X(3); %Y的三个列向量为相互正交的单位向量 Y= [X(4), X(7), X(10); X(5), X(8), X(11); X(6), X(9), X(12)]; %输出向量的初始化,必不可少 dX = zeros(12,1); % Rossler 吸引子 dx=y; dy=-k*y*((x^2)-1)/((x^2)+1)+x-x^3+B*sin(*z)+*cos(*z); dz= z;%