scilabで遊ぼう (3) [プログラム三昧]
計算が出来たので、計算結果をグラフ表示してみます。
時間軸での表示
最初は、単純な正弦波を表示させます。 ここでは、1サイクルを20点に分割して4サイクル分(8π分)表示します。 そのため、0から80まで昇順に並ぶベクトル"t"を用意します。
-->t=[0:80] t = column 1 to 7 0. 1. 2. 3. 4. 5. 6. column 8 to 13 7. 8. 9. 10. 11. 12. column 14 to 19
全部で81項目もあるので、表示しきれず、一時停止してしまいます。 こんなときには、コマンドにセミコロンを追加して表示を省略することができます。
-->t=[0:80]; -->size(t) ans = 1. 81.
"size"関数でマトリックスの大きさを確認すると、1行81列であることがわかります。 つぎは、y軸に使うベクトルを用意します。
-->y=sin(2*%pi*t/20);
”%pi"は、円周率πを表現する定数です。 ベクトルを"sin"関数に渡すとベクトルの各要素の正弦を要素とするベクトルができます。
-->scf(1);clf; -->plot(t,y);
"scf"で図番号を指定し、"clf"で図をクリアします。 そして、"plot"でx軸とy軸を指定してグラフを表示します。
同様に余弦のベクトルもグラフ表示します。
-->x=cos(2*%pi*t/20); -->scf(1);clf; -->plot(t,y,t,x);
x軸とy軸の周期を調整するといわゆるリサージュ波形を表示させることができます。
-->t=[0:2100]; -->x=cos(2*%pi*t/700); -->y=sin(2*%pi*t/300); -->scf(1);clf; -->plot(x,y);
付録 : 「scilab で遊ぼう」索引
- scilabで遊ぼう (1)
- scilabで遊ぼう (2)
- scilabで遊ぼう (3)
- scilabで遊ぼう (4)
- scilabで遊ぼう (5)
- scilabで遊ぼう (6)
- scilabで遊ぼう (7)
- scilabで遊ぼう (8)
- scilabで遊ぼう (9)
- scilabで遊ぼう (10)
- scilabで遊ぼう (11)
- scilabで遊ぼう (12)
- scilabで遊ぼう (13)
コメント 0