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

使用 filter 函数产生给定线性时不变系统的单位阶跃响应

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
$ D3 ?) o' F, E/ J! C1 d
先给出filter的大致介绍:MATLAB之filter 函数介绍(一维数字滤波器)
# g- Z, W! E3 c# e6 k  N3 M8 Q; C1 s' S+ k* p  ^1 H6 z: W
上篇博文写了产生一个系统的脉冲响应的博文,有兴趣可以看看:MATLAB之使用 impz 函数计算并画出脉冲响应) K( X4 \7 @: p1 o

& g8 Q9 h* ^7 P4 Q1 Y' q使用到的filter语法如下:! k1 `$ s/ M0 {8 E; d& P

! u- ~9 c: p6 `/ Jy = filter(b,a,x)
# d& j" _, C! N8 q  ^7 ?
% C1 ~, X! J+ ~/ i- A6 P+ e- f) B对具有分子系数向量b和分母系数向量a的系统函数对x进行滤波,得到响应y。
) B; S  g. U5 z) A+ |* j& R
9 y  R* R5 M% p' r" T- F& [  O2 J如果x为一个单位阶跃序列,则响应就是单位阶跃响应。3 v- y: ^# }/ w- \- ^3 q' @0 `, L+ P
8 Q: M( M) g" z
这就是理论上的基础思路。! L+ T: Q/ B6 k4 t5 q" s

$ J5 q5 m  Q* [$ h4 m7 T- d给出实例分析:
# A6 v- O  y6 m( O
& |: t$ H; b) `: C7 u, b- m+ ]这个系统使用差分方程来描述如下:
9 U' ], D; s8 i8 G
4 f1 h; x5 o" |. Y9 j1 Fy(n) - y(n-1) + 0.9y(n-2) = x(n)
% r; G  v$ M! n3 {& G, l& \" J. Y! |' |4 z2 ]
要求计算并画出它在 你= -20,..., 100的单位阶跃响应。
7 G2 Z$ |! E. g# S! D# R' K( Y# I+ R: R5 h% A( a. T5 T6 m/ r
写下面的MATLAB脚本之前必须声明的是下面脚本中使用了之前博文中写的一个阶跃序列的函数:
! k7 M# Z: _+ p8 t+ W+ O$ }2 K
8 ~2 d: h6 k0 p' C% Q+ F
  • function [x,n]=stepseq(n0,n1,n2);
  • % generate x(n) = u(n - n0); n1 <= n <= n2
  • %——————————————————————
  • %[x,n] = stepseq(n0, n1, n2);
  • %
  • n = [n1:n2];
  • x = [(n-n0) >= 0];
    " b1 V! c- e! X0 t4 `$ R

, J+ b, {) ?6 A- D* R4 \  _: u. H+ F; U8 M
更多相关信号的函数见博文:基本序列运算及其MATLAB的等效表示
* A2 O" C3 i2 v4 m) o& X: \
& v. g7 c) ^  h! ~下面使用MATLAB脚本实现:) ]# H1 x0 q. F# o

8 u+ Y$ w" o" y/ N# S: S3 Z* J
  • clc
  • clear
  • close all
  • b = 1;
  • a = [1,-1,0.9];
  • % n = [-20:100];
  • [x,n] = stepseq(0,-20,100);
  • stem(n,x);
  • title('step sequence');
  • xlabel('n');ylabel('x(n)');
  • ylim([0,1.2]);
  • figure;
  • s = filter(b,a,x);
  • stem(n,s);
  • title('step response');
  • xlabel('n');ylabel('s(n)');
    , m# A" s: \0 i: Y$ S1 L9 {6 l# G" i
  % q. }) ~* q, g# {" j( T% h& p

5 h' E2 T3 h6 ^ ' S& t7 r. l! y3 T( G* d2 N
3 Y7 e) G+ ]  B) k

0 W2 G* [  g% H4 t( ?0 o" G2 a3 d* W/ ~& N8 M  ^2 m
1 b9 g7 M; Z% k8 x4 m7 b
  • TA的每日心情
    奋斗
    2020-9-2 15:06
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2021-3-16 19:56 | 只看该作者
    s = filter(b,a,x)
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-10-6 03:26 , Processed in 0.125000 second(s), 26 queries , Gzip On.

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

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

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