Octave surfc プログラミング例：ガウス関数、カラーマップ、視点変更【コード付き】

Octaveプログラミングにおける surf について

Octaveプログラミングにおいて、関数 surf は、3次元の表面プロットを作成するために使用されます。これは、3つの変数（通常はx座標、y座標、z座標）の間の関係を視覚的に表現するのに役立ちます。

基本的な使い方

surf(X, Y, Z)

ここで：

• Z: 各(x, y)座標に対応する高さ（z座標）を指定する行列です。Zのサイズは、XとYによって定義されるグリッドと同じである必要があります。
• Y: 表面のy座標を指定する行列またはベクトルです。
• X: 表面のx座標を指定する行列またはベクトルです。

詳細

• 生成される表面は、これらの高さに基づいて色付けされます。デフォルトでは、高さが高いほど「暖色」（赤など）、低いほど「寒色」（青など）で表示されますが、これはカラーマップを変更することでカスタマイズできます。
• Zは、このグリッド上の各点における高さを表します。
• XとYがベクトルとして与えられた場合、surfはそれらを基にメッシュグリッドを作成します。具体的には、Xの各要素とYの各要素のすべての組み合わせからなる格子点が生成されます。

より高度な使い方

surf関数は、表面の見た目をさらに制御するためのさまざまなオプションを受け付けます。例えば：

• 軸ラベル、タイトルなどの追加
  xlabel、ylabel、zlabel、title 関数を使用して、プロットに注釈を加えることができます。
• ビューの調整
  view 関数を使用して、プロットの視点を変更できます。
• 線の表示/非表示
  EdgeColor プロパティを設定することで、表面のメッシュラインの色や表示/非表示を制御できます（例：surf(X, Y, Z, 'EdgeColor', 'none') でラインを非表示にします）。
• 照明効果の追加
  light 関数や shading コマンドを使用して、表面に陰影や光沢を加えることができます。
• カラーマップの指定
  colormap(map_name) を使用して、表面の色付けに使用するカラーマップを変更できます（例：colormap(jet)、colormap(gray)）。

例

% xとyの範囲を定義
x = -2:0.2:2;
y = -2:0.2:2;

% メッシュグリッドを作成
[X, Y] = meshgrid(x, y);

% z座標を計算（例：ガウス関数）
Z = exp(-(X.^2 + Y.^2));

% 表面プロットを作成
surf(X, Y, Z);

% 軸ラベルとタイトルを追加
xlabel('x軸');
ylabel('y軸');
zlabel('z軸');
title('ガウス関数の表面プロット');

この例では、meshgrid を使用してxとyのベクトルからグリッドを作成し、そのグリッド上の各点におけるz座標を計算しています。そして、surf 関数を使ってその3次元表面をプロットしています。

---

Octaveプログラミングにおける surfc の一般的なエラーとトラブルシューティング

surfc 関数を使用する際に遭遇しやすいエラーとその対処法を以下に示します。

引数の数の誤り

• 対処法
  関数のヘルプ (help surfc) を確認し、正しい数の引数を渡しているか確認してください。
• 原因
  surfc 関数に渡す引数の数が間違っています。通常は surfc(X, Y, Z) のように3つの行列またはベクトルが必要です。カラーデータを指定する場合は surfc(X, Y, Z, C) のように4つの引数を取ります。
• エラーメッセージの例
  error: 'surfc' called with invalid number of arguments (expected 3 or 4, got N)

引数のサイズの不一致

• 対処法
  size(X), size(Y), size(Z) を使用して各変数のサイズを確認し、矛盾がないようにデータを調整してください。meshgrid 関数を使用して X と Y のグリッドを作成し、それに対応する Z を計算することが一般的です。
• 原因
  
  • X, Y, Z が行列の場合、これらの行列のサイズが一致している必要があります。
  • X がベクトルで長さが m、Y がベクトルで長さが n の場合、Z は m x n の行列である必要があります。
• エラーメッセージの例
  error: surf: X, Y, and Z arguments must be matrices of the same size, or vectors of appropriate lengths

データ型の問題

• 対処法
  class(X), class(Y), class(Z) を使用してデータ型を確認し、必要であれば double() などで数値型に変換してください。
• 原因
  X, Y, Z のデータ型が数値型（double, singleなど）でない場合、surfc が正しく処理できないことがあります。
• エラーメッセージの例
  特に明示的なエラーが出ない場合もありますが、予期しないプロット結果になることがあります。

プロットが表示されない、または意図しない表示になる

• 対処法
  
  • figure() コマンドを使用して新しい Figure ウィンドウを作成してください。
  • 複数のプロットを重ねて表示したい場合は、hold on コマンドを使用してください。
  • データの最小値と最大値を確認し、必要に応じて軸の範囲を調整してください (xlim, ylim, zlim)。
  • colormap() 関数を使用して、適切なカラーマップを設定してください。
• 原因
  
  • 既存の Figure ウィンドウが閉じられている。
  • hold off の後に新しいプロットコマンドを実行している。
  • データの範囲が極端に大きく、プロット範囲から外れている。
  • カラーマップの設定が意図しないものになっている。

