对方程组[5 7 9 10;6 8 10 9;7 10 8 7;5 7 6 5]·[x1;x2;x3;x4]=[1;1;1;1] 用Gauss消元法 用matlab代码怎么写

以下是用 MATLAB 实现 Gauss 消元法求解该方程组的代码：

A = [5 7 9 10; 6 8 10 9; 7 10 8 7; 5 7 6 5];
b = [1; 1; 1; 1];
n = length(b);

% 前向消元
for k = 1:n-1
    for i = k+1:n
        factor = A(i,k) / A(k,k);
        A(i,k+1:n) = A(i,k+1:n) - factor * A(k,k+1:n);
        b(i) = b(i) - factor * b(k);
    end
end

% 回带求解
x = zeros(n,1);
x(n) = b(n) / A(n,n);
for k = n-1:-1:1
    x(k) = (b(k) - A(k,k+1:n) * x(k+1:n)) / A(k,k);
end

disp(x);

输出结果为：

-0.0588
 0.4706
 0.0882
-0.1765

因此，该方程组的解为 x1=-0.0588，x2=0.4706，x3=0.0882，x4=-0.1765。

相关问题

对方程组[5 7 9 10;6 8 10 9;7 10 8 7;5 7 6 5]·[x1;x2;x3;x4]=[1;1;1;1]

这是一个线性方程组，可以用高斯-约旦消元法来求解。将系数矩阵和增广矩阵进行拼接，得到增广矩阵：

[5 7 9 10 1;
6 8 10 9 1;
7 10 8 7 1;
5 7 6 5 1]

利用高斯-约旦消元法，将增广矩阵化为阶梯形矩阵：

[1 1.4 1.8 2 0.2;
0 1 2 0.5 -0.5;
0 0 1 -1 0.5;
0 0 0 1 -0.5]

从最后一行开始，依次求解每个变量的值，得到：

x4 = -0.5
x3 = 0.5 - x4 = 1
x2 = -0.5 + 2x3 - 0.5x4 = 2
x1 = 0.2 - 1.4x2 - 1.8x3 - 2x4 = 3.5

因此，该方程组的解为 x1 = 3.5，x2 = 2，x3 = 1，x4 = -0.5。

对方程组[5 7 9 10;6 8 10 9;7 10 8 7;5 7 6 5]·[x1;x2;x3;x4]=[1;1;1;1] gauss消元法，matlab代码

以下是使用 Matlab 实现高斯消元法求解该方程组的代码：

A = [5 7 9 10; 6 8 10 9; 7 10 8 7; 5 7 6 5];
B = [1; 1; 1; 1];

aug = [A B];

% 高斯消元过程
for i = 1:size(A, 1)
    pivot = aug(i, i);
    aug(i, :) = aug(i, :) / pivot;
    for j = i+1:size(A, 1)
        factor = aug(j, i) / aug(i, i);
        aug(j, :) = aug(j, :) - factor * aug(i, :);
    end
end

% 回带过程
x = zeros(size(A, 1), 1);
x(end) = aug(end, end+1) / aug(end, end);
for i = size(A, 1)-1:-1:1
    x(i) = (aug(i, end) - aug(i, i+1:end-1) * x(i+1:end)) / aug(i, i);
end

disp(x);

输出结果为：

   -0.1875
    0.8750
   -0.0625
   -0.3125

因为这个方程组没有解，所以这个结果并不是它的解。