「Octave入門」の版間の差分

Octave上で関数やコマンドのヘルプを表示するには，

とすると，キーワード（関数名など何でも良い）に関連した説明が表示されます．

Octaveを起動して開くウィンドウにコマンドを打ち込んでいけば利用はできますが，結果をファイルとして残したり，実験データなとを読み込む場合にはディレクトリを意識する必要があります．今現在自分が作業しているディレクトリをカレント・ディレクトリと呼びます．unix系のOSならすでに各自のホーム・ディレクトリで作業し始めているので，問題は少ないですが，Windowsの場合は注意が必要です．

Octave上でスカラー定数は，

のように入力します．要素は実数や複素数，後述の変数（数式）でも構いません．また，要素の区切りはスペース以外にカンマ（, ）でも大丈夫です．行列はベクトルをセミコロン（; ）で改行しながら入力するというイメージです．

Octave上では，実数，複素数，行列，ベクトルの区別なく変数に格納することができます．変数名は英文字で始まる英数字の文字列で，大文字小文字の区別をします．

のように，特定の変数だけをクリアすることもできます．

基本的に数値やコマンド・関数をタイプしてenterを押せば，入力されて実行されます．標準状態では結果がすぐ表示されますが，表示させないようにすることもできます．一番簡単な方法は行の最後にセミコロン（; ）をつけることです．

  >> a =[1 2 3];

四則演算は実数，複素数，行列，ベクトルおよびそれらを内容とする変数に対して実行可能です．ただし，行列，ベクトルの除算は以下のような意味となります．

Octaveではグラフ機能はgnuplotを呼び出して利用しているので，細かなコントロールにはgnuplotのコマンドを理解する必要があります．しかし，octave-forgeのパッケージにはMATLABとの互換性を高めるため，グラフ関連のコマンドが多く追加されていて（例えばgridなど），だいぶ使いやすくなっています．

1次元のベクトルの組を２次元グラフに表現できます．ベクトルXをx軸，Yをy軸としてグラフ化するには，plot(X, Y)とします．

線の種類や色を変えたいときは，変数の後に続けてダブルコーテーションかシングルコーテーションでくくったフォーマット文字列を書けばいいです．

この例は，色をblueとし，線ではなくて+記号でプロットします．その他の色や線類はhelpで調べて下さい．

なお，plot関数の引き数としてx軸, y軸, 線種を繰り返し用いることで複数のグラフを重ね書きすることもできます．

plot関数でグラフを書いた後，xlabel, ylabel, titile関数で軸ラベルとタイトルをつけることができます．

また，グリッド（格子）をつけるには，

プロットの説明（凡例）が右上に自動的にでますが，これはlegend関数で変更できます．不要な場合は，

で消えます．

なお，環境によってはlegend関数がうまく機能しない場合があります．その場合は，plot関数でプロットする際に，色や線種を指定するのと同時に説明文をセミコロンで囲むことで指定できます．

とすることで，x軸y軸それぞれの範囲を指定できます．なお，axis関数はグラフをプロットした後に実行しても効果がないので，後から実行した場合はもう一度プロットし直す必要があります．また，自動調整に戻す場合は，

   

とします．

plot関数のかわりにsemilogx関数を使うと片対数（x軸だけ対数）のグラフになります．両対数はloglog関数です．線種の指定などはplot関数と同じです

グラフをTeX文書などに取り込む場合は，EPSファイルとしていったん保存します．print関数を使うことでファイルに出力できます．

これまでのようにコマンドを1行1行打ち込んでいっても作業はできますが，ある程度まとまった処理や繰り返し試行する処理などは，面倒ですし，間違いも見つかりにくいです．そこで，その一連の作業をmファイルというテキストファイルとして保存しておくと，実行や修正が非常に楽になります．

mファイルはただのテキストファイルですから，メモ帳など使い慣れたエディタで作成・編集できます．保存する場所は自分のホームディレクトリ以下の適当なところとします．

ファイル名は変数名と同じく数字で始まらない英数字とし，拡張子は.mとします．実行はファイル名から.mを抜いた部分を入力すればいいです．カレント・ディレクトリにその名前のファイルがあれば実行されます．

自作のmファイルをあるディレクトリにまとめておき，カレント・ディレクトリがどこであっても，そのmファイルを実行できるようにすることができます．汎用性の高い関数やツールを作り，将来にわたり使いたい場合に便利です．

とすることで現在のpath設定が表示されます．（ゼミ開講現在，Windows環境ではなぜか表示されません．）これらに，自身のディレクトリを''追加''するには以下のようにします．

