《数字通信实验报告.》由会员分享，可在线阅读，更多相关《数字通信实验报告.（14页珍藏版）》请在金锄头文库上搜索。
1、 数字通信实验报告 武汉理工大学数字通信实验报告班级： 电子与通信工程143班 姓名： 王富贵 学号： 1049721403190 教师： 吕 锋 日期： 2015.04.3 实验一1、 实验项目基于MATLAB的离散无记忆高斯信源的失真-率函数曲线（图3.4-2）仿真；2、 实验目的(1)、理解信息率失真函数的定义与物理意义； (2)、分析离散信源在误码失真下的信息率失真函数表达式； (3)、提高综合运用所学理论知识独立分析和解决问题的能力； (4)、使用相关软件进行曲线的绘制。3、 实验内容与理论依据 实验内容：分析离散信源在误码失真下的信息率失真函数表达式，并绘制曲线图。 理论依据：信息
2、率失真函数的定义 研究在限定失真下为了恢复信源符号所必需的信息率，简称率失真理论。信源发出的符号传到信宿后，一般不能完全保持原样，而会产生失真。要避免这种失真几乎是不可能，而且也无必要，因为信宿不管是人还是机器，灵敏度总是有限的，不可能觉察无穷微小的失真。倘若在处理信源符号时允许一定限度的失真，可减小所必需的信息率，有利于传输和存储。率失真理论就是用以计算不同类型的信源在各种失真限度下所需的最小信息率。因此，这一理论是现代所有信息处理问题的理论基础。 香农首先定义了信息率失真函数R(D)，并论述了关于这个函数的基本定理。 定理指出：在允许一定失真度D的情况下，信源输出的信息传输率可压缩到R(D
3、)值，这就从理论上给出了信息传输率与允许失真之间的关系，奠定了信息率失真理论的基础。信息率失真理论是进行量化、数模转换、频带压缩和数据压缩的理论基础。离散信源：信源是信息的来源，是产生消息、时间离散的消息序列以及时间连续的消息的来源。信源输出的消息都是随机的，因此可以用概率来描述其统计特性。信源在数学上可以用随机变量、随机序列和随机过程来表示。信息是抽象的，信源则是具体的。离散平稳无记忆信源输出的符号序列是平稳随机序列，并且符号之间是无关的，即是统计独立的。同时，由于是平稳信源，每个随机变量的统计特性都相同。 信息率失真函数 对给定的一个信源随机变量，服从概率分布，为失真测度，定义信息率失真函
4、数为：其中为与其复制的互信息，min是对所有满足以下性质的条件概率取值。在给定的信源概率分布及条件概率的乘积所得的联合分布下的平均失真 在讨论信息率失真函数时，考虑到信源与信宿之间有一个无失真信道，称它为试验信道，对离散信源可记为，对限失真信源这一试验信道集合可定义为： 根据前面在互信息中已讨论过的性质： 且互信息是的上凸函数，其极限值存在且为信道容量： 这里，我们给出其对偶定义： 即互信息是的下凸函数。其极限值存在且为信息率失真函数。它还存在下列等效定义1：对给定的一个失真度D，率失真函数R（D）定义为： 对给定的一个码率R，率失真函数D（R）定义为：直观地结实，率失真函数R（D）是在信源序
5、列与复制序列的失真不超过D的条件下最小可能的码率（或信源最大可能的压缩率）。相对偶地，失真率函数D（R）是给定码率（压缩率）R条件所能达到的最小失真。至此，我们已给定R(D)函数一个初步描述。由定义，R(D)函数是在限定失真为最大允许失真为D时信源最小信息速率，它是通过改变试验信道特性（实际上是信源编码）来达到的。所以R(D)是表示不同D值时对应的理论上最小信息速率值。 然而对于不同的实际信源，存在着不同类型的信源编码，即不同的试验信道特性并可以求解出不同的信息率失真R(D)函数，它与理论上最佳的R(D)之间存在着差异，它反映了不同方式信源编码性能的优劣，这也正是R(D)函数的理论价值所在。特
6、别对于连续信源，无失真是毫无意义的，这时R(D)函数具有更大的价值。信息率失真函数的迭代计算 首先需要指出的是，达到率失真函数的条件概率及输出字母概率分布都不一定是唯一的。 具体迭代算法可以按如下步骤进行9： (1)先假定一个负数作为，选定初始转移概率组成阶初始矩阵。 (2)把选定的初始转移概率代入表达式中，得到相应的，然后用代入表达式中，得到相应的。 (3)再用代入表达式中，得到相应的，然后用代入表达式中，得到相应的。 (4)以此推类进行下去，直到与相当接近，其差别已在允许的精度范围之内，以及与相当接近，其差别也在允许的精度范围之内，则或就是这个值所对应的信息率失真函数的近似值。 (5)再选
7、定一个略大一些的负数作为值，重复以上的迭代计算过程，得到值的信息率失真函数的近似值。 (6)这种过程一直到信息率失真函数逼近于零为止，随着的选定就可得到信息率失真函数的曲线。4、 MATLAB程序框图与代码 function=RateDF(Pa,d,S)format long d=input(失真矩阵d=);Pa=input(输入概率分布 Pa=);r=input(输入信源数r=);s=input(输出信源数s=);S=input(拉式乘子S=);times=input(迭代次数times=);r,s=size(d);if(length(find(Pa10e-10) error(Not a p
8、rob.vector,component do not add up to 1!)endif(r=length(Pa) error(The parameters do not match!);endpba=;RS=;DS=;m=1;for z= 1: times Pba(1:r,1:s,1)=1/s*ones(r,s); for j=1:s Pb(j,1)=0; for i=1:r Pb(j,1)=Pb(j,1)+Pa(i)*Pba(i,j,1); end end for i=1:r temp(i)=0; for j=1:s temp(i)=temp(i)+Pb(j,1)*exp(S(m)*d
9、(i,j); end end for i=1:r for j=1:s Pba(i,j,2)=(Pb(j,1)*exp(S(m)*d(i,j)/temp(i); end D(1)=0; for i=1:r for j=1:s D(1)=D(1)+Pa(i)*Pba(i,j,1)*d(i,j); end end R(1)=0; for i=1:r for j=1:s if(Pba(i,j,1)=0) R(1)=R(1)+Pa(i)*Pba(i,j,1)*log2(Pba(i,j,1)/Pb(j,1); end end end n=2; while(1) for j=1:s Pb(j,n)=0; f
10、or i=1:r Pb(j,n)=Pb(j,n)+Pa(i)*Pba(i,j,n); end end for i=1:r temp(i)=0; for j=1:s % disp(SM:);disp(S(m); temp(i)=temp(i)+Pb(j,n)*exp(S(m)*d(i,j); end end for i=1:r for j=1:s if(temp(i)=0) Pba(i,j,n+1)=(Pb(j,n)*exp(S(m)*d(i,j)/temp(i); end end end D(n)=0; for i=1:r for j=1:s D(n)=D(n)+Pa(i)*Pba(i,j,n
11、)*d(i,j); end end R(n)=0; for i=1:r for j=1:s if(Pba(i,j,n)=0) R(n)=R(n)+Pa(i)*Pba(i,j,n)*log2(Pba(i,j,n)/Pb(j,n); end end end %disp(E1:);disp(abs(R(n)-R(n-1); %disp(E2:);disp(abs(D(n)-D(n-1); if(abs(R(n)-R(n-1)=10(-7) if(abs(D(n)-D(n-1)=10(-7) break; end end n=n+1; end S(m+1)=S(m)+0.5; if(abs(R(n)10(-7) end pba=Pba(:,:,:); RS=RS R(n); DS=DS D(n); m=m+1;
