amp118.m

来自「适合于数学建模的系统教程,介绍了系统模型及转换,系统时间响应和动态仿真等,并配合」· M 代码 · 共 17 行

%Matlab Program 1-18
disp('Create system model')
sys=tf([1 11 36 26],[1 14.6 74.96 153.7 99.65]);
sys=zpk(sys)
disp('A lalanced state - space realization')
[sysb,gram]=balreal(sys);
gram

disp('Check the states that cab removed')
disp('Reduce the model to lower order')
sysr=modred(sysb,3:4,'mod');
sysr=zpk(sysr)

%compare the original system to the reduced system
step(sys,'r-')
grid