matlab函数状态空间系统ss、能控性矩阵ctrb、矩阵的秩rank、能控标准型canon、零极点配置place、系统极点pole等函数（线性定常系统）-Toy模板网

这篇具有很好参考价值的文章主要介绍了matlab函数状态空间系统ss、能控性矩阵ctrb、矩阵的秩rank、能控标准型canon、零极点配置place、系统极点pole等函数（线性定常系统）。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

matlab函数能控性矩阵ctrb、能控标准型canon、零极点配置place

第一章，线性定常系统

ss

如果已知线性定常系统的ABCD四个矩阵，可以得到状态空间系统

其他更具体的用法请直接看帮助文档。

用法：ss(A,B,C,D)

假如

matlab函数状态空间系统ss、能控性矩阵ctrb、矩阵的秩rank、能控标准型canon、零极点配置place、系统极点pole等函数（线性定常系统）,matlab,matlab,控制,c,代码,状态空间

可以输入

A = [-1.5,-2;1,0];
B = [0.5;0];
C = [0,1];
D = 0;
sys = ss(A,B,C,D)

最后得到

ctrb和rank

判断系统是否能控，可以用能控性矩阵是否奇异进行判断。ctrb函数用来生成能控性矩阵，rank用来判断矩阵的秩

对于线性定常系统
$\dot{x}=Ax+Bu$

如果能控性矩阵
$CO=[B\ \ AB\ \ A^2B\ \ \cdots\ \ A^{n-1}B]$
的行秩=n，则意味着该系统完全能控（为什么是行秩呢，因为对于MIMO系统能控性矩阵不一定是方阵）

在matlab中，可以用ctrb(A,B)直接得到能控性矩阵

比如，我事先输入矩阵A和B

再输入函数

ctrb(A,B)

就会得到能控性矩阵

也可以直接用上一个函数ss生成的状态空间系统（和上个例子所用数据不同）

A = [-1.5,-2;1,0];
B = [0.5;0];
C = [0,1];
D = 0;
sys = ss(A,B,C,D)
Co=ctrb(sys)

最后用rank(Co)计算一下能控性矩阵的秩。

比如对于这个能控性矩阵，其行秩为1，显然系统不能控

canon

要会这个函数，首先去看一下本篇第一个函数ss

状态空间方程化为高阶微分方程的实现方法是非唯一的，我们就提出了标准型便于我们研究和交流。

标准型有很多种。matlab中的标准型是第一能控标准型

用法canon(sys,‘companion’)

比如对于这个系统
matlab函数状态空间系统ss、能控性矩阵ctrb、矩阵的秩rank、能控标准型canon、零极点配置place、系统极点pole等函数（线性定常系统）,matlab,matlab,控制,c,代码,状态空间

把ABCD全都输进去，然后

sys=ss(A,B,C,D)

得到
matlab函数状态空间系统ss、能控性矩阵ctrb、矩阵的秩rank、能控标准型canon、零极点配置place、系统极点pole等函数（线性定常系统）,matlab,matlab,控制,c,代码,状态空间

再输入canon

canon(sys,'companion')

place或acker

设计全状态反馈控制律中，我们把配置闭环系统极点的过程称之为极点配置。可以用place函数（acker函数可以实现）

先在matlab中输入矩阵A和B
再输入想要的极点 $P=[第一个\ \ 第二个\ \ 第三个\dots]$
最后

acker[A,B,P] //或者用place[A,B,P]也可以，计算精度不同

输出的结果就是反馈矩阵K或说R

例子：

A = [-1,-2;1,0];
B = [2;0];
p = [-1,-2];
K = place(A,B,p);

得到

pole

计算系统极点。也要看本章第一个函数ss

A = [-1,-2;1,0];
B = [2;0];
C = [0,1];
D = 0;
sys = ss(A,B,C,D);
pole(sys)文章来源地址https://www.toymoban.com/news/detail-687767.html

到了这里，关于matlab函数状态空间系统ss、能控性矩阵ctrb、矩阵的秩rank、能控标准型canon、零极点配置place、系统极点pole等函数（线性定常系统）的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！