基于MATLAB的prony方法，可用来做数据处理

function [ba] = prony（h nb na）
K = length（h） - 1;
M = nb; N = na;
if K <= max（MN）      % zero-pad input if necessary
    K = max（MN）+1;
    h（K+1） = 0;
end
c = h（1）;
if c==0    % avoid divide by zero
    c=1;
end
H = toeplitz（h