《数字通信实验报告资料.doc》由会员分享，可在线阅读，更多相关《数字通信实验报告资料.doc（14页珍藏版）》请在三一办公上搜索。
1、武汉理工大学数字通信实验报告班级： 信息154 姓名： 冯 超 学号： 1049731503280 教师： 吕 锋 日期： 2016.03.26 实验一1、 实验项目基于MATLAB的离散无记忆高斯信源的失真-率函数曲线仿真；2、 实验目的(1)、理解信息率失真函数的定义与物理意义； (2)、分析离散信源在误码失真下的信息率失真函数表达式； (3)、提高综合运用所学理论知识独立分析和解决问题的能力； (4)、使用相关软件进行曲线的绘制。3、 实验内容与理论依据 实验内容：分析离散信源在误码失真下的信息率失真函数表达式，并绘制曲线图。 理论依据：信息率失真函数的定义 研究在限定失真下为了恢复信源
2、符号所必需的信息率，简称率失真理论。信源发出的符号传到信宿后，一般不能完全保持原样，而会产生失真。要避免这种失真几乎是不可能，而且也无必要，因为信宿不管是人还是机器，灵敏度总是有限的，不可能觉察无穷微小的失真。倘若在处理信源符号时允许一定限度的失真，可减小所必需的信息率，有利于传输和存储。率失真理论就是用以计算不同类型的信源在各种失真限度下所需的最小信息率。因此，这一理论是现代所有信息处理问题的理论基础。 香农首先定义了信息率失真函数R(D)，并论述了关于这个函数的基本定理。 定理指出：在允许一定失真度D的情况下，信源输出的信息传输率可压缩到R(D)值，这就从理论上给出了信息传输率与允许失真之
3、间的关系，奠定了信息率失真理论的基础。信息率失真理论是进行量化、数模转换、频带压缩和数据压缩的理论基础。离散信源：信源是信息的来源，是产生消息、时间离散的消息序列以及时间连续的消息的来源。信源输出的消息都是随机的，因此可以用概率来描述其统计特性。信源在数学上可以用随机变量、随机序列和随机过程来表示。信息是抽象的，信源则是具体的。离散平稳无记忆信源输出的符号序列是平稳随机序列，并且符号之间是无关的，即是统计独立的。同时，由于是平稳信源，每个随机变量的统计特性都相同。 信息率失真函数R(D)对给定的一个信源随机变量，服从概率分布，为失真测度，定义信息率失真函数为：其中为与其复制的互信息，min是对
4、所有满足以下性质的条件概率取值。在给定的信源概率分布及条件概率的乘积所得的联合分布下的平均失真 在讨论信息率失真函数时，考虑到信源与信宿之间有一个无失真信道，称它为试验信道，对离散信源可记为，对限失真信源这一试验信道集合可定义为： 根据前面在互信息中已讨论过的性质： 且互信息是的上凸函数，其极限值存在且为信道容量： 这里，我们给出其对偶定义： 即互信息是的下凸函数。其极限值存在且为信息率失真函数。它还存在下列等效定义1：对给定的一个失真度D，率失真函数R（D）定义为： 对给定的一个码率R，率失真函数D（R）定义为：直观地结实，率失真函数R（D）是在信源序列与复制序列的失真不超过D的条件下最小可
5、能的码率（或信源最大可能的压缩率）。相对偶地，失真率函数D（R）是给定码率（压缩率）R条件所能达到的最小失真。至此，我们已给定R(D)函数一个初步描述，如图1。图1 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 prob.ve
8、ctor,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(i,j);
9、 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; for i=1
10、: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)*d(i,
11、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)
12、end pba=Pba(:,:,:); RS=RS R(n); DS=DS D(n); m=m+1; endend k,l,q=size(pba); Pba=pba(:,:,q); Rmin=min(RS); Dmax=max(DS); Smax=S(m-1);disp(输入正确，迭代结果如下：);disp(最小信息率Rmin:);disp(Rmin);disp(最大Dmax:);disp(Dmax);disp(最佳转移概率分布Pba:);disp(Pba);disp(最大拉式乘子Smax:);disp(Smax);plot(DS,RS)xlabel(允许的失真度D)ylabel(信息率失真函
13、数R(D)title(信息率失真函数R(D)的曲线图)5、 实验结果及分析某二元离散无记忆信源其失真矩阵为求该信源的，和函数。解：二元对称信源，其失真矩阵为可计算得：，根据参量表达式可求得：这里，MATLAB运行结果如下：输入：d=0 1;1 0;Pa=0.5 0.5;r=2;s=2;S=-99.5;times=100;输出：很好！输入正确，迭代结果如下：最小信息率Rmin: 0最大Dmax: 0.50000000000000最佳转移概率分布Pba: 0.62136760936669 0.37863239063331 0.37645018769302 0.62354981230698最大拉式乘
14、子Smax: 0ans = 0.62136760936669 0.37863239063331 0.37645018769302 0.62354981230698信息率失真函数曲线图如图2所示：图2 信息率失真函数曲线图6、 实验结论待定常数S就是函数的斜率，根据函数的S参量表述理论，也可得到函数曲线，最终完成信息率失真函数迭代计算过程。在求解信息率失真函数时，达到信息率失真函数时的转移概率及信宿的概率分布也不一定是唯一的。实验二一、实验项目1、设定符号错误概率为10的负5次方，基于MATLAB仿真分析无记忆调制的最佳接收机性能。二、实验目的1、通过实验，进一步巩固在课堂上面学到的理论知识, 学习并理解加性高斯白噪声信道的最佳接收机2、在实验过程中，对理论知识和公式进行理解，转化为相应的matlab程序语言，增强matlab编程的能力；三、实验内容与理论依据（1）最佳接收机评估各种无记忆调制方法的错误概率 二进制调制的错误概率设两个信号波形是和 (称为双极性信号）。在给定发送的情况下，s1(t)的错误概率是r1且要求有足够的SNR支持M的增加。相反地，M元正交信号的R/W0.693(-1.6dB),当时，可以得到尽可能小的错误概率。该信噪比条件是在如下极限情况下可靠传输所需要的最小SNR：信道带宽且相应的。
