5

我需要一个代码来省略矩阵的对角元素,例如,如果

A =

[1 2 3;
 1 2 3;
 1 2 3];

输出来:

[2 3;
 1 3;
 1 2];

我怎样才能简单地做到这一点(我知道很长,但我需要它简单)

4

3 回答 3

5

这是一个解决方案:

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
于 2012-07-29T14:22:44.263 回答
3

n请注意,通过nmatirx消除 a 的对角线后有两种可能性:

  1. 如果后果矩阵是nn-1就像你的问题一样),你可以这样做:

    A=A';
    A(1:n+1:n*n)=[];
    A=reshape(A,n-1,n)';
    
  2. 如果余波矩阵是n-1by n,你可以这样做:

    A(1:n+1:n*n)=[];
    A=reshape(A,n-1,n);
    
于 2012-07-29T14:30:07.560 回答
0

这是另一种方式

reshape(A(setdiff(1:9,1:4:9)),[3,2])
于 2012-07-29T14:26:12.510 回答