exm060418_1.m

来自「matlab mfiles实例分析」· M 代码 · 共 9 行

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])						%控制角度，为表现效果。