求助，hht画不出imf图形！

gunwong

我是菜鸟，在论坛找了几个程序，在matlab7.0下运行，但是画不出imf图像。
高手帮我看看怎么回事？万分感谢！！！

% 计算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);
[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-3)/N*(fs/2);

figure(1)
plot(z);

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

figure(3)
plot(f, bjp);

运行提示错误：
??? Function 'real' is not defined for values of class 'cell'.

Error in ==> real at 15
  [varargout{1:nargout}] = builtin('real', varargin{:});

Error in ==> hilbert at 40
  xr = real(xr);

Error in ==> hhspectrum at 78
  an(i,:)=hilbert(x(i,:)')';

Error in ==> hht1 at 23
[A, fa, tt] = hhspectrum(imf);

gunwong

??? Function 'real' is not defined for values of class 'cell'.
是不是要找一个转换程序？？
急切等大牛帮忙解决，万分感谢！！！

dzkt

程序一点问题也没有

shinhwaforever

我试试了一下,程序也没问题,是不是你的工具包里的文件有问题呢

youshen

??? Undefined command/function 'instfreq'.

Error in ==> hhspectrum at 79
  f(i,:)=instfreq(an(i,:)',tt,l)';
这个是怎么回事？。。

ChaChing

回复 youshen 的帖子

LS没有这个instfreq函数! 以下式检查下
which instfreq -all

youshen

程序没问题，应该是你的EMD工具箱没装好。。
我出的图：

youshen

回复 ChaChing 的帖子

谢谢前辈。。我又懂了一个。。

cassie920920

D:\My Dropbox\fig1.png你好，我也是HHT的新手，我使用您的程式，出來的結果也不對，不知道有哪位大大知道EMD工具箱該如何安裝?我點package_emd工具箱內的install_emd.m檔來安裝，安裝後hhspectrum結果還是不正確，但是程式卻沒有錯誤!!!

chenlu1986

回复 9 # cassie920920 的帖子

论坛里关于EMD工具箱和安装的问题的讨论很多很多，请您搜索一下。

jmsdxjxxy

ChaChing 发表于 2010-10-10 21:11

                               
登录/注册后可看大图

回复 youshen 的帖子

LS没有这个instfreq函数! 以下式检查下

??? Error using ==> hhspectrum at 40
Not enough input arguments.
我也是运行该程序，出现上面的问题，请问是怎么回事？

ChaChing

jmsdxjxxy 发表于 2013-3-31 10:53

                               
登录/注册后可看大图

??? Error using ==> hhspectrum at 40
Not enough input arguments.
我也是运行该程序，出现上面的问题 ...

Not enough input arguments.
字面上是说输入引数不足
检查下hhspectrum的输入及如何呼叫该程序

ChaChing

gunwong 发表于 2010-9-2 22:54

                               
登录/注册后可看大图

??? Function 'real' is not defined for values of class 'cell'.
是不是要找一个转换程序？？
急切等大 ...

刚发现没人针对LZ问题有直接回应, LZ也忘记分享处置方式!
若猜测没错, LZ应该是有自定义程序与built-in程序相同名称了, 而且其优先顺序又较为大
同样可以使用以下命令检查(同6F)
which real -all