1、3、信号的频域分析方法(1)连续非周期信号的傅里叶变换(FT)对于连续非周期信号(能量信号),即满足绝对可积条件时。可以采用傅里叶变换来分析信号的频域特性。频域为连续非周期信号。(2)连续周期信号的傅立叶级数(FS)对于连续周期信号,如果满足Dirichlet条件,可以采用傅立叶级数来分析信号的频域特性。当引入冲激函数时,可以用傅里叶变换来分析连续周期信号的频域。频域为离散非周期信号。(3)离散非周期信号的傅里叶变换(DTFT)频域为连续周期信号。(4)离散周期信号的傅立叶级数(DFS)频域为离散周期信号。时域的离散对应于频域的周期;时域的连续对应于频域的非周期;时域的周期对应于频域的离散;时
2、域的非周期对应于频域的连续。(5)离散傅里叶变换(DFT)将离散周期信号的傅立叶级数各取一个周期(N个采样点),可以得到如下离散傅里叶变换。即将离散傅里叶变换作周期扩展,可以得到傅立叶级数。(6)快速傅里叶变换(FFT)对于离散傅里叶变换的规律,可以进行快速傅里叶变换FFT。它包含两类:一类为Cooley和Tukey提出的N等于2的整数次幂的算法;另一类为Winograd提出的N不等于2的整数次幂的算法。这两类方法都大大降低了DFT的乘法运算的次数。4、连续信号的离散化为了便于对信号通过计算机进行处理,如FFT。需要对信号进行采样。根据采样定理,当带限信号带宽为W时,采样频率,即采样间隔时。采
3、样信号表示为:,其傅立叶变换为:,为的傅立叶变换因此通过一个带宽为W,增益为的低通滤波器可以将原信号恢复。重建公式为:。由于离散序列的离散傅立叶变换(DFT)表示为则。Matlab中只是对离散序列信号进行了DFT,因此,若要得到该离散信号的原连续信号的频谱,还需要将得到的信号乘以或除以,才能得到原连续信号的频谱。5、matlab对信号时频分析的方法(1)相关:Matlab中通过xcorr函数计算互相关。自相关:调用格式为c=xcorr(x)或c=xcorr(x,maxlags),其中x为输入序列,c为它的自相关函数估计。maxlags为x的最大时延,若缺省,函数返回值c长度是2N-1;若不缺省
4、,返回值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计算有偏互相关函数估计,unbiased计算无偏互相关函数估计, coeff序列归一化使零延迟的自相关函数为1,none为缺省情况。(2)傅立叶变换:Matlab中通过fft函数来计算离散傅立叶变换。(3)直接法求功率谱密度由公式可直接求得功率谱密度。
5、(4)利用自相关函数求功率谱密度也可以通过自相关函数和功率谱密度之间的关系来得到信号的功率谱密度图。M=512; 自相关函数可由下面语句求得:M=N;Ryy=zeros(1,M);for m=1:M for n=1:N-m+1 Ryy(m)=Ryy(m)+y(n)*y(n+m-1); end Ryy(m)=Ryy(m)/(N-m+1);endsubplot(3,2,4)plot(t,Ryy(1:M);title(信号自相关);xlabel(时间(s)axis(0 1 -1.2 1.2);三、实验步骤1、输入余弦信号,采样间隔取dt=0.001s,采样点数取N=512。将其变换到频域观察其频谱,
6、再将其反变换到时域与原信号进行比较,并绘制该信号的时域和频域波形。2、通过直接求解法求该余弦信号的功率谱密度,并绘制该信号的功率谱密度波形。3、分别用xcorr函数和自编制自相关函数的方式,求该余弦信号的自相关函数,并绘制自相关函数波形。4、利用自相关函数和功率谱密度的关系且该余弦信号的功率谱密度,并绘制该信号的功率谱密度波形。5、改变抽样间隔dt=0.002s,重复步骤1、2、3、4。6、改变采样点数N=1024,重复1、2、3、4、5。7、输入信号,重复1、2、3、4、5、6。四、实验结果1、写出完成实验步骤的程序。实验程序bt=0; %开始时间dt=0.001; %时间间隔N=1024;
7、 %傅立叶变换点数et=bt+N*dt-dt; %结束时间t=bt:dt:et; %时间域TT=et-bt; %总的时间y = cos(2*pi*80*t); %待分析的信号第二次为cos(2*pi*80*t)+2*cos(2*pi*50*t)figure(1);subplot(3,2,1)plot(t,y);信号Y =fftshift(fft(y,N);df=1/TT; %频率间隔Tf=N*df; %分析的频宽 f=-Tf/2+df:df:Tf/2; %频率域subplot(3,2,2)plot(f,abs(Y);频谱频率(Hz)axis(-100 100 0 max(abs(Y);Pyy=
8、abs(Y/Tf).2/TT; %计算功率谱密度subplot(3,2,3)plot(f,Pyy);根据定义求出的功率谱axis(-100 100 0 max(abs(Pyy);Pyy2=fftshift(fft(Ryy,M)/Tf;subplot(3,2,5)plot(f,abs(Pyy2);自编自相关函数求出的功率谱axis(-100 100 0 max(abs(Pyy2);Ryy2=xcorr(y,y,L=length(Ryy2);Pyy3=fftshift(fft(Ryy2(L-N+1:L),N)/Tf;subplot(3,2,6)plot(f,abs(Pyy3);系统自相关函数求出的功率谱axis(-100 100 0 max(abs(Pyy3);2、绘制实验步骤中要求的图形。1,n= 512点dt=0.0012,n=512 dt=0.0023,n=1024 dt=0.0014,n= 1024 dt=0.0025,n=512 dt=0.0016n=,512,dt=0,0027,n=1024 dt=0.0018,n=1024 dt=0.002分析 :由以上实验可以验证采样点数越多,采样间隔越小 越能放映真是的波形包络和频谱特征。