Filter matrix rows based on value of function in Matlab

5 vues (au cours des 30 derniers jours)
ARS
ARS le 18 Fév 2013
Commenté : yair suari le 1 Sep 2015
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

Réponse acceptée

Mark Whirdy
Mark Whirdy le 18 Fév 2013
Modifié(e) : Mark Whirdy le 18 Fév 2013
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)) , :)

Plus de réponses (1)

Azzi Abdelmalek
Azzi Abdelmalek le 18 Fév 2013
n=size(A,1);
A(arrayfun(@(x) expression(A(x,:)),(1:n)','un',0))=[]
  5 commentaires
ARS
ARS le 18 Fév 2013
Sorry if I'm misunderstanding, but I receive the error:
Error using subsindex
Function 'subsindex' is not defined for values of class 'cell'.
Azzi Abdelmalek
Azzi Abdelmalek le 18 Fév 2013
Modifié(e) : Azzi Abdelmalek le 18 Fév 2013
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.

Connectez-vous pour commenter.

Catégories

En savoir plus sur Logical dans Help Center et File Exchange

Produits

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by