请用matlab写出计算融合图像的AG（Average gradient），EI（Edge Intensity），SD（Standard deviation），SF（Spatial frequency）， （Gradient-based fusion performance）

以下是MATLAB代码实现：

% 读取融合图像
fusion_img = imread('fusion_image.jpg');

% 计算AG（Average gradient）
[Gx, Gy] = imgradient(fusion_img);
AG = mean(mean(sqrt(Gx.^2 + Gy.^2)));

% 计算EI（Edge Intensity）
EI = sum(sum(sqrt(Gx.^2 + Gy.^2)));

% 计算SD（Standard deviation）
SD = std2(fusion_img);

% 计算SF（Spatial frequency）
fft_img = fft2(fusion_img);
SF = sum(sum(abs(fft_img).^2)) / (size(fft_img,1)*size(fft_img,2));

% 输出结果
fprintf('AG: %.2f\n', AG);
fprintf('EI: %.2f\n', EI);
fprintf('SD: %.2f\n', SD);
fprintf('SF: %.2f\n', SF);

其中，AG表示平均梯度，EI表示边缘强度，SD表示标准差，SF表示空间频率。