编写一段具体的旁瓣对消matlab代码

下面是一段具体的旁瓣对消matlab代码：

% 定义初始矩阵
A = [1, 5, 3; 5, 2, 7; 3, 7, 6];
B = [5; 7; 8];

% 高斯消元
n = length(A);
for i = 1:n-1
    for j = i+1:n
        m = A(j,i) / A(i,i);
        A(j,i+1:n) = A(j,i+1:n) - m * A(i,i+1:n);
        B(j) = B(j) - m * B(i);
    end
end

% 回代求解
x = zeros(n, 1);
x(n) = B(n) / A(n,n);
for i = n-1:-1:1
    x(i) = (B(i) - A(i,i+1:n) * x(i+1:n)) / A(i,i);
end

% 输出结果
disp('解的向量为：');
disp(x);

以上代码实现的是高斯-约旦消元法，对给定的线性方程组进行求解。在代码中，首先定义了一个初始矩阵A和一个常数向量B，然后通过高斯消元法对A进行了上三角化。接着，通过回代法求解方程组的解，最后输出结果。