Mercurial > octave-nkf
diff scripts/control/hinf/dgkfdemo.m @ 4771:b8105302cfe8
[project @ 2004-02-16 17:45:50 by jwe]
author | jwe |
---|---|
date | Mon, 16 Feb 2004 17:45:50 +0000 |
parents | c487fb85b843 |
children | bdbee5282954 |
line wrap: on
line diff
--- a/scripts/control/hinf/dgkfdemo.m Mon Feb 16 16:22:43 2004 +0000 +++ b/scripts/control/hinf/dgkfdemo.m Mon Feb 16 17:45:50 2004 +0000 @@ -27,7 +27,7 @@ function dgkfdemo () save_val = page_screen_output; - page_screen_output = 1; + page_screen_output = 0; while (1) clc sel = 0; @@ -97,7 +97,7 @@ B = [0; 2]; G = eye(2) C = [1, 1]; - sys = ss2sys(A, [B, G], C); + sys = ss(A, [B, G], C); sys = syssetsignals(sys,"in", ... ["control input"; "disturbance 1"; "disturbance 2"]); sysout(sys) @@ -124,7 +124,7 @@ cmd = "pred_poles = sortcom([Er; Ee])"; run_cmd disp("Example 2: discrete-time example") - cmd1 = "Dsys = ss2sys(A, [G, B], C, [0, 0, 0], 1);"; + cmd1 = "Dsys = ss(A, [G, B], C, [0, 0, 0], 1);"; cmd2 = "[K,Q1,P1,Ee,Er] = lqg(Dsys,SigW, SigV,Q,R);"; disp("Run commands:") cmd = cmd1; @@ -148,7 +148,7 @@ cmd = "A = [0, 1; -2, -1]; B = [0; 1]; C = [1, 0]; sys_poles = eig(A)"; run_cmd disp("Put into Packed system form:") - cmd = "Asys = ss2sys(A,B,C);"; + cmd = "Asys = ss(A,B,C);"; run_cmd disp("Evaluate system 2-norm (impulse response energy):"); cmd = "AsysH2 = h2norm(Asys)"; @@ -172,7 +172,7 @@ run_cmd prompt disp("Put into system data structure form:") - cmd="Bsys = ss2sys(A,B,C);"; + cmd="Bsys = ss(A,B,C);"; run_cmd disp("Evaluate 2-norm:") cmd = "BsysH2 = h2norm(Bsys)"; @@ -212,7 +212,7 @@ disp("loop impulse response from w(t) =[w1; w2] to z(t) = [y1; y2]"); prompt disp("First: pack system:") - cmd="Asys = ss2sys(A, [B1, B2], [C1; C2], D);"; + cmd="Asys = ss(A, [B1, B2], [C1; C2], D);"; run_cmd disp("Open loop multivariable Bode plot: (will take a moment)") cmd="bode(Asys);"; @@ -250,7 +250,7 @@ cmd = "A = [0, 1; -2, -1]; B = [0; 1]; C = [1, 0]; sys_poles = eig(A)"; run_cmd disp("Pack into system format:") - cmd = "Asys = ss2sys(A,B,C);"; + cmd = "Asys = ss(A,B,C);"; run_cmd disp("The infinity norm must be computed iteratively by") disp("binary search. For this example, we select tolerance tol = 0.01, ") @@ -274,7 +274,7 @@ cmd = "A = [0, 1; 2, 1]; B = [0; 1]; C = [1, 0]; sys_poles = eig(A)"; run_cmd disp("Pack into system format:") - cmd = "Bsys = ss2sys(A,B,C);"; + cmd = "Bsys = ss(A,B,C);"; run_cmd disp("Evaluate with BsysH2 = hinfnorm(Bsys,tol,gmin,gmax)") BsysH2 = hinfnorm(Bsys,tol,gmin,gmax) @@ -301,7 +301,7 @@ D = [D11, D12; D21, D22] prompt disp("First: pack system:") - cmd="Asys = ss2sys(A, [B1, B2], [C1; C2], D);"; + cmd="Asys = ss(A, [B1, B2], [C1; C2], D);"; run_cmd prompt disp("Open loop multivariable Bode plot: (will take a moment)")