[SciPy-User] orthonormal polynomials (cont.)
nicky van foreest
vanforeest at gmail.com
Tue Jun 14 14:10:53 EDT 2011
Hi,
Without understanding the details... I recall from numerical recipes
in C that Gram Schmidt is a very risky recipe. I don't know whether
this advice also pertains to fitting polynomials, however,
Nicky
On 14 June 2011 18:58, <josef.pktd at gmail.com> wrote:
> (I'm continuing the story with orthogonal polynomial density
> estimation, and found a nice new paper
> http://www.informaworld.com/smpp/content~db=all~content=a933669464 )
>
> Last time I managed to get orthonormal polynomials out of scipy with
> weight 1, and it worked well for density estimation.
>
> Now, I would like to construct my own orthonormal polynomials for
> arbitrary weights. (The weights represent a base density around which
> we make the polynomial expansion).
>
> The reference refers to Gram-Schmidt or Emmerson recurrence.
>
> Is there a reasonably easy way to get the polynomial coefficients for
> this with numscipython?
>
> Thanks,
>
> Josef
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
More information about the SciPy-User
mailing list