rilt

Regularized Inverse Laplace Transform

Description

rilt Regularized Inverse Laplace Transform
  
[g,yfit,cfg] = rilt(t,y,s,g0,alpha)
  
Array g(s) is the Inverse Laplace Transform of the array y(t), calculated by a regularized least squares method.

This script is an emulation of S. Provencher CONTIN program, written in Fortran. See http://s-provencher.com/pages/contin.shtml.

• MATLAB release: MATLAB 7.3 (R2006b)

[c rod]

 

[Heiner Hegewald]

I have to test it more deeply.

[lu yijie]

its good

[Sean Gu]

It has been 20 minutes and matlab is still running with this program. Seems it needs to optimize the search to improve the time for regulization

[John]

I spent several days trying to use Provencher's original CONTIN program. It is many things, but user-friendly is not one of them, especially since I haven't written in Fortran since about 1987. In contrast, I got this program running in seconds, and convinced myself that it was working properly after playing with it for an hour or so. It will become a standard tool for my research group. Many thanks!

[L H]

The program is great and lets the user to avoid the trouble of dealing with the original Fortran CONTIN program. The program works great, although I have just a couple comments:
- The "s" array fed to rilt is actually a tau array (s = 1/tau).
- "rilt.m" should be called "rlm.m" because it is actually the forward transform by the majority of definitions I've ever seen. L[f(t)] = F(s) = F(1/tau) is the forward transform.

That said, I want to say thank-you to the programmer. This routine has helped me out with an aspect of my thesis.

[Xinjun Xinjun]

It is very easy to use,compared with CONTIN.

[trabelsi saber]

Hello, I want to use the program rilt to determine the inverse Laplace tronsform of the autocorrelation function g1(t)=exp(-t/tau). but I did not understand the input variables g0 and S.
sincerly yours.

[Xinjun Xinjun]

Usually the original CONTIN use the second-order difference operator,which gives the smoothest solution. If the rilt uses the second order difference operator,it will be more closer to the original CONTIN.