%%%%%約束最小二乘方濾波%%%%%
clear;close?all;clc;
gk=imread（‘E:\學習\A+數字圖像處理\作業4\ImageBlur_2.bmp‘）;
[MNK]=size（gk）;
gk=double（gk）;
g=gk（::1）*0.299+gk（::2）*0.587+gk（::3）*0.114;%Gray?=?R*0.299?+?G*0.587?+?B*0.114
%%?加入高斯噪聲
%%%構造一個高斯噪聲矩陣
mu=0;sigma=1;
noise=normrnd（musigmaMN）;

%%%?構造一個矩陣temp隨機選出20%的像素用于加噪聲
temp=zeros（MN）;%%%默認每個元素都為0，即不加噪聲
amount=round（M*N*0.2）;?%%%加噪像素總量
for?j=1:amount
????x=randi（M11）;
????y=randi（N11）;
????temp（xy）=1;%%%將選中的像素點的值置為1，對這些像素點加噪聲
end

%%%加高斯噪聲
gn=g+temp.*noise;
gnoise=gn-g;
gnn=uint8（gn）;
hold?on
imshow（gnn）;
imwrite（gnn‘E:\學習\A+數字圖像處理\作業4\ImageBlurN.bmp‘）;

G=fftshift（myft2mn（gn））;%%?二維離散傅立葉變換
%%?退化函數
T=1;?a=0.1;?b=0.1;
H=zeros（MN）;
for?u=1:M
????for?v=1:N
????????u1=u-M/2-1;%將H（uv）的（00）點移到矩陣中心
????????v1=v-N/2-1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-04-15?21:44??圖像運動模糊&約束最小二乘方濾波\
?????文件????????1913??2018-11-18?20:12??圖像運動模糊&約束最小二乘方濾波\CLSF.m
?????文件?????????612??2018-11-18?15:52??圖像運動模糊&約束最小二乘方濾波\ImageBlur.m
?????文件???????91078??2018-11-18?18:57??圖像運動模糊&約束最小二乘方濾波\ImageBlurN.bmp
?????文件??????270054??2018-11-18?15:55??圖像運動模糊&約束最小二乘方濾波\ImageBlur_2.bmp
?????文件??????270054??2018-11-18?17:40??圖像運動模糊&約束最小二乘方濾波\ImageDeblur.bmp
?????文件?????????473??2018-11-18?16:11??圖像運動模糊&約束最小二乘方濾波\myft2mn.m
?????文件?????????533??2018-11-18?16:11??圖像運動模糊&約束最小二乘方濾波\myift2mn.m
?????文件??????237170??2019-04-15?21:44??圖像運動模糊&約束最小二乘方濾波\實驗報告.pdf