李亚普洛夫稳定性演示图

最新推荐文章于 2024-07-08 17:40:17 发布

automan_huyaoge

最新推荐文章于 2024-07-08 17:40:17 发布

阅读量5.1k

点赞数 1

分类专栏：控制科学与工程 Matlab

原文链接：https://zhuanlan.zhihu.com/p/58738073

版权

控制科学与工程同时被 2 个专栏收录

179 篇文章 19 订阅

订阅专栏

Matlab

2 篇文章 0 订阅

订阅专栏

该博客展示了如何使用MATLAB绘制不稳定系统的轨迹图。通过ode45函数求解常微分方程，并利用quiver函数添加箭头指示轨迹方向。博客中探讨了不同初始条件下的系统行为，特别是当初始值接近和远离平衡点时的动态特性。

摘要由CSDN通过智能技术生成

不稳定系统轨迹图

function d=dxdt(t,x)

d=[ x(2)+x(1)*(2-x(1)^2-x(2)^2); 
    -x(1)+x(2)*(2-x(1)^2-x(2)^2) ];

figure('color','w');
hold on 
for theta=[0:20]*pi/10
    x0=3*[cos(theta);sin(theta)];%定义初始值数组
    [t,x]=ode45(@dxdt,[0:0.1:8],x0);
    plot(x(:,1),x(:,2),'linewidth',0.5)
    quiver(x(:,1),x(:,2),gradient(x(:,1)),gradient(x(:,2)),'linewidth',3.0);%增加轨迹方向箭头
end
for theta=[0:2:20]*pi/10
    x0=1e-5*[cos(theta);sin(theta)]; [t,x]=ode45(@dxdt,[0:0.2:20],x0);
    plot(x(:,1),x(:,2),'linewidth',0.5)
    quiver(x(:,1),x(:,2),gradient(x(:,1)),gradient(x(:,2)),'linewidth',1.5)
    xlabel('x1','FontSize',18,'FontWeight','bold','Color','r');
    ylabel('x2','FontSize',18,'FontWeight','bold','Color','r')

end

参考：

https://zhuanlan.zhihu.com/p/58738073