数数字信号处理第二版_实验三.docVIP

第二题 function [Xb,Fb] = downsin(a,f); n=[0:15]; Xb(n+1)=exp(-a.*n).*sin(2*pi*f.*n); F = fft(Xb); Fb=abs(F); clear all; k=0:15; [Xb,Fb]=downsin(0.1,0.0625); subplot(3,2,1); stem(k,Xb); hold; plot(k,Xb); xlabel(n);ylabel(ê±óòì?D?);text(8,0.5,a=0.1,f=0.0625); subplot(3,2,2); stem(k,Fb); hold; plot(k,Fb); xlabel(n);ylabel(·ù?μì?D?);text(10,3,a=0.1,f=0.0625); [Xb1,Fb1]=downsin(0.1,0.4375); subplot(3,2,3); stem(k,Xb1); hold; plot(k,Xb1); xlabel(n);ylabel(ê±óòì?D?);text(8,0.5,a=0.1,f=0.4375); subplot(3,2,4); stem(k,Fb1); hold; plot(k,Fb1); xlabel(n);ylabel(·ù?μì?D?);text(10,3,a=0.1,f=0.4375); [Xb2,Fb2]=downsin(0.1,0.5625); subplot(3,2,5); stem(k,Xb2); hold; plot(k,Xb2); xlabel(n);ylabel(ê±óòì?D?);text(8,0.5,a=0.1,f=0.5625); subplot(3,2,6); stem(k,Fb2); hold; plot(k,Fb2); xlabel(n);ylabel(·ù?μì?D?);text(10,3,a=0.1,f=0.5625); 当f=f1/fs=1/16=0.0625时，谱峰位置出现正确。时间长度为16Ts，即采样频率为16f1，满足采样定理。由于截取长度为周期的整数倍（一倍），此处不引入频谱泄漏。但是由于截断效应，相当于加一矩形窗，其频谱中旁瓣导致频率分量扩散，因此有一定程度的频谱泄漏。增加的部分高频分量会造成低频成分的混叠效应，但是由截断带来的高频旁瓣较小，故混叠效应较弱。 当f=f2/fs=7/16=0.4375时，fs=(16/7)f2，仍满足采样定理，谱峰位置正确，且为整数周期采样。同上，较弱的频谱泄漏及混叠效应来自于截断效应。 当f=f2/fs=9/16=0.5625时，fs=(16/9)f2，不满足采样定理，故产生混叠，此时谱峰位置不正确。混叠后的频谱由图知与第二种情况频谱相同，原因是超过1/2fs的高频分量叠加到低频。此时也产生了一定的频谱泄漏。 第四题 clear all; N=16;detf=1/16;n=[0:N-1]; x1(n+1)=sin(2*pi*0.125.*n)+cos(2*pi*(0.125+detf).*n); detf = 1/64;x2(n+1)=sin(2*pi*0.125.*n)+cos(2*pi*(0.125+detf).*n); subplot(2,2,1);stem(n,x1); xlabel(n);ylabel(ê±óòì?D?);text(6,1,N=16,detf=1/16); subplot(2,2,2);stem(n,abs(fft(x1))); xlabel(n);ylabel(·ù?μì?D?);text(6,4,N=16,detf=1/16); subplot(2,2,3);stem(n,x2); xlabel(n);ylabel(ê±óòì?D?);text(6,1,N=16,detf=1/64); subplot(2,2,4);stem(n,abs(fft(x2))); xlabel(n);ylabel(·ù?μì?D?);text(6,4,N=16,detf=1/64); N=128;detf=1/16;n=[0:N-1]; x3(n+1)=sin(2*pi*0.125.*n)+cos(2*pi*(0.125+detf).*n); detf = 1/64; x4(n+1)=sin(2*pi*0.125.*n)+cos(2*pi*(0.125+detf).*n); figure;subplot(2,2,1);stem(n,x3); xlabel(n);ylabel(ê±óòì?D?);axis([0 128 -2 2]);text(6,1.5,N=128,detf=1/16); subplot(