2） 计算 x(n) 的 DTFT 为 X (e ʲw ) ，绘出它的幅度和相位图;
9.1： for N=16:16:32 n=0:15; x1=cos(5*n*pi/16); n=0:N-1;k=n; x=zeros(1,N); x(find(n<=15))=x1; y=x*exp(-j*2*pi/N).^(n'*k); w=linspace(-2*pi,2*pi,500); y1=x*exp(-j*n'*w); figure; subplot(3,1,1);stem(n,x,'.');ylabel('x(n)'); subplot(3,1,2);plot(w/pi,abs(y1));xlabel('X pi');ylabel('X(w)'); subplot(3,1,3);stem(k,abs(y),'.');ylabel('mag X(k)'); end 9.2: n=0:11; x=[1,2,3,4,5,6,6,5,4,3,2,1]; k=n; N=length(n); y=x*exp(-j*2*pi/N).^(n'*k);%DFT w=linspace(-2*pi,2*pi,500); y1=x*exp(-j*n'*w);%DTFT figure; subplot(4,1,1);stem(n,x,'.');ylabel('x（n）'); subplot(4,2,3);stem(k,abs(y),'.');ylabel('mag X(k)'); subplot(4,2,4);stem(k,angle(y),'.');ylabel('ang X(k)'); subplot(4,1,3);plot(w/pi,abs(y1));xlabel('X pi');ylabel('X(jw)'); subplot(4,1,4);plot(w/pi,angle(y1));xlabel('X pi');ylabel('arg(jw)'); figure;%mag stem(2*k/N,abs(y),'filled'); hold on; plot(w/pi,abs(y1));xlabel('X pi');title('mag'); figure;%arg stem(2*k/N,angle(y),'filled'); hold on; plot(w/pi,angle(y1));xlabel('X pi');title('arg');
版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请联系我们举报，一经查实，本站将立刻删除。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/228128.html原文链接：https://javaforall.net
