0

我对手动 RREF 给定矩阵有疑问。所以我部分经历了它，但我不知道从这里去哪里。谢谢！

 M = [1 0 2 1 18; 
     0 -3 -2 0 -8; 
     -2 -3 0 0 -41; 
     1 0 -1 1 16];
 M2=M;
 M2(3,:) = M2(3,:)+(2*M2(1,:));
 M2(4,:) = M2(4,:)-M2(1,:);
 M3 = M2;
 M3(3,:) = M3(3,:)+M3(2,:);
 M3(3,:) = M3(3,:)-M3(2,:);
 M3(3,:) = M3(3,:)-M3(2,:);
 M3(2,:) = (-1/3)*M3(2,:)

最后我结束了

[1 0 2 1 18; 
0 1 .6667 0 2.6667; 
0 0 6 2 3; 
0 0 -3 0 -2]
4

1 回答 1

0

在您离开的地方接机：

x = [1 0 2 1 18; 0 1 .6667 0 2.6667; 0 0 6 2 3; 0 0 -3 0 -2];
x(3,:) = x(3,:) / 6;

x(1,:) = x(1,:) - 2 * x(3,:);
x(2,:) = x(2,:) - 2/3 * x(3,:);
x(4,:) = x(4,:) + 3 * x(3,:);

x(1,:) = x(1,:) - 1/3 * x(4,:);
x(2,:) = x(2,:) + 2/9 * x(4,:);
x(3,:) = x(3,:) - 1/3 * x(4,:);

与rref(x)此类似产生：

x =

    1.0000         0         0    0.0000   17.1667
         0    1.0000    0.0000         0    2.2223
         0         0    1.0000         0    0.6667
         0         0         0    1.0000   -0.5000

如果您想知道中间的所有步骤，这是一种方便的方法，但显然使用该rref函数通常更好地找到减少的行梯形形式。

请注意，如果您在 matlab 中执行所有步骤而不是复制诸如.6667.

于 2013-04-11T07:57:31.163 回答