th = (0:127)/128*2*pi;				%角度采樣點
rho=ones(size(th));					%單位半徑
x = cos(th);y = sin(th);
f = abs(fft(ones(10,1),128));		%對離散方波進行FFT變換，并取幅值。
rho=ones(size(th))+f';				%取單位圓為繪制幅頻譜的基準。
subplot(1,2,1),polar(th,rho,'r')
subplot(1,2,2),stem3(x,y,f','d','fill')	%取菱形離散桿頭，并填色。
view([-65 30])						%控制角度，為表現(xiàn)效果。