基于matlab的信道编码分析

1、-. z . I . . r . .题目：基于MATLAB的通信系统仿真信道编码对通信系统性能的影响专业： 通信工程：*校通*：201300484316目录 TOC o 1-3 h z u HYPERLINK l _Toc469583057一、引言 PAGEREF _Toc469583057 h 2HYPERLINK l _Toc469583058二、信道编码理论 PAGEREF _Toc469583058 h 2HYPERLINK l _Toc4695830592.1、信道编码的目的 PAGEREF _Toc469583059 h 2HYPERLINK l _Toc4695830602.2、

2、信道编码的实质 PAGEREF _Toc469583060 h 3HYPERLINK l _Toc4695830612.3、信道编码公式 PAGEREF _Toc469583061 h 3HYPERLINK l _Toc469583062三、线性分组码的编译码原理 PAGEREF _Toc469583062 h 3HYPERLINK l _Toc4695830633.1、线性分组码的根本概念 PAGEREF _Toc469583063 h 3HYPERLINK l _Toc4695830643.2、生成矩阵和校验矩阵 PAGEREF _Toc469583064 h 4HYPERLINK l _

3、Toc469583065四、MATLAB仿真 PAGEREF _Toc469583065 h 5HYPERLINK l _Toc4695830664.1仿真 PAGEREF _Toc469583066 h 5HYPERLINK l _Toc469583067原理说明 PAGEREF _Toc469583067 h 5HYPERLINK l _Toc469583068各子函数说明 PAGEREF _Toc469583068 h 5HYPERLINK l _Toc4695830694.2仿真源程序 PAGEREF _Toc469583069 h 5HYPERLINK l _Toc469583070

4、信道编码 PAGEREF _Toc469583070 h 5HYPERLINK l _Toc469583071信道解码 PAGEREF _Toc469583071 h 6HYPERLINK l _Toc469583072交织 PAGEREF _Toc469583072 h 6HYPERLINK l _Toc469583073解交织 PAGEREF _Toc469583073 h 7HYPERLINK l _Toc469583074信道衰落 PAGEREF _Toc469583074 h 7HYPERLINK l _Toc469583075六程序及仿真图 PAGEREF _Toc46958307

5、5 h 8HYPERLINK l _Toc4695830761、 file1：信道编码对通信系统性能的影响，有无信道编码的影响 PAGEREF _Toc469583076 h 8HYPERLINK l _Toc4695830772、file2：在周期性深衰落的信道条件下，交织对通信系统性能的影响 PAGEREF _Toc469583077 h 10HYPERLINK l _Toc4695830783、 file3：在交织条件下，不同时长的周期性深衰落对系统性能影响的比拟 PAGEREF _Toc469583078 h 13基于MATLAB的通信系统仿真信道编码对通信系统性能的影响摘要：简述信道

6、编码理论，详细说明分组码的编译原理、实现方法及检错纠错能力，用MATLAB仿真有无信道编码条件下对通信系统性能的影响及信道编码在不同信道下对通信系统性能的影响，如AWGN信道和深衰落信道。关键词：信道编码、分组码、MATLAB仿真、性能一、引言提高信息传输的有效性和可靠性始终是通信技术所追求的目标，而信道编码能够显著的提升信息传输的可靠性。1948年,信息论的奠基人C.E.Shannon在他的开创性论文通信的数学理论中,提出了著名的有噪信道编码定理.他指出：对任何信道,只要信息传输速率R不大于信道容量C, 就一定存在这样的编码方法：在采用最大似然译码时,其误码率可以任意小.该定理在理论上给出了

7、对给定信道通过编码所能到达的编码增益的上限,并指出了为到达理论极限应采用的译码方法.在信道编码定理中,香农提出了实现最正确编码的三个根本条件 ：(1 )采用随机编译码方式 ; (2 )编码长度 , 即分组的码组长度无限 ; (3)译码采用最正确的最大似然译码算法。二、信道编码理论2.1、信道编码的目的在数字通信系统中由于信道内存在加性噪声及信道传输特性不理想等容易造成码间串扰同时多用户干扰、多径传播和功率限制等也导致错误译码。为了确保系统的误比特率指标通常采用信道编码。信道编码是为了保证信息传输的可靠性、提高传输质量而设计的一种编码。它是在信息码中增加一定数量的多余码元，使码字具有一定的抗干扰

