p_yx=?input（‘信道矩陣P_YX=‘）;???%輸入信道矩陣
[rs]=size（p_yx）;
b=input（‘輸入迭代精度：‘）;????%輸入迭代精度
for?i=1:r
???p（i）=1.0/r;??????%初始概率為均勻分布
end
C0=-inf;
C=1;
n=0;
while（abs（C-C0）/C>=b）???%迭代計算
????for?j=1:s?????%計算q（j）
????????q（j）=0;
????????for?i=1:r
????????????q（j）=q（j）+p（i）*p_yx（ij）;
????????end
????end
????for?i=1:r?%計算a（i）
????????d（i）=0;
????????for?j=1:s
????????????if（p_yx（ij）==0）
????????????????d（i）=d（i）+0;
????????????else
????????????????d（i）=d（i）+p_yx（ij）*log（p_yx（ij）/q（j））;
????????????end
????????end
????????a（i）=exp（d（i））;
????end
????u=0;
????for?i=1:r?%計算u
????????u=u+p（i）*a（i）;
????end
????for?i=1:r
????????p（i）=p（i）*a（i）/u;?%重新賦值p（i）
????end
????C0=C;
????C=log2（u）;?%計算C
????n=n+1;
end
disp（‘輸入分布：‘）
p
disp（‘信道矩陣為：‘）;
disp（p_yx）;
disp（‘迭

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????930??2010-06-30?23:12??huffman.m

?????文件????????894??2010-06-12?19:17??c.m

?????文件????????771??2010-06-30?23:13??shannon.m

?????文件????????781??2010-06-30?23:12??fano.m

-----------?---------??----------?-----??----

?????????????????3376????????????????????4