h2des1.m

数字通信第四版原书的例程

% A function for use with the h-infinity Control Toolbox demonstration

%	Copyright (c) 1990-94 by The MathWorks, Inc.

format short 
% Construct the wieghting functions
sys = mksys(a,b,c,d);
[aw1,bw1,cw1,dw1] = tf2ss(den1,num1); sysw1 = mksys(aw1,bw1,cw1,dw1);
[aw2,bw2,cw2,dw2] = tf2ss(num2,den2); sysw2 = mksys(aw2,bw2,cw2,dw2);
[aw3,bw3,cw3,dw3] = tf2ss(den3,num3); sysw3 = mksys(aw3,bw3,cw3,dw3);
clear aw1 aw2 aw3 bw1 bw2 bw3 cw1 cw2 cw3 dw1 dw2 dw3
disp('Use AUGMENT to construct the augmented plant');
disp('   ');
sys_=augss(sys,sysw1,sysw2,sysw3,0);
disp('Build a balanced system using OBALREAL');
disp('   ');
[A,B1,B2,C1,C2,D11,D12,D21,D22]=branch(sys_);
[aa,bb,cc,mm,tt] = obalreal(A,[B1 B2],[C1;C2]);
A = aa; B1 = bb(:,1); B2 = bb(:,2); C1 = cc(1:3,:); C2 = cc(4,:);
sys_ = mksys(A,B1,B2,C1,C2,D11,D12,D21,D22,'tss');
disp('   ');
disp('Use H2LQG to design an H-2 controller');
disp('   ');
[ss_cp,ss_cl]=h2lqg(sys_);
[ae,be,ce,de]=branch(ss_cp);
clear A B1 B2 C1 C2 D11 D12 D21 D22 aa bb cc sysw1 sysw2 sysw3 ss_cp
disp('Done')