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

 

先给出filter的大致介绍：【 MATLAB 】filter 函数介绍（一维数字滤波器）

上篇博文写了产生一个系统的脉冲响应的博文，有兴趣可以看看：【 MATLAB 】使用 impz 函数计算并画出脉冲响应

使用到的filter语法如下：

y = filter(b,a,x)

对具有分子系数向量b和分母系数向量a的系统函数对x进行滤波，得到响应y。

如果x为一个单位阶跃序列，则响应就是单位阶跃响应。

这就是理论上的基础思路。

给出实例分析：

这个系统使用差分方程来描述如下：

y(n) - y(n-1) + 0.9y(n-2) = x(n)

要求计算并画出它在 你= -20，..., 100的单位阶跃响应。

写下面的MATLAB脚本之前必须声明的是下面脚本中使用了之前博文中写的一个阶跃序列的函数：

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];

更多相关信号的函数见博文：【 MATLAB 】基本序列运算及其MATLAB的等效表示

下面使用MATLAB脚本实现：

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)');

 

最后想说的是这篇博文的行文思路和上篇博文一致，但是简化了许多，如果不明白请务必看下上篇博文：

 

【 MATLAB 】使用 impz 函数计算并画出脉冲响应