利用MATLAB解特征方程,并画出特征根的分布,便于分析系统的稳定性

  本文主要介绍如何利用MATLAB解特征方程,并将特征根的分布画在坐标轴上,便于分析系统的稳定性

   我们知道,一旦求出系统的闭环特征根就很容易判定系统的稳定性,但是对于高阶系统,闭环特征根求起来是很困难的,但是利用MATLAB这一强大的数学工具,求起来就很方便,也很简单了,比如求解以下特征方程的特征根:
s 6 + 3 s 5 + 16 s 4 + 2 s 3 − 4 s 2 − 25 s − 60 = 0 s^6+3s^5+16s^4+2s^3-4s^2-25s-60= 0 s6+3s5+16s4+2s34s225s60=0
    1、把特征方程中的系数按照阶次从高到低放到一个向量里面去,若有缺项,该项系数应该补零,我们以上面的式子为例:
p=[ 1 3 16 2  -4  -25  -60 ];
    2、利用roots函数求解,roots函数可用于求解多项式的根,在完成上一步的情况下,只需要调用roots(p)就可以求出多项式的根。
r=roots(p)
    结果如下:
r =

  -1.4204 + 3.6565i
  -1.4204 - 3.6565i
   1.4475 + 0.0000i
  -0.1563 + 1.4343i
  -0.1563 - 1.4343i
  -1.2941 + 0.0000i
    3、为了更加直观的观察根的分布我们把它画在坐标轴中,调用real()函数可以获取复数的实部,调用imag()函数可以获取复数的虚部,再利用plot()以点的形式画出来就可以了,用×来标注
h=plot(real(r),imag(r),'rX');
    图像如下:

在这里插入图片描述

    4、为了更好地观察,我们对其进行放大,再加上网格线,用set()函数设置线宽linewidth为5,设置标记符号(此处为×)的大小markersize为15,利用set(gca)可以调整图形的刻度大小,语句grid on可以在画图的时候添加网格线,此部分代码如下:
set(h,'linewidth',5,'markersize',15);
set(gca,'fontsize',15);
grid on 
    5、如果需要的话,可以画出实轴和虚轴,代码如下(范围自主调整):
plot([0 0],[-5 5])
plot([-2 2],[0 0])
    6、完整的代码和图像如下:
p=[ 1 3 16 2  -4  -25  -60 ];
r=roots(p);
h=plot(real(r),imag(r),'rX');
set(h,'linewidth',5,'markersize',15);
set(gca,'fontsize',15);
hold on
plot([0 0],[-5 5])
plot([-2 2],[0 0])
grid on 

在这里插入图片描述
       (点击图上的点可以查看坐标值,对应特征根的实部和虚部,如下)
在这里插入图片描述

           很容易看出,系统是不稳定的
  • 29
    点赞
  • 146
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 3
    评论
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

慕羽★

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值