m1_fig.m

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

clear;fb=500;fc=2000;fs=12000;Ac=1;f1=500;f2=1000;f3=1500;f4=2000;
for i=1:30
   for k=1:100
       k
    d2=randint(1,300,2);d4=randint(1,300,4);d16=randint(1,300,16);d8=randint(1,300,8);
    s1=ask2(d2,fb,fc,fs,Ac);s2=ask4(d4,fb,fc,fs,Ac);s3=qam16(d16,fb,fc,fs,Ac);
    s4=fsk4(d4,fb,f1,f2,f3,f4,fs,Ac);s5=psk8(d8,fb,fc,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');
    Ma11(k)=M1_func(ss1);Ma12(k)=M1_func(ss2);Ma13(k)=M1_func(ss3);
    Ma14(k)=M1_func(ss4);Ma15(k)=M1_func(ss5);
    end
   M11(i)=mean(Ma11);M12(i)=mean(Ma12);M13(i)=mean(Ma13);
   M14(i)=mean(Ma14);M15(i)=mean(Ma15);
end
plot(M11,'-*')
hold on
plot(M12,'-^')
plot(M13,'-s')
plot(M14,'-d')
plot(M15,'-p')
xlabel('snr');ylabel('M1');
grid
legend('ask2','ask4','qam16','fsk4','psk4',0);