通信原理实验QPSK及汉明码纠错.doc

  1. 1、本文档共16页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验报告 ——QPSK 通信系统的Monte Carlo 仿真 马涛、米廷振 一、实验目的 1. 提高独立学习的能力 2. 培养发现问题,解决问题,分析问题的能力 3. 学习 Matlab 的使用 4. 掌握 4PSK 通信系统的Monte Carlo 仿真方法 5. 掌握 4PSK 通信系统的组成原理 6. 比较编码信号与未编码信号在随机信道中的传输,加深对纠错编码原理的理解。 二、实验原理 1、调制解调原理: 4PSK将四进制符号映射到四种相位,利用相位传递信息。 (1)调制 经变换后可发现一个相位调制信号可以看作两个正交载波Smc和Sms的矢量和,起幅度取决于在每个信号区间内的相位。数字相位调 制信号在几何上可用Smc和Sms的二维向量来表示。 在AWGN 信道中,在一个区间内接受到的带通信号可以表示为 r(t) = um (t)+ n(t) = um (t)+ nc (t)cos(2πfct)? ns (t)sin(2πfct) =Smc cos(2πfct)-Sms sin(2πfct) + nc (t)cos(2πfct)? ns (t)sin(2πfct) =rc cos(2πfct)-rs sin(2πfct) r(t)可用二维向量[rc,rs]表示 (2)解调 有两种方法: 最大投影点准则进行判决:将接收到的信号向量r 投射到M 个可能的传输信号向量{ } sm 之一上去,并选取对应与最大投影的向量; 最小距离法:分别计算信号到星座图上符号间的距离选取最小距离对应的符号。 2、 信道纠错编码 由发送端的信道编码器在信息码元序列中增加一些监督码元。这些监督码元和信码之间有一定的关系,使接收端可以利用这种关系由信道译码器来发现或纠正可能存在的错码。 汉明码: 发送端按照规定,计算出监督位S1 S2 S3的数值。S1 S2 S3 的值最终与错码位置的对应。 接收端收到每个码组后,先按式(8-2)~(8-4)计算出S1 、S2 和S3 ,再按表 8-2 判断错码情况。按上述方法构造的码称为汉明码。表8-2 中所列的(7,4)汉明码的 最小码距d0 =3,因此,这种码能纠正一个错码或检测两个错码。 三、系统框图 (一)未加信道纠错编码的4PSK 调制通信系统 (二)信道纠错编码(7,4)汉明码+4PSK 调制的通信系统 四、实验过程 (一)未加信道纠错编码的4PSK 调制通信系统 1.实验程序: 主程序(main) clc N=100;%码元长度 sigma2=1; s=source(2*N);%产生信源bit序列,每个码元由2bit构成,故二进制序列长度为2*N [nc,ns]=nr(sigma2,N);%nr(eb,snr_db,n) subplot(2,2,1);stem(s,.);axis([0,2*N,0,1.5]); [mc,ms]=psk(s);%调制 rc=mc+nc; rs=ms+ns; subplot(2,2,2);stem(rc,rs,.); r=ipsk(rc,rs);%最大投影法解调 subplot(2,2,3);stem(r,.);axis([0,2*N,0,1.5]); [pb,pB]=pe(s,r)%pb,pB分别为误比特率和误码率 rdis=ipskdis(rc,rs);%最小距离法解调 subplot(2,2,4);stem(rdis,.);axis([0,2*N,0,1.5]); [pbdis,pBdis]=pe(s,rdis) %pbdis,pBdis分别为误比特率和误码率 子程序: (1)信源(source) function s=source(N) s=rand(1,N);%产生二进制序列 for i=1:N if s(i)=0.5 s(i)=1; else s(i)=0; end end (2)噪声(nr) function[nc,ns]=nr(sigma2,n) %噪声 sigma=sqrt(sigma2); for i=1:n u=rand; z=sigma*sqrt(2*log10(1/(1-u))); u=rand; nc(i)=z*cos(2*pi*u); ns(i)=z*sin(2*pi*u); end (3)调制(psk) function [mc,ms]=psk(s) %映射关系:由2bit二进制数转成四进制m=(0,1,2,3)带入映射公式 %mc=1,ms=0----00 %mc=0,ms=1----01 %mc=-1,ms=0----10 %mc=0,ms=-1----11 N=length(s); for i=1:2:N-1 mc((

文档评论(0)

新起点 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档