求线性分组码的标准阵的matlab程序

信息论作业里有求线性分组码的标准阵的题目。手算太过繁琐，而且容易出错。在matlab里没有找到现成的函数，于是自己写了一个。

构造标准阵的步骤如下：

1.将2^k个码字排成一行，全零码字排在第一列。其余的(2^n-2^k)个非码字向量构成集合E。

2.从集合E中取一个重量最小的向量e置于下一行第一列，作为陪集首。对任意码字ci，将e+ci置于ci下方。如此构造出下一行。

3.从E中删去刚才构造出的那行中的出现的向量。

4.重复第2、3步，直到E为空。

可以看出，码空间中共2^n个向量，都作为标准阵的元素。标准阵每行有2^k个元素，共2^(n-k)行。

下面贴出我写的matlab程序。水平很菜，大家见谅。

%这个程序由线性分组(n,k)码的生成矩阵G

%求出最小汉明距离dmin

%产生一致校验矩阵H

%产生标准阵standard_mat

clear;

n = 7;

k = 4;

num = 2^n;

c_num = 2^k;

%G = [1 1 0 1 0; 0 1 0 1 0];

G = [1 0 0 0 1 1 1;

    0 1 0 0 1 0 1;

    0 0 1 0 0 1 1;

    0 0 0 1 1 1 0];

dmin = gfweight(G);

disp