维普资讯 / Y(f)：Cr(tl+Du(tl 式 中 (t)∈侬 ” 称为 状态 向量 ，而 日称 为 系 统 的阶次。A∈朋 ，BEJ ，CEj ，DE坤 分别为 常数矩阵．而系统的输入与输出个数分别为 P和 qo 有 了 A、B、C、D常数矩 阵，便可用 MATLAB命 令直观地表示出来，然后利用其提供的函数便可画 出系统 的频率响应曲线。 3 频率响应曲线绘制 实例 实例 l：传递函数表达式 c(s)= 聒10 而 可利 用如下一 组命 令绘 制系统模 型的 Bode 图、Nyquist图和 Nichols图 ，频率响应曲线如 图 1 所示 。 ／x／m =[10】；咖 =[1，2．10】； subploti 1】． (hum．den)； sub,orI222】．nyqu~t(hum．咖 )； subp~t(223】．n／cho／s(iTum．den)； 上述命令 中 subplot命令是将图形窗13分成若 干区域 。 图 1 传递函数模 型的频 率响应 曲线 实例 2：已给 出状态方程如下 ．绘制频率响应 曲线 。 膏f¨ — 5 2 0 0 0 _4 0 0 -3 2 _4 —1 -3 2 0 _4 Y etI=[1．2．3．4Ix(¨ ．∞∈[O．01．100] 输入如下 MATLAB命令 ，得到的频率响应曲 线 见图 2所示 。 W=logspace(一2，2) 应 用技术 ：[{ ．2．0．0；0．-4．0．0；-3，2，-4。-1；-3．2，0．-4]； 日=[1；2；0；1】； C=[1．2．3．4】；D：0 ssub／~otI221)．txxle(A．8．C．D．1． ； ＆幽嘲 f222}．,'tyquistlA．B． D'1， subp／otI2231． (A．B．C．D．1． ； 图 2 状 态方程模型的频率响应曲线 当然 ，我 们亦 可 在 MATLAB Editer／Debugger 窗 口输 入命令 ，以 M文 件的形式存盘 ，然后 在 MAⅡ ．AB Command Windows下 甩 Run Script菜 单 命令运行存人的 M文件 ，一样可以得到所需 的频 率响应曲线 。不过用 M 文件存取运行和修改方便 很 多 。 结 语 利用 MATl‘AB提供的函数，可 以绘制任意 LTI 模 型的 Bode图、Nyquist图和 Nichols图 ，模 型 可 以 是 连 续 的 或 离 散 的 ，也 可 以是 SISO 或 MIMO。 MAⅡ B的 出现．给线性 系统的频域分 析提供 了 一 个强有力的工具 。 参考文献 【1]薛定宇．控制系统计 算机 辅助设计——M B语言厦 应用．北京 ：清华大学出版杜 ．19 【2]程 卫国等．MATIAB5 3精要、编程爰高级应用．北京 ：机 械工业 出版社 ，2OO0 【3]B~yel,IM，Ford ，Mac ow出 ．Amultlvariable toolkitfor line with M Iu B．IEEE Control Systems maga~ne．1989 (收描 日期 ：2001一∞一10) (英文摘要 见第 28再) MODERN COMPUTERⅫ" 9 强 代 计 算 机 ^ 暮 第 一 一 一 期 r ， 维普资讯 / 应 用技术 地提高 了系统的开发效率，使开发 出的系统具有
