基于matlab的汉明码4FSK通信仿真实验报告

河海大学计算机及信息工程学院（常州）课程设计报告题目不同信道下汉明码4FSK系统仿真专业通信工程学号 0962310312 学生姓名程海粟指导教师高远目录一、实验目的 (3)二、实验器材 (3)三、实验内容及原理 (3)（一）汉明码编解码原理 (3)（二）4FSK调制解调原理 (6)（三）三种信道模型简介 (9)（四）程序调用函数介绍………………………………………10四、实验仿真效果图 (12)五、心得体会 (15)六、附录 (15)七、参考文献 (18)不同信道下汉明码的4FSK 系统仿真一、实验目的1、了解熟悉Matlab 仿真软件使用；2、掌握4进制频移键控（4FSK ）的调制与解调基本原理；3、掌握Matlab 仿真软件仿真4FSK 的系统设计；4、熟悉无线通信仿真过程及物理层仿真。

二、实验器材Matlab 仿真软件。

三、实验内容及原理（一）汉明码编解码原理1、编码原理一般来说，若汉明码长为n ，信息位数为k ，则监督位数r=n-k 。

若希望用r 个监督位构造出r 个监督关系式来指示一位错码的n 种可能位置，则要求21r n -≥或211rk r -≥++ (1)下面以（7，4）汉明码为例说明原理：设汉明码（n,k ）中k=4，为了纠正一位错码，由式（1）可知，要求监督位数r ≥3。

若取r=3,则n=k+r=7。

我们用6543210a a a a a a a 来表示这7个码元，用123s s s 的值表示3个监督关系式中的校正子，则123s s s 的值与错误码元位置的对应关系可以规定如表1所列。

表1 校正子和错码位置的关系则由表1可得监督关系式： 1654s a a a a =⊕⊕⊕ ()22653s a a a a=⊕⊕⊕()3 36430s a a a a =⊕⊕⊕()4在发送端编码时，信息位6543a a a a 的值决定于输入信号，因此它们是随机的。

监督位2a 、1a 、0a应根据信息位的取值按监督关系来确定，即监督位应使式（2）~式（4）中1s 、2s 、3s的值为0（表示编成的码组中应无错码）654265316430000a a a a a a a a a a a a ⊕⊕⊕=⎧⎪⊕⊕⊕=⎨⎪⊕⊕⊕=⎩ (5)式（5）经过移项运算，接触监督位265416530643a a a a a a a a a a a a=⊕⊕⎧⎪=⊕⊕⎨⎪=⊕⊕⎩ (6)式（5）其等价形式为：6543210111010001101010010110010a a a a a a a ⎡⎤⎢⎥⎢⎥⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎢⎥⎣⎦(7)式（6）还可以简记为0T T H A ∙=或0T A H ∙= (8)其中111010011010101011001H ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ []6543210A a a a a a aa =[]000=111011011011P ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦100010001r I ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦所以有[]r H PI = (9)式（6）等价于[][][]21065436543111110101011a a a a a a a a a a a Q⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎣⎦(10)其中Q 为P 的转置，即T Q P = (11)式（10）表示,信息位给定后，用信息位的行矩阵乘矩阵Q 就产生出监督位。

我们将Q 的左边加上一个k ×k 阶单位方阵，就构成一个矩阵G1000111010011000101010001011k G I Q ⎡⎤⎢⎥⎢⎥⎡⎤==⎣⎦⎢⎥⎢⎥⎣⎦(12)G 称为生成矩阵，因为由它可以产生整个码组，即有[][]65432106543a a a a a a a a a a a G =∙ (13)或者[]6543A a a a a G=∙ (14)式(13)即汉明码的编码原理2、解码和纠错原理当数字信号编码成汉明码形式（本文中即A ）后在信道中传输，由于信道中噪声的干扰，可能由于干扰引入差错，使得接收端收到错码，因此在接收端进行汉明码纠错，以提高通信系统的抗干扰能力及可靠性。

一般来说接收码组与A 不一定相同。

若设接收码组为一n 列的行矩阵B ，即[]6543210B b b b b b b b = (15)则发送码组和接收码组之差为B A E -= (16)E 就是传输中产生的错码行矩阵[]6543210E e e e e e e e = (17)若ei=0，表示接收码元无错误，若ei=1，则表示该接收码元有错。

式（16）可改写成B A E =+ (18)若E=0，即接收码组无错,则B A E A =+=，将它代人式（8），该是仍成立，即有0T B H ∙= (19)当接收码组有错时，E ≠0，将B 带入式（8）后，该式不一定成立。

在未超过检错能力时，式（19）不成立。

假设此时式（19）的右端为S,即T B H S ∙= ()20将 B A E =+代入式（20），可得()T T T S A E H A H E H =+=∙+∙由式（8）可知，所以T S E H =∙ ()21此处S 与前面的123s s s 有着一一对应关系，则S 能代表错码位置。

因此，纠错原理即，接收端收到码组后按式（20）计算出S,再根据表1判断错码情况，进行差错纠正。

（二）4FSK 调制解调原理1、调制原理二进制频移键控是数字信号调制的基本方式之一。

而多进制（MFSK ）的可降低信道系统信噪比的要求。

