通信系统MATLAB仿真资料Word文档格式.docx

通信系统MATLAB仿真资料Word文档格式.docx

《通信系统MATLAB仿真资料Word文档格式.docx》由会员分享，可在线阅读，更多相关《通信系统MATLAB仿真资料Word文档格式.docx（42页珍藏版）》请在冰豆网上搜索。

2.3.2.编码原理：

10

2.3.3.译码原理：

11

2.4.PSK调制与解调原理13

2.4.1.psk调制原理13

2.4.2.psk解调原理14

3.各部分程序的MATLAB仿真16

3.1.PCM编译码的MATLAB仿真16

3.2.HDB3码的MATLAB程序仿真18

3.3.循环码编译码的MATLAB实现20

3.4.PSK调制与编码的MATLAB仿真22

4.整体通信过程的MATLAB仿真24

5.结果分析35

6.小结体会36

7.参考文献37

摘要

本文是基于PCM编码的PSK信号的MATLAB仿真。

本文先分析了PSK系统的基本调制解调方法，用PCM编码将模拟信号转换为离散的数字信号，在基带中采用HDB3码编码，信道预编码采用循环码，在衰落信道中采用PSK调制技术，然后用MATLAB来仿真通信系统的整个过程。

通过仿真，观察了调制解调过程中各个环节时域和频域的波形，并结合调制原理，跟踪分析比较了各种调制方法的性能，并通过比较仿真与理论计算的性能，证明了仿真的可靠性。

通过这次的课程设计让我们加深了对课本知识的深入理解，而且也让我们更加熟练地运用MATLAB仿真软件。

关键词：

PCM编码，HDB3，循环码，PSK，MATLAB

Abstract

Asweallknown,PSKmodulationisthemostimportantpartinmoderncommunicationsystem.Also,thetechnologyofimprovingPSKmodulationisthemostessensialapproachtoimprovethefunctionofmoderncommunicationsystem.ThispaperisbasedonthePCMcoding,HDB3coding,

CycliccodingandMATLABtoahievePSK.Inthispaper,themethodofPSKmodulationareintroducedfirstly，andPCMcodingisusedtotransfertheanalogsignaltodigitalsignal.ThentheM-documentinMATLABisusedtosimulate.Throughobservingtheresultsofsimulation,thefactorsthataffectthecapabilityofthePSKandDPSKmodulationsystemandthereliabilityofthesimulationmodelsareanalyzed.Bythiscurriculumdesign,wenotonlymakeacomprehensiveunderstandingoftheknowledgeinthebooks,butalsobeingskilledinusingtheMATLABsimulationsoftware.

Keywords:

PCMcoding,HDB3coding,Cyclecoding,PSK,MATLAB

1.通信系统整体概述

通信系统通常由信源、信道和信宿构成，其中包括模拟信号的数字化，基带信号的编码与解码，信道编码与解码，数字信号的调制与解调。

一个典型的通信系统如图1.1所示：

本课设要求用MATLAB设计一简单的通信系统，要求采用PCM编译码，基带信号用HDB3码，信道编译码采用循环码，调制方式为PSK，信道为衰落信道。

2.各部分原理简述

2.1.脉冲编码调制（PCM）原理

脉冲编码调制就是把一个时间连续，取值连续的模拟信号进行取样，变换成时间离散的信号后，在进行量化、编码，变换成时间离散、取值离散的数字信号的过程。

脉冲编码调制包括抽样、量化、编码三个过程。

（1）抽样

模拟信号通常是在时间上连续的信号，一在系列离散的点上，对这种信号抽取样值称为抽样。

在理论上，抽样过程可以看作是用周期性单位冲激脉冲和此模拟信号相乘。

抽样结果得到的是一系列周期性的冲激脉冲，其面积和模拟信号的取值成正比。

由抽样定理，只要抽样脉冲的频率大于或等于信号最高频率的2倍时，信号就能被恢复出来。

所谓抽样就是对信号进行周期性的扫描，把时间上连续的信号变为时间上离散的信号。

抽样前后的信号波形图如图2.1所示：

