dtft-dftplot.m

Dsp中的Matlab

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('弧度')