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

Change language to:
English - Français - 日本語 - Português -

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

Справка Scilab >> CACSD > Model Transformation > equil1

equil1

balancing (nonnegative) pair of matrices

Syntax

[T [,siz]]=equil1(P,Q [,tol])

Arguments

P, Q

two non-negative symmetric matrices

T

nonsingular matrix

siz

vector of three integers

tol

threshold

Description

equil1 computes t such that:

P1=T*P*T' and Q1=inv(T)'*Q*inv(T) are as follows:

P1 = diag(S1,S2,0,0) and Q1 = diag(S1,0,S3,0) with S1,S2,S3 positive and diagonal matrices with respective dimensions siz=[n1,n2,n3]

tol is a threshold for rank determination in SVD

Examples

S1=rand(2,2);S1=S1*S1';
S2=rand(2,2);S2=S2*S2';
S3=rand(2,2);S3=S3*S3';
P=sysdiag(S1,S2,zeros(4,4));
Q=sysdiag(S1,zeros(2,2),S3,zeros(2,2));
X=rand(8,8);
P=X*P*X';Q=inv(X)'*Q*inv(X);
[T,siz]=equil1(P,Q);
P1=clean(T*P*T')
Q1=clean(inv(T)'*Q*inv(T))

See also

  • balreal — balanced realization
  • minreal — minimal balanced realization
  • equil — balancing of pair of symmetric matrices
  • hankelsv — Hankel singular values
Report an issue
<< equil Model Transformation hankelsv >>

Copyright (c) 2022-2023 (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:
Tue Feb 14 15:13:23 CET 2017