通信原理课程仿真：数字传输系统的最佳接收与误码分析

1.概论

   在数字通信系统中，信源是信息的来源，但信源输出的是消息(或符号）来表达所要传达的信息。若通信信道的传递函数为低通的，则称此信道为基带信道，适合于传送脉冲波形，所以将数字基带信号通过基带信道传输的系统称为数字基带传输系统。
   而在数字基带传输系统中，确定信号在受到加性高斯白噪声干扰下的最佳接收为采用匹配滤波器，使得在最佳采样时刻的信噪比最大，然后再选择合适的判决门限进行判决，可以使得误比特率最小。
   在实际数字通信系统中的信道往往是有限的，为此要进一步考虑适合于限带信道传输的限带信号的设计，而在本文中采用的为基于无码间干扰基带传输的升余弦滤波器。

2. 信道带宽无限时的单极性基带传输

2.1 系统模型

  当信道的带宽无限时，可以采用矩形波作为传输的信号，此处使用的为单极性不归零码。

---

发送/接收滤波器的冲激响应如下图所示：

---

系统框图如下所示:
 

---

1.单极性不归零码(NRZ)

Matlab代码如下：

%------------------
%系统参数设置
%-----------------
T_start=0;%开始时间
T_stop=10;%截止时间
T=T_stop-T_start;%仿真持续时间
T_sample=1/1000;%采样间隔
f_sample=1/T_sample; % 采样速率
N_sample=T/T_sample;% 采样点数
n=0:N_sample-1;
r_s=1;%transmission rate
alpha=0.25;%df=alpha*rs=25Hz
NumBits=T*r_s;%number of bits
NumCoff=40;%number of coefficients of RRC
SamplingIns=NumCoff/2;
SamplingStart=NumCoff;
%-------------------------------
%Transmitter
%-------------------------------
g_T=ones(1,f_sample/r_s); ;%transmission filter
b1=sign(rand(1,NumBits)-0.5);%原始-1/1序列
b2=sign(rand(1,NumBits));%原始全一序列
b3=(b1+b2)/2;%原始0/1序列
b4=zeros(f_sample/r_s,NumBits);
b4(1,:)=b3;%b7矩阵第一行
b6=reshape(b4,1,f_sample/r_s*NumBits);
for i=1:NumBits
    b4(:,i)=b4(1,i);
end
b5=reshape(b4,1,f_sample/r_s*NumBits);%NRZ单极性
s=conv(b5,g_T);
%--------------------------------
%AWGN channel
%--------------------------------
N_0=10^(-3);
noise_w=wgn(1,length(s),N_0*f_sample,'linear');%产生白噪声
r=s+noise_w;
%figure(1)
%plot(s)
%hold on
%plot(r)
%--------------------------------
%receiver
%--------------------------------
g_R=g_T;%transmission filter
y1=conv(r,g_R);
myeyediagram=eyediagram(y1(1:length(y1)),f_sample/r_s);
b_r=zeros(1,r_s*T);
for i=1:10
    x=1000*(i);
    if(y1(x)<=5*10^5)
        b_r(i)=0;
    else
        b_r(i)=1;
    end
end
b7=b6(1:1000:end);
b_t=sign(b7);
BER=length(find(b_t~= b_r))/NumBits
b_finalw=zeros(f_sample/r_s,r_s*T);
b_finalw(1,:)=b_r;
for i=1:NumBits
    b_finalw(:,i)=b_finalw(1,i);
end
b_final=reshape(b_finalw,1,f_sample/r_s*NumBits);
B1=0.5*erfc(sqrt(1./(2.*N_0)))
%%%%%%%%%%%%%%%%%%%%%画图啦
subplot(511);plot(b_r);axis([1,10,0,1]);
title('原始单极性序列');
subplot(512);plot(b5)
title('原始矩形脉冲序列');
subplot(513);plot(y1)
title('通过匹配滤波器后的信号');
subplot(514);plot(b_r)
title('判决出的结果');
subplot(515);plot(b_final)
title('接收端收到结果');

运行结果如下：

---

       ps:在图中我们可以看到五个波形，第一个是采样的[0,1]序列,第二个是矩形脉冲序列，第三个是匹配滤波后得到的波形，四五分别是判决出的采样点与零阶保持后的结果。
 

---

单极性码的眼图 $(eyediagram)$

---

 

BER(实际误码率)       B1(实际误码率）

---

  

2.4 仿真结果分析

 

• 无带限时，发送滤波器为$T_g=Tb$的矩形脉冲，匹配滤波器顺理成章地由相同的脉冲担任，要考虑的只是关于判决方面的问题，可以看到误码率还是很低的。
• 不归零码码的带宽是归零码的一半，原因是连续谱中$sinc$函数第一过零点多了一倍，导致整个功率谱展宽了一倍，但是离散谱不受影响，其只与信息速率有关。

3. 信道带宽受限时的双极性基带传输

3.1.双极性不归零码

3.2 误码性能分析

  数字PAM信号在限带及加性白高斯噪声信道条件下，最佳基带传输系统的误码率计算公式带宽不受限时的加性白高斯噪声信道条件下匹配滤波器最佳接收的误码率计算公式是相同的，根据下图很容易推得。
 $$P_e=P(binart0)\ast Q（\frac{A-V_t}{{\sigma }_0}）+P(binart1)\ast Q（\frac{A+V_t}{{\sigma }_0}）$$
 $\because P(binart0)=P(binart1),Q（\frac{A-V_t}{{\sigma }_0}）=Q（\frac{A+V_t}{{\sigma }_0}）.$
 $\therefore V_t=0,$
  $P_e=Q（\frac{A-V_t}{{\sigma }_0}$