基于MATLAB数字滤波器的设计

1. 线性平滑滤波器：使用MATLAB实现域平均法来抑制噪声。

2. 中值滤波：使用MATLAB实现中值滤波程序。

3.状态是指原李成昌统计滤波器：状态滤波程序是用ordfilt2函数实现的。

4、二维自适应去噪滤波器：利用Wiener2功能实现二维自适应去噪滤波器。

5. 过滤特定区域：使用 MATLAB 图像处理工具箱中提供的 roifilt2 函数过滤特定区域。

不好说，基本上是按照指标设计的。 让我举个例子。

如果汉明窗口设计满足索引：wp=; ws=;ap=1;AS=45的线性相位FIR高通滤波器。

FIR高通滤波器采用I型线性相位滤波器设计。

clear;

确定过滤器的顺序，使过滤器类型为I。

wp=;ws=;ap=1;as=45;

n=ceil(7*pi/(wp-ws));

n=mod(n+1,2)+n;

m=n-1;

汉明窗口长度 n 满足条件。

w=hamming(n)';

理想的低通截止值，计算截止值。

wc=(wp+ws)/2;

k=0:m;

hd=-(wc/pi)*sinc(wc*(;hd(;

h=hd.*w;

omega=linspace(0,pi,512);

mag=freqz(h,[1],omega);

plot(omega/pi,20*log10(abs(mag)))grid;% 绘制增益响应。

omega1=linspace(0,wp,512);

h1=freqz(h,[1],omega1);

omega2=linspace(0,ws,512);

h2=freqz(h,[1],omega2);

fprintf('ap=%.4f',-20*log10(max(abs(h1)))

fprintf('as=%.4f'， -20 x log10（max（abs（h2））） 计算 AS 和 AP

举个例子，我以前做过：

令人满意的WP=; ws=;ap=;as=45;低通线性相位FIR滤波器：

clear;

wp=;ws=;ap=;as=45;

n=ceil(;

n=mod(n+1,2)+n;

m=n-1;

w=blackman(n)';

wc=(wp+ws)/2;

k=0:m;

hd=(wc/pi)*sinc(wc*(;

h=hd.*w;

omega=linspace(0,pi,512);

mag=freqz(h,[1],omega);

plot(omega/pi,20*log10(abs(mag)))

grid;希望对您有所帮助 请随时询问