matlab画图程序..docxVIP

1.第一步，建立一个M文件，用来存贮函数，本例题以达芬方程（Duffing）为例，其中force为参数function df=dafen(t,x,flag,force)df=[x(2);force*cos(1.2*t)-x(1)^3+x(1)-0.3*x(2)];第二步，建立一个画图的M文件clearforce=0.222;options=odeset(RelTol,1e-7);%定义误差精度的，系统默认1e-3，如果改为1e-3，X将等于0?tt=2*pi/1.2 %定义步长的?[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],force);figureplot(x(2000:end,1),x(2000:end,2),-)%X=x-xxx %检验options的?%pojialaihold oni=2000:100:3000plot(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 onfor x=1:200for k=1:500 y=1-x*y*y/100; set(p,Xdata,x/100,Ydata,y); drawnowset(p,Xdata,x/100,Ydata,y); drawnowendend由于使用的是动画画图，所以只能使用屏幕截图保存，若此时点击窗口，已经画出的点会全部消失。当然，将命令改用打点画图的话就可以点击窗口保存。 说明： Xn+1＝1-k*Xn*Xn 图形表示的是：对于0到2之间的一个k值，任给X一个初值，经过上式迭代循环500次之后，得到的X的值与k的关系图。 可以看到： 1.最初k较小的时候(比如k＝0.1)，不论最初的X取什么，最后总是得到一个稳定值。 2.k增大到某个值时，不论最初X取什么，最后得到的是在两个X值之间跳跃的结果，即图像开始分裂了，有了二周期。 3.图像继续分裂，出现四周期和八周期，最后混沌。 4.图像中有上下两个分裂叉，其中上分裂叉有一条隐边界贯穿到下面。 5.更多性质......3. clearforce=0.05;options=odeset(RelTol,1e-7);%定义误差精度的，系统默认1e-3，如果改为1e-3，X将等于0 tt=2*pi/1.2 %定义步长的 [t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],force); figureplot(x(2000:end,1),x(2000:end,2),-)%X=x-xxx %检验options的 %pojialaihold oni=2000:100:3000plot(x(i,1),x(i,2),*)M-FILEfunction dx=duffing(t,x)mu=1.0;F=0.05;w=1;dx=[x(2);F*sin(w*t)+mu*x(1)-x(1).^3-0.12*((x(1)^2)-1)/((x(1)^2)+1)*x(2)+0.4*cos(t)]程序、[t,x]=ode45(@duffing,[0,2800],[0,1.5]);x1=x(:,1);x2=x(:,2);x1=mod(x1,2*pi);x1(x1pi)=x1(x1pi)-2*pi;plot(t(1:50:end),x1(1:50:end))%频闪采样图形figureh=plot(x1,x2)同宿轨ezplot(-sech(1.414*t),-sech(1.414*t)*tanh(1.414*t),[-4*pi,4*pi])Lyapunov指数图M文件：function dX = Rossler_ly(t,X)k=0.72;B=0.12;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(1.8*z)+0.12*cos(1.8*z);dz= z;% Rossler吸引子的Jacobi矩阵Jaco = [ 0, 1, 0; 1-3*x^2, -k*((x^2)-1)/((x^2)+1), (B/1.8)*cos(1.8*z)-(0.12/1.8)*sin(1.8*z); 0,