judge.m

来自「调制识别：基于决策论的通信信号调制识别实现」· M 代码 · 共 42 行

function result=judge(s,fs,fc)
%利用六个参数对输入信号的调制方式进行判决
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%门限值
t_m11=0.325;t_m12=0.195;t_m2=0.306;t_mf1=0.1815;
t_mf2=0.4185;t_mp1=0.395;t_mp2=0.404;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
h=hilbert(s);
a1=A_func(h);
M1=M_func(a1);
if M1>=t_m11
    a2=A_func(a1);M2=M_func(a2);
    if M2>=t_m2
        result=2; %代表4ASK信号
    else
        result=1; %代表2ASK信号
    end
elseif M1>=t_m12
    result=3; %代表16QAM信号
else
    f=dif(h,fs);Mf1=M_func(A_func(f));
    if Mf1>=t_mf1
        f0=f0_func(f);Mf2=M_func(A_func(f0));
        if Mf2>=t_mf2
            result=5; %代表4FSK信号
        else
            result=4; %代表2FSK信号
        end
    else
        fai=fai0_func(h,fs,fc);
        p1=A_func(fai);Mp1=M_func(p1);
        if Mp1>=t_mp1
            Mp2=M_func(A_func(p1));
            if Mp2>=t_mp2
                result=8; %代表8PSK信号
            else
                result=7; %代表4PSK信号
            end
        else
            result=6; %代表2PSK信号
        end
    end
end