Locating global max in 3D plot

2021 1 22
1 个回答
回答已采纳
更新时间：2021 1 23
8 次查看（30 天）

0 个投票

I'm working on some MATLAB code but I have encountered a problem that has hindered my progress.
Kindly would you assist me on how I may locate a global maxima value in a case where I have more than one maximas in a 3D plot.
I'm using this
Idx = find((Pmusic(:) == max(max(Pmusic(:)))));
[PmusicmaxRow,PmusicmaxCol] = ind2sub(size(Pmusic), Idx);
but it brings all the maximas but I am only interested with the indices of the global maxima.

请先登录，再回答此问题。

 采纳的回答

Star Strider
Star Strider 2021-1-22

0 个投票

Try this:
[maxval,Idx] = max(Pmusic(:));
[PmusicmaxRow,PmusicmaxCol] = ind2sub(size(Pmusic), Idx);
I obviously cannot test this with your matrix, however it worked correctrly when I tested it with my test matrix.

4 个评论
显示 2更早的评论 隐藏 2更早的评论

John Paul
John Paul 2021-1-22
Thank you so much sir for your expert advice. I am trying it but unfortunately, some local maximas are still being considered. I am attaching my code for your kind reference, the code is for determining 2 dimensional direction of arrival (DOA) using MUSIC algorithm. Generally, the estimator part is working well. I am supposed to use the estimated values that are the indices of the maximum spectrum value to calculate RMSE. Kindly, if you may spare sometime to advice me on what the problem maybe.
Star Strider
Star Strider 2021-1-23
The max function will only return the first maximum (and the index), in the event that there are more than one. However here, there apopears to be only one maximum.
When I run:
[maxval,Idx] = max(Pmusic(:));
[PmusicmaxRow,PmusicmaxCol] = ind2sub(size(Pmusic), Idx);
maxidxv = find(Pmusic == maxval);
the only maximum value is 18.7194, and the index is 8752. The find call (that I added here temporarily) also only detects the same one maximum, with the same index. The ‘aaa’ and ‘ccc’ vectors are simply repeats of the row and column indices, not new or different indices. So ‘aaa’ is a (1x10) vector of 63, and ‘ccc’ is a (1x10) vector of 48.
There are no local maxima, only one global maximum.
I did not examine all your code, only the part you have questions about, those lines 126 and 127 of the file you attached.
John Paul
John Paul 2021-1-23
Respected Sir, You are absolutely correct. This section of the code is supposed to genereate indices for the global maximas for 10 iterations and that is the reason ‘aaa’ and ‘ccc’ are (1x10) vectors. Thank you so much for your invalauble advice. You have made me move a step further with my work.
Star Strider
Star Strider 2021-1-23
As always, my pleasure!
I am always happy to help!

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Beamforming and Direction of Arrival Estimation 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by