【 MATLAB 】使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法

例题：

一个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)

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

第一种方法是博文里给出的：

【 MATLAB 】用 MATLAB 实现离散时间傅里叶变换（DTFT）的两个案例分析

第二个案例中，类比，如果我知道一个LTI系统的脉冲响应h(n)，那么我也能求出它的频率响应：

k = [0:M];
n = [n1:n2];
X = x * (exp(-j * pi/M)).^(n'*k);

博文里面由具体的推荐，上述程序中X就是频率响应，也就是x的DTFT，如果x换成h，则X可以换成H。

第二种方法是，通过差分方程直接求出系统的频率响应，求解的方法是通过向量化的方法。

设某一LTI系统的差分方程表示为：

$y(n)+ \sum_{l=1}^{N}a_ly(n-l)= \sum_{m = 0}^{M}b_mx(n-m)$

可以用一种简单的矩阵向量乘法来完成。如果在 $[0,\pi]$ 的 $k=0,1,...,K$ 个等分频率上求 $H(e^{jw})$ ,那么

注意，上述的b以及a向量都是行向量。

先给出第一种方法的脚本：

clc
clear
close allb = [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 allb = [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');

从图可以看出这确实是一个低通滤波器。

【 MATLAB 】使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法相关推荐

python求列表的中位数_python中获取中位数的两种方法
普通方法: 对列表进行排序,然后根据长度为奇数或者偶数的不同情况计算中位数 def huahua(x): length = len(x) print(length) x.sort() print(x) ...
python创建树结构、求深度_Python实现二叉树的最小深度的两种方法
找到给定二叉树的最小深度最小深度是从根节点到最近叶子节点的最短路径上的节点数量注意:叶子节点没有子树 Example: Given binary tree [3,9,20,null,null,15 ...
求1+2！+3！+...20！的两种方法
方法一单层循环的方式 <script>var sum = 0;var fac = 1;//原理:用sum 存储和 fac 用于计算阶乘for(var i = 1; i <= 20; ...
C++求整数的个（十）位数的两种方法
个位数一.将这个整数转换为字符串,截取字符串的最后一位即为整数的个位数 #include<iostream> #include <sstream> using namespa ...
求100到999中的水仙花数（两种方法，c语言）
#include<stdio.h> int main() {int a,b,c;for(a=1;a<=9;a++)for(b=0;b<=9;b++)for(c=0;c<= ...
matlab 求矩阵秩,求矩阵秩的两种方法及MATLAB的应用
摘要: 高等代数是一门逻辑思维比较强和理论知识比较深的学科, 它具有丰富的数学知识, 涉及许多重要的数学思想, 其在数学领域的应用很广泛, 如行列式.矩阵的相关计算和求解线性方程组的解方面的应 ...
用matlab计算稳态误差,利用Matlab求稳态误差的两种方法.
利用Matlab求稳态误差的两种方法摘要:稳态误差是系统控制精度或抗扰动能力的一种度量,它是稳态性能的一个重要指标.本文介绍利用Matlab的控制系统工具箱和Simulink工具箱求取系统误差稳态的 ...
matlab中给信号添加高斯白噪声的两种方法，awgn计算过程,randn函数
y=awgn(x,snr,px_dBW) 给信号x添加噪声功率为某个值的高斯白噪声. snr为信噪比,单位dB. px_dBW为信号x的指定功率(注意,是指定功率,而不是x本身的功率),单位dBW. ...
Matlab对含噪声图像的滤波操作_两种噪声_三种滤波器_两种方法
注释很重要 Matlab对含噪声图像的滤波操作. 噪声: 高斯噪声(正态分布) 均匀噪声用到的滤波器: 高斯滤波器盒型滤波器中值滤波器用到的两种方法: 直接conv2 fft2 %%C1 fi ...

【 MATLAB 】使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法

【 MATLAB 】使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法相关推荐

最新文章

热门文章