integrold.m

来自「一种新的时频分析方法的matlab源程序。」· M 代码 · 共 21 行

function [xx]=integr(x,dt)

%   [xx]=integr(x,dt)
%   [m,n]=size(x)
%   defult x(:,1) to use least
%   column operation
%   Processes only one column!
%   NOTE: sum should be initialized to 0 (not xx) JM

[m,n]=size(x) ;

xx(1,1)=dt*x(1,1)/2 ;
sum=xx(1,1) ;
for jj=2:m ;
   sum=sum+(x(jj-1,1)+x(jj,1))*dt/2    ;
   xx(jj,1)=sum ;
end