ex731.m

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

syms s a n
   G1=1/(s+10);G2=(s+a)/(s+1);                % 、 的符号表达式
 Go=G1*G2;Gb=Go/(1+Go);         % ：开环传递函数， ：闭环传递函数
Go=simplify(Go),Gb=simplify(Gb)           %化简
[num,den]=numden(Gb),                  %读出 的分子和分母的符号表达式
p=solve(den)                            %求出系统闭环的极点
z=solve(num)                            %求出系统闭环的零点
n=length(p);
for i=1:n
a(i)=solve(p(i))                         %求出系统闭环临界稳定极点的a值
end