课程作业记录2：自然映射和格雷映射下的QPSK与16QAM的Matlab仿真

一次作业记录，出的图全为SNR vs BER。

1.QPSK

1.1 自然映射
自然映射下的QPSK如图所示：

BER计算代码如下：

%本代码为自然映射下的QPSK调制
clear all;                                            
close all;  

%生成10000个01/11/00/10的随机数列作为初始信号  
r=randint(1,20000);

%分成两列，第一列代表实部，第二列代表虚部
s=reshape(r,10000,2);
rs=s;
k=1;

%按星座图对信号进行调制
for i=1:1:10000
    if (s(i,1)==1 && s(i,2)==0)
        s(i,1)=-sqrt(2)/2;
        s(i,2)=-sqrt(2)/2;
    end
    if (s(i,1)==0 && s(i,2)==1)
        s(i,1)=-sqrt(2)/2;
        s(i,2)=sqrt(2)/2;
    end
    if (s(i,1)==0 && s(i,2)==0)
        s(i,1)=sqrt(2)/2;
        s(i,2)=sqrt(2)/2;
    end
    if (s(i,1)==1 && s(i,2)==1)
        s(i,1)=sqrt(2)/2;
        s(i,2)=-sqrt(2)/2;
    end
end

for snrdb=1:1:10
   
    %给信号加上功率为snr（db）的高斯白噪声信号
    x=awgn(s,snrdb,'measured');
    
    %解调
    for j=1:1:10000
        if x(j,1)>=0 && x(j,2)>