projsl - Linear system projection (original) (raw)

Please note that the recommended version of Scilab is 2026.0.1. This page might be outdated.
See the recommended documentation of this function

Scilab help >> CACSD > projsl

Calling Sequence

Arguments

sl,slp

syslin lists

Q,M

matrices (projection factorization)

Description

slp= projected model of sl where Q*M is the full rank factorization of the projection.

If (A,B,C,D) is the representation of sl, the projected model is given by (M*A*Q,M*B,C*Q,D).

Usually, the projection Q*M is obtained as the spectral projection of an appropriate auxiliary matrix W e.g. W = product of (weighted) gramians or product of Riccati equations.

Examples

rand('seed',0);sl=ssrand(2,2,5);[A,B,C,D]=abcd(sl);poles=spec(A) [Q,M]=pbig(A,0,'c');
slred=projsl(sl,Q,M);spec(slred('A')) sl('D')=rand(2,2);
trzeros(sl)
wi=inv(sl);
[q,m]=psmall(wi('A'),2,'d');
slred2=projsl(sl,q,m); trzeros(slred2)

A=diag([-1,-2,-3]); sl=syslin('c',A,rand(3,2),rand(2,3));[nk2,W]=hankelsv(sl) [Q,M]=pbig(W,nk2(2)-%eps,'c');
slr=projsl(sl,Q,M);
hankelsv(slr)

See Also