File: sysdiag.sci

package info (click to toggle)
scilab 4.0-12
  • links: PTS
  • area: non-free
  • in suites: etch, etch-m68k
  • size: 100,640 kB
  • ctags: 57,333
  • sloc: ansic: 377,889; fortran: 242,862; xml: 179,819; tcl: 42,062; sh: 10,593; ml: 9,441; makefile: 4,377; cpp: 1,354; java: 621; csh: 260; yacc: 247; perl: 130; lex: 126; asm: 72; lisp: 30
file content (27 lines) | stat: -rw-r--r-- 634 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
27
function r=sysdiag(varargin)
//Returns the block-diagonal system made with subsystems put in the main
//diagonal
// Syntax:
// r=sysdiag(a1,a2,...,an)
//
// ai    : subsystems (i.e. gains, or linear systems in state-space or
//                     transfer form)
//Remark:
//  At most 17 arguments...
//Example
// s=poly(0,'s')
// sysdiag(rand(2,2),1/(s+1),[1/(s-1);1/((s-2)*(s-3))])
// sysdiag(tf2ss(1/s),1/(s+1),[1/(s-1);1/((s-2)*(s-3))])
// Copyright INRIA

//!
r=varargin(1);
[m1,n1]=size(r);
for k=2:size(varargin)
  ak=varargin(k)
  [mk,nk]=size(ak);
  r=[r,0*ones(m1,nk);0*ones(mk,n1),ak]
  m1=m1+mk
  n1=n1+nk
end
endfunction