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

求助，hht画不出imf图形！

我是菜鸟，在论坛找了几个程序，在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);
= 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);

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

Error in ==> real at 15
= 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
= hhspectrum(imf);

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

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

dzkt 发表于 2010-9-4 20:50

程序一点问题也没有

shinhwaforever 发表于 2010-10-10 16:35

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

youshen 发表于 2010-10-10 20:13

??? Undefined command/function 'instfreq'.

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

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

回复 youshen 的帖子

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

youshen 发表于 2010-10-11 10:16

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

youshen 发表于 2010-10-11 10:18

回复 ChaChing 的帖子

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

cassie920920 发表于 2011-9-13 19:09

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

chenlu1986 发表于 2011-9-15 08:49

回复 9 # cassie920920 的帖子

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

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

ChaChing 发表于 2010-10-10 21:11 static/image/common/back.gif
回复 youshen 的帖子

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


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

ChaChing 发表于 2013-3-31 11:14

jmsdxjxxy 发表于 2013-3-31 10:53 static/image/common/back.gif
??? Error using ==> hhspectrum at 40
Not enough input arguments.
我也是运行该程序，出现上面的问题 ...

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

ChaChing 发表于 2013-3-31 11:25

gunwong 发表于 2010-9-2 22:54 static/image/common/back.gif
??? Function 'real' is not defined for values of class 'cell'.
是不是要找一个转换程序？？
急切等大 ...

刚发现没人针对LZ问题有直接回应, LZ也忘记分享处置方式!
若猜测没错, LZ应该是有自定义程序与built-in程序相同名称了, 而且其优先顺序又较为大
同样可以使用以下命令检查(同6F)
which real -all
页: [1]
查看完整版本: 求助，hht画不出imf图形！