diff scripts/control/mb.m @ 3228:dbcc24961c44

[project @ 1998-12-09 18:42:12 by jwe]
author jwe
date Wed, 09 Dec 1998 18:42:13 +0000
parents ba1c7cdc6090
children 28aba52a2368
line wrap: on
line diff
--- a/scripts/control/mb.m
+++ b/scripts/control/mb.m
@@ -1,4 +1,7 @@
-# $Revision: 1.1 $
+# I think that this m-file can be deleted
+# a.s.hodel@eng.auburn.edu - 4 Dec. 1998
+
+# $Revision: 2.0.0.0 $
 
 Ap = [0 1;1960 0];
 Bp = [0;-6261];
@@ -6,8 +9,8 @@
 Dp = 0;
 
 Gp = ss2sys(Ap,Bp,Cp,Dp,0,2,0,[],"delta_i","delta_y");
-Gp = syschnames(Gp,"st",1,"delta_x1");
-Gp = syschnames(Gp,"st",2,"delta_x2");
+Gp = syssetsignals(Gp,"st","delta_x1",1);
+Gp = syssetsignals(Gp,"st","delta_x2",2);
 
 Ak = [-20 1;-22160 -200];
 Bk = [-20;-2160];
@@ -15,8 +18,8 @@
 Dk = 0;
 
 Gk = ss2sys(Ak,Bk,Ck,Dk,0,2,0,[],"y","i");
-Gk = syschnames(Gk,"st",1,"x1");
-Gk = syschnames(Gk,"st",2,"x2");
+Gk = syssetsignals(Gk,"st","x1",1);
+Gk = syssetsignals(Gk,"st","x2",2);
 
 Gc = sysgroup(Gp,Gk);
 
@@ -34,7 +37,10 @@
 # Gc = sysprune(Gc,1,1);
 
 is_stable(Gc)
-eig(Gc.a)
+Gca = sys2ss(Gc);
+eig(Gca)
 
-Acl = [Gp.a, -Gp.b*Gk.c; Gk.b*Gp.c, Gk.a]
+[Gpa,Gpb,Gpc,Gpd] = sys2ss(Gp);
+[Gka,Gkb,Gkc,Gkd] = sys2ss(Gk);
+Acl = [Gpa, -Gpb*Gkc; Gkb*Gpc, Gka]
 eig(Acl)