exa6_3.m

来自「MATLAB 的一些编程,基本的程序,为出学MATLAB的初学者」· M 代码 · 共 24 行

function [ht,gt,t]=exa6_3(b,a)
%求系统冲激响应和阶跃响应
if length(b)>length(a)
    error('此系统为物理不可实现系统!');
end
[r,p]=residue(b,a);
t=0:0.1:10;
ht=zeros(size(t));
for ii=1:length(r);
    ht=ht+r(ii)*exp(p(ii)*t);
end
figure;plot(t,ht);grid
hold on;

a=[a 0];
[r,p]=residue(b,a);
t=0:0.1:10;
gt=zeros(size(t));
for ii=1:length(r);
    gt=gt+r(ii)*exp(p(ii)*t);
end
plot(t,gt,'r-');legend('ht','gt');