mf2_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=2250;
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=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');
         Maf21(k)=Mf2_func(ss1,fs,Ac);Maf22(k)=Mf2_func(ss2,fs,Ac);Maf23(k)=Mf2_func(ss3,fs,Ac);
    end
    Mf21(i)=mean(Maf21);Mf22(i)=mean(Maf22);Mf23(i)=mean(Maf23);
end
plot(Mf21,'-*')
hold on
plot(Mf22,'-^')
plot(Mf23,'-s')
xlabel('snr');ylabel('Mf2');
grid
legend('fsk2','fsk4','fsk8',0);