ex631.m

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

A=[-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0];
B=[1;0;0;0];C=[1 7 24 24];D=[0];
disp('原系统的极点为');p=eig(A)'
P=[-30;-1.2;-2.4+sqrt(-16);-2.4-sqrt(-16)];                      %P为所期望的极点
K=place(A,B,P)                         %求状态反馈增益 ；或用acker(A,B,P)命令
disp('配置后系统的极点为');p=eig(A-B*K)'
disp('极点配置后的闭环系统为')
sysnew=ss(A-B*K,B,C,D)                                  %（6.3.1）式
step(sysnew/dcgain(sysnew))                                %绘制时间响应曲线