生成用于彩色模糊图像复原的维纳滤波、逆滤波、lucy-richardson的matlab代码

以下是使用MATLAB实现的彩色模糊图像复原的维纳滤波、逆滤波、lucy-richardson的代码：

function img_restore = weiner_filter(img, psf, snr)
    [H, W, ~] = size(img);
    psf = psf / sum(psf(:));
    psf_pad = zeros(H, W);
    psf_pad(1:size(psf, 1), 1:size(psf, 2)) = psf;
    psf_shift = fftshift(psf_pad);
    img_fft = fft2(img);
    psf_fft = fft2(psf_shift);
    psf_fft_abs = abs(psf_fft) .^ 2;
    snr = 10 ^ (snr / 10);
    weiner_filter = conj(psf_fft) ./ (psf_fft_abs + snr ./ img_fft);
    img_restore = real(ifft2(img_fft .* weiner_filter));
    img_restore = uint8(max(min(img_restore, 255), 0));
end

function img_restore = inverse_filter(img, psf)
    [H, W, ~] = size(img);
    psf = psf / sum(psf(:));
    psf_pad = zeros(H, W);
    psf_pad(1:size(psf, 1), 1:size(psf, 2)) = psf;
    psf_shift = fftshift(psf_pad);
    img_fft = fft2(img);
    psf_fft = fft2(psf_shift);
    inverse_filter = conj(psf_fft) ./ (abs(psf_fft) .^ 2 + eps);
    img_restore = real(ifft2(img_fft .* inverse_filter));
    img_restore = uint8(max(min(img_restore, 255), 0));
end

function img_restore = lucy_richardson(img, psf, n_iter)
    [H, W, ~] = size(img);
    psf = psf / sum(psf(:));
    psf_pad = zeros(H, W);
    psf_pad(1:size(psf, 1), 1:size(psf, 2)) = psf;
    psf_shift = fftshift(psf_pad);
    img_restore = ones(size(img));
    for i = 1:n_iter
        error = real(ifft2(fft2(img_restore) .* conj(fft2(psf_shift))) - img);
        error = ifftshift(error);
        error_fft = fft2(error);
        img_restore = real(ifft2(conj(fft2(psf_shift)) .* error_fft));
        img_restore = max(min(img_restore, 255), 0);
    end
    img_restore = uint8(img_restore);
end

其中，`weiner_filter`函数实现了维纳滤波，`inverse_filter`函数实现了逆滤波，`lucy_richardson`函数实现了逆滤波lucy-richardson，`img`为输入的彩色模糊图像，`psf`为点扩散函数，`snr`为信噪比，`n_iter`为迭代次数。函数返回复原后的彩色图像。注意，这里使用了`uint8`类型，因此输出的图像像素值范围在[0, 255]之间。