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