1、c)三角波序列d)反三角波序列上机实验内容:(1)观察高斯序列的时域和幅频特性,固定信号中参数,改变的值,使分别等于2、4、8,观察他们的时域和幅频特性,了解当取不同值时,对信号的时域和幅频特性的影响;固定,使分别等于8、13、随着q值的增大,时域信号幅值变化缓慢,频域信号频谱泄露程度减小。随着p的增大,时域信号幅值不变,会在时间轴移位。(2)观察衰减正弦序列的时域和幅频特性,检查普峰出现的位置是否正确,注意频谱的形状,绘出幅频特性曲线,改变分别等于0.4375和0.5625,观察这两种情况下,频谱的形状和普峰出现的位置,有无混叠和泄漏现象?说明产生现象的原因。解答: n=0:1:15; xn
2、=exp(-0.1*n).*sin(2*pi*0.0625*n); subplot(1,2,1);stem(n,xn);xlabel(t/T);ylabel(x(n) xk1=fft(xn);xk1=abs(xk1); subplot(1,2,2);stem(n,xk1);kX(k) xn=exp(-0.1*n).*sin(2*pi*0.4375*n); xn=exp(-0.1*n).*sin(2*pi*0.5625*n);(3)观察三角波和反三角波的时域和幅频特性,用点FFT分析信号序列和的幅频特性,观察两者的序列形状和频谱曲线有什么异同?绘出两序列及其幅频特性曲线。在末尾补零,用点FFT分
3、析这两个信号的幅频特性,观察幅频特性发生了什么变化?两种情况下的FFT频谱还有相同之处吗?这些变化说明了什么? for n=0:3xcn(n+1)=n;end; for n=4:7xcn(n+1)=8-n; xcnxcn = 0 1 2 3 4 3 2 17;stem(n,xcn); xk1=fft(xcn);xdn(n+1)=4-n;xdn(n+1)=n-4; xdnxdn = 4 3 2 1 0 1 2 3stem(n,xdn); xk1=fft(xdn); xcn=xcn,zeros(1,24);31; xdn=xdn,zeros(1,24);时,的幅频特性相同,在点FFT分析这两个信号
4、的幅频特性时,它们还有相同之处,即当取4的整数倍时对应幅值相等。分析:点FFT分析信号的幅频特性:由上两式可知,当k2=4k1时,两个信号的对应频率幅值相等,即对信号末尾补零加长整数个周期可以对原信号达到细化频谱的作用。(4)一个连续时间信号含两个频率分量,经采样得已知分别为1/16和1/64,观察其频谱;当不变,其结果有何不同,为什么? x1n=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1/16)*n); xk1=fft(x1n);subplot(1,2,1);legend(f=1/16 x2n=sin(2*pi*0.125*n)+cos(2*pi*(0.125+1
5、/64)*n); xk2=fft(x2n);xk2=abs(xk2);subplot(1,2,2);stem(n,xk2);f=1/64127; stem(n,xk1); stem(n,xk2);由于离散傅里叶变换的选频性质:不等于整数时,则信号频谱会发生泄漏。(5)用FFT分别计算()和)的16点循环卷积和线性卷积。 xan=exp(-(n-8).2/2); xbn=exp(-0.1*n).*sin(2*pi*0.0625*n); subplot(4,1,1);stem(n,xan);nxa(n) subplot(4,1,2);stem(n,xbn);xb(n) xak=fft(xan);x
6、bk=fft(xbn);x1k=xak.*xbk; x1n=ifft(x1k);subplot(4,1,3);stem(n,x1n);x1(n)循环卷积 x2n=conv(xan,xbn); m=0:length(x2n)-1;subplot(4,1,4);stem(m,x2n);x2(n)线性卷积(6)产生一512点的随机序列,并用做线性卷积,观察卷积前后频谱的变化。要求将分成8段,分别采用重叠相加法和重叠保留法。在编辑调试窗中编写程序:function yy=xeni(N2,xen,i)for n=N2*i:N2*(i+1)-1 xeni(n-N2*i+1)=xen(n+1);endyy=
7、xeni;将上述文件存盘,文件名为xeni.m。function yy=xenni(N1,N2,xen,i)N1+N2*(i+1)-2将上述文件存盘,文件名为xenni.m。function t=shiftmm(a,n)m=length(n);for i=1:a; for j=m+i-1:-1:1 n(j+1)=n(j); end;a n(i)=0;t=n;将上述文件存盘,文件名为shiftmm.m。退回到指令窗: xcn=0 1 2 3 4 3 2 1;xen=rand(1,512); qqqqq=conv(xcn,xen); stem(0:518,qqqqq);幅度 N1=length(x
8、cn);N2=length(xen)/8; xcn=xcn zeros(1,N2-1); xck=fft(xcn); for i=1:8xenii=xeni(N2,xen,i-1);xenii=xenii zeros(1,N1-1);xeki=fft(xenii);yki=xck.*xeki;yni=ifft(yki);y(i,:)=yni; for i=0:for j=0:i*N2-1ynii(i+1,0+1:i*N2-1+1)=0;for j=i*N2:N1+(i+1)*N2-2ynii(i+1,i*N2+1:N1+(i+1)*N2-2+1)=y(i+1,:for j=N1+(i+1)*N
9、2-1:N1+8*N2-2ynii(i+1,N1+(i+1)*N2-1+1:N1+8*N2-2+1)=0; yn=zeros(1,N1+8*N2-1);yn=yn+ynii(i,:N1+8*N2-2; stem(n,yn);重叠相加法 xen21=shiftmm(N1-1,xen);xen2i(i,:)=xenni(N1,N2,xen21,i-1);xek2i=fft(xen2i(i,:);yk2i=xck.*xek2i;yn2i=ifft(yk2i);y2(i,:)=yn2i; y2(:,1:N1-1)=; n2=0:8*N2-1; stem(n2,y2(1,:) y2(2,:) y2(3,
10、:) y2(4,:) y2(5,:) y2(6,:) y2(7,:) y2(8,:);重叠保留法(7)用FFT分别计算)的16点循环相关和线性相关,问一共有多少种结果,它们之间有何异同点。1)求线性相关 k=length(xbn); xan1=xan zeros(1,k-1); xbn1=xbn zeros(1,k-1); xak=fft(xan1); xbk=fft(xbn1); rm=real(ifft(conj(xak).*xbk); rm1=rm(k+1:2*k-1) rm(1:k); m=(-k+1):(k-1); stem(m,rm1);线性相关2)求循环相关 xbk=fft(xbn); stem(n,rm);循环相关(8)用FFT分别计算)的自相关函数。 k=length(xan); xak=fft(xan,2*k); rm=real(ifft(conj(xak).*xak); rm=rm(k+2:2*k) rm(1: stem(m,rm);m(2) xbk=fft(xbn,2*k); rm=real(ifft(conj(xbk).*xbk);