基于MATLAB的调制解调及信道编译码仿真_第1页
基于MATLAB的调制解调及信道编译码仿真_第2页
基于MATLAB的调制解调及信道编译码仿真_第3页
基于MATLAB的调制解调及信道编译码仿真_第4页
基于MATLAB的调制解调及信道编译码仿真_第5页
已阅读5页,还剩13页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、基于MATLAB的调制解调与信道编译码仿真摘要:随着信息时代的步伐,通信技术得到了全面的发展,信息技术已成为了21世纪最强大的国际化动力。在通信技术中,信息的调制、解调和误码纠错都占有重要的地位。MATLAB作为一款功能强大的数学工具软件,在通信领域中得到了很广泛的应用。本文基于MATLAB对信号进行模拟仿真设计,实现对二进制相移键控、循环码的纠错仿真、BPSK的调制解调等进行仿真设计。关键字:MATLAB、调制解调、2PSK、BPSK、重复码。一 、二进制和四进制相移键控调制仿真设计1.1 二进制相移键控(2PSK)原理 相移键控是利用载波的相位变化来传递数字信息的,而振幅和频率保持不变。在

2、2PSK中常用0和分别表示二进制“0”和“1”,2PSK的信号时域表达式为: 表示为第n个符号的绝对相位,因此上式可改写为: 由于表示信号的两种码形完全相同,极性相反,故2PSK信号一般表示一个双极性全占空矩形脉冲序列与一个正弦载波相乘。=s(t)coswct其中: s(t)=g(t-nTs);这里,g(t)为脉宽Ts的单个矩形脉冲;的统计特征为 =1, & 概率为P-1, & 概率为1-P即发送二进制“0”时(),取0相位;发送二进制符号“1”时(),取相位。这种以载波的不同相位直接法去表示相应二进制数字信号的调制方式,称为二进制绝对相移方式。 2PSK信号时间波形2PSK信

3、号的调制原理如下图所示,与2ASK的产生方法相比,只是对s(t)的要求不同,在2ASK中s(t)是单极性的,而在2PSK中s(t)是双极性的基带信号。(a) 模拟调制方法(b) 键控法 2PSK信号的解调通常采用相干解调法,解调原理如下原理框图所示,在相干解调中,如何得到与接受的2PSK信号同频同相的相干载是关键问题,后续进一步介绍。 2PSK信号的调解原理框图 1.2 二进制差分相移键控2PSK相干调解时,由于载波恢复中相位有0、模糊性,导致调解过程出现”反相工作”现象,恢复出的信号“1”和“0”倒置,从而使2PSK难以实际应用。二进制差分相移键控(2DPSK)能够克服次缺点。2DPSK是利

4、用前后码元的载波相位变化传递数字信息,假设为当前码元与前一码元的载波相位差,定义为一种数字信息与之间的关系为 =0 表示数字信息“0”1 表示数字信息“1”与是可将一组二进制数字信息与其对应的2DPSK信号的载波相位关系示例如下:二进制数字信息: 1 1 0 1 0 0 1 1 02DPSK信号或相位:(0) 0 0 0 或 :()0 0 0 0 0 0相应的2DPSK的典型波形如下图: =0 表示数字信息“1” 表示数字信息“0”(a)绝对码 1 1 0 1 0(c)2DPSK(b)相对码 0 1 0 0 1 1 对于相同的基带数字信息序列,由于序列初始码元的参考相位不同,2DPSK信号的位

5、置可以不同。2DPSK信号的相位并不直接代表基带信号,而前后码元相对相位的差才唯一决定信息的符号。 1.3 四进制幅度调制原理振幅键控时利用载波的振幅变化来传递数字信息,而其频率和初始相位保持不变。4ASK信号的一般表达式为e4ASK=nang(t-nTB)cosct如果其中是0的四进制信号,只要让载波信号与多进制信号通过乘法器即可调制完成。如果两个电平都不是0,只要让载波信号的振幅固定,通过乘法器与多进制信号相乘。而解调的两种方法包括非相干解调和相干解调,其原理图如下:三1.3 程序流程图:2PSK调制流程图二、二进制对称信道下的重复码编译码及错误概率统计仿真设计2.1 循环码的编码循环码是

6、采用循环移位的特性界定的一种线性分组码,除了具有线性分组码的一般特性之外,还具有循环性,是一种无权码。循环码的最大特点就是码字循环特性,所谓循环特性是指循环码中任何一组循环码循环移位后,所得到的码组仍是循环码。若(an-1 an-2 ,a1an)为一组循环码组 则(an-2 an-3 ,a0 an-1)、(an-3 an-4 ,an-1 an-2)、,还是循环码组,也就是说,不论是左移还是右移,也不论移多少位,仍然循环的循环码组。例如(3,1)重复码两个码字是(0,0,0)和(1,1,1)。 编码 for j=1:L for k=1:n code=code;data(j,1);end end其

