卷积编码译码通信原理课程设计MATLAB样本.doc

一． 卷积码背景
Elias于1955年提出卷积码，1967年Vterbi提出最大似然译码法，广泛运用于当代通信。
二． 编码原理
1. 解析法
a. 离散卷积法（用于定义）
涉及：connection pictorial，connection vector，connection polynomial
如下图：设输入信息序列为u=(u0,u1,u2……),相应输c出序列为（c10,c11,c12,c13…）,c2=(c20,c21,c22,c23…),则编码方程可写为：c1=u*g1,c2=u*g2,c=(c1,c2),其中g1,g2表达编码器两个脉冲冲激响应，即编码可由输入信息序列U和编码器两个冲激响应卷积得到，故如此称卷积码，其中脉冲冲激响应应是指当输入为(1000…)时所观测到得两个输出序列值，编码器m=3级寄存器冲激响应可持续到k=m+1=4位。
g1=（1011），g2=（1111）

如下图：

普通来说，g1=（g10,g11,g12…g1m） g2=(g20,g21,g22…g2m)编码输出序列合并为一种输出序列为c=(c10,c20,c11,c21…)
b. 生成矩阵（用于理论分析）

c. 码多项式（用于工程）
g1=(1011)=1+x^2+x^3
g2=(1111)=1+x+x^2+x^3
输入序列u=（10111）1+x^2+x^3+x^4
c1=ug1=1+x^7=(10000001)
c2=ug2=(11011101)

2.卷积图形表达法（state diagram，treediagram，trellis diagram）
a．状态图：
编码器下一时刻输出取决于取决于编码器当前状态及下一时刻输入，而当前状态取决于各寄存器存储状态内容，因而总状态数达2^mk个。
如图：

状态表为：
ui
ui-1
si
0
0
a
1
0
b
0
1
c
1
1
d
有如下状态图：

b．树图（将状态图准时间展开）
设初始状态s0=00为树根，对每个时刻也许输入进行分支，分数级数L表达，L=0时，u0=0向上，u0=1向下，依次向后无限延伸，分支上数字表达相应输出，a，b，c，d表达状态。（长处：时序关系清晰）

d．网格图（研究卷积码最大似然译码维特比算法工具）
纵坐标表状态，横坐标表时间。
从图中咱们可看出编码过程中使用K-1个冲洗比特使得编码器恢复初始状态是有必要。

三、卷积码译码
Viterbi译码（最大似然译码）
先验概率条件下，后验概率最大者似然函数最大，最值MAP即最大似然(ML)译码。而最大对数似然函数即计算最小汉明距，如此，比较接受序列和发送序列汉明距，选出最小汉明距序列作为最佳译码即可。
维特比译码思想：
译码过程选取整个集合途径中使使然函数最大途径，若在某节点上发现某条途径不也许获得最大对数似然函数则放弃该途径，然后在剩余“幸存”途径中重新选取直至最后一级，从而大大减轻译码工作量。
长处：
1.途径度量可加性，满足Bellmax最优化原理。
2.局部判决及时除去大量非最优途径节约运算。

译码环节：
1. 从L=m时刻开始，计算每一状态单个途径某些度量值，并存储幸存途径及度量值。
2. L增长，L=m+1，将进入某一状态时分支度量值与前一时段幸存途径度量值相加，选取最大度量值途径进入，存储新幸存途径及度量值，删除所存其她途径。
3.若L节2，否则停止。
基于局部最优化运算等效于整体最优化，且在译码工程中，到达每一状态有两条也许途径，通过ACS工程，择取最佳途径即是最佳译码。

四、卷积劫难性错误（Catastrophic error）
有限数量码无差错引起无限数量已译码数据比特差错，浮现劫难传播充要条件是：生成多项式有共同因子。例如：g1=1+x；
g2=1+x^2；则编码器会引起catastrophic error propagation。

五．Simulink仿真

蒙特卡罗（Monte Carlo，MC）是一种以概率记录理论为指引一类通过仿真来分析科学现象重要办法，使用随机数（或伪随机数）来解决诸多计算问题体现出了很强灵活性和易实现性，该办法重要分三步：构造或描述概率过程；在概率过程中随机抽样；建立各种预计量并给出近似解。MC办法被广泛运用于通信领域,通信中人们最关怀重要问题之一就是系统误码率性能。如果把码元与否出错看作服从0-1分布随机变量,那么误码率正好是该随机变量数学盼望,而对随机变量数学盼望进行预计,正是MC办法中最为典型问题。下面内容正是基于MC而建立起得仿真模型，通过实验现象咱们将会看到卷积编码译码过程和系统性能分析，涉及劫难性性编码现象。

原理方框图：

Simulink模型

重要模块参数设立：
1. 伯努力信号发生器

2. 卷积编码器

3. BPSK调制

4. 加性高斯噪声信道

5. BPSK解调

6. Viterbi译码

Matlab调试程序：
clear;
clc;
x=-10:10；
y=x；
h=1:3;
colordef black
A=[poly2trellis(3,[3 5]),poly2trellis(3,[7 5]),poly2trellis(3,[4 5 7])];
for j=1:3
CHOICE=A(j);
for i=1:length(x)
SNR=x(i);
sim('liuzekun');
y(i)=mean(BitErrorRats);
end
h=semilogy(x,y);
if j==1
set(h,'color','b');
elseif j==2
set(h,'color','r');
elseif j==3
set(h,'color','y');
end
hold on
end

title('卷积编码译码');
legend('劫难编码器','1/2编码器','1/3编码器')
xlabel('信道信噪比');
ylabel('误码率');

六．实验成果及分析

编码器劫难性编码系统成果

对比三个仿真成果可知：劫难性编后果十分严重，在实际应用中必要编码，此外1/3编码器比1/2编码器性能要好。

七、参照资料
1.樊昌信.通信原理(第6版).国防工业出版社
2.周炯盘.通信原理（第3版）.北京邮电出版社
3.黄载禄.通信原理.科学出版社
4. 斯克拉.数字通信.电子工业出版社