Graphics

極座標でのコンターの表示

極座標系で定義されたデータのコンターを表示することができます。たとえば、極座標でグリッドを設定し、座標をCartesian座標に変換します。

[th,r] = meshgrid((0:5:360)*pi/180,0:.05:1);
[X,Y] = pol2cart(th,r);

つぎに、単位円の内部で複素行列Zを生成します。

Z = X+i*Y;

X、Y、Zは、円の内部の点です。

関数のサーフェスを作成し、表示します。

f = (Z.^4-1).^(1/4);
surf(X,Y,abs(f))

つぎのステートメントを使って、サーフェスの下に単位円を表示します。

hold on
surf(X,Y,zeros(size(X)))
hold off

グラフのラベル付け

つぎのステートメントによりプロットにラベルを付加します。

xlabel('Real','FontSize',14);
ylabel('Imaginary','FontSize',14);
zlabel('abs(f)','FontSize',14);

Cartesian座標でのコンター

また、以下のステートメントによりCartesian座標で、このサーフェスのコンターを表示し、x-軸とy-軸にラベルを付けます。

contour(X,Y,abs(f),30)
axis equal
xlabel('Real','FontSize',14);
ylabel('Imaginary','FontSize',14);

極座標でのコンター

コンターは極座標軸内に表示することもできます。関数polarを使って、極座標軸を作成した後、関数polarで設定したラインを削除します。

h = polar([0 2*pi], [0 1]);
delete(h)

holdをonにして、極座標グリッド上にコンターを表示します。

hold on
contour(X,Y,abs(f),30)

コンターのオフセットの変更

対話型プロット