利用伪随机数产生水印，用DFT将水印嵌入彩色图像的相位区域的matlab代码

以下是利用伪随机数产生水印，用DFT将水印嵌入彩色图像的相位区域的matlab代码：

% 读取彩色图像
img = imread('color_image.jpg');
imshow(img);
title('Original Image');

% 将彩色图像转换为灰度图像
gray_img = rgb2gray(img);
imshow(gray_img);
title('Gray Image');

% 生成水印
watermark = randi([0 1], [size(gray_img, 1) size(gray_img, 2)]);
imshow(watermark);
title('Watermark');

% 对水印进行DFT
dft_watermark = fft2(watermark);

% 对图像进行DFT
dft_img = fft2(double(gray_img));

% 将水印嵌入图像的相位区域
dft_img_phase = exp(1i * angle(dft_img));
dft_img_watermark = abs(dft_img) .* dft_img_phase .* exp(1i * angle(dft_watermark));
img_watermarked = real(ifft2(dft_img_watermark));

% 显示嵌入水印后的图像
imshow(uint8(img_watermarked));
title('Watermarked Image');

其中，`randi([0 1], [size(gray_img, 1) size(gray_img, 2)])`用于生成与灰度图像大小相同的随机01矩阵作为水印，`fft2()`用于进行二维DFT，`exp(1i * angle(dft_img))`用于获取图像的相位信息，`abs(dft_img)`用于获取图像的幅度信息，`.*`用于进行矩阵对应元素相乘，`real(ifft2(dft_img_watermark))`用于进行二维IDFT，最终得到嵌入水印后的图像。