我需要一个代码来省略矩阵的对角元素,例如,如果
A =
[1 2 3;
1 2 3;
1 2 3];
输出来:
[2 3;
1 3;
1 2];
我怎样才能简单地做到这一点(我知道很长,但我需要它简单)
这是一个解决方案:
Alower = tril(A, -1);
Aupper = triu(A, 1);
result = Alower(:, 1:end-1) + Aupper(:, 2:end)
演示:
> A = [1 2 3; 1 2 3; 1 2 3]
A =
1 2 3
1 2 3
1 2 3
> tril(A, -1)(1:end, 1:end-1) + triu(A, 1)(1:end, 2:end)
ans =
2 3
1 3
1 2
n
请注意,通过n
matirx消除 a 的对角线后有两种可能性:
如果后果矩阵是n
(n-1
就像你的问题一样),你可以这样做:
A=A';
A(1:n+1:n*n)=[];
A=reshape(A,n-1,n)';
如果余波矩阵是n-1
by n
,你可以这样做:
A(1:n+1:n*n)=[];
A=reshape(A,n-1,n);
这是另一种方式
reshape(A(setdiff(1:9,1:4:9)),[3,2])