按照你之前的提示，>> untitled5 对于此运算，数组的大小不兼容。 出错 untitled5 (第 28 行) if size(medianFilter3x3_gaosi) ~= size(originalImage)，报错如何解决？

如果`medianFilter3x3_gaosi`的大小与`originalImage`不匹配，你需要调整滤波后的图像大小使其与原图像尺寸一致。这通常通过使用`imresize`函数实现，确保新生成的图像在行和列上都与原图像有相同的像素数量。以下是修复代码的一段示例：

% 如果滤波后的图像大小与原图像不符
if size(medianFilter3x3_gaosi) ~= size(originalImage)
    % 使用imresize函数调整滤波后的图像尺寸
    medianFilter3x3_gaosi = imresize(medianFilter3x3_gaosi, size(originalImage));
end

% 现在你可以安全地计算PSNR
psnr_value = psnr(medianFilter3x3_gaosi, originalImage);
fprintf('3x3 中值滤波后的PSNR: %.2f dB\n', psnr_value);

确保每次你想比较的图像都经过这样的校准步骤，然后才能正确计算PSNR值。

相关问题

对于此运算，数组的大小不兼容。 出错 Untitled1 (第 25 行) P_rec = P' * H' .* AF;

这个错误一般是因为数组的维度不匹配导致的。在这里，P'、H'和AF的维度可能不兼容，导致无法进行 .*（点乘）运算。你可以使用 size() 函数检查一下这三个数组的维度，看看它们是否匹配。如果维度不匹配，你需要对其中一个或多个数组进行调整，以使它们的维度相同。你可以使用 reshape() 函数或者其他数组操作函数来实现这个目标。

对于此运算，数组的大小不兼容。 出错 untitled11 (第 78 行) error(i) = sum(c~=a); % 统计误码数

这个错误通常是由于向量长度不匹配引起的。在这种情况下，您需要确保在进行向量运算时，所有的向量长度必须相同。

在您的代码中，`a` 和 `c` 的长度应该都是 `i`，但是在某些情况下，它们可能会因为计算错误而变得不同。您可以在计算 `error` 和 `BER` 之前添加以下代码，以确保向量长度相同：

a = a(1:length(c));

这将截取 `a` 的长度，使其与 `c` 的长度相同。这样，您就可以在不出现数组大小不兼容的情况下计算误码率了。

修改后的代码如下：

SNR = 0:0.1:10; % 信噪比范围
error = zeros(1,length(SNR)); % 误码率数组
for i = 1:length(SNR)
    fsk_noisy = awgn(e_fsk, SNR(i), 'measured');  % 加入噪声
    subplot(414); plot(t, fsk_noisy);  % 绘制加入噪声后的时域波形
    title('受干扰时域波形');
    %解调
    demod1=fsk_noisy.*s1;
    demod2=fsk_noisy.*s2;
    [b,a]=butter(6,2*fm/j,'low');
    df1=filter(b,a,demod1);%低通滤波
    df2=filter(b,a,demod2);%低通滤波
    z1=df1(j/i/2:j/i:j);
    z2=df2(j/i/2:j/i:j);

    %判决
    for k=1:i
        if z1(k)>z2(k)
            c(k)=1;
        else
            c(k)=0;
        end
    end

    %比较
    a = a(1:length(c)); % 确保向量长度相同
    error(i) = sum(c~=a); % 统计误码数
    BER(i) = error(i)/i; % 计算误码率
end

%绘制误码率曲线
figure(3);
semilogy(SNR, BER);
xlabel('信噪比(dB)');
ylabel('误码率');
title('2FSK误码率曲线');

希望这可以解决您的问题！