图2.1抽样前后的信号波形图

但由于实际的电路不可能产生理想的冲激脉冲，故抽样的结果是平顶脉冲，如图2.2所示。

图2.2抽样示意图

所谓量化就是把抽样得到的瞬时值将其幅度离散，即用一组规定的电平，把瞬时抽样值用最接近的电平值表示。

量化分为均匀量化和非均匀量化。

均匀量化时，量化间隔是确定的，量化器的平均输出信号量噪比随量化电平数M的增大而提高，故这种均匀量化器对于小输入信号很不利。

在非均匀量化时，量化间隔是随信号抽样值的不同而变化的。

信号抽样值小时，量化间隔也小，信号抽样值大时，量化间隔也变大。

故在实际应用中常采用非均匀量化。

设模拟信号的抽样值为m（KT），其中T是抽样周期，k为整数。

此抽样值仍然是一个取值连续的变量。

若仅用N个二进制数字码元来代表此抽样值的大小，则N个二进制码元只能代表M=

个不同的抽样值。

这样，共有M个离散电平，它们称为量化电平。

用这M个量化电平表示连续抽样值的方法称为量化。

在原理上，量化过程可以认为是在一个量化器中完成的。

在实际中，量化过程常是和后续的编码过程结合在一起完成的，不一定存在量化器。

非均匀量化的实现方法通常是在进行量化之前，先将信号抽样值压缩再进行均匀量化。

通常使用的压缩器中，大多采用对数式压缩。

广泛采用的两种对数压缩律是

压缩律和A压缩律。

（3）编码

所谓编码就是用一组二进制码组来表示每一个固定电平的量化值。

然而，实际上量化是在编码过程中完成的，故编码过程也称为模数变换（A/D）。

最常用的编码是用二进制的符号，例如“0”和“1”，表示此离散数值。

在用电路实现时，最常用的是一种方案称为逐次比较法编码，其基本原理方框图如图2.3所示。

此图示出的是一个3位编码器。

编码器的输入信号抽样脉冲值在0和7.5之间。

它将输入的信号模拟抽样脉冲编成3位二进制编码c1c2c3.

图2.3编码电路

如图2.3所示，输入信号抽样脉冲电流Is由保持电路短时间保持，并各几个称为权值电流的标准电流Iw逐次比较。

每比较一次，得出1位二进制码。

权值电流Iw是在电路中预先产生的。

Iw的个数决定于编码的位数，现在共有三个不同的Iw值。

因为表示量化值的二进制码有3位，即c1c2c3。

它们能够表示8个十进制数，从0到7，如表2.1所列。

表2.1编码表

量化值

c1

c2

c3

4

5

6

7

在15折线法中采用的折叠码有9位。

其中第一位c1表示量化值的极性正负。

后面的8位分为段落码和段内码两部分，用于表示量化值的绝对值。

其中第2到4位（c2c3c4）是段落码，共计3位，可以表示8种斜率的段落；

其他5位（c5~c9）为段内码，可以表示每一段落内的32种量化电平。

段内码代表的32个量化电平是均匀划分的，所以，这8位共能表示

种量化值。

在进行PCM编码时，我国大陆、欧洲各国以及国际间互连时采用A压缩律及相应的13折线法，北美、日本和韩国等少数国家和地区采用u压缩律及15折线法。

A律压缩先将原信号进行非均匀量化，再进行均匀量化。

对小信号进行细致编码，对大信号进行粗略编码，结果是小信号的量噪比高，大信号的量噪比低。

PCM码由极性码、段落码、段内码组成。

极性码位于最高位，若抽样值为正，则极性码为“１”，若抽样值为负，则极性码为“０”。

段落码对应于A律编码非均匀量化区间，其编码方法见表2.2。

　　　　　　　　　　　　　　表　2.2　段落码的编码

段落序号

段落码（C２C３C４）

量化范围（量化单位）

８

１１１

１０２４—２０４８

７

１１０

５１２——１０２４

６

１０１

２５６———５１２

５

１００

１２８———２５６

４

０１１

６４————１２８

３

０１０

３２—————６４

２

