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