[SciPy-User] [OT] Transform (i.e., Fourier, Laplace, etc.) methods in Prob. & Stats.

S. Chris Colbert sccolbert at gmail.com
Thu Nov 26 17:45:49 EST 2009


i dont know if it will be of any use to you, but:

Laplace transforms and their inversions are also used extensively in control 
theory. 

I wrote some python for the numerical inversion awhile back (before i was any 
good in python, be warned!) 

There are two different inversion methods in the attached file: the method of 
riemann sums is the faster of the two, and here is the reference from which I 
made my implementation:

http://books.google.com/books?id=CmX1aHur7jcC&pg=PA410&lpg=PA410&dq=%27%27%27This+algorithm+is+proposed+by+Tzou,
+Ozisik+and+Chiffelle+%281994%29%27%27%27&source=bl&ots=NSiw3tKRvG&sig=cqfa_ka_baPbcnhoSA9Gcxo8Vj8&hl=en&ei=YgQPS-3dKJ2qmwPO6LncBQ&sa=X&oi=book_result&ct=result&resnum=1&ved=0CAgQ6AEwAA#v=onepage&q=%27%27%27This%20algorithm%20is%20proposed%20by%20Tzou%2C%20Ozisik%20and%20Chiffelle%20%281994%29%27%27%27&f=false

You can probably throw out the stehfest method. I was using it originally as 
it was faster, but then I managed to vectorize the riemann method. I had no 
luck vectorizing the stehfest method. 

Cheers, 

Chris

> Sounds good, thanks again!
> 
> DG
> 
> On Thu, Nov 26, 2009 at 12:34 PM, nicky van foreest 
<vanforeest at gmail.com>wrote:
> > > Chung, K. L., 2000. "A Course In Probability Theory, 2nd Ed." Academic.
> > > looks like a really good general reference, Nicky - I assume this is
> > > the "Chung" to which you were referring?  Thanks!!!
> >
> > That is the one indeed.
> >
> > For z transforms you might like generatingfunctionoly (or something
> > like this). Search on Herbert Wilf, On his website you can find a very
> > nice book on the uses of z transforms. The first chapter is very
> > accessible.
> >
> > HOpe this helps.
> >
> > NIcky
> > _______________________________________________
> > SciPy-User mailing list
> > SciPy-User at scipy.org
> > http://mail.scipy.org/mailman/listinfo/scipy-user
> 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inverselaplace.py
Type: text/x-python
Size: 4864 bytes
Desc: not available
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20091126/eece7333/attachment.py>


More information about the SciPy-User mailing list