invlap.m

来自「Programs for the book Advanced Engineeri」· M 代码 · 共 12 行

% INVLAP.M Factors of Inverse Laplace transform of 
%  (s^2+4s+2)/s(s+1)(s+2) using conv and residue
%    
% Compute the coefficients of denominator polynomial
Q1=[1 0];            % s
Q2=[1 1];            % s+1
QD1=conv(Q1,Q2);
Q3=[1 2];            % s+2
Q=conv(QD1,Q3)       % Coefficients of denominator
% Compute partial fraction coefficients
P=[1 4 2]            % Coefficients of numerator
[K,S,p]=residue(P,Q) % K coefficients, S roots of Q