(完整版)16QAM理論誤碼率與實(shí)際誤碼率MATLAB仿真程序

1、%16QAM%等效基帶調(diào)制clc;clearall;closeall;%關(guān)閉所有窗口，比如之前畫的圖nsymbol=10000;%每種信噪比下符號數(shù)的發(fā)送符號數(shù)，10八6個數(shù)errornumber=zeros(1,10);errorratio=zeros(1,10);data=randi(0,1,1,4*nsymbol);%產(chǎn)生1行，4*nsymbol列均勻分布的隨機(jī)數(shù)data_mod1=zeros(1,nsymbol);%各參量初始化data_mod2=zeros(1,nsymbol);data_demod1=zeros(1,nsymbol);data_demod2=zeros(1,nsym

2、bol);data_demod3=zeros(1,nsymbol);data_demod4=zeros(1,nsymbol);s_receive1=zeros(1,nsymbol);s_receive2=zeros(1,nsymbol);fori=1:nsymbol%格雷碼調(diào)制到星座圖的映射ifdata(4*i-3)=0&data(4*i-2)=0data_mod1(i)=3;ifdata(4*i-1)=0&data(4*i)=0data_mod2(i)=3;elseifdata(4*i-1)=0&data(4*i)=1data_mod2(i)=1;elseifdata(4*i-1)=1&da

3、ta(4*i)=1data_mod2(i)=-1;elseifdata(4*i-1)=1&data(4*i)=0data_mod2(i)=-3;endelseifdata(4*i-3)=0&data(4*i-2)=1data_mod1(i)=1;ifdata(4*i-1)=0&data(4*i)=0data_mod2(i)=3;elseifdata(4*i-1)=0&data(4*i)=1data_mod2(i)=1;elseifdata(4*i-1)=1&data(4*i)=1data_mod2(i)=-1;elseifdata(4*i-1)=1&data(4*i)=0data_mod2(i

4、)=-3;endelseifdata(4*i-3)=1&data(4*i-2)=1data_mod1(i)=-1;ifdata(4*i-1)=0&data(4*i)=0data_mod2(i)=3;elseifdata(4*i-1)=0&data(4*i)=1data_mod2(i)=1;elseifdata(4*i-1)=1&data(4*i)=1data_mod2(i)=-1;elseifdata(4*i-1)=1&data(4*i)=0data_mod2(i)=-3;endelseifdata(4*i-3)=1&data(4*i-2)=0data_mod1(i)=-3;ifdata(4*

5、i-1)=0&data(4*i)=0data_mod2(i)=3;elseifdata(4*i-1)=0&data(4*i)=1data_mod2(i)=1;elseifdata(4*i-1)=1&data(4*i)=1data_mod2(i)=-1;elseifdata(4*i-1)=1&data(4*i)=0data_mod2(i)=-3;endendendSNR_dB=1:10;%信噪比dB形式SNR=104(SNR_dB/10);%信噪比轉(zhuǎn)化為線性值forloop=1:10sigma=sqrt(5/SNR(loop);s_receive1=data_mod1+sigma*randn(1

6、,length(data_mod1);s_receive2=data_mod2+sigma*randn(1,length(data_mod2);%?fork=1:nsymbolifs_receive1(k)2data_demod1(k)=0;data_demod2(k)=0;ifs_receive2(k)2data_demod3(k)=0;data_demod4(k)=0;elseifs_receive2(k)0&s_receive2(k)-2&s_receive2(k)0data_demod3(k)=1;data_demod4(k)=1;elseifs_receive2(k)0&s_rece

7、ive1(k)2data_demod3(k)=0;data_demod4(k)=0;elseifs_receive2(k)0&s_receive2(k)-2&s_receive2(k)0data_demod3(k)=1;data_demod4(k)=1;elseifs_receive2(k)-2data_demod3(k)=1;data_demod4(k)=0;endelseifs_receive1(k)-2data_demod1(k)=1;data_demod2(k)=1;ifs_receive2(k)2data_demod3(k)=0;data_demod4(k)=0;elseifs_re

8、ceive2(k)0&s_receive2(k)-2&s_receive2(k)0data_demod3(k)=1;data_demod4(k)=1;elseifs_receive2(k)-2data_demod3(k)=1;data_demod4(k)=0;endelseifs_receive1(k)2data_demod3(k)=0;data_demod4(k)=0;elseifs_receive2(k)0&s_receive2(k)-2&s_receive2(k)0data_demod3(k)=1;data_demod4(k)=1;elseifs_receive2(k)-2data_de

9、mod3(k)=1;data_demod4(k)=0;endendif(data_demod1(k)=data(k*4-3)|(data_demod2(k)=data(k*4-2)|(data_demod3(k)=data(k*4-1)|(data_demod4(k)=data(k*4)errornumber(loop)=errornumber(loop)+1;endenderrorratio(loop)=errornumber(loop)/nsymbol;endt=0.75*erfc(sqrt(3*SNR/30);ser_theory=1-(1-t)42;%理論誤碼率，注意Q函數(shù)和誤差函數(shù)的對應(yīng)