频谱分析实验报告

《频谱分析实验报告》由会员分享，可在线阅读，更多相关《频谱分析实验报告（14页珍藏版）》请在装配图网上搜索。

1、基于的频谱分析一、实验目的1频谱分析的基本概念及应用；2常用信号的频谱：余弦信号、方波信号、三角波信号3.应用DFT进行频谱分析，MATLAB编程仿真(1) 整周期采样、非整周期采样一频谱泄漏；(2) 填零运算一高密度频谱；(3) 非整周期采样时，加窗改进；4. 应用频谱分析仪测量常用信号的频谱。二、实验原理1、利用DFT计算连续信号的频谱分析过程X(jG)X(ej)X(ej)X(e)aNN，2兀k/N从图中看出，这是一次近似过程，首先，用离散采样信号的DTFTX(ej)来近似连续信号x(t)的傅立叶变换X(jG)，其次，将x(n)截短，这一过程aa等效于用一矩形序列RN(n)与x(n)相乘，

2、其DTFT为X(ej)，X(ej)*R(ej)NN最后，再对截短的信号作DFT2、频谱分析的近似过程出现问题及解决办法(1) 混叠混叠现象对连续信号x(t)进行数字处理前，要进行采样x(nT)=x(t)为8(t-nT)an，采样序列的频谱是连续信号频谱的周期延拓，周期为fs，如采样率过低，不满足采样定理，fsv2fh，则导致频谱混迭，使一个周期内的谱对原信号谱产生失真，无法恢复原信号，进一步的数字处理失去依据。另外，干扰信号的频率高于折叠频率f/2的也会产生频谱混叠。s对于干扰信号的频率高于f/2产生频谱混叠的解决办法：s在采样前加一抗混叠滤波器，滤除高于折叠频率f/2的频率成分。s(2) 泄

3、漏 产生原因处理实际信号序列x(n)时，一般总要将它截断为一有限长序列，长为N点，相当于乘以一个矩形窗w(n)=RN(n)。矩形窗函数，其频谱有主瓣，也有许多副瓣，窗口越大，主瓣越窄，当窗口趋于无穷大时，就是一个冲击函数。我们知道，时域的乘积对应频域的卷积，所以，加窗后的频谱实际是原信号频谱与矩形窗函数频谱的卷积，卷积的结果使频谱延伸到了主瓣以外，且一直延伸到无穷。当窗口无穷大时，与冲击函数的卷积才是其本身，这时无畸变，否则就有畸变。例如，信号为ej0nT，是一单线谱，但当加窗后，线谱与抽样函数进行卷积，原来在Q0处的一根谱线变成了以Q0为中心的，形状为抽样函数的谱线序列，原来在一个周期(Qs

4、)内只有一个频率上有非零值，而现在一个周期内几乎所有频率上都有非零值，即x(ejt)的频率成份从Q0处“泄漏”到其它频率处去了。考虑各采样频率周期间频谱“泄漏”后的互相串漏，卷积后还有频谱混迭现象产生。 解决办法：增加截短的长度N，以减小泄漏。(3)栅栏效应栅栏效应的产生N点DFT是在频率区间0，2n上对信号频谱进行N点等间隔采样，得到的是若干个离散的频谱点X(k)，且它们限制在基频的整数倍上，这就好像在栅栏的一边通过缝隙看另一边的景象一样，只能在离散点处看到真实的景象，其余部分频谱成分被遮挡，所以称之为栅栏效应。减小栅栏效应方法：尾部补零，使谱线变密，增加频域采样点数，原来漏掉的某些频谱分量

5、就可能被检测出来。(4)DFT的分辨率填补零值可以改变对DTFT的采样密度，人们常常有一种误解，认为补零可以提高DFT的频率分辨率。事实上我们通常规定DFT的频率分辨率为f/n，s这里的N是指信号x(n)的有效长度，而不是补零的长度。不同长度的x(n)其DTFT的结果是不同的；而相同长度的x(n)尽管补零的长度不同其DTFT的结果应是相同的，他们的DFT只是反映了对相同的DTFT采用了不同的采样密度。参数选择的一般原则：(1) 若已知信号的最高频率f，为防止混叠，选定采样频率f2f；maxsmax(2) 根据频率分辩率Af，确定所需DFT的长度n，f/Af；s(3) f和N确定以后，即可确定相

