【MATLAB】 求矩阵中椭圆形范围内最大值及其行列号

---

前言

在矩阵中给定椭圆形范围内，求最大矩阵值，及其在矩阵中的行列号
注：矩阵中先行号后列号，图形中先x（列表示）后y（行表示）

计算函数

输入椭圆圆心坐标、长半轴、短半轴（均为整数）及矩阵M
输出椭圆内最大值及其行列号

注：本文计算的椭圆方程为焦点在x轴的版本

代码

求最值的代码如下（示例）：
新建m文件：

function[height, xmax, ymax] = TuoMax(ra, rb, x0, y0, M)
% b^2*x^2+a^2*y^2=a^2*b^2
% index1-行号-y; index2-列号-x
D = M(y0-rb:y0+rb,x0-ra:x0+ra);
while 1
    [maxcolumn,indexcolumn] = max(D);
    [Dmax,index2] = max(maxcolumn);
    index1 = indexcolumn(index2);
    ymax = index1+y0-rb-1;
    xmax = index2+x0-ra-1;
    veri = (xmax-x0)^2*rb^2+(ymax-y0)^2*ra^2-ra^2*rb^2;
    if veri<=0
        height = Dmax;
        break
    else
        D(index1,index2)=0;
        fprintf('错误结果 %f\n',xmax,ymax);
    end
end
end

测试及调用时在命令行输入（示例）：

clear;clc; 

M = rand(2402,2402);
M = M.*4000+5000;
ra=100;rb=20;x0=1500;y0=2000;
[height, xmax, ymax] = TuoMax(ra, rb, x0, y0, M)

---

欢迎在评论区讨论或补充说明。