matlab筛选矩阵列,【转载】[Matlab]在矩阵中寻找满足条件的元素

参见张志涌精通matlab

>> A&＃61;randn(3,4)

A &＃61;

-0.5883    0.1139   -0.0956   -1.3362

2.1832    1.0668   -0.8323    0.7143

-0.1364    0.0593    0.2944    1.6236

以逻辑数组表达满足条件的元素&＃xff1a;

>>L&＃61;(A>0)

L &＃61;

0     1     0     0

1     1     0     1

0     1     1     1

%寻找符合条件元素位置(返回为一列&＃xff0c;并以单下标形式表示位置&＃xff0c;故转置为一行)

>>s&＃61;(find(L))&＃39;或者

>> s&＃61;(find(A>0))&＃39;

s &＃61;

2     4     5     6     9    11    12

可以对满足条件的元素集体赋值

>> A(s)&＃61;0

A &＃61;

-0.5883         0   -0.0956   -1.3362

0         0   -0.8323         0

-0.1364         0         0         0

如果要以双下标形式返回符合条件的元素位置&＃xff1a;

>> [ii,jj]&＃61;find(A&＃61;&＃61;0);

>> disp(ii&＃39;),disp(jj&＃39;)

2     1     2     3     3     2     3

1     2     2     2     3     4     4

最后&＃xff0c;与arrayfun综合运用的一个例子&＃xff0c;问题见于simwe forum&＃xff0c;想了下可以采用下面的两种方式来完成&＃xff1a;

请问如何查找矩阵中相同向量的位置&＃xff1f;

比如矩阵A&＃xff1d;【1&＃xff0c;2&＃xff0c;3&＃xff0c;1&＃xff1b;

2&＃xff0c;3&＃xff0c;4&＃xff0c;2&＃xff1b;

3&＃xff0c;3&＃xff0c;8&＃xff0c;3】

向量为b &＃xff1d;【1&＃xff1b;2&＃xff1b;3】

如何能不用循环找出b在A矩阵列中的位置&＃xff0c;loc&＃xff1d;【1 4】呢&＃xff1f;

方法一&＃xff1a;

A&＃61;[1,2,3,1;2,3,4,2;3,3,8,3];

B&＃61;[1;2;3];

C&＃61;find((ismember(A&＃39;,B&＃39;,&＃39;rows&＃39;))&＃39;);

方法二&＃xff1a;

A&＃61;[1,2,3,1;2,3,4,2;3,3,8,3]

B&＃61;[1;2;3];

[M,N]&＃61;size(A);

N&＃61;1:N;

P&＃61;find(arrayfun(&＃64;(n) isequal(A(:,n),B),N));

版上回复的其它方法&＃xff1a;

方法1&＃xff1a;find(all(bsxfun(&＃64;eq,A,B)))

方法2&＃xff1a;find(arrayfun(&＃64;(n) all(A(:,n)&＃61;&＃61;B),1:4))

方法3&＃xff1a; find(all(~bsxfun(&＃64;minus,A,B)))

方法4&＃xff1a; find(arrayfun(&＃64;(n) isequal(A(:,n),B) ,1:4))

方法5&＃xff1a;find(B&＃39;*A&＃61;&＃61;sum(B.^2))

bsxfun和arrayfun函数非常有用&＃xff0c;建议楼主参考帮助文档了解其用法。方法5是用了点数学方法&＃xff0c;有时候我们多分析下问题会得到意想不到的结果。

也可以不那么复杂&＃xff0c;我只会些简单的函数命令(其实和你的方法3一个思路)

%% 矩阵中寻找向量

clc

clear all

A&＃61;[1 2 3;1 2 3;2 3 4;3 4 5]&＃39;

B&＃61;repmat([1:3]&＃39;,1,4)

D&＃61;find(~all(B-A,1))

%B &＃61; all(A, dim) tests along the dimension of A specified by scalar dim.

%dim&＃61;1&＃xff1a;每一列都为非0&＃xff0c;该列取1&＃xff0c;否则取0&＃xff0c;缺省&＃xff1b;dim&＃61;2&＃xff1a;行&＃xff1b;dim&＃61;3&＃xff1a;页&＃xff0c;以此类推