Scilab Website | Contribute with GitLab | Mailing list archives | ATOMS toolboxes
Scilab Online Help
2025.0.0 - Русский


specfact

spectral factor

Syntax

[W0, L] = specfact(A, B, C, D)

Description

Given a spectral density matrix phi(s):

-1                   -1
R + C*(s*I-A) * B  +  B'*(-s*I-A') * C'  with R=D+D' > 0

specfact computes W0 and L such that W(s)=W0+L*(s*I-A)^-1*B is a spectral factor of of PHI(s), i.e.

phi(s)=W'(-s)*W(s)

Examples

A=diag([-1,-2]);B=[1;1];C=[1,1];D=1;s=poly(0,'s');
W1=syslin('c',A,B,C,D);
phi=gtild(W1,'c')+W1;
phis=clean(ss2tf(phi))
clean(phis-horner(phis,-s)');   //check this is 0...
[A,B,C,D]=abcd(W1);
[W0,L]=specfact(A,B,C,D);
W=syslin('c',A,B,L,W0)
Ws=ss2tf(W);
horner(Ws,-s)*Ws

See also

  • gtild — tilde operation
  • sfact — discrete time spectral factorization
  • fspecg — stable factorization of continuous time dynamical systems
Report an issue
<< rowinout Linear System Factorization sysfact >>

Copyright (c) 2022-2024 (Dassault Systèmes)
Copyright (c) 2017-2022 (ESI Group)
Copyright (c) 2011-2017 (Scilab Enterprises)
Copyright (c) 1989-2012 (INRIA)
Copyright (c) 1989-2007 (ENPC)
with contributors
Last updated:
Thu Oct 24 11:18:32 CEST 2024