function feature_vector=five_features(Amplititude,Phase,Frequency)
feat
%提取零中心歸一化瞬時幅度之譜密度的最大值garma_max

    
   amp_mean=mean(Amplititude_2ask);
   amp_cn=Amplititude_2ask/amp_mean-1;
   garma_max1=max(abs(fft(amp_cn.^2)))/size(Amplititude_2ask,1);
   garma_max_2ask=garma_max1;
   feature_vector(k,1)=garma_max_2ask;
%提取零中心歸一化瞬時幅度絕對值的標(biāo)準(zhǔn)偏差delta_aa

    amp_mean=mean(Amplititude_2ask);
    amp_cn=Amplititude_2ask/amp_mean-1;
    amp_cn_mean_square=(sum(abs(amp_cn))/size(amp_cn,1)).^2;
    amp_cn_square_mean=sum((amp_cn).^2)/size(amp_cn,1);
    delta_aa_2ask=sqrt(amp_cn_square_mean-amp_cn_mean_square);
    stochastic(k,2)=delta_aa_2ask;

%提取零中心非弱信號段瞬時相位非線性分量絕對值的標(biāo)準(zhǔn)偏差delta_ap

   
    phase_mean=mean(Phase_2ask);
    phase_nl=Phase_2ask-phase_mean;
 
    amp_mean=mean(Amplititude_2ask);
    amp_n=Amplititude_2ask/amp_mean;
    nCount=0;
    for j=1:size(Amplititude_2ask)
        if amp_n(j)>0.9
            nCount=nCount+1;
        else phase_nl(j)=0;
        end
    end
    phase_nl_mean_square=(sum(abs(phase_nl))/nCount).^2;
    phase_nl_square_mean=sum((phase_nl).^2)/nCount;
    delta_ap_2ask=abs(sqrt(phase_nl_mean_square-phase_nl_square_mean));
    stochastic(k,3)=delta_ap_2ask;



%提取零中心非弱信號段瞬時相位非線性分量的標(biāo)準(zhǔn)偏差delta_dp

    phase_mean=mean(Phase_2ask);
    phase_nl=Phase_2ask-phase_mean;
  
    amp_mean=mean(Amplititude_2ask);
    amp_n=Amplititude_2ask/amp_mean;
    nCount=0;
    for j=1:size(Amplititude_2ask)
        if amp_n(j)>0.9
            nCount=nCount+1;
        else phase_nl(j)=0;
        end
    end
    phase_nl_mean_square=(sum(phase_nl)/nCount).^2;
    phase_nl_square_mean=sum((phase_nl).^2)/nCount;
    delta_dp_2ask=abs(sqrt(phase_nl_mean_square-phase_nl_square_mean));
    stochastic(k,4)=delta_dp_2ask;


%提取零中心歸一化非弱信號段瞬時頻率絕對值的標(biāo)準(zhǔn)偏差delta_af

    fre_mean=mean(Frequency_2ask);
    fre_nl=Frequency_2ask-fre_mean;
 
    amp_mean=mean(Amplititude_2ask);
    amp_n=Amplititude_2ask/amp_mean;
    nCount=0;
    for j=1:size(Amplititude_2ask)
        if amp_n(j)>1
            nCount=nCount+1;
        else fre_nl(j)=0;
        end
    end
    fre_nl_mean_square=(sum(abs(fre_nl))/nCount).^2;
    fre_nl_square_mean=sum((fre_nl).^2)/nCount;
    delta_af_2ask=abs(sqrt(fre_nl_mean_square-fre_nl_square_mean));
    stochastic(k,5)=delta_af_2ask;
    
end