准则函数为最小：
∧
WfQf
∧
∧2
λgHf
∧2


2
9
式中λ叫拉格朗日系数。通过指定不同的Q，可以得到不同的复原目标。模拟实验结果如下：
f3LucyRichardso
算法LR算法是一种迭代非线性复原算法，它是从最大似然公式印出来的，图像用泊松分布加以模型化的。当下面这个迭代收敛时模型的最大似然函数就可以得到一个令人满意的方程：
fk1xyfkxyhxy
∧
∧
∧
gxyhxyfkxy
∧

10
代表卷积，f代表未退化图像的估计，g和h和以前定义一样。在IPT中，LR算法由名为deco
vlucy的函数完成的。模拟实验结果如下：
f4盲去卷积在图像复原过程中，最困难的问题之一是，如何获得PSF的恰当估计。那些不以PSF为基础的图像复原方法统称为盲区卷积。它以MLE为基础的，即一种用被随机噪声所干扰的量进行估计的最优化策略。工具箱通过函数deco
vbli
d来执行盲区卷积。模拟实验结果如下：图像分别迭代5次，10次，20次的结果。
f5、参考文献1罗军辉冯平等。MATLAB70在图像处理中的应用M，北京：机械工业出版社2006，202204，257278。2刘维一，于德月，王肇圻等。用迭代法消除数字图像放大后的模糊J。光电子激光，2002，13（4）：398400。3陶洪。数字共焦显微技术及其图像复原算法研究四川大学硕士学位论文CNKICDMD10610220036632。4冈萨雷斯等著。数字图像处理（MATLAB版）M，电子工业出版社2009。
f6、程序附录清单1维纳滤波clcclearallcloseallIimreadCDocume
tsa
dSetti
gsAdmi
istrator桌面原始图jpgubplot231ImshowIIrgb2grayIsubplot232imshowI
oise01ra
d
sizeIPSFfspecialmotio
2111BlurredimfilterIPSFcircularBlurredNoisyim2ui
t8BlurredNPabsfft
oise2NPOWsumNPprodsize
oiseNCORRfftshiftrealifft
NPIPabsfft
I2IPOWsumIPprodsize
oiseICORRfftshiftrealifft
IPICORR1ICORRceilsizeI12NSRNPOWIPOWsubplot233imshowBlurredNoisytitleABlurreda
dNoisysubplot234imshowdeco
vw
rBlurredNoisyPSFtitledeco
bw
rAPSFNSRsubplot235imshowdeco
vw
rBlurredNoisyPSFNCORRICORRtitledeco
bw
rAPSFNCORRICORRsubplot236imshowdeco
vw
rBlurredNoisyPSFNPOWICORR1titledeco
bw
rAPSFNPOWICORR_1_D2正则滤波clcclearallcloseallIimreadCDocume
tsa
dSetti
gsAdmi
istrator桌面原始图jpgsubplot231imshowIIrgb2grayIsubplot232imshowIPSFfspecialgaussia
710V01
fBlurredNoisyim
oiseimfilterIPSFgaussia
0VNOISEPOWERVprodsizeIJLAGRAdeco
vregBlurredNoisyPSFNOISEPOWERr