if2phase.m

matlab工具箱,可以为一般的时频分析的用户提供很好的学习教材.里面的代码详细周全!注释清楚!

function phi=if2phase(iflaw)
%IF2PHASE Generate the phase from the instantaneous frequency

%	O. Lemoine - February 1996.
% experimental
	 
N=length(iflaw);
phi=zeros(N,1);

%phi=2*pi*cumsum(iflaw);
for k=1:N,
 t=1:k;
 phi(k)=2*pi*integ(iflaw(t),t);
end