five_features.asv

通信系统中用于设别调制方式的5个特征参数的提取

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;
%提取零中心归一化瞬时幅度绝对值的标准偏差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;

%提取零中心非弱信号段瞬时相位非线性分量绝对值的标准偏差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;



%提取零中心非弱信号段瞬时相位非线性分量的标准偏差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;


%提取零中心归一化非弱信号段瞬时频率绝对值的标准偏差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