Filter matrix rows based on value of function in Matlab
6 次查看(过去 30 天)
显示 更早的评论
This seems like a simple question but I have been unable to find an answer anywhere. If I have a Matlab matrix "A" consisting of an arbitrary number of rows, how would I filter these rows based on the value of some function "f" (the argument of which is a row vector)? In other words, how would I keep only the rows of matrix A for which f is true? I tried
A(f(A(:)), :)
but to no success. Any help would be greatly appreciated.
So for example, if I have the matrix:
A =
1 1
1 2
1 3
1 4
2 1
2 2
2 3
2 4
3 1
3 2
3 3
3 4
4 1
4 2
4 3
4 4
And a function f that returns true only for inputs [1 2], [1 3], [2 4] and [3 4], my desired output would be:
ans =
1 2
1 3
2 4
3 4
采纳的回答
Mark Whirdy
2013-2-18
编辑:Mark Whirdy
2013-2-18
Shown below using an anonymous function for neatness, but the trick is just to make sure that your function (anonymous, nested, sub, or separate) returns a LOGICAL VECTOR (or a vector of linear indexes)
fn = @(x)(x>1);
a = randn(20,1);
fn(a); % this would return a logical vector if you ran it
myfilteredVector = a(fn(a));
If I've misunderstood, could you give me an example with an input & a required output.
------
So the bones of the required operation is actually a row-wise matching? I do this by transforming each row into a single unique element by 10*A(:,1) + A(:,2), as below. Then ismember() with its mex engine can then be put to work on matching individual element. There are other matching/indexing approaches though but the primary aim is to have the function outputting a logical vector.
A = [...
1, 1;
1, 2;
1, 3;
1, 4;
2, 1;
2, 2;
2, 3;
2, 4;
3, 1;
3, 2;
3, 3;
3, 4;
4, 1;
4, 2;
4, 3;
4, 4;];
fn = @(x)(ismember(x,[12;13;24;34]));
% fn(10*A(:,1) + A(:,2)) would return a logical vector
A(fn(10*A(:,1) + A(:,2)),:); % logical vector selects rows, & we want every column so (logic, :)
You can of course also specify the matchset a=[12;13;24;34] as an argument as below.
fn = @(x,a)(ismember(x,a));
a = [12;13;24;34];
A(fn(10*A(:,1) + A(:,2)),:);
To be honest, your function is kind of unnecessary really (unless there's a reason outside of the example you've given), as you can just use ismember with the unique encoder 10*A(:,1)+A(:,2).
a = [1,2;1,3;2,4;3,4];
A( ismember(10*A(:,1)+A(:,2), 10*a(:,1)+a(:,2)) , :)
更多回答(1 个)
Azzi Abdelmalek
2013-2-18
n=size(A,1);
A(arrayfun(@(x) expression(A(x,:)),(1:n)','un',0))=[]
5 个评论
Azzi Abdelmalek
2013-2-18
编辑:Azzi Abdelmalek
2013-2-18
Try
n=size(A,1)
A(cell2mat(arrayfun(@(x) f(A(x,:)),(1:n)','un',0)),:)=[]
Also this is not a general solution, each expression has to be treated individually.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!