信息科学基础实验报告[2]

信息科学基础验证性实验报告其失真矩阵为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=input(拉式乘子 S=)times=input(迭代次数 times= )r,s=size(d);if (len gth(fi nd(Pa10e-10)error(Not a prob.vector,comp onent do not add up to r)endif (r=le ngth(Pa);error( The parameters do not 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=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)=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(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(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