(41) 実数 xの n乗根を知りたい。
1. やりたいこと 急に 1001/3 すなわち 100の三乗根 を求めなければならなくなった。 ぱっと値が出てこないので、octaveを使って求めることにした。 ちなみに、二乗根は sqrt で算出する。 octave… 続きを読む »
1. やりたいこと 急に 1001/3 すなわち 100の三乗根 を求めなければならなくなった。 ぱっと値が出てこないので、octaveを使って求めることにした。 ちなみに、二乗根は sqrt で算出する。 octave… 続きを読む »
octaveのコマンドを入力すると、直後にその結果が表示される。 octave:1> a=rand(5) a = 0.948277 0.879068 0.984755 0.911834 0.582972 0.986024… 続きを読む »
始点(x,y)=(1,2)から、XY方向成分(3,2)のベクトルと、 始点(x,y)=(3,5)から、XY方向成分(-1,-4)のベクトルを quiver で描いてみる。 (1) まずは1本ずつ描く場合 octave:1… 続きを読む »
こんなCSVファイルを用意した。 [user@dog-server]$ cat test.csv 10,1 20,6 30,18 40,11 50,14 [user@dog-server$ このCSVファイルを csvr… 続きを読む »
rand()は、一様分布の乱数生成、 randn() は、正規分布に従った乱数生成をしてくれる。 (1) まずは rand() で 200×200個の乱数を生成し、ヒストグラムでその分布を見てみる。 octav… 続きを読む »
eye()を使うと任意のサイズの 単位行列 を作ることができる。 octave:1> eye(2) ans = Diagonal Matrix 1 0 0 1 octave:2> eye(3) ans = D… 続きを読む »
ones()は、任意の次数、要素数で、全要素の値が1の配列を作ることができる。 octave:1> ones(1) ans = 1 octave:2> ones(2) ans = 1 1 1 1 octave… 続きを読む »
zeros()は、任意の次数、要素数で、全要素の値が0の配列を作ることができる。 octave:6> zeros(1) ans = 0 octave:7> zeros(2) ans = 0 0 0 0 oct… 続きを読む »
numel( A ) は、指定した変数Aの全要素数を返す。 octave:1> a=[1 2 3;4 5 6;7 8 9] a = 1 2 3 4 5 6 7 8 9 octave:2> numel(a) a… 続きを読む »
size( A ) は、配列Aの各次元の要素数を見ることができる。 octave:5> a=[1 2 3;4 5 6] a = 1 2 3 4 5 6 octave:6> size(a) ans = 2 3 … 続きを読む »