Appendix B, matlab code for den Iseger numerical inverse Laplace transform. In the function DenIseger(Lf, t), Lf is a 'model', e.g. Lf = 1/(s^2+1). We obtain sin(t). But what if Lf is a vector? I failed to find a matlab function like 'fit' that would work with complex arguments. 'Fit' could give me a model to evaluate in ft = real(Lf(s)).
Article Numerical Inversion of Laplace Transform for Time Resolved T...
The algorithm developped by Peter Den Iseger needs a model.
You may try this:
Fo = fit(s, real(Lf), FT); % FT is the fittype
L = @(s) Fo(s);
[ft] = Heavy_DenIseger(L, t);
Maybe there's a probleme for using complex numbers with the fit object. You can obtain the coefficients of your fitting law with the function coeffvalues
Yes, I know about this problem. That is the reason why I asked my question :)
- Badges
- Science topic
- Similar topics
- Mathematics
- Statistics