ex8_23.m

华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 1. 安装MATLAB软件。安装时

% ex8_23
% 计算单位矩阵对角化的预对角化网络
G=tf({1.02 -0.52;-0.54 1.04},{[11.76 1],[10.1 1];[10.4 1],[2.6 1]});
syms k11 k12 k21 k22 s s11 s12 s21 s22
G11=G(1,1);G12=G(1,2);
G21=G(2,1);G22=G(2,2);
m1=G11*G22;m2=G12*G21;
mn=conv(m1.num{1},m2.den{1})-conv(m2.num{1},m1.den{1});
md=conv(m1.den{1},m2.den{1});
kp11n=deconv(conv(G22.num{1},md),G22.den{1});
kp11d=deconv(conv(G22.den{1},mn), G22.den{1});
kp11=tf(kp11n,kp11d)
kp12n=deconv(conv(G12.num{1},md), G12.den{1});
kp12d=deconv(conv(G12.den{1},mn), G12.den{1});
kp12=tf(-kp12n,kp12d)
kp21n=deconv(conv(G21.num{1},md), G21.den{1});
kp21d=deconv(conv(G21.den{1},mn), G21.den{1});
kp21=tf(-kp21n,kp21d)
kp22n=deconv(conv(G11.num{1},md), G11.den{1});
kp22d=deconv(conv(G11.den{1},mn), G11.den{1});
kp22=tf(kp22n,kp22d)