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

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

*

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

      2015/04/26

連立方程式

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

次のような連立一次方程式

連立方程式を解く

を行列を用いて表示することにします。

SPONSORED LINK

A, x, b を次のように

連立方程式を解く

定めると、連立方程式は

Ax = b

と表されます。行列 xは

x=A(-1)b

として求めることができます。これを、scilabで計算する場合には、

x=A\b

のように表します。ただし、Windowsの場合には「バックスラッシュ」が「¥ 」で表示されます。
x=A\b は、x=inv(A)*b と同じです。

この行列は A,x,bで表すと次のようになります。

連立方程式を解く

1. A=[ 1 -2 2; 2 -4 5; 1 -3 -1];
2. b=[2 1 7]’;
3. x=A\b

–>  A=[ 1 -2 2; 2 -4 5; 1 -3 -1];
–> b=[2 1 7]’;
–> x=A\b
x =
16.
4.
– 3.

(2)x=A\bとy=inv(A)*bの比較

1. A=[ 1 3 2; 2 5 8; 3 9 7];
2. b=[5 12 12]’;
3. x=A\b
4. y=inv(A)*b

–>A=[ 1 3 2; 2 5 8; 3 9 7];
–> b=[5 12 12]’;
–> x=A\b
x =
53.
– 14.
– 3.
–>y=inv(A)*b
y =
53.
– 14.
– 3.

同じ結果が得られます。
実際に計算して画面は、図のようになります。

連立方程式の計算結果

大きなサイズの連立方程式を解く場合には、工夫が必要です。

scilabの良いところは、数式レベルで記述できることです。物事を考える場合イメージしやすいし、簡単な計算で結果を把握できます。


 - 数値計算, 行列計算 ,

        

SPONSORED LINK

SPONSORED LINK

  関連記事

行列
行列演算の基礎計算 |scilab入門

行列演算はScilabの基本となっています。行列を用いることにより数式レベルで扱 …

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

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

関数のオンライン定義deff
関数のオンライン定義 deff を使ってみる!

関数のオンライン定義を用いると、数式だけを呼び出して計算ができます。 関数 fu …

no image
逆行列 inv|scilab入門

逆行列は、方程式を解くときに使います。 逆行列 inv を使用して解く場合につい …

行列の計算
行列の和の計算|scilab入門

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

scilab 絶対値
絶対値の計算 abs|scilab入門

絶対値の計算:abs の計算の紹介です。  5 – 9 の計算は   …

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

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

線形システムの定義
多項式による定義 poly

多項式の計算について。 Silabでは多項式と有理式を定義できます。また、有理式 …

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

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

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

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