|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
各位老师:
你们好!
我是一个初学EMD的菜鸟,应用xray老师的程序,对一调制信号进行分析,但是还是有些问题不太明白,希望各位老师给予指导!
%%%%%程序
clear;clc;close all;
fs=100;T=10;N=T*fs;
t=0:1/fs:T; % 采样率100HZ
z=(1+sin(2*pi*t)).*sin(2*pi*10*t);%调幅信号
imf = emd(z); %%%%%%%%%%EMD
[A, fa, tt] = hhspectrum(imf);
if size(imf,1) > 1
[A,fa,tt] = hhspectrum(imf(1:end-1, :));
else
[A,fa,tt] = hhspectrum(imf);
end
[E, tt1] = toimage(A,fa,tt,length(tt));
for k = 1:size(E,1)
bjp(k) = sum(E(k,:))*1/fs*1/T;
end
f = (0:N-2)/N*(fs/2);
figure(1)
plot(t,z);
figure(2)
imagesc(tt1,[0,0.5*fs],E);
set(gca,'YDir','normal')
% 使用灰度图显示
% colormap(flipud(gray))
figure(3)
plot(f, bjp);
分析结果如图:
我有几个问题不明白:
1。EMD分解后,各阶的IMF是存在什么位置的?工作空间里的矩阵A,和fa都是9行的矩阵,不知道哪个才是?
2。边际谱中只包含10HZ的频率成分,而我的主要目的是获得调制信号中的包络信号的频率(1HZ),应该怎么获得?
3。应用EMD分解,我是否可以直接从分解中的IMF中提取1HZ的频率信号,再进行FFT变化,从而,获得1HZ的频率?
以上是我的问题,希望各位老师给予指导!!!
|
|