后续陆续更新:传递函数讲解、工作点线性化、系统辨识、LMS
一、从函数数学角度看互质分解
互质,代数学术语。指几个整数或几个多项式的一种关系。
从整数来看,最大公约数是1的几个整数称为互质的整数.
从多项式来看,除常数(零次多项式)外无其他公因式的几个多项式称为互质多项式。
两个多项式 f ( x ) f(x) f(x) 与 g ( x ) g(x) g(x) 互质的充分必要条件是存在多项式 u ( x ) u(x) u(x) 与 v ( x ) v(x) v(x) ,使 f ( x ) u ( x ) + g ( x ) v ( x ) = 1 f(x)u(x)+g(x)v(x)=1 f(x)u(x)+g(x)v(x)=1。
从矩阵角度来看,若两个矩阵 M M M 和 N N N 有相同的列数且存在矩阵 X r X_r Xr 和 Y r Y_r Yr 使得 [ X r Y r ] [ M N ] = X r M + Y r N = I \begin{bmatrix} X_r&Y_r \end{bmatrix}\begin{bmatrix} M\\N \end{bmatrix}=X_rM+Y_rN=I [XrYr][MN]=XrM+YrN=I,则称矩阵 M 、 N M、N M、N 右互质。
同理,若两个矩阵 M ~ \tilde M M~ 和 N ~ \tilde N N~ 有相同的行数且存在矩阵 X l X_l Xl 和 Y l Y_l Yl 使得 [ M ~ N ~ ] [ X l Y l ] = M ~ X l + N ~ Y l = I \begin{bmatrix}\tilde M&\tilde N \end{bmatrix}\begin{bmatrix} X_l\\Y_l \end{bmatrix}=\tilde MX_l+\tilde N Y_l=I [M~N~][XlYl]=M~Xl+N~Yl=I,则称矩阵 M 、 N M、N M、N 左互质。
以上三种是数学中常用的关于互质的定义。
二、从控制系统角度看互质分解
1.传递函数视角
如果 G ( s ) = N ( s ) D ( s ) = a 0 + ⋯ + a n s n b 0 + ⋯ + b m s m G(s)=\frac{N(s)}{D(s)}=\frac{a_0+\cdots+a_ns^n}{b_0+\cdots+b_ms^m} G(s)=D(s)N(s)=b0+⋯+bmsma0+⋯+ansn 满足 G ( ∞ ) = 0 G(\infty)=0 G(∞)=0,则称其为有理传递函数。
G ( s ) G(s) G(s) 的一个右互质分解(RCF)为 G ( s ) = N ( s ) M − 1 ( s ) G(s)=N(s)M^{-1}(s) G(s)=N(s)M−1(s),一个左互质分解(LCF)
为 G ( s ) = M ~ − 1 ( s ) N ~ ( s ) G(s)=\tilde M^{-1}(s)\tilde N(s) G(s)=M~−1(s)N~(s),且存在 X r , Y r , X l , Y l X_r,Y_r,X_l,Y_l Xr,Yr,Xl,Yl 使得 [ X r Y r − N ~ M ~ ] [ M − Y l N X l ] = I \begin{bmatrix} X_r & Y_r \\ -\tilde N&\tilde M \end{bmatrix}\begin{bmatrix} M & -Y_l \\ N& X_l \end{bmatrix}=I [Xr−N~YrM~][MN−YlXl]=I,则称 G ( s ) G(s) G(s) 有双互质分解。
上述左右互质分解皆有Matlab相关函数可进行调用。
2.状态方程视角
下讲解以状态方程描述的控制系统左右互质分解实现。
G ( s ) = [ A B C D ] = ( A , B , C , D ) G(s)=\begin{bmatrix} A & B \\ C&D \end{bmatrix}=(A,B,C,D) G(s)=[ACBD]=(A,B,C,D)为传递函数的一个最小状态空间实现,令 F F F 和 L L L 使得 A + B F A+BF A+BF 和 A − L C A-LC A−LC 都稳定(特征值小于0),且定义如下,
左互质分解相关传函状态空间实现为
M
~
=
(
A
−
L
C
,
−
L
,
C
,
I
)
N
~
=
(
A
−
L
C
,
B
−
L
D
,
C
,
D
)
X
l
=
(
A
+
B
F
,
L
,
C
+
D
F
,
I
)
Y
l
=
(
A
+
B
F
,
−
L
,
F
,
0
)
\tilde M=(A-LC,-L,C,I)\\ \tilde N=(A-LC,B-LD,C,D)\\ X_l=(A+BF,L,C+DF,I)\\ Y_l=(A+BF,-L,F,0)
M~=(A−LC,−L,C,I)N~=(A−LC,B−LD,C,D)Xl=(A+BF,L,C+DF,I)Yl=(A+BF,−L,F,0)
右互质分解相关传函状态空间实现为
M
=
(
A
+
B
F
,
B
,
F
,
I
)
N
=
(
A
+
B
F
,
B
,
C
+
D
F
,
D
)
X
r
=
(
A
−
L
C
,
−
B
+
L
D
,
F
,
I
)
Y
r
=
(
A
−
L
C
,
−
L
,
F
,
0
)
M=(A+BF,B,F,I)\\ N=(A+BF,B,C+DF,D)\\ X_r=(A-LC,-B+LD,F,I)\\ Y_r=(A-LC,-L,F,0)
M=(A+BF,B,F,I)N=(A+BF,B,C+DF,D)Xr=(A−LC,−B+LD,F,I)Yr=(A−LC,−L,F,0)
针对一般以状态方程描述的控制系统,常以如上状态空间实现加以仿真。