matlab-信号的时频分析.doc

《matlab-信号的时频分析.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、信号的时频分析实验一、实验目的1、掌握Matlab对信号时频分析方法。2、掌握能量信号、周期性功率信号和非周期性功率信号的概念。3、掌握能量信号和功率信号的截断信号的时频域特性。4、掌握相关函数的概念及与功率谱的关系。二、实验原理1、能量信号的时频分析（1）能量信号：能量有限的信号，满足。如时间受限信号。（2）能量信号的频谱密度能量的频谱密度为：的逆变换为原信号：也可以表示为：，（3）能量谱密度根据帕塞瓦定理（能量守恒），可以知道因此，可以将看成是信号的能量谱密度，表示能量随频率的分布。（4）能量信号的相关函数能量信号的自相关函数定义：能量信号的互相关函数定义：（5）能量

2、信号的相关函数和能量谱密度的关系由此可看到，能量信号的自相关与其能谱密度是一对傅立叶变换对。2、功率信号的时频分析（1）功率信号：如果信号的能量无穷大，但其功率存在，则称该信号为功率信号，（2）功率信号的频谱函数设为周期性功率信号，为周期，则有式中，为复数，表示为式中，是频率为的分量的振幅；是频率为的分量的相位。由傅立叶级数知识可知，上式为周期信号s(t)的傅里叶级数的系数，即周期信号s(t)可以表示为：（3）功率谱密度令的截短信号为，则为能量信号。有定义功率谱密度为：得到信号功率：若功率信号为周期信号，则：，其中，由于，则（4）功率信号的相关函数功率信号的自相关函数定义

3、：功率信号的互相关函数定义：（5）功率信号的相关函数和功率谱密度的关系相关函数和功率谱密度互为傅立叶变换。3、信号的频域分析方法（1）连续非周期信号的傅里叶变换（FT）对于连续非周期信号（能量信号），即满足绝对可积条件时。可以采用傅里叶变换来分析信号的频域特性。频域为连续非周期信号。（2）连续周期信号的傅立叶级数（FS）对于连续周期信号，如果满足Dirichlet条件，可以采用傅立叶级数来分析信号的频域特性。当引入冲激函数时，可以用傅里叶变换来分析连续周期信号的频域。频域为离散非周期信号。（3）离散非周期信号的傅里叶变换（DTFT）频域为连续周期信号。（4）离散周期信号的

4、傅立叶级数（DFS）频域为离散周期信号。时域的离散对应于频域的周期；时域的连续对应于频域的非周期；时域的周期对应于频域的离散；时域的非周期对应于频域的连续。（5）离散傅里叶变换（DFT）将离散周期信号的傅立叶级数各取一个周期（N个采样点），可以得到如下离散傅里叶变换。即将离散傅里叶变换作周期扩展，可以得到傅立叶级数。（6）快速傅里叶变换（FFT）对于离散傅里叶变换的规律，可以进行快速傅里叶变换FFT。它包含两类：一类为Cooley和Tukey提出的N等于2的整数次幂的算法；另一类为Winograd提出的N不等于2的整数次幂的算法。这两类方法都大大降低了DFT的乘法运算的次

5、数。4、连续信号的离散化为了便于对信号通过计算机进行处理，如FFT。需要对信号进行采样。根据采样定理，当带限信号带宽为W时，采样频率，即采样间隔时。采样信号表示为：，其傅立叶变换为：，为的傅立叶变换因此通过一个带宽为W，增益为的低通滤波器可以将原信号恢复。重建公式为：。由于离散序列的离散傅立叶变换(DFT)表示为则。Matlab中只是对离散序列信号进行了DFT，因此，若要得到该离散信号的原连续信号的频谱，还需要将得到的信号乘以或除以，才能得到原连续信号的频谱。5、matlab对信号时频分析的方法（1）相关：Matlab中通过xcorr函数计算互相关。自相关：调用格式为c=

6、xcorr(x)或c=xcorr(x,maxlags)，其中x为输入序列，c为它的自相关函数估计。maxlags为x的最大时延，若缺省，函数返回值c长度是2N-1；若不缺省，返回值c长度是2*maxlags+1。互相关：调用格式有四种c=xcorr(x,y)；c=xcorr(x,y,’option’)；c=xcorr(x,y,maxlags,’option’)；[c,lags]=xcorr(x,y,maxlags,’option’)。其中x,y为两个独立的信号序列，长度相同。maxlags意义同上。option为选择项：’bised’计算有偏互相关函数估计，'unbias

7、ed'计算无偏互相关函数估计，‘coeff’序列归一化使零延迟的自相关函数为1，‘none’为缺省情况。（2）傅立叶变换：Matlab中通过fft函数来计算离散傅立叶变换。（3）直接法求功率谱密度由公式可直接求得功率谱密度。（4）利用自相关函数求功率谱密度也可以通过自相关函数和功率谱密度之间的关系来得到信号的功率谱密度图。M=512;自相关函数可由下面语句求得：M=N;Ryy=zeros(1,M);form=1:Mforn=1:N-m+1Ryy(m)=Ryy(m)+y(n)*y(n+m-1);endRyy(m)=Ryy(m)/