基于MATLAB的布莱克曼窗FIR数字低通滤波器设计程序.pdf

基于MATLAB的布莱克曼窗FIR数字低通滤波器设计程序 本文主要介绍了基于MATLAB的布莱克曼窗FIR数字低通滤波器设计程序。该程序使用MATLAB语言编写，实现了基于布莱克曼窗函数设计的FIR数字低通滤波器。 一、程序概述 该程序的主要目的是设计一个基于布莱克曼窗函数的FIR数字低通滤波器。该滤波器的设计主要通过以下步骤实现： 1. 定义了滤波器的参数，包括通带频率Wp、阻带频率Ws、过渡带宽度tr_width和滤波器长度N。 2. 然后，使用ideal_hp1函数计算理想低通率的单位冲击响应hd。 3. 接下来，使用blackman函数计算布莱克曼窗函数w_ham。 4. 然后，计算实际的单位冲击响应h，通过将理想脉冲响应hd与窗函数w_ham相乘。 5. 使用freqz_m2函数计算实际滤波器的幅度响应和相位响应。 二、理想低通率的单位冲击响应 ideal_hp1函数用于计算理想低通率的单位冲击响应hd。该函数的实现方式是使用sin函数计算单位冲击响应。 hd = [sin(pi*m) - sin(Wc*m)]./(pi*m) 其中，m是时域采样点，Wc是理想低通滤波器的截止频率。 三、布莱克曼窗函数 blackman函数用于计算布莱克曼窗函数w_ham。该函数的实现方式是使用黑曼窗函数公式： w_ham = blackman(N) 其中，N是滤波器长度。 四、实际的单位冲击响应 实际的单位冲击响应h是通过将理想脉冲响应hd与窗函数w_ham相乘计算得到。 h = hd.*w_ham 五、实际滤波器的幅度响应和相位响应 freqz_m2函数用于计算实际滤波器的幅度响应和相位响应。该函数的实现方式是使用MATLAB的freqz函数计算滤波器的幅度响应和相位响应。 db = 20*log10((mag+eps)/max(mag)) pha = angle(H) 其中，mag是滤波器的幅度响应，H是滤波器的频率响应。 六、滤波器的时域和频域性能 程序还计算了滤波器的时域和频域性能，包括滤波器的时域响应、频域响应和幅度响应。 本文介绍了基于MATLAB的布莱克曼窗FIR数字低通滤波器设计程序的设计和实现过程。该程序可以用于设计和实现FIR数字低通滤波器，满足实际应用中的需求。