8、能力。2.2、信道编码的实质 信道编码的实质就是在信息码中增加一定数量的多余码元称为监视码元，使它们满足一定的约束关系，这样由信息码元和监视码元共同组成一个由信道传输的码字。举例而言，欲传输k位信息，经过编码得到长为n(nk)的码字，则增加了 n - k = r 位多余码元，我们定义 R = k / n 为编码效率。【2】2.3、 信道编码公式 令信息速率为fb，经过编码以后的速率为ft，定义：Rfb/ft为编码率。则对于任何一个信道，总存在一个截止速率R0，只要RR0，总可以到达：BERCR2-nR0，其中CR为*个常数，n为编码的约束长度。 对于等概二进码、AWGN信道，有：三、线性分组码

9、的编译码原理3.1、 线性分组码的根本概念一个n ，k线性分组码， 是把信息划成k个码元为一段(称为信息组)， 通过编码器变成长为n个 码元的一组， 作为n， k线性分组码的一个码字。 假设每位码元的取值有q种(q为素数幂)， 则共有qk个码字。 n长的数组共有qn组， 在二进制情况下， 有2n个数组。 显然， qn个n维数组(n重)组成一个GF(q)上的n维线性空间。 如果qk(或2k)个码字集合构成了一个k维线性子空间， 则称它是一个n，k线性分组码。即将k维k重信息空间的元素线性映射到n维n重矢量空间接收矢量/收码 的k维n重子空间码空间。如下列图为7,3码3.2、生成矩阵和校验矩阵生成

10、矩阵：G称为生成矩阵，因为可以用它产生整个码组A，即有生成矩阵的性质：具有IkQ形式的生成矩阵称为典型生成矩阵。由典型生成矩阵得出的码组A中，信息位的位置不变，监视位附加于其后。这种形式的码组称为系统码。 矩阵G的各行也必须是线性无关的。如果已有k个线性无关的码组，则可以将其用来作为生成矩阵G，并由它生成其余码组。【3】监视矩阵：监视矩阵可用来校验和纠错。四、MATLAB仿真4.1仿真4.1.1原理说明采用模块化编程，力求把每个功能独立成各个模块，让程序更清晰。首先介绍各个子程序及其实现的根本功能。此次还采用了交织与解交织，其主要作用是，一个实际上的突发信道，经过发送端交织和接收端解交织的信息

11、处理后，就完全等效成一个随机独立过失信道，信道交织编码实际上就是一类信道改造技术，它将一个突发信道改造成一个随机独立过失信道。它本身并不具备信道编码检、纠错功能，起到信号预处理的作用4.1.2各子函数说明运行环境为Matlab7.0版本通信过程的每个模块写成子程序函数：Channelcoding 为信道编码函数Channeldecoding 为信道解码纠错子函数Interwaving 为交织子函数Deinterwaving 为解交织子函数addfade为向信道参加衰落参数的子函数awgn 为库函数，向信源加高斯白噪声pskmod 为库函数，用于信号调制，输出为复数pskdemod 为库函数，用

12、于信号解调4.2仿真源程序4.2.1信道编码%信道编码子函数，sym为编码码流，G为生成矩阵，k为编码方式的长度，如7,4码的4function bitcoded=channelcoding(sym,G,k)A=vec2mat(sym,k);U=A*G;U=mod(U,2);bitcoded=reshape(U,1,);4.2.2信道解码function bitdecoded=channeldecoding(recode,Etab,Smatri*,H,n,k)% 前向纠错函数，实现纠错功能% bidecoded为纠错后返回的比特流% recode为输入的比特流% E为错误图样表，S为对应的伴随

