dtft-dftplot.m

来自「Dsp中的Matlab」· M 代码 · 共 12 行

function[y]=dtft-dftplot(w,Xw,k,Xk)
magXw=abs(Xw);angXw=angle(Xw);
magXk=abs(Xk);angXk=angle(Xk);
subplot(2,2,1);plot(w/pi,magXw);grid
xlabel('以pi为单位的频率');title('幅度部分');ylabel('幅度')
subplot(2,2,3);plot(w/pi,angXw);grid
xlabel('以pi为单位的频率');title('相角部分');ylabel('弧度')
subplot(2,2,2);stem(k,magXk);grid
xlabel('以pi为单位的频率');title('幅度部分');ylabel('幅度')
subplot(2,2,4);stem(k,angXk);grid
xlabel('以pi为单位的频率');title('相角部分');ylabel('弧度')