用matlab写dct算法水印嵌入

以下是一个简单的Matlab程序，演示了如何实现DCT算法水印嵌入：

% 原始图像文件名
filename = 'lena.bmp';
% 水印文件名
watermark = 'watermark.bmp';
% 嵌入后的图像文件名
outputfile = 'watermarked.bmp';

% 读入原始图像和水印图像
I = imread(filename);
w = imread(watermark);

% 将原始图像转换为灰度图像
if size(I,3) == 3
    I = rgb2gray(I);
end

% 将水印图像转换为二值图像
w = im2bw(w, 0.5);

% 将图像分成8x8的块，进行DCT变换
dctI = dct2(I);
dctW = dct2(w);

% 将水印嵌入DCT系数中
alpha = 0.1; % 嵌入强度
dctI(1:8,1:8) = dctI(1:8,1:8) + alpha * dctW;

% 将嵌入水印后的DCT系数进行反变换，得到嵌入水印后的图像
Iw = idct2(dctI);

% 将结果保存到文件中
imwrite(Iw, outputfile);

这个程序首先读入原始图像和水印图像，将原始图像转换为灰度图像，将水印图像转换为二值图像。然后将原始图像分成8x8的块，进行DCT变换，将水印嵌入DCT系数中，再进行反变换，得到嵌入水印后的图像，最后将结果保存到文件中。

需要注意的是，嵌入强度alpha的取值应该根据实际情况进行调整，过大的嵌入强度会导致图像失真，过小的嵌入强度则可能导致水印难以检测。