abc198478 发表于 2009-2-21 11:47

求助:使用G. rilling编的EMD程序处理信号的具体例子

各位大虾:
      小弟近来要用G.rilling的程序来处理信号,但是程序里面有好多细节问题一时难以看懂,所以求助已经会使用该程序的高手分享一下所用的例子,帮助我更好地理解使用G.rilling写的程序!先谢了!

zhangnan3509 发表于 2009-2-21 13:26

本身那个程序包里就有例子啊！很清楚的，楼主好好看看啊

abc198478 发表于 2009-2-21 17:13

我看了,但是copy上去运行不出来,出现错误很多!

zhangnan3509 发表于 2009-2-21 21:27

是什么问题啊具体说说看

大鹏之举 发表于 2009-3-3 13:19

我把程序打开运行一下也有错误呀？请高手解释一下！

我把程序打开运行一下也有错误呀？请高手解释一下！附件是我下载的的M文件和出现的错误。

anfengping 发表于 2010-1-25 16:00

同问

我也是这个问题，麻烦各位大牛给解释一下。非常感谢！！！

zhjinctgu 发表于 2011-12-23 21:56

% 计算2FSK信号的HHT谱和边际谱
% 作者：xray    2007.11
clear

rand('seed', 0);

T = 0.05;                % 仿真时间
R = 500;               % 码速率
fd = 1000;               % 载波频差
fc = 2000;               % 载波频率
fs = 20000;            % 采样率
samp = fs/R;             % 每个码元上的采样点数
N = T*fs;
n = 1:N;
x = randint(1, R*T, 2);
y = fskmod(x, 2, fd, samp, fs);
y = y .* exp(i*2*pi*fc/fs*n);
y = real(y);
% z = awgn(y, 20, 'measured');
z = y;

imf = emd(z);
= hhspectrum(imf);
if size(imf,1) > 1
    = hhspectrum(imf(1:end-1, :));
else
    = hhspectrum(imf);
end
= 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-3)/N*(fs/2);

figure(1)
plot(z);

figure(2)
imagesc(tt1,,E);
set(gca,'YDir','normal')
% 使用灰度图显示
% colormap(flipud(gray))

figure(3)
plot(f, bjp);
页: [1]
查看完整版本: 求助:使用G. rilling编的EMD程序处理信号的具体例子