上篇博文写了产生一个系统的脉冲响应的博文,有兴趣可以看看:【 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脚本实现:
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)');
最后想说的是这篇博文的行文思路和上篇博文一致,但是简化了许多,如果不明白请务必看下上篇博文:
本文同步分享在 博客“李锐博恩”(CSDN)。
如有侵权,请联系 support@oschina.cn 删除。
本文参与“OSC源创计划”,欢迎正在阅读的你也加入,一起分享。