找回密码
 注册
关于网站域名变更的通知
查看: 711|回复: 1
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-5-8 14:41 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
例题:
( t, o) E1 h( {* y3 Z; a& E
  M7 x- h/ i  l0 ?  @一个3阶低通滤波器由下面差分方程描述:
- d8 S4 e2 X8 o" W$ k9 c, v8 [. i) c; i4 b5 U" Q9 O6 ~
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)9 u# r3 O  U" B7 O% P9 v

, p8 _" D" r" i4 @" I画出这个滤波器的幅度和相位响应,并验证它是一个低通滤波器。8 P/ B9 R3 Q* {3 {0 C) F% z* ]
, P: y2 r. ~- {0 Q3 m" K7 d
第一种方法是博文里给出的:用 MATLAB 实现离散时间傅里叶变换(DTFT)的两个案例分析+ G9 F5 W2 T* y2 B% A7 {

/ `$ W  R% b; |第二个案例中,类比,如果我知道一个LTI系统的脉冲响应h(n),那么我也能求出它的频率响应:3 Z, Y9 e( i' j
- V2 V4 F. l( }; i$ Z% ]+ X+ \
  • k = [0:M];
  • n = [n1:n2];
  • X = x * (exp(-j * pi/M)).^(n'*k);

  • , N, Z. l0 ]+ U/ Z

2 K: W- I9 ?! ~
1 ~5 z6 F. n  h博文里面由具体的推荐,上述程序中X就是频率响应,也就是x的DTFT,如果x换成h,则X可以换成H。
+ c& [! S7 x' ]# `( \. C3 ?: ^* [/ ~# r4 H. J5 n6 l  ^" U
第二种方法是,通过差分方程直接求出系统的频率响应,求解的方法是通过向量化的方法。7 m0 x% y5 D+ N' S( F! D

/ a( n- ^7 v8 w设某一LTI系统的差分方程表示为:
0 |" Q, i" n+ y- G. r" m
, l; v, u4 B! [; ~0 ?
% Z2 b" W9 s& }6 V2 v, x
/ h0 a) q0 ]& U* }可以用一种简单的矩阵向量乘法来完成。如果在 个等分频率上求 ,那么/ ~, h9 `( |! @9 X$ z
) d$ h" o  l+ \! s
+ L. m/ o# \3 c$ e

* C1 y- U2 A* E* G8 w3 U注意,上述的b以及a向量都是行向量。8 v2 i4 H' C" ^9 Z6 o9 `) v, ]
7 Z9 k$ \) q6 m& O- D3 }4 U
先给出第一种方法的脚本:
4 {3 e! c  _- s! Q; w& D, }' D0 p8 T( ^' X0 Y) x, v& y  I: G2 [# H
  • clc
  • 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');
    9 y. l4 J2 B# }1 A" _  T9 Z/ J- {

# o+ y, e( V* A& L
' t; K4 w$ P$ W9 q 5 Z: A" ~* s& h' L6 H3 P
& I) ?  y# R% h" l

  X) G9 z8 V# z+ h1 w' Q这种方法的思路是通过差分方程可以得到有理传递函数或者频率响应的分子和分母系数,通过impz函数得到脉冲响应,之后由脉冲响应h(n)得到频率响应。
5 r( ^0 W: |/ @3 z$ D0 X* A* x
6 C* s: Z6 L8 A; h! d6 n; V" ^
  • 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');
    . q) |6 `, ]2 y$ F

9 U8 P- D, g6 U0 ]+ ~& o 0 ?- m& c' D/ Y3 Z

; g) u% D& v1 k: |8 S
% P' q+ h8 S2 Y: c( @2 R从图可以看出这确实是一个低通滤波器。/ ?* ^8 e! a- ^

/ a1 {+ d6 |  [, w/ }! Q, F8 |5 `* y& f  e( R0 j; c$ }/ g. H+ Q/ [7 [" g6 r

该用户从未签到

2#
发表于 2021-5-8 16:45 | 只看该作者
使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-10-29 22:15 , Processed in 0.187500 second(s), 26 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表