求助：disp_hhs.m
经过研究 ，我的程序终于没有错误了，但是边际谱出来了，可是HHT谱却没有显示。（只有运行出了，图谱版但是上面什么也没有：没有坐标轴，也没有坐标）请求高手指点，能够简要的帮我解释一下disp_hhs.m是什么原理画谱的，
function disp_hhs(varargin)
error(nargchk(1,3,nargin));
fs = 0;
inf = -20;
im = varargin{1};
t = 1:size(im,2);
switch nargin
case 1
%raf
case 2
if isscalar(varargin{2})
inf = varargin{2};
else
t = varargin{2};
end
case 3
if isvector(varargin{2})
t = varargin{2};
inf = varargin{3};
else
inf = varargin{2};
fs = varargin{3};
end
case 4
t = varargin{2};
inf = varargin{3};
fs = varargin{4};
end
if isempty(inf)
inf = -20;
end
if inf > 0
inf = -inf;
elseif inf == 0
error('inf must be nonzero')
end
M=max(max(im));
warning off
im = 10*log10(im/M);
warning on
figure
if fs == 0
imagesc(t,,im,);
ylabel('normalized frequency')
else
imagesc(t/fs,,im,);
ylabel('frequency')
end
set(gca,'YDir','normal')
xlabel('time')
title('Hilbert-Huang spectrum')
研究了一下，可是发现基本没看明白：:@(
小妹在此万分感谢！:@)
附加：和实际品频率和采样点个数有关吗？
和实际品频率和采样点个数有关吗？我分析的实验测得的数据，在读取数据仿真是设定其
Ns=length(y2);%采样点的长度 （y2为采样数据）
fs=1000; % Sampling frequenc
t=0:(1/fs):((Ns-1)/fs); % time series
freq=*fs/Ns; % frequency axis
请求论坛的朋友们给予帮助，谢谢了:@)
页:
[1]