deo3s.m

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

function f=DEO3S(x)
%  x 必须为行数据
%  Discrete Energy Operator 3 Samples.
%  三点对称差分能量算子
%  使用传递函数法
%  H(z)=z(1+2*z^-1+z^-2)/4;
Px=Phid(x);

Ns=length(Px);
w=2*pi*(-Ns/2:Ns/2)/Ns;
w=[w(1:Ns/2),w(Ns/2+1:Ns)];%去0点
z=exp(i*w);
Hz=z.*( 1+2*z.^-1+z.^-2 )/4;

Xz=fft(Px,Ns);
Xz=[Xz(Ns/2+1:Ns),Xz(1:Ns/2)];%重新排列

Yw=Hz.*Xz;
 Yw=[Yw(Ns/2+1:Ns),Yw(1:Ns/2)];
 f=real(ifft(Yw));