2FSK 信号的产生方法主要有两种：采用模拟调频电路实现；采用键控法来实现，即在二进制基带脉冲序列的控制下通过开关电路对两个不同的独立频率源进行选通，使其在每个码元期间输出f1和f2两个载波之一。

频移键控是利用载波的频率变化来传递信息的。

在2FSK中，载波的频率随二进制基带信号在f1和f2两个频率点间变化。

同理4FSK中基带脉冲序列四个码元（00 01 10 11）可用f1,f2,f3,f4四个载波之一；本文讨论4FSK是通过并联输入两位基带信号，两位二进制来表示四进制的频移键控。

2FSK键控法调频原理图如下：图 1.1(1) 2FSK调制原理4FSK可通过基带信号（00,01,10,11）并联传输0或1来分别用f1,f2,f3,f4四个载频表示，两路基带信号作为控制选通选通开关，1路选通开关发送0时选通载频f1, 发送0时选通载频f2, 1路选通开关发送0时选通载频f3, 送1时选通载频f4。

两路不同载频通过相加器得到已调信号发送出去。

4FSK键控法调频原理图（图1.1）如下：图 3.1.1(2) 4FSK调制原理2、解调原理4FSK信号的相干解调法原理框图如图1.2所示。

其原理是：4FSK信号先经过带通滤波器去除调制信号频带以外的在信道中混入的噪声，此后该信号分为四路，每路信号与相应载波相乘，再经过低通滤波器去除高频成分，得到包含基带信号的低频信号，将其送入抽样判决器中进行抽样判决，抽样判决器的输出分别得到两路原基带信号表示四进制得到原始码元。

图 1.2 4FSK解调原理（三）三种信道模型简介1、瑞利信道瑞利衰落信道（Rayleigh fading channel）是一种无线电信号传播环境的统计模型。

这种模型假设信号通过无线信道之后，其信号幅度是随机的，即“衰落”，并且其包络服从瑞利分布。

这一信道模型能够描述由电离层和对流层反射的短波信道，以及建筑物密集的城市环境。

12瑞利衰落只适用于从发射机到接收机不存在直射信号（LoS，Line of Sight）的情况，否则应使用莱斯衰落信道作为信道模型。

瑞利分布是一个均值为0，方差为σ2的平稳窄带高斯过程，其包络的一维分布是瑞利分布。

其表达式及概率密度如图所示。

瑞利分布函数瑞利分布是最常见的用于描述平坦衰落信号接收包络或独立多径分量接受包络统计时变特性的一种分布类型。

两个正交高斯噪声信号之和的包络服从瑞利分布。

瑞利衰落能有效描述存在能够大量散射无线电信号的障碍物的无线传播环境。

若传播环境中存在足够多的散射，则冲激信号到达接收机后表现为大量统计独立的随机变量的叠加，根据中心极限定理，则这一无线信道的冲激响应将是一个高斯过程。

如果这一散射信道中不存在主要的信号分量，通常这一条件是指不存在直射信号（LoS），则这一过程的均值为0，且相位服从0 到2π的均匀分布。

即，信道响应的能量或包络服从瑞利分布。

若信道中存在一主要分量，例如直射信号（LoS），则信道响应的包络服从莱斯分布，对应的信道模型为莱斯衰落信道。

通常将信道增益以等效基带信号表示，即用一复数表示信道的幅度和相位特性。

由此瑞利衰落即可由这一复数表示，它的实部和虚部服从于零均值的独立同分布高斯过程。

2、多径信道一种扩频通信系统的多径解扩方法：将来自n个信道的数据分别经过模/数转换电路进行量化，然后把经过量化后的数据保存在各对应的数据存储器中，之后根据多径的峰值和多径延时的值，从各数据存储器的不同位置依序读出n个信道的多径数据后，进行该多径解扩：首先是预解扩；其次是二次解扩累加；再者是多次解扩累加；最后输出n个信道的多径解扩数据。

3、高斯信道高斯信道是一个射频通信信道，其包含了各种频率的特定噪声频谱密度的的特征，从而导致了信道中错误的任意分布。

常指加权高斯白噪声(AWGN)信道。

这种噪声假设为在整个信道带宽下功率谱密度(PDF)为常数，并且振幅符合高斯概率分布。

信号经过高斯信道会受到加性干扰的影响。

（四）程序调用函数介绍MATLAB中提供了汉明码的编解码、4FSK调制解调函数和不同通信信道模型函数，本程序直接调用进行编程。

1、encode函数功能：编码函数语法：code=encode(msg,N,K,’hamming’)说明：该函数对二进制信息msg进行汉明编码，K为信息位长度，N为码字长度。

msg是一个K列矩阵。

2、decode函数功能：译码函数语法：rcvcode=decode(code,N,K,’hamming’)说明：该函数对接受码字进行译码，恢复出原始信息，译码参数及方式必须和编码时采用的完全相同。

3、dmod函数功能：调制函数语法：fskout=dmod(signal,Fc,Fd,Fs,'fsk',M,df)说明：该函数可以对输入的数字信号进行4FSK调制，Fs为系统的采样频率，Fc 为载波频率，Fd为码元速率，df为频差。

4、ddemod函数功能：调制函数语法：y=ddemod(signal,Fc,Fd,Fs,'fsk',M,df)说明：该函数可以对输入的4FSK信号进行解调，Fs为系统的采样频率，Fc为载波频率，Fd为码元速率，df为频差。