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

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)

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

你可能感兴趣的:(matlab,矩阵,开发语言)