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
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA -
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
function [gm]=fspecg(g)
// gm=fspecg(g) : stable factorization.
// g and gm are continuous-time linear systems in state-space form.
// gtild(g)*g = gtild(gm)*gm
//with gm stable.
// Imaginary-axis poles are forbidden.
//
//!
[a,b,c,d]=g(2:5),g=[];
a=-a;
b=-b;
h=[-a',c'*c;
0*eye(a),a];
x=ric_desc(h);h=[]
gm=syslin('c',-a'+c'*c*x,-c',b'-d'*c*x,d')';
endfunction
|
