用MAtlab写一个双边滤波算法代码

可以这样写：

function output = bilateral_filter(input, sigma_d, sigma_r)
% bilateral_filter - Bilateral filtering of a grayscale image
%
% input - input image
% sigma_d - spatial standard deviation
% sigma_r - range standard deviation

% Get the size of the input image
[m, n] = size(input);

% Create