信号处理仿真：滤波器设计与仿真_（13）.MATLAB在滤波器设计中的应用.docxVIP

PAGE1

PAGE1

MATLAB在滤波器设计中的应用

1.滤波器设计的基本概念

在信号处理中，滤波器是一种用于提取或抑制信号中特定频率成分的系统。滤波器设计的目标是根据特定的应用需求，构建一个能够满足这些需求的滤波器。滤波器可以分为两大类：模拟滤波器和数字滤波器。模拟滤波器通常用于处理连续时间信号，而数字滤波器则用于处理离散时间信号。

1.1模拟滤波器

模拟滤波器包括低通滤波器、高通滤波器、带通滤波器和带阻滤波器。这些滤波器的设计通常基于电路理论，使用电阻、电容和电感等元件来实现。模拟滤波器的设计方法有多种，常见的有巴特沃斯滤波器、切比雪夫滤波器和椭圆滤波器等。

1.2数字滤波器

数字滤波器的设计则主要基于离散数学和计算机算法。数字滤波器可以分为无限脉冲响应（IIR）滤波器和有限脉冲响应（FIR）滤波器。IIR滤波器的特点是其脉冲响应是无限的，而FIR滤波器的脉冲响应是有限的。FIR滤波器通常具有线性相位特性，适用于对相位要求较高的应用。

2.MATLAB中的滤波器设计工具

MATLAB提供了丰富的工具和函数来设计和仿真滤波器，其中最常用的有FilterDesignandAnalysisTool（fdatool）、fir1、butter、cheby1、cheby2和ellip等函数。

2.1FilterDesignandAnalysisTool(fdatool)

fdatool是一个图形用户界面（GUI）工具，用于设计和分析滤波器。通过fdatool，用户可以方便地选择滤波器类型、设置滤波器参数、查看滤波器的频率响应和相位响应等。

2.1.1使用fdatool设计滤波器

打开fdatool

fdatool

选择滤波器类型在fdatool界面中，选择“DesignMethod”和“FilterType”。例如，选择“FIR”和“Lowpass”。

设置滤波器参数输入滤波器的频率规格，如通带边缘频率、阻带边缘频率、通带最大衰减和阻带最小衰减等。

设计滤波器点击“DesignFilter”按钮，fdatool会生成滤波器的系数，并显示其频率响应和相位响应。

导出滤波器设计完成后，可以将滤波器导出到MATLAB工作区，以便在后续的仿真中使用。

2.2使用MATLAB函数设计滤波器

MATLAB提供了多个函数用于设计滤波器，这些函数可以根据不同的滤波器类型和需求进行选择。

2.2.1设计FIR滤波器

FIR滤波器的设计可以使用fir1函数。fir1函数基于窗函数法，可以设计低通、高通、带通和带阻滤波器。

2.2.1.1低通FIR滤波器设计示例

%设计一个低通FIR滤波器

%采样频率为1000Hz，通带边缘频率为100Hz，窗函数选择汉宁窗

%采样频率

fs=1000;

%通带边缘频率

fc=100;

%滤波器阶数

N=50;

%汉宁窗

window=hann(N+1);

%设计低通滤波器

b=fir1(N,fc/(fs/2),window);

%计算滤波器的频率响应

[H,f]=freqz(b,1,512,fs);

%绘制频率响应

figure;

plot(f,20*log10(abs(H)));

xlabel(频率(Hz));

ylabel(幅度(dB));

title(低通FIR滤波器的频率响应);

gridon;

2.2.1.2高通FIR滤波器设计示例

%设计一个高通FIR滤波器

%采样频率为1000Hz，阻带边缘频率为150Hz，窗函数选择矩形窗

%采样频率

fs=1000;

%阻带边缘频率

fc=150;

%滤波器阶数

N=50;

%矩形窗

window=rectwin(N+1);

%设计高通滤波器

b=fir1(N,fc/(fs/2),high,window);

%计算滤波器的频率响应

[H,f]=freqz(b,1,512,fs);

%绘制频率响应

figure;

plot(f,20*log10(abs(H)));

xlabel(频率(Hz));

ylabel(幅度(dB));

title(高通FIR滤波器的频率响应);

gridon;

2.2.2设计IIR滤波器

IIR滤波器的设计可以使用butter、cheby1、cheby2和ellip函数。这些函数分别用于设计巴特沃斯滤波器、切比雪夫I型滤波器、切比雪夫II型滤波器和椭圆滤波器。

2.2.2.1巴特沃斯低通IIR滤波器设计示例

%设计一个