自动对焦求傅里叶谱的matlab代码

以下是一个简单的自动对焦求傅里叶谱的 Matlab 代码示例：

% 加载图像
img = imread('test.jpg');

% 设定搜索范围
search_range = 50;

% 设定对焦窗口大小
window_size = 50;

% 获取图像大小
[height, width] = size(img);

% 初始化自相关矩阵
auto_corr = zeros(height, width);

% 计算自相关矩阵
for i = 1:height
    for j = 1:width
        auto_corr(i,j) = sum(sum((img - img(i,j)).^2));
    end
end

% 初始化傅里叶谱矩阵
fft_auto_corr = fft2(auto_corr);

% 计算傅里叶谱的幅值
fft_auto_corr = abs(fft_auto_corr);

% 初始化最优位置和最优值
best_pos = [0, 0];
best_val = 0;

% 在搜索范围内寻找最优位置
for i = -search_range:search_range
    for j = -search_range:search_range
        % 计算对焦窗口
        x1 = max(1, j);
        x2 = min(width, j + window_size);
        y1 = max(1, i);
        y2 = min(height, i + window_size);
        window = fft_auto_corr(y1:y2, x1:x2);
        
        % 计算窗口内的平均值
        window_val = mean(window(:));
        
        % 如果窗口内的平均值比当前最优值更优，则更新最优值和最优位置
        if window_val > best_val
            best_pos = [j, i];
            best_val = window_val;
        end
    end
end

% 输出最优位置和最优值
fprintf('Best focus position: (%d, %d)\n', best_pos(1), best_pos(2));
fprintf('Best focus value: %f\n', best_val);

请注意，这只是一个简单的示例代码，实际应用中可能需要进行更多的优化和改进。