等高線が期待通りに表示されない

• 対処法
  contour 関数や contourc 関数を別途使用して等高線を細かく制御するか、surfc のオプションで等高線のレベルを調整できるか確認してください（Octaveのバージョンによって利用可能なオプションが異なる場合があります）。
• 原因
  
  • Z データの変化が小さすぎる、または大きすぎるため、デフォルトの等高線の間隔では適切な表示にならない。

メモリ不足

• 対処法
  データのサイズを減らす、またはより多くのメモリを搭載した環境で実行してください。
• 原因
  X, Y, Z のサイズが非常に大きく、システムのメモリを使い果たしている。
• エラーメッセージの例
  error: out of memory or dimension too large for Octave's index type

1. エラーメッセージをよく読む
   エラーメッセージは問題の原因を示唆していることが多いです。
2. ヘルプを参照する
   help surfc コマンドを実行して、関数の正しい使い方や利用可能なオプションを確認してください。
3. 簡単なデータで試す
   小さなサイズの単純なデータで surfc が正しく動作するかどうかを確認し、問題の切り分けを行います。
4. 変数の内容を確認する
   disp(X), disp(Y), disp(Z) などを使用して、変数の値やサイズが期待通りであるかを確認してください。
5. プロット関連のコマンドを確認する
   figure, hold on/off, subplot, axis, xlabel, ylabel, zlabel, title, colormap など、プロットに関連する他のコマンドの設定も確認してください。

---

基本的な surfc の使用例

この例では、基本的なガウス関数の表面を surfc でプロットします。

% xとyの範囲を定義
x = -2:0.2:2;
y = -2:0.2:2;

% メッシュグリッドを作成
[X, Y] = meshgrid(x, y);

% z座標を計算（ガウス関数）
Z = exp(-(X.^2 + Y.^2));

% 表面と等高線をプロット
figure; % 新しい Figure ウィンドウを作成
surfc(X, Y, Z);

% 軸ラベルとタイトルを追加
xlabel('x軸');
ylabel('y軸');
zlabel('z軸');
title('ガウス関数の表面プロットと等高線');

% グリッド線を追加
grid on;

解説

1. x = -2:0.2:2; と y = -2:0.2:2;: x軸とy軸の値をそれぞれ-2から2まで0.2刻みで定義したベクトルを作成します。
2. [X, Y] = meshgrid(x, y);: meshgrid 関数は、ベクトル x と y を基に、3次元プロットで使用する格子点（メッシュグリッド）を作成します。X と Y は同じサイズの行列となり、それぞれの要素が対応する点のx座標とy座標を表します。
3. Z = exp(-(X.^2 + Y.^2));: 各格子点 (X(i, j), Y(i, j)) におけるz座標（高さ）を計算します。ここではガウス関数を使用しています。.^2 は要素ごとのべき乗を表します。
4. figure;: 新しい描画ウィンドウを開きます。
5. surfc(X, Y, Z);: X, Y, Z のデータを用いて3次元の表面プロットを作成し、底面に等高線を追加します。
6. xlabel('x軸');, ylabel('y軸');, zlabel('z軸');, title('ガウス関数の表面プロットと等高線');: プロットの各軸のラベルとタイトルを設定します。
7. grid on;: プロットにグリッド線を表示します。

カラーマップの変更例

colormap 関数を使って、表面の色付けを変更できます。

% (上記の例と同じデータ生成部分)

figure;
surfc(X, Y, Z);
colormap(jet); % カラーマップを 'jet' に設定
xlabel('x軸');
ylabel('y軸');
zlabel('z軸');
title('ガウス関数の表面プロットと等高線 (jetカラーマップ)');
grid on;

解説

• colormap(jet);: colormap 関数に 'jet' という文字列を渡すことで、カラーマップをジェットスケールに変更します。Octaveには様々な組み込みのカラーマップがあります（例：gray, hot, cool, spring, summer, autumn, winter など）。

等高線のプロパティを変更する例

surfc で描画される等高線のプロパティを直接制御する標準的な方法は、contour 関数を hold on で重ねて使用することです。

% (上記の例と同じデータ生成部分)

figure;
surf(X, Y, Z); % まず表面のみをプロット
hold on;
contour(x, y, Z, 'LevelStep', 0.1, 'LineColor', 'k'); % 等高線を重ねてプロット
hold off;

xlabel('x軸');
ylabel('y軸');
zlabel('z軸');
title('ガウス関数の表面プロットとカスタマイズされた等高線');
grid on;

解説

1. surf(X, Y, Z);: まず surf 関数で表面のみをプロットします。
2. hold on;: これ以降のプロットコマンドを現在の Figure 上に重ねて描画します。
3. contour(x, y, Z, 'LevelStep', 0.1, 'LineColor', 'k');: contour 関数を使って等高線をプロットします。
   • x, y, Z: 表面プロットと同じデータを使用します。
   • 'LevelStep', 0.1: 等高線の間隔を 0.1 に設定します。
   • 'LineColor', 'k': 等高線の色を黒（'k'）に設定します。
