一个关于hht_toolbox_20040808(Huang)的问题啊
求助一下hhtpm.m函数的注释问题啊？function = hhtpm(data,nyy,t0,t1,ni)
%
% =hhtp_m(data,ny,t0,t1,ni): Hilbert spectrum -- phase-image.
% ny: the frequency resolution.
% t0:
the true strat time.
% t1: the true end time.
% ni: vertical smouth points (Ex. ni=5)
%The outputs are:
% nt: a 2-D matrix of the hilbert spectrum
% t:
the time axis
% f: the frequency axis
%NOTE: (1) On the screen it will appear 'max-frequency' and 'min-frequency'
% (2) contour(t,f,nt) will plot the contour of the hilbert spectrum.
% Z. SHEN 07-2-1995 Initial
%
D. Xiang 03-26-2002Modify
% At The Johns Hopkins University.
= size(data); %read the dimensions
%-----Hilbert Transform --------------------!
dt=(t1-t0)/(npt-1);
data=hilbert(data);
a=angle(exp(j*pi)*data);
omg=diff(data)./data(1:npt-1,:)/(2*pi*dt);
%---- force the negative omg to be zero, DX. ---!
omg=(omg+abs(omg))/2.;
% add the last row to omg so that it has the same dimention
omg=;
clear data
%----- get local frequency -----------------!
wmx=max(max(omg))
wmn=min(min(omg))
dw=wmx-wmn;
if wmn<0.
error('Error: negative frequency appears!');
end
clear p;
%----- Construct the ploting matrix --------!
h1=zeros(npt,nyy+1);
p=round(nyy*(omg-wmn)/dw)+1;
for j1=1:npt
for i1=1:knb
ii1=p(j1,i1);
h1(j1,ii1)=h1(j1,ii1)+a(j1,i1);
if (ii1 > 1)&(ii1 <= nyy)
h1(j1,ii1-1)=h1(j1,ii1);
h1(j1,ii1+1)=h1(j1,ii1);
end
end
end
%---- the results ------------------!
w=linspace(wmn,wmx,nyy+1)';
xs=linspace(t0,t1,npt)';
h=flipud(rot90(abs(h1)));
(data,nyy,t0,t1,ni)这几个参数怎么设置呢？
用过这个工具箱的大虾们指导一下啊 ！！！
页:
[1]