pfe.m

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

function [K,S,p] = pfe
% PFE.M
%  Partial fraction expansion of polynomial P/Q
% CALL: [K,S,p]= pfe
% INPUTS:  From Keyboard as array [a0 a1 ... an]
%	   P Numerator polyomial coefficients in decending order
%          Q Denominator polynomial coefficients in decending order
% OUTPUTS: K Column vector of coefficients in the expansion
%          S Column vector of corresponding roots of the denominator 
%            in the expansion
%          p Coefficients of polynomial if P/Q is not proper fraction	al
%
P=input('Coefficients of numerator in decending order = ')
Q=input('Coefficients of denominator in decending order = ')
[K,S,p]=residue(P,Q);	

% 
% Note: residue can also reconstruct P/Q if three arguments are
%        given. See command nargin.