4. hold off;: これ以降のプロットコマンドは新しい Figure に描画されます。

異なる視点からのプロット例

view 関数を使うと、3次元プロットの視点を変更できます。

% (上記の例と同じデータ生成部分)

figure;
surfc(X, Y, Z);
view(30, 45); % 仰角 30度、方位角 45度から見た図
xlabel('x軸');
ylabel('y軸');
zlabel('z軸');
title('異なる視点からのガウス関数の表面プロットと等高線');
grid on;

• view(30, 45);: view(az, el) は、方位角 (azimuth) と仰角 (elevation) を度単位で指定して視点を設定します。

---

surf と contour の組み合わせ

最も一般的で柔軟性の高い代替方法は、surf 関数で3次元表面をプロットし、その上に contour 関数で等高線を重ねて描画する方法です。

% (データの準備は surfc の例と同じ)

figure;
surf(X, Y, Z); % 3次元表面をプロット
hold on; % 複数のプロットを重ねて描画するモードをオン

% 底面に等高線をプロット
contour(x, y, Z, 'LevelStep', 0.1, 'LineColor', 'k', 'LineWidth', 1.5);

% 必要に応じて、他の等高線（例えば表面上）を追加することも可能
% contour3(x, y, Z, 10, 'LineColor', 'r'); % 3次元空間に等高線をプロット

hold off; % 描画モードをオフ

xlabel('x軸');
ylabel('y軸');
zlabel('z軸');
title('surf と contour の組み合わせ');
grid on;

利点

• 表面と等高線のプロパティを個別に制御
  それぞれの関数のオプションを使って、色、透明度、照明などを別々に調整できます。
• 等高線の位置の制御
  底面だけでなく、任意のz平面や3次元空間に等高線を描画できます（contour3 関数）。
• 等高線のカスタマイズ
  contour 関数は、等高線のレベル、線の色、線の太さ、ラベルの表示など、surfc よりも詳細なカスタマイズが可能です。

欠点

• 正確に重ね合わせるためには、データの x, y, Z の整合性を保つ必要があります。
• surfc に比べてコードが少し長くなる可能性があります。

pcolor (または imagesc) と contour の組み合わせ

表面を直接3次元で描画する代わりに、pcolor 関数（擬似カラープロット）や imagesc 関数を使って、2次元の色付きグリッドとして表面の高さ情報を表現し、その上に contour で等高線を重ねる方法です。

% (データの準備は surfc の例と同じ)

figure;
pcolor(X, Y, Z);
shading interp; % 色を滑らかに補間

hold on;
contour(x, y, Z, 'LevelStep', 0.1, 'LineColor', 'w'); % 等高線を重ねてプロット
hold off;

xlabel('x軸');
ylabel('y軸');
title('pcolor と contour の組み合わせ');
colorbar; % カラーバーを表示
axis equal; % x軸とy軸のスケールを等しくする

利点

• パフォーマンス
  大規模なデータの場合、surf よりも高速に描画できる可能性があります。
• 2次元的な表現
  3次元的な視覚化よりも、特定の上から見た高さの情報に焦点を当てたい場合に有効です。

欠点

• 視点の変更は基本的にできません。
• 3次元的な形状の把握が surf ほど直感的ではありません。

スライスプロット (slice)

特定の平面でデータを切り取り、その断面を表示する方法です。等高線と直接的な代替ではありませんが、3次元データの内部構造を理解するのに役立ちます。

% 3次元のボリュームデータ (例として簡単な関数) を作成
[x, y, z] = meshgrid(-2:0.2:2, -2:0.2:2, -2:0.2:2);
V = x .* exp(-(x.^2 + y.^2 + z.^2));

% スライスする平面を定義
xs = [-1, 0, 1];
ys = 0;
zs = []; % 全てのz値でスライス

figure;
slice(x, y, z, V, xs, ys, zs);
colorbar;
xlabel('x'); ylabel('y'); zlabel('z');
title('スライスプロット');

利点

• 特定の断面におけるデータの値の変化を詳細に確認できます。
• 3次元データの内部構造を視覚化できます。

欠点

• 表面全体の形状や等高線とは異なる情報を提供します。

isosurface と patch の組み合わせ

特定の等値面を抽出して3次元形状として表示する方法です。これも等高線の直接的な代替ではありませんが、特定の高さを持つ領域を視覚化するのに役立ちます。

% (3次元のボリュームデータ V は上記と同じ)

level = 0.1; % 表示する等値面のレベル

fv = isosurface(x, y, z, V, level);
figure;
patch(fv, 'FaceColor', 'red', 'EdgeColor', 'none');
view(3); axis vis3d;
xlabel('x'); ylabel('y'); zlabel('z');
title(['等値面 (レベル = ', num2str(level), ')']);
lighting gouraud;

利点

• 特定のデータ値を持つ3次元的な領域を明確に視覚化できます。

• 表面全体の形状や連続的な高さの変化を示すのには適していません。