保留和识别数组中满足 MATLAB 阈值的元素

Question

我有一个由C_{ij} 表示的矩阵（大约 342 x 342），我想识别满足条件C_{ij} > rho 的所有索引 i,j，其中 rho 是某个固定值。我正在使用 MATLAB

例如，如果我将矩阵 C_{ij} 设为：

C = 1     0.7    0.8
    0.7   1      0.5
    0.8   0.5    1

而 rho = 0.6 那么满足条件 C_{ij} > 0.6 的索引是 i,j = 1,2 为 C_{11}=C_{22}=1 和 C_{12}=C_{21}=0.7 但请注意，i,j=3 不满足这个条件，因为虽然 C_{13}=C_{31}=0.8, C_{23}=C_{32}=0.5

我不确定在 MATLAB 中如何/最好的方法是？

Answer 1

C = rand(342,342);
rho = 0.6;
res = C(C>rho); %// contains all values that are above the threshold
[row,col,val] = find(res); %// returns the indices.

row 会包含你的 i，col 你的 j，但我会提防 using i or j as a variable。 val 包含相应的值，但您可以省略它以释放 RAM，因为您也可以使用 C(row,col) 来获取值。

如果你的矩阵是对称的，你可以使用它：

[rowt,colt,val] = find(triu(C)>rho); %// find only in the upper triangle
row = [rowt;colt]; %// flip rows and columns to obtain all results
col = [colt;rowt]; %// flip rows and columns to obtain all results

供您编辑：

IdxR = find(diff(unique(rows))==1,1,'first');
IdxC = find(diff(unique(columns))==1,1,'first');
Result = C(1:IdxR,1:IdxC);