Octave の使い方その3
数値解析
□ 連立一次方程式の解
x + 2*y = 2, 3*x + 4*y = 3 の解の計算
octave: > A = [1, 2; 3, 4]
octave: > b = [2; 3]
octave: > A\b
最小二乗の意味で最適な解を求めることもできる.
□ 多項式
多項式 x^3 - 8*x^2 -3*x + 9 の表現.
octave: > p = [1 -8 -3 9]
零点の計算(求根)
octave: > r = roots(p)
零点から多項式を決定
octave: > poly(r)
特性多項式(|sI-A|)の計算
octave: > A = [1, 2, 3; 5, 4, 6; 2, 1, 3]
octave: > poly(A)
多項式の積
s^2 + 2*s + 1 と s^3 + 2*s^2 + 2*s + 1 の積
octave: > conv([1, 2, 1], [1, 2, 2, 1])
□ 最小二乗近似
x = 1 y = 4
2 5
3 7
4 9
5 12
に対する近似.
octave: > x = [1, 2, 3, 4, 5]
octave: > y = [4, 5, 7, 9, 12]
直線近似
octave: > polyfit(x, y, 1)
二次関数近似
octave: > polyfit(x, y, 2)
□ ファイルの利用
data.dat という名前で
1 4
2 5
3 7
4 9
5 12
というデータが保存されているとする.
octave: > load data.dat
octave: > data
octave: > data(:,1)
octave: > data(:,2)
octave: > polyfit(data(:,1), data(:,2), 1)
octave: > sum(data(:,2))
octave: > mean(data(:,2))
octave: > plot(data(:,1), data(:,2), 'x')
octave: > plot(data(:,1), data(:,2))
octave: > closeplot
Octave
naniwa@rbt.his.u-fukui.ac.jp