invprony.m

这是一个 York大学开发的FDTD代码。很好用的工具

function fp=invprony(C,Sp,Ts,Nts,Nte)
% fp=invprony(C,Sp,Ts,Nts,Nte)
% Determine the representation of a time response from the poles discovered 
% the method of prony which respresents the time series
% as a set of exponential basis functions such that:
%    f(N*Ts)=C0*exp(Sp0*N*Ts) + C1*exp(Sp1*N*Ts) + ... + CNp-1*exp(SpNp-1*N*Ts) 
%  J. F. Dawson 27 Oct 1999

Np=size(C,1)

mu=[exp((Nts-1)*Sp.*Ts).'];
for cnt=Nts:Nte-1
    mu=[mu;exp(cnt*Sp.*Ts).'];
end
fp=mu*C;