VandePol非线性方程(极限环)特性分析.docVIP

Van de Pol非线性方程（极限环）特性分析 【实验目的】 1、了解Van de Pol所代表的物理现象。 2、熟悉相平面法分析问题，并理解极限环。 【实验设备与软件】 MATLAB/Simulink软件 【实验原理】 Van de Vol方程及其物理现象 Van de Pol方程为 式中，阻尼有正负值（非线性阻尼）。能量保守的自由振动只是一种理想化模型，在实际系统中总会由于阻尼而振动衰减下来。强迫的周期力可以使振动不衰减而维持着，但这是靠外来的力量。自激振动第一次阐明了一个系统怎样依靠内部的常能量源而维持振动。 【实验内容、方法、过程与分析】 %求著名的Van Der Pol 方程x+( x^2 ?1)x + x = 0的数值解并绘制其时间响应 %曲线和状态轨迹图 %1.演化为状态方程 %令x1 = x,x2= x, 把x+( x^2 ?1)x + x = 0写成状态方程x1=(1-x2^2)*x1-x2,x2=x1 xdot=zeros(2,1);%使xdot 成为二元零向量采用列向量以便被matlab 其他指令调用 xdot(1)=(1-x(2)^2)*x(1)-x(2); xdot(2)=x(1); Command Window t0=0;tf=20;x0=[0,0.25];[t,x]=ode45(vdp,[t0,tf],x0);subplot(1,2,1),plot(t,x(:,1),:b,t,x(:,2),-r), legend(速度,位移),subplot(1,2,2),plot(x(:,1),x(:,2)) 奇点及极限环： 【实验结论与总结】