高斯白噪声加噪声MATLAB

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，强度为3的高斯白噪声
N=imabsdiff（n img）;%图像减法运算提取高斯白噪声
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;