信号时频域分析方法。实验.docVIP

通过本实验掌握基本信号的时频域分析方法。 实验仪器（软、硬件）： 1、 计算机 1台 2、 Matlab软件 1套 3、 激光打印机 1台 实验步骤 1、 在Matlab中产生不同的信号，其中主要包括正弦信号、方波、冲激信号、随机噪声、矩形窗函数、三角波等； 2、 对产生的信号进行Fourier级数展开、Fourier变换； 3、 产生一个由正弦信号和随机信号叠加的混合信号，并对其进行进行FFT计算； 4、 应用不同窗函数对一正弦信号进行采样，其中包括矩形窗、Hamming窗、Hanning窗。比较不同窗函数采样得到的结果。 实验结果 一、 单个信号 1. 正弦信号及FFT： y=2*sin(2*pi*80*t) 2. 随机信号及FFT： y=randn(size(t)); 3. 方波信号及FFT： y=square(2*pi*10*t); 4. 锯齿波信号及FFT： y=sawtooth(2*pi*10*t); 二、 复合信号 1. 三个正弦信号叠加及FFT： y=sin(2*pi*30*t)+2*sin(2*pi*40*t)+3*sin(2*pi*50*t); 分析：由频谱可以看出各个分量的频率及幅度，幅值大的在时域上占的比重大 2. 正弦信号叠加随机信号及FFT： y=2*sin(2*pi*50*t)+randn(size(t)); 分析：随机信号每个频率都有但是比重都不大，时域上正弦信号的趋势还在，频域上除50hz的以外其他幅度很小。 3. 正弦信号叠加方波信号接FFT： y=2*sin(2*pi*50*t)+square(2*pi*10*t); 分析：时域上正弦信号随方波的叠加跳跃性波动，频域上可以看出频率50的正弦波幅值为2，占主导。 三、 信号加窗 1. 正弦信号加矩形窗： 程序： fs=1000; det=1/fs; t=0:det:1.5; L=length(t); w=zeros(1,length(t)); window_width = 1*fs; w(1:window_width)=rectwin(window_width); y=sin(2*pi*50*t)+sin(2*pi*51*t); y=y.*w; subplot(2,1,1); plot(t,y); subplot(2,1,2); NFFT = 2^nextpow2(L); Y = fft(y,NFFT)/L; f = fs/2*linspace(0,1,NFFT/2+1); plot(f,2*abs(Y(1:NFFT/2+1))) 结果图： 2. 正弦信号加hann窗： 程序： 将w(1:window_width)=rectwin(window_width); 改为w(1:window_width)=hann(window_width); 结果图： 3. 正弦信号加hamming窗： 程序： 将w(1:window_width)=hann(window_width); 改为w(1:window_width)=hamming(window_width); 结果图： 分析：加窗时如果窗的长度大于信号长度则信号后面补0，若窗长度小于信号长度则后面的数据丢失。矩形窗的优点是主瓣比较集中，缺点是旁瓣较高；汉宁窗主瓣加宽并降低，旁瓣则显著减小，但汉宁窗主瓣加宽，相当于分析带宽加宽，频率分辨力下降；而海明窗是改进了的汉宁窗，主瓣比较集中，旁瓣则显著减小。 实验一，混叠现象的时域与频域表现 考虑频率分别为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