function?f=DEO3S（x）
%??x?必須為行數(shù)據(jù)
%??Discrete?Energy?Operator?3?Samples.
%??三點對稱差分能量算子
%??使用傳遞函數(shù)法
%??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（PxNs）;
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））;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????413??2006-02-25?16:17??Teager\DEO3S.m

?????文件????????444??2008-03-28?10:05??Teager\DESA.m

?????目錄??????????0??2014-12-25?21:46??Teager

-----------?---------??----------?-----??----

??????????????????857????????????????????3