ex732.m

来自「详细介绍了Mathlab的使用方法与技巧」· M 代码 · 共 9 行

syms s
   G=1/(s*(s+2));H=1/(s+1);                 % 、 的符号表达式
   Gb=G/(1+G*H);Gb=simplify(Gb)          %闭环传递函数 的符号表达式
   [num,den]=numden(Gb)                  %读出 的分子和分母的符号表达式
z=solve(num);p=solve(den);               %z：零点；p：极点
z=numeric(z),p=numeric(p)               %将零点和极点数字化
num=sym2poly(num);den=sym2poly(den);   %将符号多项式转化为多项式的系数向量
step(num,den)                          %绘出系统的单位阶跃响应图