2022年通信原理实验报告含matlab程序 2.pdf

《2022年通信原理实验报告含matlab程序 2.pdf》由会员分享，可在线阅读，更多相关《2022年通信原理实验报告含matlab程序 2.pdf（17页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、通信原理实验报告实验二 QPSK通信系统的Monte Carlo仿真一、实验目的1.提高独立学习的能力2.培养发现问题，解决问题，分析问题的能力3.学习 Matlab 的使用4.掌握 4PSK通信系统的Monte Carlo 仿真方法5.掌握 4PSK通信系统的组成原理6.比较编码信号与未编码信号在随机信道中的传输，加深对纠错编码原理的理解。二、实验原理1 调制解调原理一组 M 载波相位调制信号波形的一般表示式为：1,.,1 ,0),22cos(MmMmtftAgtucTm式中tgT是发送滤波器的脉冲形状，决定了传输信号的频谱特性，A 是信号的幅度，PSK信号对所有的m 都具有相同的能量，即s

2、mmdttu2其中的s代表每个传输符号的能量。当tgT是一个矩形脉冲时，定义为TtTtgT0,2这时在符号区间Tt0内传输的信号波形可以表示为（用sA）1,.,1 ,0),22cos(2MmMmtfTtucsm把式中的余弦函数的相角看成两个相角的和，可以将上式表示为tststfMmtgtfMmtgtumsmccTscTsm212sin2sin2cos2cos这里MmsMmssmssmc2sin,2cos，而t1和t2是两个正交基函数，定义为tftgttftgtcTcT2sin,2cos21适当地将tgT归一化，就可以将这两个基函数的能量归一化到1。这样一个相位调制信号可以看作两个正交载波，起幅

3、度取决于在每个信号区间内的相位，因此，数字相位调制信号在几何上可用mcs和mss的二维向量来表示，即MmMmsssm2sin,2cos名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 17 页 - - - - - - - - - 在 AWGN信道中，在一个区间内接受到的带通信号可以表示为tftntftntutntutrcsccmm2sin2cos其中tnc和tns是加性噪声的同相分量和正交分量，将这个接收信号与给出的t1和t2作相关，两个相关器的输出可以表示为：sscsmn

4、MmnMmnsr2sin,2cos这 两个 正交的 噪声 分量tnc和tns是 零均值， 互不 相关的 高斯 随机过 程， 这样200022NnEnEnnEnEnEscscsc，方差为。最佳检测器将接收到的信号向量r 投射到 M 个可能的传输信号向量ms之一上去，并选取对应与最大投影的向量，从而得到相关准则为mmsrsrC ,。由于二相相位调制与二进制PAM 是相同的，所以差错概率为022NQPbpsk式中b是每比特能量。4PSK可以看作是两个在正交载波上的二相相位调制系统，所以1 个比特的差错概率与二相相位调制是一致的。2 信道纠错编码在随机信道中， 错码的出现是随机的，且错码之间是统计独立

5、的。例如，由高斯白噪声引起的错码就具有这种性质。因此，当信道中加性干扰主要是这种噪声时，就称这种信道为随机信道。 由于信息码元序列是一种随机序列，接收端是无法预知的，也无法识别其中有无错码。为了解决这个问题， 可以由发送端的信道编码器在信息码元序列中增加一些监督码元。这些监督码元和信码之间有一定的关系，使接收端可以利用这种关系由信道译码器来发现或纠正可能存在的错码。在信息码元序列中加入监督码元就称为差错控制编码，有时也称为纠错编码。不同的编码方法有不同的检错或纠错能力。有的编码就只能检错不能纠错。那么，为了纠正一位错码，在分组码中最少要增加多少监督位才行呢？编码效率能否提高呢？从这种思想出发进

6、行研究，便导致汉明码的诞生。汉明码是一种能够纠正一位错码且编码效率较高的线性分组码。下面我们介绍汉明码的构造原理。一般说来，若码长为n，信息位数为k，则监督位数rn-k 。如果希望用r 个监督位构造出 r 个监督关系式来指示一位错码的n 种可能位置，则要求2r - 1 n 或 2r k + r + 1 （81）下面我们通过一个例子来说明如何具体构造这些监督关系式。设分组码 (n，k)中 k4，为了纠正一位错码，由式（81）可知，要求监督位数r 3 。若取 r=3，则 n= k + r =7。我们用 650 表示这 7 个码元，用S1 、S2 、S3 表示三个监督关系式中的校正子，则S1 S2

7、S3 的值与错码位置的对应关系可以规定如表18-1 所列。表 8-1 S1 S2 S3错码位置S1 S2 S3错码位置名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 17 页 - - - - - - - - - 001 010 100 011 0 1 2 3101 110 111 000 45 6 无错由表中规定可见，仅当一错码位置在2 、4 、5 或 6 时，校正子S1 为 1；否则 S1 为0。这就意味着2 、4 、5 和 6四个码元构成偶数监督关系S1 65 4 2

8、 （82）同理， 1 、3 、5 和 6构成偶数监督关系S2 65 3 1 （83）以及 0 、3 、4 和 6构成偶数监督关系S3 64 3 0 （84）在发送端编码时，信息位6 、5 、4 和 3 的值决定于输入信号，因此它们是随机的。监督位 2 、 1 和 0 应根据信息位的取值按监督关系来确定，即监督位应使上三式中S1 、 S2 和S3 的值为零（表示变成的码组中应无错码）000034613562456185）由上式经移项运算，解出监督位340635164562（86）给定信息位后，可直接按上式算出监督位，其结果如表8-2 所列。表8-2 信息位监督位信息位监督位65432106543

9、210名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 17 页 - - - - - - - - - 0000 0001 0010 0011 0100 0101 0110 0111 000 011 101 110 110 101 011 000 1000 1001 1010 1011 1100 1101 1110 1111 111 100 010 001 001 010 100 111 接收端收到每个码组后，先按式（82）（ 84）计算出S1 、S2 和 S3 ，再按表8-

10、2判断错码情况。按上述方法构造的码称为汉明码。表8-2 中所列的（ 7，4）汉明码的最小码距 d0 3，因此，这种码能纠正一个错码或检测两个错码。3、两实验的系统框图如下：（一）未加信道纠错编码的4PSK调制通信系统（二）信道纠错编码（7，4）汉明码 +4PSK调制的通信系统AWGN信道输入数据信道编码(汉明码 ) 4-PSK 调制4-PSK 解调信道解码(汉明码 ) 输出数据名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 17 页 - - - - - - - - - 三

11、、实验内容（一）未加信道纠错编码的4PSK 调制通信系统1、实验要求1)根据系统组成框图用Matlab 编写程序，并在每条语句后进行注释2)最大投影点准则进行判决a, 计算噪声方差2分别为 0、0.1、0.5、1.0 时的符号差错概率和比特差错概率；b, 画出在每种2时，在检测器输入端1000 个接收到的信号加噪声的样本（星座图）；c, 分别画出数据点为1000、5000、10000、100000 时的 Monte Carlo 仿真误比特率曲线和理论误比特率曲线，比较差别，分析数据点的数量对仿真结果的影响；3)将检测器的判决准则改为最小距离法（星座图上符号间的距离），比较与上面结果的区别。4)

