mp1_fig.m

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

clear;fb=500;fc=2000;fs=12000;Ac=1;
for i=1:30
    for k=1:100
       d2=randint(1,200,2);d4=randint(1,200,4);d8=randint(1,200,8);
       s3=psk2(d2,fb,fc,fs,Ac);s4=psk4(d4,fb,fc,fs,Ac);s5=psk8(d8,fb,fc,fs,Ac);
       ss3=awgn(s3,i,'measured');ss4=awgn(s4,i,'measured');ss5=awgn(s5,i,'measured');
       Map13(k)=Mp1_func(ss3,fs,fc);Map14(k)=Mp1_func(ss4,fs,fc);Map15(k)=Mp1_func(ss5,fs,fc);
    end
    Mp13(i)=mean(Map13);Mp14(i)=mean(Map14);Mp15(i)=mean(Map15);
end
plot(Mp13,'-*')
hold on
plot(Mp14,'-^')
plot(Mp15,'-s')
xlabel('snr');ylabel('Mp1');
grid
legend('psk2','psk4','psk8',0);