在控制系统中巧解Diophantine 方程——MATLAB 语言应用.pdfVIP

在控制系统中巧解 Diophantine 方程 ——MATLAB 语言应用 Find The Solutions of Diophantine Equation Smartly With MATLAB 胡巧会 刘贺平 孟庆元 摘要 本文就控制系统中常遇到的Diophantine 方 （2 ）式左边的矩阵称为Sylvester 矩阵，也称 程求解问题，为得到一个可以方便嵌入极点配置控 为结式。Sylvester 矩阵的一个重要的性质是如果 制、最小方差控制、广义最小方差控制、广义预测 A(z −1 ), B (z −1 ) 两个多项式是互质的，则由 控制的仿真程序中的子程序, 介绍了一种基于 A(z −1 ), B (z −1 ) 两个多项式系数组成的Sylvester 矩 MATLAB 语言的简洁的求解方法，这个程序兼容性 阵是满质的。（2 ）式 Sylvester 矩阵中由A(z −1 ) 的 好，而且长度短，可以大大提高运行速度，极为方 便。本文还给出了程序清单。 系数组成的行数为n+n +1,列数为n +1 ，由B (z −1 ) x x 关键词 Diophantine 方程 MATLAB 控制系统 的系数组成的行数为m+n +1 ，列数为n +1 。因此 y y Abstract As far as Diophantine equation in the Sylvester 矩阵为方阵的条件是 control system is concerned, this paper develop a new and smart solution applying MATLAB and the max{ n+n +1, m+n +1}= n +1+n +1 （3 ） x y x y program list is given also. T (z −1 ) 次数应满足 Key words Diophantine equation MATLAB n ≤max{n+ n ,m+ n } （4 ） t x y Control system 于是根据Sylvester 矩阵满秩对A(z −1 ), B (z −1 ) 互质 性的要求，以及公式（3 ）和（4 ）得出确定 1 引言 X (z −1 ), Y (z −1 ) 最小次数解唯一存在的条件是： 在控制器的设计中，经常遇到Diophantine