matlab绘制hpof振荡器的极限环和输出结果

虚心学习四足步态中.....................

先贴个公式

其实简单说就是，用matlab解微分方程组

函数里的代码

function FF = f( t,x )%UNTITLED 此处显示有关此函数的摘要
%   此处显示详细说明
a=100; 
B=0.75;
u=1;
W=2*pi;

r1=x(1)^2+x(2)^2;

FF=[a*(u-r1)*x(1)-W*x(2);
    a*(u-r1)*x(2)+W*x(1);
     ];

end

脚本里的代码

[t,y] = ode45(@hopf,[0 5],[0;0.1]);

subplot(2,1,1)
plot(y(:,1),y(:,2))


subplot(2,1,2)
plot(t,y(:,1),t,y(:,2))

结果如图

后续慢慢更新