mf1_fig.m

调制识别：基于决策论的通信信号调制识别实现

clear;fb=500;fc=2000;fs=12000;Ac=1;f1=500;f2=1000;f3=1500;f4=2000;
f5=750;f6=1250;f7=1750;f8=3250;
for i=1:30
    for k=1:100
        d2=randint(1,300,2);d4=randint(1,300,4);d8=randint(1,300,8);
        s1=fsk2(d2,fb,f1,f4,fs,Ac);s2=fsk4(d4,fb,f1,f2,f3,f4,fs,Ac);
        s3=psk2(d2,fb,fc,fs,Ac);s4=psk4(d4,fb,fc,fs,Ac);
        s5=fsk8(d8,fb,f1,f2,f3,f4,f5,f6,f7,f8,fs,Ac);
        ss1=awgn(s1,i,'measured');ss2=awgn(s2,i,'measured');
        ss3=awgn(s3,i,'measured');ss4=awgn(s4,i,'measured');
        ss5=awgn(s5,i,'measured');
        Maf11(k)=Mf1_func(ss1,fs,Ac);Maf12(k)=Mf1_func(ss2,fs,Ac);
        Maf13(k)=Mf1_func(ss3,fs,Ac);Maf14(k)=Mf1_func(ss4,fs,Ac);
        Maf15(k)=Mf1_func(ss5,fs,Ac);
    end
    Mf11(i)=mean(Maf11);Mf12(i)=mean(Maf12);Mf13(i)=mean(Maf13);
    Mf14(i)=mean(Maf14);Mf15(i)=mean(Maf15);
end
plot(Mf11,'-*')
hold on
plot(Mf12,'-^')
plot(Mf13,'-s')
plot(Mf14,'-d')
plot(Mf15,'-p')
xlabel('snr');ylabel('Mf1');
grid
legend('fsk2','fsk4','psk2','psk4','fsk8',0);