diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /104/CH5/EX5.13/5_13.sce | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '104/CH5/EX5.13/5_13.sce')
| -rwxr-xr-x | 104/CH5/EX5.13/5_13.sce | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/104/CH5/EX5.13/5_13.sce b/104/CH5/EX5.13/5_13.sce new file mode 100755 index 000000000..02f1f03b0 --- /dev/null +++ b/104/CH5/EX5.13/5_13.sce @@ -0,0 +1,24 @@ +//OCF form
+A=[1 2 1;0 1 3;1 1 1]
+B=[1;0;1]
+C=[1 1 0]
+D=0
+[row,col]=size(A)
+c=s*eye(row,col)-A
+x=det(c)
+r=coeff(x)
+M=[r(1,2) r(1,3) 1;r(1,3) 1 0;1 0 0]
+V=[C;C*A;C*A^2]
+disp(V,"observability matrix=")
+if (det(V)==0) then
+ printf("system cannot be transformed into ocf form")
+else
+ printf("system can be transformed into ocf form")
+end
+Q=inv(M*V)
+disp(Q,"Q=")
+Aocf=inv(Q)*A*Q
+Cocf=C*Q
+B=inv(Q)*B
+disp(Aocf,"Aocf=")
+disp(Cocf,"Cocf=")
\ No newline at end of file |