下面求各信号序列之各基本模式分量IMF之瞬时振幅与瞬时频率。

 

 （1）对TW486012各IMF分量求取瞬时振幅与瞬时频率。

 

[A_imf_TW486012,f_imf_TW486012]=hhspectrum(imf_TW486012);

%A_imf_TW486012是TW486012各个IMF分量的瞬时振幅组成的矩阵,f_imf_TW486012是TW486012各个IMF分量对应的瞬时频率向组成的矩阵。

 

    运行，A_imf_TW486012、f_imf_TW486012分别得到一个15*58318的矩阵。

前14行分别是IMF分量瞬时振幅与瞬时频率，最后一行分别是对应残余项的瞬时振幅与瞬时频率。58318是时间长度，比原始信号序列少两个采样点。

 

%A_imf_TW486012各IMF分量瞬时振幅图

for k=1:14
subplot(5,3,k)
plot(A_imf_TW486012(k,:))
end

 

    运行，得：

    图14-1 A_imf_TW486012

  

%f_imf_TW486012各IMF分量瞬时频率图
for k=1:14
subplot(5,3,k)
plot(f_imf_TW486012(k,:))
end

 

    运行，得：

    图14-2 f_imf_TW486012

  

（2）对GY291712各IMF分量求取瞬时振幅与瞬时频率。

 

[A_imf_GY291712,f_imf_GY291712]=hhspectrum(imf_GY291712);

%A_imf_GY291712是GY291712各个IMF分量的瞬时振幅组成的矩阵,f_imf_GY291712是GY291712各个IMF分量对应的瞬时频率组成的矩阵。

 

    运行，A_imf_GY291712、f_imf_GY291712分别得到一个12*35002的矩阵。

前11行分别是IMF分量瞬时振幅与瞬时频率，最后一行分别是对应残余项的瞬时振幅与瞬时频率。35002是时间长度，比原始信号序列少两个采样点。

 

%A_imf_GY291712各IMF分量瞬时振幅图

for k=1:11

subplot(4,3,k)
plot(A_imf_GY291712(k,:))
end

 

    运行，得：

    图14-3 A_imf_GY291712 

 

%f_imf_GY291712各IMF分量瞬时频率图
for k=1:11

subplot(4,3,k)
plot(f_imf_GY291712(k,:))
end

 

    运行，得：

    图14-4 f_imf_GY291712

 

（3）对DY291712各IMF分量求取瞬时振幅与瞬时频率。

 

[A_imf_DY291712,f_imf_DY291712]=hhspectrum(imf_DY291712);

%A_imf_DY291712是DY291712各个IMF分量的瞬时振幅组成的矩阵,f_imf_DY291712是DY291712各个IMF分量对应的瞬时频率组成的矩阵。

 

    运行，A_imf_DY291712、f_imf_DY291712分别得到一个13*35002的矩阵。

前12行分别是IMF分量瞬时振幅与瞬时频率，最后一行分别是对应残余项的瞬时振幅与瞬时频率。35002是时间长度，比原始信号序列少两个采样点。

 

%%A_imf_DY291712各IMF分量瞬时振幅图

for k=1:12
subplot(4,3,k)
plot(A_imf_DY291712(k,:))
end

 

    运行，得：

    图14-5 A_imf_DY291712

 

%%f_imf_DY291712各IMF分量瞬时频率图
for k=1:12

subplot(4,3,k)
plot(f_imf_DY291712(k,:))
end

 

    运行，得：

    图14-6 f_imf_DY291712

  

（4）对JY291712各IMF分量求取瞬时振幅与瞬时频率。

 

[A_imf_JY291712,f_imf_JY291712]=hhspectrum(imf_JY291712);

%A_imf_JY291712是JY291712各个IMF分量的瞬时振幅组成的矩阵,f_imf_JY291712是JY291712各个IMF分量对应的瞬时频率组成的矩阵。

 

    运行，A_imf_JY291712、f_imf_JY291712分别得到一个14*35002的矩阵。

