Index of maximum values in accumarray
3 次查看(过去 30 天)
显示 更早的评论
Hello all,
I'm working with wind data that is measured every hour, the data is in the form of a matrix of 3 columns, the first is the date, the second is wind speed and the third is wind direction. I have created a code using accumarray in order to get the maximum daily value of wind speed. However I want to be able to get the corresponding wind direction for those maximum values, so for this, I need the indexes of the Maximum values.
I have tried the suggestion from here given by Walter. However, the resulting indexing shows a different result from the one I need (numerical example below):
My code is the following:
max_val = accumarray(c,wind(:,2),[],@nanmax);
max_index= accumarray(c,wind(:,2),[],@max_and_idx);
where:
function idx = max_and_idx(x)
[~,idx] = max(x);
end
the answer I get is:
max_val =[4.5 5.1 3.6 2.5 .....]
max_index=[16 15 15 12 .....]
I checked my data and 4.5 corresponds to position #15 in day 1, 5.1 corresponds to position #16 on day 2 and so on. What I want is the actual position of these numbers within the entire wind matrix. How can I fix this problem?
0 个评论
采纳的回答
Ameer Hamza
2018-5-20
编辑:Ameer Hamza
2018-5-22
A better approach might be to use the splitapply to divide the data based on days. Then the better will be better organized and you can easily process for required results. For example
splitData = splitapply(@(x) {x}, wind(:,2:3), findgroups(c));
now you can search splitData to find the max value and the corresponding direction,
result = cellfun(@(x) x(x(:,1)==max(x(:,1)), :), splitData, , 'UniformOutput', 0);
the result will contain the maximum daily value and wind direction.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Matrices and Arrays 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!