6、应模拟信号的时间长度工,n/f，NT(这ss里T是采样周期)。三、实验内容及结果1、余弦波频谱分析对余弦函数x(t)=cos(2*25t)进行频谱分析，采样频率为1000Hz,对其进行整周期采样，非整周期采样，频谱泄露及其改善方法。其程序代码和运行结果如下。程序代码：整周期采样整周期采样点，占八、加长矩形窗长度矩形窗加长到改窗运行结果：on;拒形窗仙长到1西2、方波信号频谱对方波函数X(t)=square(2*25t)进行频谱分析，采样频率为1000Hz，对其进行整周期采样，非整周期采样，频谱泄露及其改善方法。其程序代码和运行结果如下。程序代码：fs=1000;t=0:1/fs:1.5;%xx

7、=square(2*pi*25*t);subplot(321),plot(t,x),axis(00.2-1.21.2);title(x(t)=square(2*pi*25*t);gridon;%x(n)nx=0:99;x0=square(2*pi*1/25*nx);subplot(322),stem(nx,x0);%整周期采样N1=100;x1=x0(1:N1);X1=fft(x1,N1);magX1=abs(X1(1:N1/2+1)*2/N1;w1=0:N1/2*2*pi/N1;subplot(323);stem(w1/pi,magX1);title(整周期采样)；grid%265点DFT,

8、highresolutionspectrumN2=256;X=fft(x0,N2);magX=abs(X(1:N2/2+1)*2/N2;w=0:N2/2*2*pi/N2;subplot(324);stem(w/pi,magX);title(256点DFT,highresolutionspectrum);gridon;%加长矩形窗长度mx=0:149;xm=square(2*pi*1/25*mx);Nm=256;Y=fft(xm,Nm);magY=abs(Y(1:Nm/2+1)*2/Nm;wm=0:Nm/2*2*pi/Nm;subplot(325);stem(wm/pi,magY);title(

9、矩形窗加长到);gridon;%改窗N=256;wn=hamming(length(x0);hn=x0.*wn;H=fft(hn,N);magH=abs(H(1:N/2+1)*2/N;on;wH=0:N/2*2*pi/N;subplot(326);stem(wH/pi,magH);title(力口窗)；grid运行结果：x(t)=square(2*pi*25*t)2256点DFT,highresolutionspectrum1-：00.5越hmrnming窗00.51矩環窗加怅到1药0.400.53、锯齿波信号频谱对方波函数x(t)=sawtooth(2*25t)进行频谱分析，采样频率为100

10、0Hz，对其进行整周期采样，非整周期采样，频谱泄露及其改善方法。其程序代码和运行结果如下。程序代码：fs=1000;t=0:1/fs:1.5;%xx=sawtooth(2*pi*25*t);subplot(321),plot(t,x),axis(00.2-1.21.2);title(x(t)=sawtooth(2*pi*25*t);gridon;%x(n)nx=0:99;x0=sawtooth(2*pi*1/25*nx);subplot(322),stem(nx,x0);title(x(n);gridon;%整周期采样on;on;N1=100;x1=x0(1:N1);X1=fft(x1,N1)

11、;magX1=abs(X1(1:N1/2+1)*2/N1;w1=0:N1/2*2*pi/N1;subplot(323);stem(w1/pi,magX1);title(整周期采样)；grid%265点DFT,highresolutionspectrumN2=256;X=fft(x0,N2);magX=abs(X(1:N2/2+1)*2/N2;w=0:N2/2*2*pi/N2;subplot(324);stem(w/pi,magX);title(256点DFT,highresolutionspectrum);gridon;%加长矩形窗长度mx=0:149;xm=sawtooth(2*pi*1/2

12、5*mx);Nm=256;Y=fft(xm,Nm);magY=abs(Y(1:Nm/2+1)*2/Nm;wm=0:Nm/2*2*pi/Nm;subplot(325);stem(wm/pi,magY);title(矩形窗加长到);gridon;N=256;%改窗wn=hamming(length(x0);hn=x0.*wn;H=fft(hn,N);magH=abs(H(1:N/2+1)*2/N;wH=0:N/2*2*pi/N;subplot(326);stem(wH/pi,magH);title(力口窗)；grid运行结果：title(x(n);gridon;%整周期采样title(x(n);g

13、ridon;%整周期采样x(t)=sawtooth(2ffpi*25*t)4、三角波信号频谱对方波函数x(t)=sawtooth(2*25t,0.5)进行频谱分析，采样频率为1000Hz,对其进行整周期采样，非整周期采样，频谱泄露及其改善方法。其程序代码和运行结果如下。程序代码：fs=1000;t=0:1/fs:1.5;%xx=sawtooth(2*pi*25*t,0.5);subplot(321),plot(t,x),axis(00.2-1.21.2);title(x(t)=sawtooth(2*pi*25*t,0.5);gridon;%x(n)nx=0:99;x0=sawtooth(2*p

14、i*1/25*nx,0.5);subplot(322),stem(nx,x0);N1=100;x1=x0(1:N1);X1=fft(x1,N1);magX1=abs(X1(1:N1/2+1)*2/N1;w1=0:N1/2*2*pi/N1;subplot(323);stem(w1/pi,magX1);title(整周期采样)；grid%265点DFT,highresolutionspectrumN2=256;X=fft(x0,N2);magX=abs(X(1:N2/2+1)*2/N2;w=0:N2/2*2*pi/N2;subplot(324);stem(w/pi,magX);title(256点

15、DFT,highresolutionspectrum);gridon;%加长矩形窗长度mx=0:149;xm=sawtooth(2*pi*1/25*mx,0.5);Nm=256;Y=fft(xm,Nm);magY=abs(Y(1:Nm/2+1)*2/Nm;wm=0:Nm/2*2*pi/Nm;subplot(325);stem(wm/pi,magY);title(矩形窗加长到);gridon;%改窗N=256;wn=hamming(length(x0);hn=x0.*wn;H=fft(hn,N);magH=abs(H(1:N/2+1)*2/N;wH=0:N/2*2*pi/N;subplot(32

16、6);stem(wH/pi,magH);title(力口窗)；gridon;on;运行结果：title(x(n);gridon;%整周期采样5、频率分辨率频率分辨率是指所用算法能把信号中两个靠得很近的谱峰保持分开的能力。若两个信号的谱峰靠的很近，很难分开。为了解决这个问题，可以进行填零运算，但是填零运算产生的是高密度频谱，不能很好的解决。另一个方法是加长原始信号长度，这样便是产生高分辨率频谱，便可以分开很靠近的两个谱峰。对函数nx=0:99;x=cos(0.48*pi*nx)+cos(0.52*pi*nx)进行频谱分析，对其进行DTFT,10点DFT,100点DFT，填零运算，高密度频谱和高分

17、辨率频谱。其程序代码和运行结果如下。程序代码：clear;clc;%x(n)nx=0:99;x=cos(0.48*pi*nx)+cos(0.52*pi*nx);subplot(321);stem(nx,x);title(x(n)=cos(0.48*pi*n)+cos(0.52*pi*n);gridon;%X(ejw)w=0:200*pi/200;X=dtft(x,nx,w);magX=abs(X);subplot(322);plot(w/pi,magX);title(X(ejw);gridon;%10paDFTN1=10;x1=x(1:N1);X1=fft(x1,N1);magX1=abs(X

18、1(1:N1/2+1)*2/N1;w1=0:N1/2*2*pi/N1;subplot(323);stem(w1/pi,magXl);title(10paDFT);gridon;%100paDFT-highresolutionspectrumN2=100;X=dft(x,N2);magX=abs(X(1:N2/2+1)*2/N2;w=0:N2/2*2*pi/N2;subplot(324);stem(w/pi,magX);title(100paDFT-highresolutionspectrum);gridon;%ii90$00,100yaDFT,high-densityspectrumx2=x(1:N1)zeros(1,90);N2=100;n2=0:N2-1;subplot(325);stem(n2,x2);title(x(n),0=n=9+90zeros);gridon;X2=fft(x2,N2);magX2=abs(X2(1:N2/2+1);w2=0:N2/2*2*pi/N2;subplot(326);stem(w2/pi,magX2);title(highdensityspectrum);gridon;运行结果：