polyadd.m

来自「Mathematical Methods by Moor n Stiling.」· M 代码 · 共 19 行

function p = polyadd(a,b)
%
% Add the polynomials p=a+b
% 
% function p = polyadd(a,b)
%
% a,b = polynomial
%
% p = polynomial sum.

% Copyright 1999 by Todd K. Moon

k = length(b) - length(a);
p = [zeros(1,k) a] + [zeros(1,-k) b];
if(all(p==0))
  p = 0;
else
  p = p(find((p==0)==0):end);   % get rid of leading zeros
end