7、中L为信息序列的码长,n为重复码字的长度。这种编码方式时间开销较大。2.2循环码的译码 decode=reshape(demodcode,n,L); decode=decode; weight=sun(decode); weight=weight;weight(weight<n/2)=0;weight(weight>n/2);dedata=weight;2.3循环码的纠错当传输信息比特b=1,得到码字c=b*G=(1,1,1);若信道中存在噪声干扰,导致接收端收到一个错误的码元,如v=(1,0,1) 此时根据多判决法可以判断发送的信息比特为1。同理,若b=0,收到码字v=(1,0,

8、0),可以判断出发送的信息bite为0。通过这样的方式实现纠错。有纠错定理:若分组码有最小距离d,那么该码的最大检错数tmd和最大纠错数tmc满足:(1) 在检错模式时,有tmd=d-1(2) 在纠错模式下,tmc=(d-1)/2(3) 在混合检错模式下有 Tmc+tmdd-1并同时又tmc<tmd 2.4 程序流程图如下:续下页循环码的编码译码流程图2.5 程序脚本语言:n=1;N=100000;for p=0.05:0.01:0.1 a=randi(0 1,1,N); A=bsc(a,p); B=bsc(a,p); C=bsc(a,p); D=bsc(a,p); for i=1:N

9、X(i)=A(i)+B(i)+C(i)+D(i); end for i=1:100000 if X(i)<=2 e(i)=0; else e(i)=1; end end K=find(e=a); E(n)=length(K)/N; n=n+1;endp=0.05:0.01:0.1;plot(p,E);2.6 调试结果 :三、BPSK的调制与解调3.1 BPSK 简介:均匀分布在给定频带上的高斯噪声,概率分布正态分布。高斯噪声是一种随机噪声,在任选瞬时中任取n个,其值按n个变数的高斯概率定律分布。如果一个噪声,它的幅度分布服从高斯分布,而它的功率谱密度又是均匀分布的,则称它为高斯白噪声。如

10、热噪声,散粒噪声等都是高斯白噪声。 加性高斯白噪声(AWGN)从统计上而言是随机无线噪声,其特点是其通信信道上的信号分布在很宽的频带范围内。高斯白噪声的概念:“白”指功率谱恒定;高斯指幅度取各种值时的概率p(x)是高斯函数。加性高斯白噪声在通信领域中指的是一种各频谱分量服从均匀分布(即白噪声),且幅度服从高斯分布的噪声信号。因其可加性、幅度服从高斯分布且为白噪声的一种而得名。该噪声信号为一种便于分析的理想噪声信号,实际的噪声信号往往只在某一频段内可以用高斯白噪声的特性来进行近似处理。3.2程序设计流程图: 3.3程序脚本语言 :用AWGN产生高斯白噪声   clear,cl

11、c;   N=0:1000;   fs=1024;   t=N./fs;      y=3*sin(2*pi*t);   x=wgn(1,1001,2);   i=y+x;   subplot(3,1,1),plot(x);   subplot(3,1,2),plot(y);   subplot(3,1,3),plot(i);3.4调试结果 高斯白噪声 运行结果四、重复码和BPS K调制解调4.1

12、译码简介:译码器分两种,一种是检错译码,一种是纠错译码,对于检错译码,记译码器输出为y,则译码器输出y=(r,s),r为当前接受向量,s是r是否具有差错的标志;对于纠错译码,其工作区又分为译码成功和译码失败两种状态。译码器能够在达到译码码字差错概率最小的条件下输出一个确切的码字C,即y=C。译码失败是指译码器不能输出一个正确的码字C,通常此时的译码器输出y与检测译码相同。伴随式纠错译码是通用译码方法。(1) 按最可能出现的2r个差错图案e,计算相应的伴随式s,并构建伴随式-差错检测图案表s,e;(2) 对接受向量r计算伴随式s(3) 查(s,e)表的e;(4) 纠错计算C=r-e;一个二源数据

13、传输系统,在不同观察统计点上的误码率是不同的,尽管一个码字的差错检测并不一定导致全部消息的差错,但是对大多数纠错码,在统计上总可以认为消息与码字之间的映射,不改变码字差错导致在整个码长内比特差错的均匀分布特性,这样在统计意义上有pbmpbd若考虑码字是等概率发送的,则错误译码的误码率pbe为 pbe=1nj=dnjpwe(j)其中,pwe(j)是发送全0码字并错为j重码字的概率。另一方面,一个码字差错冰然有至少2t+1位码字比特错,每个码字平均有(2t+1)k/n位消息比特错,所以pbe与pwe有如下渐进关系kpbe(2t+1)knpwe对于译码失败造成的误码率为pbf,有1ni=t+1nni

14、-Diipi1-pn-i因此译码后总的误码率为pbd, pbd=pbe+pbf 二元数据传传输系统 Pbc是信道的比特差错概率。4.2 流程图 4.3 程序代码 clear all ;clc;t=0:0.01:1; L=length(t);N=10000; ES_NO=2; s=sin(2*pi*t);A=randi(1,N);b=zeros(1,5);C=zeros(1,N);for k=1:N for i=1:5 b(i)=A(k); end x=b; y=1-2*x; z=; for i=1:5 if y(i)=1 z=z s; else z=z -s; end end r=awgn(z,ES_NO-10*log10(L/2)-10*log10(5),'measured'); for u=1:5 a=r(u*L-(L-1): u*L).*s; d(u) = sum(a); end Y=sign(d); for j=1:5 if Y(j)=-1 C(k)=C(k)+1; end end if C(k)>=3 C(k)=1; else C(k)=0; end endnumber=length(find(A=C

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论