请用matlab写出求融合图像交叉熵，熵，互信息，峰值信噪比的程序

以下是一些示例代码，用于在MATLAB中计算融合图像的交叉熵、熵、互信息和峰值信噪比。

1. 计算交叉熵：

% 假设有两幅图像 A 和 B，它们的融合图像为 C
C = imread('fused_image.png');
A = imread('image1.png');
B = imread('image2.png');

% 计算交叉熵
H_AB = -sum(sum(C.*log(A+B)))/numel(C);

2. 计算熵：

% 计算熵
H_C = -sum(sum(C.*log(C)))/numel(C);

3. 计算互信息：

% 计算联合概率
joint_prob = histcounts2(A(:),B(:))/numel(A);

% 计算边缘概率
p_A = histcounts(A(:))/numel(A);
p_B = histcounts(B(:))/numel(B);

% 计算互信息
H_A = -sum(p_A.*log(p_A));
H_B = -sum(p_B.*log(p_B));
I_AB = H_A + H_B - H_AB;

4. 计算峰值信噪比：

% 计算信号的能量
signal_energy = sum(sum(C.^2));

% 计算噪声的能量
noise_energy = sum(sum((A-C).^2)) + sum(sum((B-C).^2));

% 计算峰值信噪比
PSNR = 10*log10(signal_energy/noise_energy);

请注意，这些示例代码仅作为参考。实际应用中，您可能需要根据自己的数据格式和计算要求进行修改。