Hi Matthias,
I understand that you want to find values in the “experimental” dataset that fall within a particular tolerance of the elements in the “simulated” dataset. You are using the ismembertol function to do the same:
>>CommonValues = ismembertol(experimental,simulated,tolerance)
This function will return a logical array “CommonValues” the same size as “experimental”. Since your “experimental” array is a column vector, the return value will also be a column vector. The location of “1”s in this vector will correspond to elements in “experimental” array that are within tolerance of elements in “simulated”.
If you need the position of matched elements in the “experimental” array, you can use the “find” function:
>>expRows = find(CommonValues==1);
This will give you a vector “expRows” containing the locations of elements “experimental” which are within a tolerance of elements in “simulated” (These will be the row numbers from the “experimental” as it is a column vector).
You can also get the actual elements in “experimental” using the below command:
>> expValues = experimental(CommonValues);
This will give you a vector with elements from “experimental” that matched to elements in “simulated”.
If you need to get the index location in “simulated” for each element in “experimental” that is a member of “simulated”, you can specify an array in the output parameters:
>>[ CommonValues, locSimulated] = ismembertol(experimental,simulated,tolerance);
“locSimulated” will be an array the same size as “experimental” (unless the “ByRows” option is specified) and will contain indices to the elements in “simulated” that are found in “experimental” (within tolerance). It will contain a 0 whenever an element in “experimental” is not a member of “simulated”. Note that if there are multiple matches in “simulated” for a single entry in “experimental”, only the first matched location in “simulated” would be returned. The contents will be the indices of “simulated” as a linearly indexed array. Please refer to the below link for details on linear indexing:
These linear indices can be converted to subscripts (for example, rows and columns in the case of “simulated” being a 2D matrix) using the ind2sub function.
You can get the elements in “simulated” using the above indices by the following command:
>> simulated(locSimulated (locSimulated ~=0))
This will first get all the indices in “simulated” that are not equal to 0 (that is, the ones that are within tolerance) and will return values in “simulated” based on the non-zero positions.
Hope that helps!
-Shruti
