用matlab编写一段dct变换提取水印的程序

以下是一个简单的MATLAB程序，用于在DCT变换后提取图像水印：

% 读取图像和水印
img = imread('image.jpg');
watermark = imread('watermark.jpg');

% 将水印调整为与图像大小相同
watermark = imresize(watermark, size(img));

% 将图像和水印转换为灰度图像
img_gray = rgb2gray(img);
wm_gray = rgb2gray(watermark);

% 对图像和水印进行DCT变换
img_dct = dct2(img_gray);
wm_dct = dct2(wm_gray);

% 将水印嵌入到DCT系数的LSB中
alpha = 0.1; % 水印强度
wm_dct_lsb = bitget(wm_dct, 1); % 获取水印的LSB
img_dct_lsb = img_dct; % 复制DCT系数
img_dct_lsb(1:size(wm_dct,1), 1:size(wm_dct,2)) = img_dct_lsb(1:size(wm_dct,1), 1:size(wm_dct,2)) + alpha * wm_dct_lsb; % 将水印嵌入到DCT系数的LSB中

% 反DCT变换
img_wm = idct2(img_dct_lsb);

% 显示图像和提取的水印
figure;
subplot(1,2,1), imshow(img), title('Original Image');
subplot(1,2,2), imshow(img_wm), title('Watermarked Image');

该程序将图像和水印转换为灰度图像，并对它们进行DCT变换。然后，将水印嵌入到DCT系数的LSB中，使用一个名为alpha的参数来控制水印的强度。最后，使用反DCT变换生成带有嵌入水印的图像，并显示原始图像和提取的水印。