matlab - 在matlab中为布尔矩阵求解Ax=b的高斯消元法

Question

为极长的帖子道歉。我有以下代码，我想使用 matlab 和高斯消除来计算 Ax=b。当我将输出与实际预期输出进行比较时，输出部分正确。我不确定我实际上哪里出错了。也包括代码和示例矩阵。

[ grow, gcol ] = size(generator_matrix);
for i = 1 : grow - 1
    if ( generator_matrix(i,i) == 0 )
        t =  min ( find ( generator_matrix(i+1:grow,i) ~= 0 ) + i );
        temp = generator_matrix(i,:);
        generator_matrix(i,:) = generator_matrix(t,:);
        generator_matrix(t,:) = temp;
        tb = encoded_data(i,:);
        encoded_data(i,:)=encoded_data(t,:);
        encoded_data(t,:) = tb;
    end;
    for j = i+1 : grow
        m = -generator_matrix(j,i) / generator_matrix(i,i);
        generator_matrix(j,i) = 0;
        generator_matrix(j, i+1:grow) = generator_matrix(j, i+1:grow) + m * generator_matrix(i, i+1:grow);
        encoded_data(j,:) = mod(encoded_data(j,:) + m * encoded_data(i,:),2);
    end;
end;


x(grow,:) = encoded_data(grow,:) / generator_matrix(grow, grow);
for i = grow - 1 : -1 : 1
    x(i,:) = mod(( encoded_data(i,:) - sum ( x(i+1:grow) .* generator_matrix(i, i+1:grow) ) ),2) / generator_matrix(i,i);
end;

样本值为：生成矩阵的值为

    [ 1 1 1 1          
    0 1 0 0
    0 0 1 0
    0 0 0 1]    

编码数据的值

 [  1 1 1 1 0
    0 0 0 0 1
    0 1 1 1 0
    0 0 1 0 0 ]

x 的值应该是（但我目前得到的输出不是下面的矩阵）

     [1 1 1 1 0
      0 0 0 0 1
      0 1 1 1 0
      0 0 1 0 0 ]

任何建议或提示都会非常有帮助。

我也尝试为更大的矩阵执行代码 20x20，但我没有得到正确的结果。巴维亚