测试信号分析与处理实验.docVIP

第四节 实验 实验一，混叠现象的时域与频域表现 考虑频率分别为3Hz，7Hz，13Hz 的三个余弦信号，即：g1(t)=cos(6πt), g2(t)=cos(14πt), g3(t)=cos(26πt),当采样频率为10Hz 时，即采样间隔为0.1秒，则产生的序列分别为：g1[n]=cos(0.6πn), g2[n]=cos(1.4πn), g3[n]=cos(2.6πn) 对g2[n]，g3[n] 稍加变换可得： g2[n]=cos(1.4πn)＝cos((2π-0.6π)n)= cos(0.6πn) g3[n]=cos(2.6πn)= cos((2π+0.6π)n)=cos(0.6πn) 因此它们在时域表现为一个序列。实际上，当给定频率为(10k±3Hz)的余弦信号，且采样频率为10Hz 时，均表现为g1[n]=cos(0.6πn)的序列。以下的matlab 程序画出三序列的时域和频域图（图9―32），非常直观地说明了混叠现象。 n=1:300; t=(n-1)*1/300; g1=cos(6*pi*t); g2=cos(14*pi*t); g3=cos(26*pi*t); plot(t,g1,t,g2,t,g3); k=1:100; s=k*0.1; q1=cos(6*pi*s); q2=cos(14*pi*s); q3=cos(26*pi*s); hold on; plot(s(1:10),q1(1:10),bd); figure subplot(2,2,1);plot(k/10,abs(fft(q1))) subplot(2,2,2);plot(k/10,abs(fft(q2))) subplot(2,2,3);plot(k/10,abs(fft(q3))) 如果将采样频率改为30Hz，则三信号采样后不会发生频率混叠，可运行以下的程序，观察序列的频谱。 k=1:300; q=cos(6*pi*k/30); q1=cos(14*pi*k/30); q2=cos(26*pi*k/30); subplot(2,2,1);plot(k/10,abs(fft(q))) subplot(2,2,2);plot(k/10,abs(fft(q1))) subplot(2,2,3);plot(k/10,abs(fft(q2))) 问题与实践：保证采样后的信号不发生混叠的条件是什么？若信号的最高频率为17Hz，采样频率为30Hz，问是否会发生频率混叠？混叠成频率为多少Hz的信号？编程验证你的想法。 实验二 泄漏效应 由于泄漏效应的复杂性，下面的实验演示单一频率正弦信号由于截断引起的泄漏。首先考察频率为10Hz的正弦信号，采样频率为64Hz，对32点的采样序列进行DFT(FFT),由于是整周期截取，所以不会产生泄漏，程序和运行结果如下(图9-34)： N=input(type in the length of dft=);%采样点数 t=input(type in the sampling period=);%采样时间间隔 freq=input(type in the sinusoid frequency =);%信号频率 k=0:N-1; f=sin(2*pi*freq*k*(1/t)); F=fft(f); stem(k,abs(F));grid; xlabel(k);ylabel(|x(k)|); 如果采样频率仍为64Hz，采样点数仍为32点，将信号频率改为11Hz，由于非整周期截取，发生了明显的泄漏，如图9-35所示： 问题与实践：如何选取采样点数、采样频率，才能保证频率为11Hz 的正弦信号不会发生泄漏现象？选取采样点数为30，采样频率为66Hz，运行程序看一看是否有泄漏发生。 实验三 FFT(DFT)参数的选取 本实验考察有不同频率的两个正弦信号加上随机信号组成的信号的傅立叶变换，通过改变序列的采样点数、DFT长度以及信号的频率，观察谱分辨率以及泄漏对谱估计的影响。 N=input(signal length=);%采样点数 R=input(DFT length=);%DFT的变换长度 fr=input(type in the sinusoid frequences=);%输入正弦信号的频率 y=hamming(N);%可采用其他的窗函数,进行比较 n=0:N-1; x=cos(2*pi*n*fr(1))+cos(2*pi*n*fr(2))+randn(1,N)/5; x1=x.*y; Fx=fft(x1,R); k=0:R-1; subplot(2,1,1) stem(k/R,abs(Fx));grid 当fr=[0.22 0.34]，N=16时，分别取R=16，32，64，128，256，可以看出随R的增大，