数字通信实验-武理.docxVIP

武汉理工大学 数字通信》课程实验报告 学生班级： 信研 18 学生姓名： 陈 学生学号： 104972180 2 任课教师： 韩 日 期： 2019.06.20 数字通信课程实验报告 数字通信课程实验报告 PAGE PAGE # 实验一 一、实验项目 基于 MATLAB 的离散无记忆高斯信源的失真 -率函数曲线仿真 二、实验目的 1、了解信号失真 -率函数的定义、性质及其物理性质； 2、掌握离散无记忆高斯信源的率失真函数计算方法； 3、学会使用 MATLAB 实现离散无记忆高斯信源的失真 -率函数的计算方法； 4、提高综合运用所学理论知识独立分析和解决问题的能力。 三、实验内容与理论依据 1、实验内容 使用 MATLAB 工具仿真离散无记忆高斯信源的失真 -率函数曲线，并对其 失真率函数曲线的特点性质进行分析。 2、基础理论 利用采样定理，模拟信源的输出可转化成一个等效的离散时间抽样序列， 然后对样本幅度进行量化和编码。信号值幅度量化带来了数据的压缩，但同时 引入了某些波形失真或信号保真度的损失。 假设有一个无记忆的信源， 其连续幅值输出 X的 PDF是 p(x)，量化后幅值 符号集是 X ，每符号的失真度是 d x,x ，那么，为了以小于等于 D 的失真来 表示无记忆信源的输出 X ，每信源符号需要的最低比特率称之为率失真函数 R(D)，定义为 R(D) minp(x?|x):Edd(x|X?) D I(X;X?) 式中，I X;X 是互信息，一般，当D增大时 R(D)减小，反之 D减小时 R(D) 增大。从率失真定义可以看到， R(D)取决于 p(x) 的统计特性和失真量度，两者 中不管哪个有变化都会改变 R(D) 。这里还要提到一点， 对于很多信源的统计特 性和失真量度而言，并不存在封闭形式的率失真函数 R(D)。 香农第三定理指出，只要 RR(D)，一个无记忆信源 X 能以码率 R 编码而 失真不会超过 D，相反若 RR(D)，则任何编码的失真度都将超过 D。显然，对 于一个给定的失真度， 任何一个信源的率失真函数 R(D)代表了信源速率的下边 界。这就从给出了信息传输率与允许失真之间的关系，奠定了信息率失真理论 的基础 在讨论信息率失真函数时，考虑到信源与信宿之间有一个无失真信道，称 它为试验信道， 对离散信源可记为 Pji ，对限失真信源这一试验信道集合可定义 为： PD Pji : D d i j pi Pjidij . 根据前面在互信息中已讨论过的性质： I(U;V) I(pi;Pji ) 且互信息是 pi 的上凸函数，其极限值存在且为信道容量： C maxI(pi;Pji ) pi 这里，我们给出其对偶定义： R(D) PminP I(U;V) PminP I(pi;Pji)Pji PD Pji PD 即互信息是 Pji 的下凸函数。其极限值存在且为信息率失真函数。 由定义， R(D) 函数是在限定失真为最大允许失真为 D 时信源最小信息速率， 它是通过改 变试验信道特性(实际上是信源编码)来达到的。 (1)平方误差失真下高斯信源的率失真函数 常用模型之一是连续幅度、无记忆信源的高斯信源模型。这种模型具有高 斯统计特性，以平方误差 d(x,x) (x x)2为失真量度，其率失真函数已知， 为 1log 1log 2 Rg(D) 2log D 0 D 2 D2 式中， δ2 是信源的方差。注意， 式中， δ2 是信源的方差。注意， Rg(D)与信源均值 E[X] 无关。其曲线如 下图所示。如果将上式中 D 和R 的函数关系式倒一 下就可以用 R 来表示 D， Dg(R) 2 2R 2 这个函数称为离散时间无记忆高斯信源的失真 率函数。如果用分贝表示，可得 2 10log Dg (R) 6R 10log 2 R(d)的上边界，零均值、有限方差 δ2(用均方 误差失真量度）、连续幅度的无记忆信源，其率失真函数的上边界为 R(D) 1 R(D) 12 log 2 D 0 D x2 （2）一般信源信息率失真函数的迭代计算 首先需要指出的是，达到率失真函数的条件概率 P vj |ui 及输出字母概率 分布 q vj 都不一定是唯一的。 具体迭代算法可以按如下步骤进行： 1 先假定一个负数作为 S1 ，选定初始转移概率 p1 bj /ai 1 组成 rs r s 阶初始矩阵。 r ②把选定的初始转移概率 p1 bj /ai 代入表达式 p bj p ai p bj /ai i1 中，得到相应的p1 bj 中，得到相应的 p1 bj ，然后用 p1 bj 代入表达式 p bj /ai Sd ai ,bj p bj eSd ai,bj s Sd ai ,bj p bj e i j j1 中，得到相应的 p2