００１

１６—————３２

１

０００

０——————１６

段内码就是对段落码的每个区间在进行均匀量化，若采用自然二进制码的编码方式，编码方法见表2.3。

　　　　　　　　　　　　　　表　2.3　段内码的编码

量化间隔

段内码（Ｃ５Ｃ６Ｃ７Ｃ８）

段内码（C５C６C７C８）

１５

１１１１

０１１１

１４

１１１０

０１１０

１３

１１０１

０１０１

１２

１１００

０１００

１１

１０１１

００１１

１０

１０１０

００１０

９

１００１

０００１

１０００

０

００００

由此可见，脉冲编码调制方式就是一种传递模拟信号的数字通信方式。

PCM编码的框图如图2.4所示：

图2.4PCM编码的系统框图

将PCM脉冲流解调为原模拟信号的过程即为PCM译码过程，PCM译码过程是PCM编码过程的逆过程。

其过程为：

（1）按照编码表将二进制的脉冲流还原为对应的十进制数值。

（2）对照量化表将

（1）中的十进制的数值转化为对应的量化电平值。

（3）将

（2）中得到的平顶脉冲通过一低通滤波器，得到还原的模拟信号。

2.2.HDB3码原理

HDB3（3ndOrderHighDensityBipolar）码的全称是三阶高密度双极性码。

它是AMI码的一种改进型，改进的目的是为了保持AMI码的优点而克服其缺点，使连零的个数不超过三个。

下面先对AMI码作简要介绍。

AMI（AlternativeMarkInversion）码的全称是传号交替反转码，其编码的规则是将消息码的“1”（传号）交替的变化为“+1”和“-1”，而“0”（空号）保持不变。

AMI码的优点是没有直流成分，且高低频分量少，能量集中频率为1/2的码速处，其编译码电路简单。

但当信源中出现长串的“0”符号时，信号的电平长时间不跳变，造成提取定时信号的困难。

解决长串“0”码问题的有效方法之一是采用HDB3码。

HDB3码的编码规则是：

（1）检查消息码中“0”的个数。

当“0”的数目小于等于3时，HDB3码与AMI码一样，+1与-1交替。

（2）当连“0”数目超过3时，将每4个连“0”化作一小节，定义为B00V，称为破坏节，其中V称为破坏脉冲，B称为调节脉冲。

（3）V与前一个相邻的非“0”脉冲的极性相同（这就破坏了极性交替的规则），且相邻的V码之间的极性必须交替。

V的取值为+1或-1。

（4）B的取值可选0，+1，-1，以使V同时满足（3）中要求。

（5）V码后面的传号极性也要交替。

例如：

消息码：

110000000011

AMI码：

-1100000000-11

HDB3码：

-11-B00-V+B00+V-11

HDB3码的编码虽然比较复杂，但译码却比较简单。

从上述的编码规则可以看出，每一个破坏脉冲V总是与前一个非“0”脉冲同极性（包括B在内）。

这就是说，从收到的符号序列中可以容易的找到破坏点V，于是断定V符号及其前面的三个符号必须是连“0”符号，从而恢复出四个连“0”码，再将所有的-1变成+1后便得到原消息码。

HDB3码除了具有AMI码的优点外，同时还将连“0”码限制在三个以内，使得接收时能保证定时信息的提取。

因此，HDB3码是目前应用最为广泛的码型，A律PCM四次群以下的接口码型均为HDB3码。

2.3.循环码编码原理

2.3.1．RS循环码编译码原理与特点

设C使某（n,k）线性分组码的码字集合，如果对任

，它的循环移位

也属于C，则称该（n,k）

码为循环码。

该码在结构上有另外的限制，即一个码字任意循环移位的结果仍是一个有效码字。

其特点是：

（1）可以用反馈移位寄存器很容易实现编码和伴随式的计算；

（2）由于循环码有很多固有的代数结构，从而可以找到各种简单使用的译码办法。

如果一个（n,k）线性码具有以下的属性，则称为循环码：

如果n元组

是子空间S的一个码字，则经过循环移位得到的

也同样是S中的一个码字；