13、式表% H为监视矩阵，n,k为码的类型，如7，4码，n=7，k=4 row=length(recode)/n; %行数 E=zeros(row,n); %错误图样 RM=zeros(row,n); %纠错之后的矩阵 R=vec2mat(recode,n); S=R*H; %伴随矩阵 S=mod(S,2); for i=1:row for j=1:2(n-k) %查表纠错 if(S(i,:)=Smatri*(j,:) E(i,:)=Etab(j,:); RM(i,:)=R(i,:)+E(i,:); RM(i,:)=mod(RM(i,:),2); break; end end end bitdec

14、oded=reshape(RM,1,); %转化为比特流4.2.3交织function retbit=interweaving(bitstream,row,col)%功能：实现对输入比特的交积% retbit为交积后返回的比特流向量% bitstream 为需要交积的比特流向量% row 和 col为交积器的行和列，% 通过改变col就可以改变交积深度 retbit=zeros(1,length(bitstream); bitarr=vec2mat(bitstream,row); bitarr=bitarr; for i=1:length(bitstream)/(row*col) temp=b

15、itarr(:,(i-1)*col+1):i*col); retbit(1,(i-1)*(row*col)+1):(i*(row*col)=reshape(temp,1,); end4.2.4解交织function retbits=deinterweaving(bitstream,row,col)%功能：实现对输入比特的解交积%rebits为解交积后返回的比特流% bitstream输入的比特流%row 和 col为交积器的行和列，通过改变col就可以改变交积器的长度 retbits=zeros(1,length(bitstream); bitarr=vec2mat(bitstream,col

16、); for i=1:length(bitstream)/(row*col) temp=bitarr(i-1)*row+1:i*row,:); retbits(1,(i-1)*row*col+1:i*row*col)=reshape(temp,1,); end4.2.5信道衰落function code=addfade(modcode,Tf,isperiod,isfade)%功能：向传输序列modcode叠加衰落性信道的衰落参数k(t)%code为参加衰减参数之后返回的序列。% modcode为调制之后的序列% Tf 为衰落时间，以ms为单位，小于10ms，% Tf=1，表示衰落1ms% is

17、period 周期衰落和一次性衰落的标志，% isperiod=1表示周期性衰落，0表示一次性衰落% isfade表示是否存在衰落，1存在，0不存在衰落直接返回modcodeif(isfade=1) if(isperiod=1) %周期性衰落 for k=1:length(modcode)/(100*Tf) a=(k-1)*100*Tf+31; b=(k-1)*100*Tf+30+10*Tf; modcode(1,a:b)=0.1*modcode(1,a:b); end else %一次衰落 a=31; b=30+10*Tf; modcode(1,a:b)=0.1*modcode(1,a:b)

18、; end code=modcode;else code=modcode;end六 程序及仿真图1、 file1：信道编码对通信系统性能的影响，有无信道编码的影响执行时间：Elapsed time is 355.023518 seconds.ticclc%功能：有无信道编码性能比拟M=2; %进制b=log2(M) ; %每符号比特数n=128*10000 ; %符号数G=1 1 1 1 0 0 0;1 0 1 0 1 0 0;0 1 1 0 0 1 0;1 1 0 0 0 0 1 ; %生成矩阵H=1 0 0 1 1 0 1;0 1 0 1 0 1 1;0 0 1 1 1 1 0 ; %监视

19、矩阵Etab= 0 0 0 0 0 0 0;0 0 0 0 0 0 1; %错误图样 0 0 0 0 0 1 0;0 0 0 0 1 0 0; 0 0 0 1 0 0 0;0 0 1 0 0 0 0; 0 1 0 0 0 0 0;1 0 0 0 0 0 0;Smatri*=Etab*H; %对应的伴随式sym=randint(n,1,M);sym=de2bi(sym,left-msb); %模拟信源编码bitcoded=channelcoding(sym,G,4); %信道编码，7,4码modbit=pskmod(bitcoded,M);%在传输序列modbit参加AWGN噪声Snr=0:0.

20、2:15; %噪声为0到15dL=length(snr)ser=zeros(1,L);ser2=zeros(1,L);for k=1:L y=awgn(modbit,10*log10(b)+snr(k),measured); zsym=pskdemod(y,M); %复数解调 zbit=de2bi(zsym,left-msb); recode=reshape(zbit,1,); Rstream=recode; err=(Rstream=bitcoded); errnum=sum(err); ser(k)=log10(errnum/length(bitcoded); %纠错 bitdecoded

21、=channeldecoding(Rstream,Etab,Smatri*,H,7,4); err=(bitdecoded=bitcoded); errbits=sum(err); ser2(k)=log10(errbits/(length(bitcoded);endplot(snr,ser,b-*)hold onplot(snr,ser2,r-o)grid onlegend(没有信道编码,信道编码);*label(Eb/No(dB);ylabel(SER);title(2PSK有无信道编码性能比拟);toc %Elapsed time is 278.288819 seconds.程序运行时间

22、结论：由上图在较高信噪比的时候编码增益很明显大大提高了误码性能，但是在低信噪被的时候编码增益较小甚至可能是负值，则是因为编码后信息比特能量扩散到多个编码比特上，每个编码比特能量有所降低，如果信噪比低的话，编码冗余带来的性能增益可能弥补不了编码比特的能量的降低，因此信道*噪比的波动会显著降低编码性能2、file2：在周期性深衰落的信道条件下，交织对通信系统性能的影响ticclc%功能：有无信道编码性能比拟M=2; %进制b=log2(M) ; %每符号比特数n=128*10000 ; %符号数interrow=8;intercol=10;%交积矩阵的行和列G=1 1 1 1 0 0 0;1 0

23、1 0 1 0 0;0 1 1 0 0 1 0;1 1 0 0 0 0 1 ; %生成矩阵H=1 0 0 1 1 0 1;0 1 0 1 0 1 1;0 0 1 1 1 1 0 ; %监视矩阵Etab= 0 0 0 0 0 0 0;0 0 0 0 0 0 1; %错误图样 0 0 0 0 0 1 0;0 0 0 0 1 0 0; 0 0 0 1 0 0 0;0 0 1 0 0 0 0; 0 1 0 0 0 0 0;1 0 0 0 0 0 0;Smatri*=Etab*H; %对应的伴随式sym=randint(n,1,M);sym=de2bi(sym,left-msb); %模拟信源编码bit

24、coded=channelcoding(sym,G,4); %信道编码，7,4码interv=interweaving(bitcoded,interrow,intercol); %交积向量modbit=pskmod(bitcoded,M);modbit2=pskmod(interv,M);%向传输序列modcode叠加衰落性信道的衰落参数k(t)modbitfade=addfade(modbit,1,1,1);modbitfade2=addfade(modbit2,1,1,1);%1ms周期性衰落modbitfade3=addfade(modbit2,2,1,1);%衰落时长2ms%在传输序列

25、modbit参加AWGN噪声snr=0:0.2:25; %噪声为0到25dL=length(snr)ser=zeros(1,L);ser2=zeros(1,L);for k=1:L y=awgn(modbitfade,10*log10(b)+snr(k),measured); y2=awgn(modbitfade2,10*log10(b)+snr(k),measured); y3=awgn(modbitfade3,10*log10(b)+snr(k),measured); zsym=pskdemod(y,M); %复数解调 zsym2=pskdemod(y2,M); zsym3=pskdemo

26、d(y3,M); zbit=de2bi(zsym,left-msb); zbit2=de2bi(zsym2,left-msb); zbit3=de2bi(zsym3,left-msb); recode=reshape(zbit,1,); recode2=reshape(zbit2,1,); recode3=reshape(zbit3,1,); deinterv=deinterweaving(recode2,interrow,intercol);%解交积向量 deinterv3=deinterweaving(recode3,interrow,intercol); Rstream=recode;

27、Rstream2=deinterv; Rstream3=deinterv3; %纠错 bitdecoded=channeldecoding(Rstream,Etab,Smatri*,H,7,4); bitdecoded2=channeldecoding(Rstream2,Etab,Smatri*,H,7,4); bitdecoded3=channeldecoding(Rstream3,Etab,Smatri*,H,7,4); err=(bitdecoded=bitcoded); errbits=sum(err); ser(k)=log10(errbits/(length(bitcoded); err2=(bitdecoded2=bitcoded); errbits2=sum(err2); ser2(k)=log10(errbits2/(length(bitcoded); err3=(bitdecoded3=bitcoded); errbits3=sum(err3); ser3(k)=log10