1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
function [S,Series]=sysfact(Sys,Gain,flag)
//For Sys=[A,B,C,D],
//if flag=='post' returns S = [A+B*Gain, B , Gain, I]
// and Series = minss(Sys*S)
//if flag=='pre' returns S=[A+Gain*C, Gain , C, I]
// and Series = minss(S*Sys)
if flag=='post' then
S=syslin(Sys('dt'),Sys('A')+Sys('B')*Gain, Sys('B'), Gain, eye(Gain*Sys('B')));
//Series=minss(Sys*S);
Series= syslin('c',Sys('A')+Sys('B')*Gain, Sys('B'), Sys('C')+Sys('D')*Gain, Sys('D'));
end
if flag=='pre' then
S=syslin(Sys('dt'),Sys('A')+Gain*Sys('C'), Gain, Sys('C'), eye(Sys('C')*Gain));
//Series=minss(S*Sys);
Series=syslin('c',Sys('A')+Gain*Sys('C'), Sys('B')+Gain*Sys('D'), Sys('C'), Sys('D'));
end
//Example:
//Sys=ssrand(3,2,4);Sys('D')=rand(3,2);
//S=sysfact(Sys,lqr(Sys),'post');
//ww=minss(Sys*S);
//ss2tf(gtild(ww)*ww),Sys('D')'*Sys('D')
//Syst=Sys';
//S1=sysfact(Syst,lqe(Syst),'pre');
//ww=minss(S1*Syst);
//ss2tf(ww*gtild(ww)),Syst('D')*Syst('D')'
//Sys=ssrand(2,3,4);
//[X,d,F,U,k,Z]=abinv(Sys);
//ss2tf(Sys*Z)
//ss2tf(Sys*sysfact(Sys,F,'post')*U)
endfunction
|
