Mercurial > octave-nkf

comparison scripts/control/hinf/hinfsyn_chk.m @ 7131:a184bc985c37

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-11-08 15:55:02 by jwe]
author jwe
date Thu, 08 Nov 2007 15:55:02 +0000
parents f084ba47812b
children
comparison
equal deleted inserted replaced
7130:5eeb46c784d7 7131:a184bc985c37
86 86
87 Pc = Pf = []; 87 Pc = Pf = [];
88 88
89 ## Construct the two Hamiltonians 89 ## Construct the two Hamiltonians
90 g2 = 1/(g*g); 90 g2 = 1/(g*g);
91 Hc = [ A , g2*B1*B1' - B2*B2'; -C1'*C1 , -A']; 91 Hc = [A, g2*B1*B1' - B2*B2'; -C1'*C1, -A'];
92 Hf = [ A' , g2*C1'*C1 - C2'*C2; -B1*B1' , -A]; 92 Hf = [A', g2*C1'*C1 - C2'*C2; -B1*B1', -A];
93 93
94 ## check if Hc, Hf are in dom(Ric) 94 ## check if Hc, Hf are in dom(Ric)
95 Hcminval = min (abs (real (eig (Hc)))); 95 Hcminval = min (abs (real (eig (Hc))));
96 Hfminval = min (abs (real (eig (Hf)))); 96 Hfminval = min (abs (real (eig (Hf))));
97 if (Hcminval < ptol); 97 if (Hcminval < ptol);