数値解析と制御のためのScilab入門

数値解析と制御のためのScilab入門では、行列、微分方程式、制御などの基礎的な計算を紹介しています。また、計算結果をScilabを用いたグラフィック表示による可視化について紹介しています。

*

行列の成分を取り出す tril|scilab入門

   

scilab2

行列の成分を取り出す場合には、diag、tril、triuを用います。
・diagは行列の対角の要素を取り出します。
・trilは行列の対角の要素と下の要素を取り出します。
・triuは行列の対角の要素と上の要素を取り出します。

SPONSORED LINK

次のような行列の対角成分、下の成分,上の成分を取り出してみたいと思います。
行列

の対角要素、対角要素の上の要素、対角要素の下の要素を取り出すことにします。

1. A=[4 1 5; 3 -1 1; 2 -1 -3];
2. D=diag(A)
3. L=tril(A)
4. U=triu(A)

 

1~4までの数字は打ち込まないで下さい。
->A=[4 1 5; 3 -1 1; 2 -1 -3];
–>D=diag(A)
D =
4.
– 1.
– 3.
–>L=tril(A)
L =
4. 0. 0.
3. – 1. 0.
2. – 1. – 3.

–>U=triu(A)
U =
4. 1. 5.
0. – 1. 1.
0. 0. – 3.

実行するとこのようになります。

行列 tril

方程式を解く場合にLU分解の手法を取り入れているので
ざっくばらんに理解しておくといいです。


 - 数値計算, 行列計算 , , ,

        

SPONSORED LINK

SPONSORED LINK

  関連記事

連立方程式
連立一次方程式を解いてみる|scilab入門

逆行列 inv と (¥)を用いて、簡単なの連立一次方程式を解いてみました。 次 …

no image
行列の積|scilab入門

今回は行列の積について実際に計算してみました。 簡単な計算でも実際計算してみない …

多項式
多項式が有利式の場合の計算

フリーのシュミレーションソフトscilabを用いて色々な計算ができます。 常にバ …

gyoretu2-13 単位行列
ゼロ行列や定数行列を作る |scilab入門

ゼロ行列や定数行列を作ることにします。  ゼロ行列:zeros  定数行列:on …

2次元プロットにグリッドを追加
2次元プロットにグリッドを追加:xgrid|グラフィック

グラフのデータを比較したいとか、見やすくしたい場合には、グリッドを入れて表示させ …

線形システムの定義
多項式を用いた線形システムの定義 syslin

多項式による伝達関数や線形すステムの定義は、syslinを用いて定義できます。 …

行列の計算
行列の差の計算

行列の引き算についてです。 2つの行列の和について考えします。例として3行3列の …

ifによる分岐と関係演算子・論理演算子
ifによる分岐と関係演算子・論理演算子

繰り返し処理に if による分岐と関係演算子・論理演算子 を用いることにします。 …

転置行列
転置行列の表示 |scilab

転置行列の表示について 行列 Bの転置行列は、B にt を付けて表します。行の要 …

3次スプライン補間 :splin
3次スプライン補間 splin について!

計測したデータを補間して表示させたい場合には、3次スプライン補間 splinを用 …