desa.m

来自「能量算子的matlab实现程序及其详细算法」· M 代码 · 共 14 行

function  f=DESA(x)
%能量分离算法求解|a(n)|
% 算法详细描述： open TheoryofDEO.doc
%                     2Ψ[x(n)]
%|a(n)|=-------------------------------------
%         sqrt(  Ψ[ x(n+1)-x(n-1) ]  )
U=2*DEO3S(x);
%L=abs(   real(  sqrt(   DEO3S(  circshift(x',-1)'-circshift(x',1)'  )   )  )   );
L=abs(   (  sqrt(   DEO3S(  circshift(x',-1)'-circshift(x',1)'  )   )  )   );

f=abs(U./L);%应去除分母为0的点
A=find(f==inf);
f(A)=0;