11/06/2006

9.8.3極座標指令 polar、ezpolar

極座標指令則是在一般二維座標中產生極座標之繪圖。其語法為:


polar(theta, rho, LineSpec)

其第一項輸入參數theta為水平角度,以弧度表示,rho則為向量。若rho隨theta而變化時,則其在極座標軸上亦會產生變化。下面之八葉扇即為其應用之一例:

theta = linspace(0, 2*pi);
polar(theta, cos(4*theta)); % 進行極座標繪圖


與polar對應,ezpolar亦可作極座標函數之繪圖,只是其所用的參數係以字串轉換為函數之內容,而自變數之範圍則可自訂,或以[0, 2*pi]為預設範圍。其語法如下:

ezpolar('function', [a, b])


例如:
ezpolar('sin(4*t)*cos(3*t)')  




以另一字串參數取代原方程式字串之方法,亦可使用。例如:

r = '100/(100+(t-1/2*pi)^8)*(2-sin(7*t)-1/2*cos(30*t))';
ezpolar(r,[-pi/2,3*pi/2])



利用inline轉換函數,亦可直接呼叫其手把,繪出其函數圖:

h=inline('sin(t)/t');
ezpolar(h,[-6*pi 6*pi])