if2phase.m

来自「一套关于时域分析的matlab工具箱软件集,包括源码书籍」· M 代码 · 共 14 行

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