matlab最小二乘法滤波

 clear all;
I=imread（‘peppers.png‘）;
I=I（125+[1:256]1:256:）;    %图像剪切
subplot（231）;imshow（I）;
xlabel（‘（a）原始图像‘）;
%模拟运动模糊和噪声
LEN=11;                      %设置长度
THETA=5;                    %设置角度
PSF=fspecial（‘motion‘LENTHETA）;  %生成滤波器
blurred=imfilter（IPSF‘conv‘）;   %图像卷积运算
subplot（232）;imshow（blurred）;
xlabel（‘（b）模糊图像‘）;
%添加高斯噪声
V=0.02;
BlurredNoise=imnoise（blurred‘gaussian‘0V）;   %添加高斯白噪声
subplot（233）;imshow（BlurredNoise）;
xlabel（‘（c）带高斯噪声模糊图像‘）
%使带高斯噪声模糊图像复原
NP=V*prod（size（I））;        %噪声功率
[J1LAGRA]=deconvreg（BlurredNoisePSFNP）; %图像复原
subplot（234）;imshow（J1）;
xlabel（‘（d）使用真实NP复原‘）
J2=deconvreg（BlurredNoisePSFNP*1.3）;   %采用放大1.3倍的NP
subplot（235）;imshow（J2）;
xlabel（‘（e）采用放大1.3倍的NP复原‘）
%使用小于真实NP值的复原
J3=deconvreg（BlurredNoisePSFN