matlab生成m元n次扩展信源

在信息论与编码的时候，我们经常使用二元n次扩展信源，做一些实验的时候，需要快速生成一个n次扩展信源，下面我们利用matlab来实现n次扩展信源的生成： 首先我们来生成它的编码，这里按照自然顺序生成每一行，没有按照格雷编码。

code=zeros(2^n,n); p=zeros(n,1); for i=0:2^n-1 c=dec2bin(i,n); for j=1:n code(i+1,j)=str2num(c(j)); end end

接着我们计算每行编码的概率：

pp=1; c=code(i+1,:); for k=1:n if c(k)==0 pp=pp*P(1); else pp=pp*P(2); end end p(i+1)=pp;

汇总的matlab代码如下：

function [code,p]=n_kuo(P,n) %二元信源的n次扩展信源产生器 code=zeros(2^n,n); p=zeros(n,1); for i=0:2^n-1 c=dec2bin(i,n); for j=1:n code(i+1,j)=str2num(c(j)); end pp=1; c=code(i+1,:); for k=1:n if c(k)==0 pp=pp*P(1); else pp=pp*P(2); end end p(i+1)=pp; end end

我们以0-0.7 1-0.3的信源概率，4次扩展为例，结果如下： 接着，我们也做了m元信源做n次扩展的结果，代码见资源： m元信源做n次扩展的matlab实现.rar