请给我一段基于区域能量融合的灰度可见光图像与红外灰度图像融合的matlab代码，窗口区域设置为3×大小

好的，以下是基于区能量融合的灰度可见光图像与红外灰度图像融合的 MATLAB 代码：

% 读取可见光图像和红外图像
visibleImg = imread('visible.jpg');
infraredImg = imread('infrared.jpg');

% 将图像转换为灰度图像
visibleGray = rgb2gray(visibleImg);
infraredGray = rgb2gray(infraredImg);

% 设置窗口大小
windowSize = 3;

% 计算区能量
visibleEnergy = energyFunc(visibleGray, windowSize);
infraredEnergy = energyFunc(infraredGray, windowSize);

% 计算融合图像
alpha = 0.5;
fusionImg = alpha * visibleGray + (1 - alpha) * infraredGray;

% 显示融合图像
imshow(fusionImg);

% 定义区能量计算函数
function energy = energyFunc(img, windowSize)
    [m, n] = size(img);
    energy = zeros(m, n);
    halfWindow = floor(windowSize/2);
    for i = halfWindow+1:m-halfWindow
        for j = halfWindow+1:n-halfWindow
            % 计算当前窗口中像素点的灰度值平均值
            window = img(i-halfWindow:i+halfWindow, j-halfWindow:j+halfWindow);
            meanValue = mean(window(:));
            % 计算当前窗口中像素点的能量
            energy(i, j) = sum(sum((window - meanValue).^2));
        end
    end
end

以上代码中，`energyFunc` 函数用于计算图像的区能量，`fusionImg` 变量用于存储融合后的图像，`alpha` 变量用于控制可见光图像和红外图像的权重。你可以根据自己的需要进行调整。