写一段在图像DCT域嵌入图像水印的matlab程序

以下是一个简单的 MATLAB 程序，用于在图像的 DCT 域嵌入图像水印：

% 读取原始图像
img = imread('original_image.jpg');

% 将图像转换为灰度图像
gray_img = rgb2gray(img);

% 对灰度图像进行 DCT 变换
dct_img = dct2(double(gray_img));

% 设置水印信息
watermark = [1 0 1 0 1 1 0 1];

% 定义嵌入水印的位置
embed_positions = [2, 7, 10, 15, 20, 25, 30, 35];

% 嵌入水印信息
for i = 1:length(embed_positions)
    [row, col] = ind2sub(size(dct_img), embed_positions(i));
    dct_img(row, col) = dct_img(row, col) + watermark(i);
end

% 进行逆 DCT 变换
watermarked_img = idct2(dct_img);

% 显示原始图像和嵌入水印后的图像
subplot(1,2,1), imshow(gray_img), title('Original Image');
subplot(1,2,2), imshow(uint8(watermarked_img)), title('Watermarked Image');

% 保存嵌入水印后的图像
imwrite(uint8(watermarked_img), 'watermarked_image.jpg');

请确保将程序中的 `original_image.jpg` 替换为你想要嵌入水印的图像文件，并将 `watermarked_image.jpg` 替换为保存嵌入水印后图像的文件名。在程序中，我们使用了一个简单的水印信息 `[1 0 1 0 1 1 0 1]`，以及一个嵌入位置数组 `embed_positions`，你可以根据需要修改这些参数。

需要注意的是，这只是一个简单的示例程序，并不考虑一些实际应用中的细节和保护措施。在实际应用中，你可能需要更复杂的算法和技术来提高水印的鲁棒性和透明性，以及进行水印的检测和提取。