この例では，Documentsディレクトリの下のm_filesディレクトリにmファイルを入れておけば，カレントディレクトリがどこであっても実行できるようになります．なお，間違って

とすると，pathがm_filesだけになり，自作のコマンド以外は全く実行できなくなりますので，注意です．（一度Octaveを終了して再度起動すれば元に戻ります）

Octave上では，システムは伝達関数表現（分子分母それぞれの係数を降べきの順に並べた2つのベクトルの組），状態空間表現（4つの行列の組），システム行列表現（1つの行列で表現）のいずれかで表現できます．しかし，後述する解析や設計の関数はシステム行列表現を引数としてとるので，その他の表現で与えられた物をシステム行列表現に変換する必要があります．変換を行う関数はtf()です（古いバージョンのoctaveの場合はtf2sys）．

代表例としてボード線図を書いてみます．関数はbodeを利用します．

一旦Octaveを終了して次回続きをやりたいときなど，変数の内容をとっておきたいことがあります．スクリプトファイルとして作業手順を保存することはできますが，時間のかかる計算が必要な場合その度に計算し直すのも面倒です．

  >> a = [1 2 3];

  >> B = [1 2 3; 4 5 6; 7 8 9];

とすれば，ファイルalldata.datに，現在のすべての変数の内容がアスキー形式（テキスト形式）で書き出されます．メモ帳など適当なテキストエディタで内容を確認できるはずです．なお，バイナリ形式で保存する場合は-binaryというオプションをつけます．テキスト形式にくらべて精度が高く，ファイルサイズも小さくなります．

のように-forceオプションをつけます．

すべての変数を保存する必要がない場合や，一つのファイルに保存してある複数の変数のうちのいくつかだけを復元したい場合もあります．そのような場合，saveやloadコマンドでファイル名の後に変数名を羅列します．

できます．なお，先と同様にloadコマンドではファイルの形式は自動的に判別されるので，オプションは通常不要です．

実験結果をOctaveで処理したりグラフ化する場合，実験用のプログラム（通常はC言語などで自作するでしょう）側はOctaveに渡したいデータをテキスト形式のファイルとして保存するのが簡単です．

例えば，次のような内容のテキストファイルresult.txtがカレントディレクトリにあるとします．

  >> y = result(:,3);

まず，書き出したい変数を結合して一つの変数にしておきます．このときサイズ（行数）が揃っていないとエラーになりますので，揃えるか，一緒にすることをあきらめて別のファイルにするか決めます．シミュレーション結果などは，時刻ベクトル，入力ベクトル，出力ベクトルなどが大抵同じ行数でしょうから，まとめて1つの行列にできるでしょう．

  >> outfile = [t u y];

保存されたファイルの最初の5行は，日付や変数名，サイズなどの情報で先頭文字が#になっています．受け取る側のアプリケーションでその行を無視するように設定できれば，そのまま使用できますが，そうでない場合はテキストエディタで削除してから読み込ませる必要あります．

スクリプト実行中にキーボードから数値や文字列を入力させるにはinput命令を使います．同時にメッセージとして表示させる文字列を指定することができます．

画面への出力はC言語と同じくprintf命令を使います．

条件判定にはC言語と同じくif命令が利用できます．ただし，C言語のように括弧で条件式は囲みません．また，中括弧（{ }）は使用せず，endifでif命令の範囲を指定します．さらに，C言語ではelseとifを分けて書きますが，Octaveではelseifと続けて書きます．

  if a > 0

  elseif a > -5

  else

while命令は条件が成り立っている間，命令群を繰り返します．やはり中括弧は用いずに，命令群の終わりはendwhileで指定します．なお，中断は\break命令です．

  while 条件

  　命令群

  sum = 0;

  while  sum < 100

  　sum = sum + a;

for命令はC言語のforを限定的にした感じで，インデックスとする変数の初期値，増分，終値を指定して命令群を繰り返します．命令群の終わりはendfor，中断はbreakなのはwhileと同様です．なお，増分を省略すると1ととられます．

  　if sum >100

  　　break

初期値:増分:終値という書き方は，ベクトルの作成にも使えますので覚えておくと便利です．

スクリプト内でユーザ独自の関数を定義したり，別ファイルで定義したりできます．

  % 関数の説明

  　命令群

関数の終わりを示すためにendfunction命令が必要です．また，function命令直後のコメント文はhelpで表示されるので，簡単な説明を入れておくと良いでしょう．

なお，C言語の関数と同様に，関数内で変数はローカル変数となるので，変数名がかぶっても内容は別物となることに注意が必要です．

  c  = a - b;

  err = c' * c;