信号的时频域分析实验MATLAB

信号的时频域分析信号的时频域分析信号的时频域分析信号的时频域分析

clear all;
bt=0;                           %开始时间
dt=0.001;                       %时间间隔
N=512;                          %傅立叶变换点数
et=bt+N*dt-dt;                  %结束时间
t=bt:dt:et;                     %时间域
TT=et-bt;                        %总的时间
y = cos（2*pi*80*t）;             %待分析的信号
figure（1）;
subplot（321）
plot（ty）;
title（‘信号‘）;
xlabel（‘时间（s）‘）;    
axis（[0 0.5 -1.2 1.2]）; 

Y =fftshift（fft（yN））;
df=1/TT;                         %频率间隔
Tf=N*df;                         %分析的频宽                         
f=-Tf/2+df:df:Tf/2;             %频率域
subplot（322）
plot（fabs（Y））;
title（‘频谱‘）;
xlabel（‘频率（Hz）‘）;
axis（[-100 100 0 max（abs（Y））]）;

Pyy=abs（Y/Tf）.^2/TT;   %计算功率谱密度
subplot（323）
plot（fPyy）;
title（‘根据定义求出的功率谱‘）;
xlabel（‘频率（Hz）‘）;
axis（[-100 100 0 max（abs（Pyy））]）;

M=N;
Ryy=zeros（1M）;
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）;
end
subplot（324）
plot（tRyy（1:M））;
title（‘信号自相关‘）;
xlabel（‘时间（s）‘）;      
axis（[0 0.5 -1 1]）;

Pyy2=fftshift（fft（RyyM））/Tf;
subplot（325）
plot（fabs（Pyy2））;
title（‘自编自相关函数求出的功率谱‘）;
xlabel（‘频率（Hz）‘）;
axis（[-100 100 0 max（abs（Pyy2））]）;

Ryy2=xcorr（yy‘unbiased‘）;
L=length（Ryy2）;
Pyy3=fftshift（fft（Ryy2（L-N+1:L）N））/Tf;
subplot（326）
plot（fabs（Pyy3））;
title（‘系统自相关函数求出的功率谱‘）;
xlabel（‘频率（Hz）‘）;
axis（[-100 100 0 max（abs（Pyy3））]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1498  2009-11-03 09:28  实验一 信号的时频分析实验\t1.m

     文件       1494  2009-11-03 09:28  实验一 信号的时频分析实验\t2.m

     文件       1495  2009-11-03 09:28  实验一 信号的时频分析实验\t3.m

     文件       1495  2009-11-03 09:28  实验一 信号的时频分析实验\t4.m

     文件       1504  2009-11-03 09:08  实验一 信号的时频分析实验\t5.m

     文件       1507  2009-11-03 09:21  实验一 信号的时频分析实验\t6.m

     文件       1508  2009-11-03 09:21  实验一 信号的时频分析实验\t7.m

     文件       1509  2009-11-03 09:21  实验一 信号的时频分析实验\t8.m

     文件     732160  2009-11-03 09:35  实验一 信号的时频分析实验\实验一 信号的时频分析实验.doc

     目录          0  2010-03-02 21:51  实验一 信号的时频分析实验

----------- ---------  ---------- -----  ----

               744170                    10