或者，一般来说，经过j次循环移位后得到的

也是S中的一个码字。

RS码的编码系统是建立在比特组基础上的，即字节，而不是单个的0和1，因此它是非二进制BCH码，这使得它处理突发错误的能力特别强。

码长：

信息段：

（t为纠错符号数）

监督段：

最小码段：

最小距离为d的本原RS码的生成多项式为：

g（x）=（x-α）（x-α2）（x-α3）…（x-αd-2）

信息元多项式为：

：

m（x）=m0+m1x+m2x2+…+mk-1xk-1

循环码特点有：

（1）循环码是线性分组码的一种，所以它具有线性分组的码的一般特性，且具有循环性，纠错能力强。

（2）循环码是一种无权码，循环码编排的特点为相邻的两个数码之间符合卡诺中的邻接条件，即相邻数码间只有一位码元不同，因此它具有一个很好的优点是它满足邻接条件，没有瞬时错误（在数码变换过程中，在速度上会有快有慢，中间经过其他一些数码形式，即为瞬时错误）。

（3）码字的循环特性，循环码中任一许用码经过牡环移位后，所得到的码组仍然是许用码组。

循环码中的生产多项式g（x）的次数为n-k次，且g（x）的最低位不为“0”。

由信息码构成信息多项式

，其中最高幂次为k-1；

用

乘以信息多项式m（x），得到的

，最高幂次为n-1，该过程相当于把信息码（

，

，……，

）移位到了码字德前k个信息位，其后是r个全为零的监督位；

用g（x）除

得到余式r（x）,其次数必小于g（x）的次数，即小于（n-k），将此r（x）加于信息位后做监督位，即将r（x）于

相加，得到的多项式必为一码多项式。

综上所述，循环码的编码步骤为：

（1）有信息码构成信息多项式m（x）=mk-1xk-1+``````m0，其中高幂次为k-1。

（2）用xn-k乘上信息多项式m（x），得最高幂次为n-1，做移位。

（3）用g（x）除xn-km（x）得到余式r（x）。

（4）用xn-km（x）加余式r（x）得到循环码。

编码过程流程图如图2.5所示：

图2.5循环码编码流程图

由接收到的y（x）计算伴了随式s（x）。

根据伴随式s（x）找出对应的估值错误图样。

计算c^（x）=y（x）+e^（x），得估计码字。

若c^（x）=c（x），则译码正确，否则错误。

由于g（x）的次数为n-k次，g（x）除E（x）后得余式（即伴随式）的最高次数为n-k-1次，故S（x）共有2n-k个可能的表达式，每一个表达式对应一个错误格式。

可以知道（7,4）循环码的S（x）共有2（7-4）=8个可能的表达式，可根据错误图样表来纠正（7,4）循环码中的一位错误。

解码过程流程图如图2.6示：

否

存储c（x）

图2.6循环码解码流程图

2.4.PSK调制与解调原理

2.4.1.psk调制原理

2PSK信号用载波相位的变化来表征被传输信息的状态，通常规定0相位载波和π相位载波分别表示传“1”和传“0”。

设二进制单极性码为an,其对应的双极性二进制码为bn,则2PSK信号的一般时域信号可以表示为：

S2psk（t）=[

bng（t-nTs）]cosωct

式中

bn=-1（当an=0时，概率为P）

bn=1（当an=1时，概率为1-P）

则时域信号可以变为:

S2psk（t）=[

g（t-nTs）]cos（ωct+π）当an=0时

g（t-nTs）]cos（ωct+0）当an=1时

由此可知2PSK信号是一种双边带信号,2PSK信号的带宽为:

B2PSK=（ƒc+Rs）-（ƒc-Rs）=2Rs

式中Rs为码元速率。

值得注意的是，2PSK码元序列的波形与载频和码元持续时间之间的关系有关。

当一个码元中包含有整数个载波周期时，在相邻码元的边界处波形是不连续的，或者说相位是不连续的。

当一个码元中包含的载波周期数比整数个周期多半个周期时，则相位连续。

当载波的初始相位差90度时，即余弦波改为正弦波时，结果类似。

