function?cha=ninefeature（xmtafcGfs）

fd=1000;%symbol?freqency
hi=hilbert（xm）;
sph=2*pi*fc*ta;
samp=abs（hi）;
spha=unwrap（angle（hi））-sph;
at=1;
%first?rmax
ma=mean（samp）;
sampn=samp/ma;
sampcn=sampn-1;
rmax=max（abs（fft（sampcn））.^2/length（sampcn））;
%second?aap
ph=mean（spha）;
phnl=spha-ph;
ind=find（sampn>at）;
c=length（ind）;
for?m=1:c
????phse（m）=phnl（ind（m））;
end;
aap=sqrt（mean（phse.^2）-mean（abs（phse））^2）;
%third?adp
adp=sqrt（mean（phse.^2）-mean（phse）^2）;
%four?p
fcn=round（fc*G/fs-1）;
xfft=fft（xm）;
spl=xfft（1:fcn）;
spu=xfft（（fcn+1+1）:（2*fcn+1））;
pl=sum（abs（spl）.^2）;
pu=sum（abs（spu）.^2）;
p=（pl-pu）/（pl+pu）;
%five?aaa
aaa=sqrt（mean（sampcn.^2）-mean（abs（sampcn））^2）;
%six?aaf
freq（1）=0;
for?k=1:lengt