%功能:本程序实现信息率失真函数R(D)的迭代计算算法
%说明:可以按照教科书上以前后两次的R的误差来停止迭代,但是我做出的结果图在D=0.5以后多出个尾巴,故采用比较D的方法。
%时间:2010年1月6日
%作者、版权:本程序由+++编写完成,您可以任意使用或修改,无需通知作者本人,我们的口号是“分享才是王道”!
clc
clear all;
p=[0.5,0.2,0.3];
d=[0 1 1;1 0 1;1 1 0];
[row,col]=size(d);
err_liufuliang=0.0001;
S=-100;
k=0;
while S<-0.01
k=k+1;
D1=100;
D2=10;
%R0=100;
%R=10;
pji=ones(row,col)./col;
while abs(D2-D1)>err_liufuliang %abs(R-R0)>err_liufuliang
D1=D2;
%R0=R;
D2=0;
R=0;
qj=p*pji;
for i=1:row
fm=0;
for j=1:col
fm=fm+qj(j)*exp(S*d(i,j));
end
for j=1:col
pji(i,j)=qj(j)*exp(S*d(i,j))/fm;
end
end
for i=1:row
for j=1:col
D2=D2+p(i)*pji(i,j)*d(i,j);
R=R+p(i)*pji(i,j)*log2(pji(i,j)/qj(j));
end
end
end
a(1,k)=D2;
a(2,k)=R;
[k a(1,k) a(2,k)]
S=S*2/3;
end
plot(a(1,1:k),a(2,1:k),'.r');
xlabel('D');
ylabel('R(D)');
RD.rar_RD_rd 算法_rd函数_rd算法
版权申诉
149 浏览量
2022-07-15
03:13:24
上传
评论
收藏 796B RAR 举报
JonSco
- 粉丝: 72
- 资源: 1万+
最新资源
- keil2 + proteus + 8051.exe
- 1961ee27df03bd4595d28e24b00dde4e_744c805f7e4fb4d40fa3f695bfbab035_8(1).c
- mediapipe-0.9.0.1-cp37-cp37m-win-amd64.whl.zip
- windows注册表编辑工具
- mediapipe-0.9.0.1-cp37-cp37m-win-amd64.whl.zip
- 校园通行码预约管理系统20240522075502
- 车类型数据集6250张VOC+YOLO格式.zip
- The PyTorch implementation of STGCN.STGCN-main.zip
- 092300108.cpp
- 车类型数据集6000张VOC+YOLO格式.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