相图的函数分析与分型线.doc

第三讲 相图的函数分析与分型线 3.1 能量函数与保守系统 对于无阻尼的，质量是m，受力为的质点，其运动方程为 可以设，那么化成方程组 对于方程的两边进行计算 ， = 这里函数称为能量函数，包含两部分，一是动能，二是势能，那么方程的含义就是能量保持不变，那么我们称方程表示的系统为保守系统，是一个能量守恒的系统。 例子：双势井的Duffing方程 利用第二章判断平衡点（0,0）是鞍点，（-1,0）和（1,0）是中心。 能量函数 [x,y]=meshgrid(-2:0.1:2); z=0.5*y.^2-0.5*x.^2+0.25*x.^4; mesh(x,y,z) 等位线 比如对应着（-1,0）和（1,0） 是两个封闭的曲线，分布围绕（-1,0）和（1,0）的轨迹。 是一个封闭的曲线。 function out1 = duffexam(t,y) out1 = [y(2); y(1)-y(1).^3]; options = odeset(RelTol,1e-4,AbsTol,[1e-4 1e-4]); [t,y] = ode45(duffexam,[0 10],[2; 0],options); plot(y(:,1),y(:,2)) hold on [t,y] = ode45(duffexam,[0 10],[-0.5; 0],options);%多次取不同的初始向量 这