对于此运算，数组的大小不兼容。 出错 untitled (第 18 行) Fw = 4 * (F / 4) + w / 64; % 添加水印

非常抱歉，我在代码中犯了一个错误。由于数组大小不兼容，导致无法执行水印添加操作。请使用以下更新后的代码：

% 读取原始灰度图像
image = imread('原始图像.jpg');
image = rgb2gray(image);

% 读取水印图像
watermark = imread('水印图像.jpg');
watermark = rgb2gray(watermark);
watermark = imresize(watermark, size(image)); % 调整水印图像大小与原始图像相同

% 将原始图像和水印图像转换为double类型
image = im2double(image);
watermark = im2double(watermark);

% 添加水印
F = fft2(image); % 对原始图像进行傅里叶变换
w = fft2(watermark); % 对水印进行傅里叶变换
Fw = F + w / 64; % 添加水印
watermarked_image = real(ifft2(Fw)); % 对添加水印后的图像进行逆傅里叶变换

% 提取水印
extracted_watermark = (Fw - F) * 64; % 提取水印

% 高质量JPEG压缩和解压缩
imwrite(watermarked_image, '加过水印的图像.jpg', 'Quality', 100);
compressed_image = imread('加过水印的图像.jpg');
decompressed_image = im2double(imread('加过水印的图像.jpg'));

% 从压缩后的图像中提取水印
F_compressed = fft2(decompressed_image);
extracted_watermark_compressed = (Fw - F_compressed) * 64;

% 显示结果
figure;
subplot(2,2,1), imshow(watermarked_image), title('加过水印的图像');
subplot(2,2,2), imshow(extracted_watermark), title('提取的水印');
subplot(2,2,3), imshow(compressed_image), title('经过JPEG压缩和解压缩后的图像');
subplot(2,2,4), imshow(extracted_watermark_compressed), title('从JPEG压缩后的图像中提取的水印');

这次更新后的代码应该可以正常运行了。如果您还有任何问题，请随时提问。非常抱歉给您带来了困扰！