|||
下面求各信号序列之各基本模式分量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)
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-1 07:30
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社