matlab低通滤波器库函数代码_如何在Matlab中实现低通巴特沃斯滤波器？

From this answer, I know how to create a High-pass Butterworth filter.

From this video, I know that, lowpasskernel = 1 - highpasskernel.

So, I created the following Low-pass Butterworth Filter,

function [out, kernel] = butterworth_lp(I, Dl, n)

height = size(I, 1);

width = size(I, 2);

[u, v] = meshgrid(-floor(width/2):floor(width/2)-1,-floor(height/2):floor(height/2)-1);

% lp_kernel = 1 - hp_kernel

kernel = 1 - butter_hp_kernel(u, v, Dl, n);

% fft the image

I_fft_shifted = fftshift(fft2(double(I)));

% apply lowpass filter

I_fft_shift_filtered = I_fft_shifted .* kernel;

% inverse FFT, get real components

out = real(ifft2(ifftshift(I_fft_shift_filtered)));

% normalize and cast

out = (out - min(out(:))) / (max(out(:)) - min(o