信号检测实验四参考答案一.docx

实验四、误码率的分析1、假定QPSK信号具有如下的星座图模式：1) 试利用二元数字通信系统的误码性能精确推导QPSK系统的误码率。2) 假设系统中的噪声是高斯白噪声，利用matlab对该系统在工作区间（信噪比010dB）下的误码性能进行仿真。答： 表示每比特信号的能量，表示每个符号信号的能量 对于二元数字通信系统：二进制相移键控BPSK： 对于QPSK系统：其比特错误率等于BPSK的比特错误率，而误码率则是 当性噪比较高时，MATLAB 程序代码如下代码如下：clear;clc;SINR_dB=0:0.2:10;%信噪比m=length(SINR_dB);number_of_errors=zeros(1,m);%用于保存Qpsk调制下每一个信噪比下对应的误码率for i=1:1:m %4PSK调制,也叫QPSKM=4; k=log2(M); n=1200; %比特序列长度 samp=1; %过采样率 x=randint(n,1); %生成随机二进制比特流 x4=reshape(x,k,length(x)/k); %将原始的二进制比特序列每2个一组xsym=bi2de(x4.,left-msb); %将矩阵转化为相应的16进制信号序列 y=modulate(modem.pskmod(M),xsym); %用16QAM调制器对信号进行调制 yn=awgn(y,SINR_dB(i),measured); %加入高斯白噪声 yd=demodulate(modem.pskdemod(M),yn); %此时解调出来的是16进制信号 z=de2bi(yd,left-msb); %转化为对应的二进制比特流 z=reshape(z.,numel(z),1); number_of_errors(i),bit_error_rate2(i)=biterr(x,z); endfiguresemilogy(SINR_dB,number_of_errors,blue*-);xlabel(信噪比);ylabel(误码率);title(QPSK系统的误码率);legend(QPSK);仿真结果如下：2、若信号为16PSK，重做1。答：（1）有题意设其误码率为：且有，其中，当性噪比较高时，2）MATLAB仿真代码如下：clear;clc;SINR_dB=0:0.2:10;%信噪比m=length(SINR_dB);number_of_errors=zeros(1,m);%用于保存16psk调制下每一个信噪比下对应的误码率for i=1:1:m %采用16psk调制M=16; k=log2(M); n=1200; %比特序列长度 samp=1; %过采样率 x=randint(n,1); %生成随机二进制比特流 x4=reshape(x,k,length(x)/k); %将原始的二进制比特序列每四个一组xsym=bi2de(x4.,left-msb); %将矩阵转化为相应的16进制信号序列 y=modulate(modem.pskmod(M),xsym); %用16QAM调制器对信号进行调制 yn=awgn(y,SINR_dB(i); %加入高斯白噪声 yd=demodulate(modem.pskdemod(M),yn); %此时解调出来的是16进制信号 z=de2bi(yd,left-msb); %转化为对应的二进制比特流 z=reshape(z.,numel(z),1); number_of_errors(i),bit_error_rate(i)=biterr(x,z); endfiguresemilogy(SINR_dB,number_of_errors,red*-); xlabel(信噪比);ylabel(误码率);title(16PSK系统的误码率);legend(16PSK);仿真结果如下：3、若信号为16QAM，且具有如下的星座图，重做1。答：（2）MATLAB仿真代码如下：clear;clc;SINR_dB=0:0.2:10;%信噪比m=length(SINR_dB);number_of_errors=zeros(1,m);%用于保存16QAM调制下每一个信噪比下对应的误码率for i=1:1:m%采用16QAM调制M=16; k=log2(M); n=1200; %比特序列长度 samp=1; %过采样率 x=randint(n,1); %生成随机二进制比特流 x4=reshape(x,k,length(x)/k); %将原始的二进制比特序列每四个一组xsym=bi2de(x4.,left-msb); %将矩阵转化为相应的16进制信号序列 y=modulate(modem.qammod(M),xsym); %用16QAM调制器对信号进行调制 yn=awgn(y,SINR_dB(i),measured); %加入高斯白噪声 yd=demodulate(modem.qamdemod(M),yn); %此时解调出来的是16进制信号 z=de2bi(yd,left-msb); %转化为对应的二进制比特流 z=reshape(z.,numel(z),1); number_of_errors(i),bit_error_rate(i)=biterr(x,z); endfiguresemilogy(SINR_dB,number_of_errors,red*-); xlabel(信噪比);ylabel(误码率);title(16QAM系统的误码率);legend(16QAM);仿真结果