Matlab设计数字滤波器入门

最新推荐文章于 2024-10-28 00:04:27 发布

山无忧

最新推荐文章于 2024-10-28 00:04:27 发布

阅读量766

点赞数 1

分类专栏：数字信号处理文章标签： matlab 数字信号处理

本文链接：https://blog.csdn.net/weixin_44316365/article/details/104525436

版权

数字信号处理专栏收录该内容

5 篇文章 0 订阅

订阅专栏

该博客介绍了如何使用MATLAB代码分析3阶低通滤波器的幅度和相位响应，以验证其低通特性。通过差分方程获取滤波器的系数，然后利用impz函数和直接计算方法分别绘制了幅度和相位响应曲线，展示了滤波器在不同频率下的行为。

摘要由CSDN通过智能技术生成

一个3阶低通滤波器由下面差分方程描述：

y(n) = 0.0181 x(n) + 0.0543 x(n-1) + 0.0543 x(n-2) + 0.0181 x(n-3) + 1.76 y(n-1) - 1.1829 y(n-2) + 0.2781 y(n-3)

画出这个滤波器的幅度和相位响应，并验证它是一个低通滤波器。
代码1：

clear
close all
 
b = [0.0181,0.0543,0.0543,0.0181];
a = [1.0000,-1.7600,1.1829,-0.2781];
[h,t]=impz(b,a);
k = [0:500];
w = (pi/500)*k;
t = t';
h = h';
H = h * ( exp(-j*pi/500) ).^(t'*k);
magH = abs(H);
angH = angle(H);
subplot(2,1,1);
plot(w/pi,magH);
title('Magnitude part');
subplot(2,1,2);
plot(w/pi,angH);
title('Angle part');

在这里插入图片描述

这种方法的思路是通过差分方程可以得到有理传递函数或者频率响应的分子和分母系数，通过impz函数得到脉冲响应，之后由脉冲响应h(n)得到频率响应
代码二：

clc
clear
close all
 
b = [0.0181,0.0543,0.0543,0.0181];
a = [1.0000,-1.7600,1.1829,-0.2781];
m = 0:length(b)-1;
l = 0:length(a)-1;
k = 0:500;
w = (pi/500)*k;
nume = b * exp(-j * m' * w);
den = a * exp(-j * l' * w);
H = nume ./ den;
magH = abs(H);
angH = angle(H);
 
subplot(2,1,1);
plot(w/pi,magH);
title('Magnitude Response');
 
subplot(2,1,2);
plot(w/pi,angH);
title('Phase Response');