Scilab Website | Contribute with GitLab | Mailing list archives | ATOMS toolboxes
Scilab Online Help
2024.1.0 - 日本語


penlaur

行列ペンシルのローラン係数

呼び出し手順

[Si,Pi,Di,order]=penlaur(Fs)

[Si,Pi,Di,order]=penlaur(E,A)

引数

Fs

標準ペンシル s*E-A

E, A

実数正方行列

Si,Pi,Di

実数正方行列

order

整数

説明

penlaur は,無限大における (s*E-A)^-1の第一ローラン係数を 計算します.

s = 無限大 において, (s*E-A)^-1 = ... + Si/s - Pi - s*Di + ...

order = 特異点の次数 (order=index-1).

行列ペンシル Fs=s*E-A は可逆である必要があります.

添字0のペンシルの場合, Pi, Di,... は 0,および Si=inv(E) となります.

添字1のペンシル(order=0)の場合, Di =0 となります.

より大きな添字のペンシルの場合, 項 -s^2 Di(2), -s^3 Di(3),... は次のように指定されます:

Di(2)=Di*A*Di, Di(3)=Di*A*Di*A*Di (最大 Di(order)).

注意

実験的なバージョン: so*E-Aの条件数が悪い場合に 問題を発生します

F=randpencil([],[1,2],[1,2,3],[]);
F=rand(6,6)*F*rand(6,6);[E,A]=pen2ea(F);
[Si,Pi,Di]=penlaur(F);
[Bfs,Bis,chis]=glever(F);
norm(coeff(Bis,1)-Di,1)

参照

  • glever — 行列ペンシルの逆
  • pencan — 行列ペンシルの正準形
  • rowshuff — シャッフルアルゴリズム
Report an issue
<< pencan Matrix Pencil quaskro >>

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:
Mon Jun 17 17:54:17 CEST 2024