《信息科学基础实验报告[2].doc》由会员分享，可在线阅读，更多相关《信息科学基础实验报告[2].doc（4页珍藏版）》请在三一文库上搜索。
1、信息科学基础验证性实验报告其失真矩阵为f 0D=I11，求该信源的0丿Dmax，Dmin 和 R(D)函数。专业信息与计算科学班级112班姓名日饶燕学号 55实验日期2013年12月6星期五成绩评定一、题目率失真函数的迭代运算Eg.某二元离散无记忆信源yP(X)J=广01/21、1/2二、程序fun cti on =RateDF(Pa,d,S)format longd=input('失真矩阵 d=')Pa=input( '输入概率分布Pa=')r=input( '输入信源数r=') s=input('输岀信源数s=')S=inpu
2、t('拉式乘子 S=')times=input('迭代次数 times=' )r,s=size(d);if (len gth(fi nd(Pa<=0)=0)error(end'Not a prob.vector,comp onent do not add up to r if (abs(sum(Pa)-1)>10e-10)error('Not a prob.vector,comp onent do not add up to r)endif (r=le ngth(Pa);error( 'The parameters do no
3、t match!' end pba=; RS=; DS=; m=1;for z=1:timesPba(1:r,1:s,1)=1/s*o nes(r,s); for j=1:sPb(j,1)=0;for i=1:rPb(j,1)=Pb(j,1)+Pa(i)*Pb(i,j,1);endendfor i=1:rtemp(i)=0;for j=1:stemp(i)=temp(i)+Pb(j,1)*exp(S(m)*d(i,j);endendfor i=1:rfor j=1:sPba(i,j,2)=(Pb(j,1)*exp(S(m)*d(i,j)/temp(i);endD(1)=0;for i=
4、1:rfor j=1:sD(1)=D(1)+Pa(i)*Pba(i,j,1)*d(i,j);endendR(1)=0;for i=1:rfor j=1:sif (Pba(i,j,1)=0)R(1)=R(1)+Pa(i)*Pba(i,j,1)*log2(Pba(i,j,1)/Pb(j,1);endendendn=2;while (1)for j=1:sPb(j,n)=0;for i=1:rPb(j,n)=Pb(j,n)+Pa(i)*Pba(i,j,n);endendfor i=1:rtemp(i)=0;for j=1:s%disp('SM');disp(S(m); temp(i)
5、=temp(i)+Pb(j,n)*exp(S(m)*d(i,j); endendfor i=1:rfor j=1:sif (temp(i)=0)Pba(i,j,n+1)=(Pb(j,n)*exp(S(m)*d(i,j)/temp(i);endendendD(n)=0;for i=1:rfor j=1:sD(n)=D(n)+Pa(i)*Pba(i,j,n)*d(i,j);endendR(n)=0;for i=1:rfor j=1:sif (Pba(i,j,n)=0)R(n)=R(n)+Pa(i)*Pba(i,j,n)*log2(Pba(i,j,n)/Pb(j,n);endendend%disp(
6、'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)v=10*7)break ;endendn=n+1;endS(m+1)=S(m)+0.5;if (abs(R( n)<10*7)endpba=Pba(:,:,:);RS=RS R(n);DS=DS D(n);m=m+1;endendk,1,q=size(pba);Pba=pba(:,:,q);Rmin=min(RS);Dmax=max
7、(DS);Smax=S(m-1);disp( ' 输入正确，迭代结果如下： ' );disp( ' 最小信息率 Rmin:' );disp(Rmin);disp( ' 最大 Dmax:' );disp(Dmax);disp( ' 最佳转移概率分布 Pba:' );disp(Pba);disp( ' 最大拉氏乘子 Smax:' );disp(Smax);三、运行结果失真矩阵 d=0 1;1 0; 输入概率分布 Pa=0.5 0.5; 输入信源数 r=2;输出信源数 s=2;拉式乘子 S=-99.5; 迭代次数 times=100; 输入正确，迭代结果如下： 最小信息率 Rmin:0最大 Dmax:'0.000000最佳转移概率分布 Pba:0.936669 0.0633310.769302 0.230698最大拉氏乘子 Smax:0ans =0.936669 0.0633310.769302 0.230698