12、分析整个通信系统各组成部分的原理、分析实验结果，总结概括。2、实验程序（1）主程序1 N=1000;%数据点数为1000 xulie=subserial(N);%生成序列m=psktiaozhi(xulie);%4psk 调制第一步，算出码元所对应的m for j=1:N/ 2 %两路信号a(j)=cos(pi*m(j)/ 2); b(j)=sin(pi*m(j)/ 2); end; t=1; G=char(*r,*g,*b,*y);%！ ！ ！ ！ ！ ！ ！ ！ ！ ！ ！ ！for R=0,0.1,0.5,1.0 for i=1:N/2%生成噪声nc(i),ns(i)=gnguass(0

13、,sqrt(R);%生成噪声序列end; for j=1:N/ 2 %两路信号加噪声Rc(j)=a(j)+nc(j); Rs(j)=b(j)+ns(j); end; subplot(2,2,t) hold on; for i=1:N/2 plot(Rs(i),Rc(i),G(m(i)+1,:);axis(-4 4 -4 4);% 根据不同的m（i）值可以判断某一码元的来源，名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 17 页 - - - - - - - - - 从而赋

14、予其相应颜色end hold off; title( 星座图 ); s1=zuixiaojuli(Rc,Rs)%最小距离法shuchuxulie1=pskjietiao(s1);%4psk 解调errorate1juli=byterror(xulie,shuchuxulie1,t)%比特错误errorate2juli=fuhaoerror(xulie,shuchuxulie1,t)%符号错误s2=zuidatouying(Rc,Rs);% 最大投影法shuchuxulie2=pskjietiao(s2);%4psk 解调errorate1touying=byterror(xulie,shuch

15、uxulie2,t)%比特错误errorate2touying=fuhaoerror(xulie,shuchuxulie2,t)%符号错误t=t+1; end; （2）主程序2 N0=0:0.1:2; Eb=1;%平均码元能量为1 Pb=0.5*erfc(sqrt(Eb./N0);% 理论差错概率曲线公式semilogy(N0,Pb);%画出 N0 与 Pb（差错概率）的关系曲线，理论为蓝色，1000 点红色， 5000点黑色， 10000 点绿色， 100000 点黄色hold on; k=1;%控制大循环，每次循环改变点数G=char(rd-,kd-,gd-,yd-); for N=100

16、0,5000%,10000,100000 xulie=subserial(N);%生成序列m=psktiaozhi(xulie);%4psk 调制第一步，算出码元所对应的m for j=1:N/ 2 %两路信号a(j)=cos(pi*m(j)/ 2); b(j)=sin(pi*m(j)/ 2); end; t=1;%控制循环，每次循环改变N 值for RR=0:0.1:2 %不同 N 值R=sqrt(RR/2);%不同 N 值所对应的不同方差值for i=1:N/2%生成噪声nc(i),ns(i)=gnguass(0,R); 名师资料总结 - - -精品资料欢迎下载 - - - - - - -

17、 - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 17 页 - - - - - - - - - end; for j=1:N/ 2 %两路信号加噪声Rc(j)=a(j)+nc(j); Rs(j)=b(j)+ns(j); end; s=zuidatouying(Rc,Rs);%最大投影法shuchuxulie=pskjietiao(s);%4psk 解调byteerro(t)=0;% 计算比特错误，计数变量清零for i=1:N if (xulie(i)=shuchuxulie(i)% 比较输出序列与原序列，若相同无操作，若不同计数变量增一b

18、yteerro(t)=byteerro(t); else byteerro(t)=byteerro(t)+1; end; end; errorate(t)=byteerro(t)/N;% 算出误比特率t=t+1; end; RR=0:0.1:2; semilogy(RR,errorate,G(k);%画图显示误比特率与N 值的函数关系k=k+1; end; （3）高斯噪声产生子函数function gsrv1,gsrv2=gnguass(m,sgma) %本函数生成以高斯分布生成两相互正交的随机数%m 为期望， sgma 为方差if nargin=0 m=0; sgma=1; elseif n

19、argin=1 sgma=m; m=0; end%以上几句的含义是，若未输入参量，则期望为0，方差为 1，即正态分布；如输入一个参量，则该参量为sgma，期望 0；u=rand;%u 为随机产生的一个数，在01 中服从均匀分布z=sgma*sqrt(2*log10(1 /(1-u); u=rand; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 17 页 - - - - - - - - - gsrv1=m+z*cos(2*pi*u); gsrv2=m+z*sin(2*pi*u);% 通过上述四行公式将均匀分布映射成高斯分布（4）序列产生子函数function source=subserial(N) %本函数的功能是产生一个长度为N 的随机 0、1 序列%N为序列长度，source 为生成序列source=zeros(1,N); for i=1:N; temp=rand;%产生随机数，在01 上服从均匀分布if (temp序列的转换。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 17 页 - - - - - - - - -