clear?all
close?all
clc
[fnpnfi]=uigetfile（‘*.jpg‘‘選擇原始圖像‘）;
I=imread（[pn?fn]）;
img=I;
img=rgb2gray?（img）;%灰度化
img=double?（img）;
a=input（‘請輸入密鑰：‘）;
c=1/10^abs（9-a）;
figure;
subplot（221）;
imshow?（uint8?（img））;
title（‘原始的圖像‘）;
[pq]=size（img）;
%生成密鑰
n=awgn（img?4?3‘linear‘）;%在img中添加信噪比為4，強(qiáng)度為3的高斯白噪聲
N=imabsdiff（n?img）;%圖像減法運(yùn)算提取高斯白噪聲
subplot（222）;
imshow（double（N））;
title（‘密匙‘）;
[LoHiLo_R?Hi_R]=wfilters（‘haar‘）;%生成haar小波的濾波器
[ca?ch?cv?cd]=dwt2?（img?Lo?Hi）;%圖像dwt分解
y=[ca?ch;cv?cd];
r=N（:599）;
N=[Nr];
Y=y+c*abs?（y）.*N;