通过 a(i,j)=-1+2max(i,j) 和 b(j)=sum(j to n) a(i,j) 定义一个 nxn 矩阵 A。然后在这个系统上测试 Naive-Gauss (Matlab)

Question

使用 Naive Gauss 测试系统很容易，但我无法定义 A 和 b。我现在的代码如下：

A=zeros(n);


for i=1:n

    for j=1:n
        smax=max(i,j);
        A(i,j)=-1+2*smax;
    end
end

for j=1:n

    sum=0;
    sum=sum+A(i,j);
    b(j)=sum;
end

Answer 1

b 的定义不清楚。 i 在 sum 中的值是多少？

假设总和扩展到所有i：

n = 5; %// matrix size. Set as needed
[ii jj] = ndgrid(1:n);
A = -1+2*max(ii,jj);
b = sum(fliplr(cumsum(fliplr(A),2)));