File: dcf.sci

package info (click to toggle)
scilab 2.4-1
  • links: PTS
  • area: non-free
  • in suites: potato, slink
  • size: 55,196 kB
  • ctags: 38,019
  • sloc: ansic: 231,970; fortran: 148,976; tcl: 7,099; makefile: 4,585; sh: 2,978; csh: 154; cpp: 101; asm: 39; sed: 5
file content (26 lines) | stat: -rw-r--r-- 834 bytes parent folder | download | duplicates (2) 
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
 
function [n,m,x,y,nt,mt,xt,yt]=dcf(g,polf,polc,tol)
//[n,m,x,y,nt,mt,xt,yt]=dcf(g,[polf,polc,[tol]]) returns eight
//stable systems (n,m,x,y,nt,mt,xt,yt) for the doubly coprime factorization
//
//                 !xt  -yt! ! m  y !
//                 !       !*!      ! = eye
//                 !-nt  mt! ! n  x !
// G must be stabilizable and detectable.
// See copfac for a description of parameters.
//!
// Copyright INRIA
[lhs,rhs]=argn(0),
       n1=contr(g(2),g(3)),n2=contr(g(2)',g(4)'),
       select rhs,
        case 1 then
        polc=-ones(1,n1),polf=-ones(1,n2),tol=1000*%eps,
        case 2 then
        tol=polf,polc=-ones(1,n1),polf=-ones(1,n2),
        case 3 then tol=1000*%eps,
       end,
       [n,m,xt,yt]=copfac(g,polf,polc,tol),
       [nt,mt,x,y]=copfac(g',polc,polf,tol),
       nt=nt',mt=mt',x=x',y=y',