前13行分别是IMF分量瞬时振幅与瞬时频率，最后一行分别是对应残余项的瞬时振幅与瞬时频率。35002是时间长度，比原始信号序列少两个采样点。

 

%%A_imf_JY291712各IMF分量瞬时振幅图

for k=1:13

subplot(5,3,k)
plot(A_imf_JY291712(k,:))
end

 

    运行，得：

    图14-7 A_imf_JY291712

 

%%f_imf_JY291712各IMF分量瞬时频率图
for k=1:13

subplot(5,3,k)
plot(f_imf_JY291712(k,:))
end

 

    运行，得：

    图14-8 f_imf_JY291712

 

（5）对CY291712各IMF分量求取瞬时振幅与瞬时频率。

 

[A_imf_CY291712,f_imf_CY291712]=hhspectrum(imf_CY291712);

%A_imf_CY291712是CY291712各个IMF分量的瞬时振幅组成的矩阵,f_imf_CY291712是CY291712各个IMF分量对应的瞬时频率组成的矩阵。

  

    运行，A_imf_CY291712、f_imf_CY291712分别得到一个14*35002的矩阵。

前13行分别是IMF分量瞬时振幅与瞬时频率，最后一行分别是对应残余项的瞬时振幅与瞬时频率。35002是时间长度，比原始信号序列少两个采样点。

 

%%A_imf_CY291712各IMF分量瞬时振幅图

for k=1:13

subplot(5,3,k)
plot(A_imf_CY291712(k,:))
end

 

    运行，得：

    图14-9 A_imf_CY291712

%%f_imf_CY291712各IMF分量瞬时频率图
for k=1:13

subplot(5,3,k)
plot(f_imf_CY291712(k,:))
end

 

    运行，得：

    图14-10 f_imf_CY291712

 

（6）对MB291712各IMF分量求取瞬时振幅与瞬时频率。

 

[A_imf_MB291712,f_imf_MB291712]=hhspectrum(imf_MB291712);

%A_imf_MB291712是MB291712各个IMF分量的瞬时振幅组成的矩阵,f_imf_MB291712是MB291712各个IMF分量对应的瞬时频率组成的矩阵。

 

    运行，A_imf_MB291712、f_imf_MB291712分别得到一个14*35002的矩阵。

前13行分别是IMF分量瞬时振幅与瞬时频率，最后一行分别是对应残余项的瞬时振幅与瞬时频率。35002是时间长度，比原始信号序列少两个采样点。

 

%%A_imf_MB291712各IMF分量瞬时振幅图

for k=1:13

subplot(5,3,k)
plot(A_imf_MB291712(k,:))
end

 

    运行，得：

    图14-11 A_imf_MB291712

  

%%f_imf_MB291712各IMF分量瞬时频率图
for k=1:13

subplot(5,3,k)
plot(f_imf_MB291712(k,:))
end

 

    运行，得：

    图14-12 f_imf_MB291712

  

    以上各信号序列各IMF分量都有瞬时振幅与瞬时频率两列向量， 可以用多种方法画出它们的时间－瞬时频率－振幅三维图，即所谓的时频图来。但数据太多，画在一张图上将会是乱糟糟、黑乎乎的一团，什么东西也看不出来。下面随意选取MB291712的第5个IMF分量的一小段数据，画一张三维时频图。

 

t30001_30500=[30001:30500];

stem3(t30001_30500,f_imf_MB291712(5,30001:30500),A_imf_MB291712(5,30001:30500))

 

    运行，得：

    图14-13 MB291712的第5个IMF分量的第30001～30500组数据时频图

 

    将上图旋转一下，得：

    图14-14 图14-13之旋转图

 

    图14-15 图14-13之旋转图

 

 

（本文首发于：http://blog.sina.com.cn/s/blog_6ad0d3de0100zggw.html

首发时间：2011-09-14 16:07:03）