以上说明，相邻码元的相位是否连续与相邻码元的初始相位是否相同不可混为一谈。

只有当一个码元中包含有整数个载波周期时，相邻码元边界处的相位跳变才是由调制引起的相位变化[16]。

2PSK信号的产生方法主要有两种。

第一种叫相乘法，是用二进制基带不归零矩形脉冲信号与载波相乘，得到相位反相的两种码元。

第二种方法叫选择法，是用此基带信号控制一个开关电路，以选择输入信号，开关电路的输入信号是相位相差

的同频载波。

这两种方法的复杂程度差不多，并且都可以用数字信号处理器实现。

2PSK信号的调制原理框图如图2.7所示

图2.72PSK信号的调制原理框图

2.4.2.psk解调原理

2PSK信号的解调方法是相干解调法。

由于PSK信号本身就是利用相位传递信息的，所以在接收端必须利用信号的相位信息来解调信号。

下图2-3中给出了一种2PSK信号相干接收设备的原理框图。

图中经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。

判决器是按极性来判决的。

即正抽样值判为1，负抽样值判为0.

2PSK信号相干解调各点时间波形如图2.4所示.当恢复的相干载波产生180°

倒相时,解调出的数字基带信号将与发送的数字基带信号正好是相反,解调器输出数字基带信号全部出错.

图2.8相干解调各点波形

这种现象通常称为"

倒π"

现象.由于在2PSK信号的载波恢复过程中存在着180°

的相位模糊,所以2PSK信号的相干解调存在随机的"

现象,从而使得2PSK方式在实际中很少采用.

2PSK信号的调制原理框图如图2.9所示：

图2.92PSK信号的解调原理框图

说明：

由于PSK信号的功率谱中无载波分量，所以必须采用相干解调的方式。

在相干解调中，如何得到同频同相的本地载波是个关键问题。

只有对PSK信号进行非线性变换，才能产生载波分量。

2PSK信号经过带通滤波器得到有用信号，经相乘器与本地载波相乘再经过低通滤波器得到低频信号v（t），再经抽样判决得到基带信号。

3.各部分程序的MATLAB仿真

3.1.PCM编译码的MATLAB仿真

用MATLAB进行PCM编码的程序如下所示：

clear;

clc;

N=8;

%编码位数为8位

dt=1/8;

%抽样时间间隔

endtime=4;

t=0:

dt:

endtime;

%时间范围设定

%抽样

fori=1:

endtime/dt

x（i）=sin（dt*i*pi）;

end

a=-1;

b=1;

%量化范围设定

dv=（b-a）/2^7;

%量化间隔

128

m（i）=a+i*dv;

%量化值

128-1;

q（i）=（m（i）+m（i+1））/2;

%量化中值

minm=min（m）;

%量化

endtime/dt;

forj=1:

128-1

if（（x（i）>

=m（j））&

&

（x（i）<

=m（j+1）））

xq（i）=q（j）;

xc（i）=j-64+1;

elseifx（i）==0

xq（i）=0;

xc（i）=0;

elseifx（i）<

minm

xq（i）=minm;

xc（i）=-64+1;

end

stem（xc）;

gridon;

ifxc（i）<

=0%负极性最高位为0

xcc（i）=abs（xc（i））;

elseifxc（i）>

0%正极性最高位为1

xcc（i）=（xc（i））+128;

xpcm=dec2bin（xcc,8）;

%将量化值转化为0，1代码

抽样图如图3.1示：

图3.1原信号的抽样图样

程序运行后得到的PCM序列为：

>

xpcm

xpcm=

10011001

10101110

10111100

11000000

10000001

00011000

00101101

00111011

00111111

00000000

3.2.HDB3码的MATLAB程序仿真

%HDB3码的编译码过程

%输入为xpcm

32

ff（8*（i-1）+1:

8*i）=xpcm（i,:

）;

end%将得到的PCM码转换为一维序列

AMI=[];

%定义AMI码

k=1;

%将PCM码转化为AMI码

endtime/dt*8

ifff（i）==dec2bin（