写一个用于激光光束质量的 matlab程序

编写一个用于激光光束质量分析的MATLAB程序，通常涉及到计算光束的几个关键参数，如M2系数（光束发散度）、稳定性（波前误差）以及光强分布等。以下是一个简单的示例，演示如何创建一个基础的函数来评估光束的质量：

function [M2, StrehlRatio] = laserBeamQuality(Ix, Iy)
% M2系数计算
[height, width] = size(Ix); % 获取图像尺寸
sigma_x = std(Ix(:)); % 横向光斑的标准偏差
sigma_y = std(Iy(:)); % 纵向光斑的标准偏差
M2 = sqrt((sigma_x^2 + sigma_y^2) / (2 * mean(Ix(:))^2));

% Strehl Ratio（斯特列尔比）计算
if max(Ix(:)) > 0 % 如果有非零光强度
    background = mean(Ix(:)) - median(Ix(:)); % 背景值
    peakIntensity = max(max(Ix)); % 光峰强度
    StrehlRatio = peakIntensity / (pi * (sigma_x^2 + sigma_y^2));
else
    StrehlRatio = NaN; % 如果光强度全为零，无法计算
end
end

% 示例使用
Ix_data = randn(512, 512); % 创建一个模拟的激光光斑图像
Iy_data = randn(512, 512);
[M2_result, StrehlRatio_result] = laserBeamQuality(Ix_data, Iy_data);

% 输出结果
disp(['M2 Coefficient: ', num2str(M2_result)]);
disp(['Strehl Ratio: ', num2str(StrehlRatio_result)]);

%