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

1、精選優(yōu)質(zhì)文檔-傾情為你奉上%16QAM%等效基帶調(diào)制clc;clear all;close all;%關(guān)閉所有窗口，比如之前畫(huà)的圖nsymbol = 10000;%每種信噪比下符號(hào)數(shù)的發(fā)送符號(hào)數(shù)，106個(gè)數(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,nsymb

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

3、(4*i-1)=0 & data(4*i)=1 data_mod2(i) = 1; elseif data(4*i-1)=1 & data(4*i)=1 data_mod2(i) = -1; elseif data(4*i-1)=1 & data(4*i)=0 data_mod2(i) = -3; end elseif data(4*i-3)=0 & data(4*i-2)=1 data_mod1(i) = 1; if data(4*i-1)=0 & data(4*i)=0 data_mod2(i) = 3; elseif data(4*i-1)=0 & data(4*i)=1 data_mo

4、d2(i) = 1; elseif data(4*i-1)=1 & data(4*i)=1 data_mod2(i) = -1; elseif data(4*i-1)=1 & data(4*i)=0 data_mod2(i) = -3; end elseif data(4*i-3)=1 & data(4*i-2)=1 data_mod1(i) = -1; if data(4*i-1)=0 & data(4*i)=0 data_mod2(i) = 3; elseif data(4*i-1)=0 & data(4*i)=1 data_mod2(i) = 1; elseif data(4*i-1)=

5、1 & data(4*i)=1 data_mod2(i) = -1; elseif data(4*i-1)=1 & data(4*i)=0 data_mod2(i) = -3; end elseif data(4*i-3)=1 & data(4*i-2)=0 data_mod1(i) = -3; if data(4*i-1)=0 & data(4*i)=0 data_mod2(i) = 3; elseif data(4*i-1)=0 & data(4*i)=1 data_mod2(i) = 1; elseif data(4*i-1)=1 & data(4*i)=1 data_mod2(i) =

6、 -1; elseif data(4*i-1)=1 & data(4*i)=0 data_mod2(i) = -3; end end endSNR_dB = 1:10;%信噪比dB形式SNR = 10.(SNR_dB/10);%信噪比轉(zhuǎn)化為線性值for loop= 1:10 sigma = sqrt(5/SNR(loop); s_receive1 = data_mod1 + sigma * randn(1,length(data_mod1); s_receive2 = data_mod2 + sigma * randn(1,length(data_mod2);%? for k = 1:nsym

7、bol if s_receive1(k)2 data_demod1(k) = 0; data_demod2(k) = 0; if s_receive2(k)2 data_demod3(k) = 0; data_demod4(k) = 0; elseif s_receive2(k)0 & s_receive2(k)-2 & s_receive2(k)0 data_demod3(k) = 1; data_demod4(k) = 1; elseif s_receive2(k)0 & s_receive1(k)2 data_demod3(k) = 0; data_demod4(k) = 0; else

8、if s_receive2(k)0 & s_receive2(k)-2 &s_receive2(k)0 data_demod3(k) = 1; data_demod4(k) = 1; elseif s_receive2(k)-2 data_demod3(k) = 1; data_demod4(k) = 0; end elseif s_receive1(k)-2 data_demod1(k) = 1; data_demod2(k) = 1; if s_receive2(k)2 data_demod3(k) = 0; data_demod4(k) = 0; elseif s_receive2(k)

9、0 & s_receive2(k)-2 & s_receive2(k)0 data_demod3(k) = 1; data_demod4(k) = 1; elseif s_receive2(k)-2 data_demod3(k) = 1; data_demod4(k) = 0; end elseif s_receive1(k)2 data_demod3(k) = 0; data_demod4(k) = 0; elseif s_receive2(k)0 & s_receive2(k)-2 & s_receive2(k)0 data_demod3(k) = 1; data_demod4(k) =

10、1; elseif s_receive2(k)-2 data_demod3(k) = 1; data_demod4(k) = 0; end end if (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; end end errorratio(loop) = errornumber(loop)/nsymbol;endt = 0.75*erfc(sqrt(3*SNR/30);ser_theory = 1-(1-t).2;%理論誤碼率，注意Q函