在MATLAB中实现滤波程序，可以使用内置的滤波函数，如`filter`、`filtfilt`、`conv`等，也可以使用自定义的滤波器设计函数，如`butter`、`firpm`等。以下是几种常见滤波器的MATLAB实现方法：
1. 线性平滑滤波器
```matlab
% 读取图像
I = imread（'c4.jpg'）；
% 转换为灰度图像
flag = isrgb（I）；
if flag == true
I = rgb2gray（I）；
else
I = I；
end
% 添加椒盐噪声
I1 = imnoise（I, 'salt & pepper', 0.02）；
% 3x3模板平滑滤波
k1 = filter2（fspecial（'average', 3）, I1） / 255；
% 5x5模板平滑滤波
k2 = filter2（fspecial（'average', 5）, I1） / 255；
% 7x7模板平滑滤波
k3 = filter2（fspecial（'average', 7）, I1） / 255；
% 9x9模板平滑滤波
k4 = filter2（fspecial（'average', 9）, I1） / 255；
% 显示滤波结果
subplot（233）； imshow（k1）；
subplot（234）； imshow（k2）； title（'5x5 模板平滑滤波'）；
```
2. 中值滤波器
```matlab
% 读取图像
I = imread（'c4.jpg'）；
% 转换为灰度图像
flag = isrgb（I）；
if flag == true
I = rgb2gray（I）；
else
I = I；
end
% 中值滤波
J = medfilt2（I, [3 3]）；
% 显示原图和处理后的图像
figure；
subplot（1, 2, 1）； imshow（I）； title（'原图'）；
subplot（1, 2, 2）； imshow（J）； title（'处理后'）；
```
3. 巴特沃斯低通滤波器
```matlab
% 设置技术指标
Fs = 1000； % 采样频率
fp = 200； % 通带截止频率
Rp = 3； % 通带衰减
fs = 500； % 阻带截止频率
Rs = 30； % 阻带衰减
% 设计滤波器
d = fdesign.lowpass（'N,F3dB', 5, fp, Fs）；
Hd = design（d, 'butter'）；
% 应用滤波器
x = cos（2 * pi * 60 * linspace（0, 1, Fs）） + sin（2 * pi * 120 * linspace（0, 1, Fs）） + randn（size（linspace（0, 1, Fs）））；
y = filter（Hd, x）；
```
4. 使用`filter`函数进行滤波
```matlab
% 定义滤波器系数
a = [1.0000； 0； 0.4860； 0； 0.0177]；
b = [0.0940； 0.3759； 0.5639； 0.3759； 0.0940]；
% 原始信号
OriginalSignal = [1； 0.5； 0.8； 0.5； 1]；
% 应用滤波器
SignalFilter = filter（b, a, OriginalSignal）；
```
5. 使用`filtfilt`函数进行双向滤波
```matlab
% 定义滤波器系数
b = [0.0940； 0.3759； 0.5639； 0.3759； 0.0940]；
a = [1.0000； 0； 0.4860； 0； 0.0177]；
% 原始信号
x = [1； 0.5； 0.8； 0.5； 1]；
% 应用双向滤波
y = filtfilt（b, a, x）；
```
6. 使用`conv`函数进行卷积滤波