ctrb

能控性矩阵

折叠所有页面

语法

有限公司= ctrb (A, B) 有限公司= ctrb(系统)

有限公司= ctrb (A, B)返回可控性矩阵:

$C o ＝［ \begin{matrix} B & 一个 B & {一个}^{2} B & ．．． & {一个}^{n - 1} B \end{matrix} ］$

在哪里一个是一个n——- - - - - -n矩阵,B是一个n——- - - - - -米矩阵,有限公司有n行和纳米列。

有限公司= ctrb(系统)计算状态空间LTI对象的能控矩阵sys．该语法等价于:

有限公司= ctrb (sys.A sys.B);

系统是可控的，如果有限公司已满秩n．

打开生活的脚本

定义一个和B矩阵。

A = [1 1;4 2];B = [1 -1;1 1];

计算可控性矩阵。

有限公司= ctrb (A, B);

确定不可控状态的数量。

unco = length(A) - rank(Co)

不可思议的= 1

不可控状态表明有限公司没有完全排名2。因此，这个系统是不可控制的。

可控矩阵的秩估计是病态的;也就是说，它对舍入误差和数据中的误差非常敏感。从这个简单的例子可以看出这一点。

$一个＝ \begin{matrix} ［ \begin{matrix} 1 & δ \\ 0 & 1 \end{matrix} ］， & B ＝［ \begin{array}{l} 1 \\ δ \end{array} ］ \end{matrix}$

这一对是可控的，如果<年代pan class="inlineequation"> $δ \neq 0$ 但是,如果<年代pan class="inlineequation"> $δ < \sqrt{e p 年代}$ ,在那里每股收益是相对机器精度。ctrb (A, B)返回

$［ \begin{matrix} B & 一个 B \end{matrix} ］＝［ \begin{matrix} 1 & 1 \\ δ & δ \end{matrix} ］$

这不是完全的等级。对于这样的情况，最好使用ctrbf．

ctrbf|<年代pan itemscope itemtype="//www.tatmou.com/help/schema/MathWorksDocPage/SeeAlso" itemprop="seealso">obsv

之前介绍过的R2006a

下载代码示例