22用MATLAB语言求下列系统的状态方程（1）状态方程模型参数编写matlab程序如下
um172424de
110355024ABCDtf2ss
umde
2零极点增益编写程序
um172424de
110355024ZPKtf2zp
umde
得到结果Z27306285312730628531i15388P4321K13部分分式形式编写程序
um172424de
110355024RPHresidue
umde
得到结果R40000600002000010000P40000300002000010000H（2）解：）（1）传递函数模型参数：编写程序A225512505225425125025025051251125175025075；B4220C0202D0
umde
ss2tfABCD得到结果
um040000140000220000150000de
10000400006250052500225002零极点增益模型参数A225512505225425125025025051251125175025075；B4220C0202D0ZPKss2zpABCD（3）部分分式形式的模型参数：A225512505225425125025025051251125175025075；B4220C0202D0
umde
ss2tfABCDRPHresidue
umde
261解：m文件程序为h01disp函数的数值解为显示‘’中间的文字dispy同上y1fort0h1my
fdispy显示y的当前值ymmhe
d保存文件da2m在matalb命令行中键入da得到结果函数的数值解为y109000081000729006561059050531404783043050387403487（2）另建一个m文件求解ye在t∈01的数值程序为h01disp函数的离散时刻解为dispyfort0h1yexptdispye
d保存文件da3m在matalb命令行中键入da3函数的离散时刻解为y10904808187074080670306065054880496604493040660367925解：1编辑m文件程序h01disp四阶龙格库塔方法求解函数数值解为dispyy1fort0h1dispyk1yk2yk1h2k3yk2h2k4yk3hyyk12k22k3k4h6e
d保存文件q5m在matlab命令行里键入q5得到结果四阶龙格库塔方法求解函数数值解为y109048081870740806703060650548804966044930406603679
27a10b146c1341635dco
vab；eco
vdce10000800003199007521000f0005100gefg10000800003199008021001000000
f以上是计算闭环传递函数的特征多项式以上是计算闭环传递函数的特征多项式以上是计算闭环传递函数的特征多项式以上是计算闭环传递函数的特征多项式prootsg计算特征多项式的根计算特征多项式的根计算特征多项式的根计算特征多项式的根，，就是闭环，就是闭环就是闭环就是闭环，传递函数的极点传递函数的极点传r