Orthogonal polynomials/functions on the interval [0,1] but with same weight as Gegenbauer polynomials - MathOverflow most recent 30 from http://mathoverflow.net 2013-05-24T02:54:39Z http://mathoverflow.net/feeds/question/81553 http://www.creativecommons.org/licenses/by-nc/2.5/rdf http://mathoverflow.net/questions/81553/orthogonal-polynomials-functions-on-the-interval-0-1-but-with-same-weight-as-ge Orthogonal polynomials/functions on the interval [0,1] but with same weight as Gegenbauer polynomials Marcus P S 2011-11-21T21:50:48Z 2013-01-01T19:30:29Z <p>I am looking for an othogonal basis of functions over the interval $[0,1]$ with weight function $(1-x^2)^{\alpha-1/2}$. Gegenbauer polynomials are frustratingly close to what I need, but they are defined over the interval $[-1,1]$, and a change of variables ends up changing the weight function.</p> <p>None of the orthogonal polynomial families I have looked at (Chebyshev, Gegenbauer, Legendre, Laguerre, Jacobi, Hermite) have this property. </p> <p>Does anyone know of a family that does? Suggestions for references that may point the way would also be very helpful.</p> <p>Thanks!</p> http://mathoverflow.net/questions/81553/orthogonal-polynomials-functions-on-the-interval-0-1-but-with-same-weight-as-ge/81555#81555 Answer by robot for Orthogonal polynomials/functions on the interval [0,1] but with same weight as Gegenbauer polynomials robot 2011-11-21T23:09:57Z 2011-11-21T23:09:57Z <p>If your polynomial is related to symmetries in one way or another it may be a <a href="http://en.wikipedia.org/wiki/Matrix_coefficient" rel="nofollow">matrix coefficient</a> of a Lie group representation.</p> http://mathoverflow.net/questions/81553/orthogonal-polynomials-functions-on-the-interval-0-1-but-with-same-weight-as-ge/81574#81574 Answer by J. M. for Orthogonal polynomials/functions on the interval [0,1] but with same weight as Gegenbauer polynomials J. M. 2011-11-22T04:32:23Z 2011-11-22T13:12:06Z <p>There is <a href="http://www.cs.kuleuven.be/publicaties/rapporten/tw/TW534.pdf" rel="nofollow">this paper</a> and <a href="http://dx.doi.org/10.1016/j.cam.2011.10.006" rel="nofollow">this paper</a> which treat the special case of "half-range Chebyshev polynomials" (both kinds, corresponding to the weights $\dfrac1{\sqrt{1-x^2}}$ and $\sqrt{1-x^2}$ over $[0,1]$) to deal with Fourier expansions of nonperiodic functions. I have a feeling that half-range Gegenbauer polynomials have been treated before, and I'll try to see what I can dig up.</p> <p>In the meantime, one can use the Stieltjes procedure to build up the recursion relations for these half range Gegenbauers. Letting</p> <p>$$\langle f(x),g(x) \rangle^{(\alpha)}=\int_0^1 (1-t^2)^{\alpha-1/2} f(t)g(t)\mathrm dt$$</p> <p>be the associated inner product, the Stieltjes procedure for generating monic orthogonal polynomials $\phi_k(x)$ uses the formulae</p> <p><code>\begin{align*}b_k&amp;=\frac{\langle x\phi_k(x),\phi_k(x)\rangle^{(\alpha)}}{\langle\phi_k(x),\phi_k(x)\rangle^{(\alpha)}}\\ c_k&amp;=\frac{\langle\phi_k(x),\phi_k(x)\rangle^{(\alpha)}}{\langle\phi_{k-1}(x),\phi_{k-1}(x)\rangle^{(\alpha)}}\end{align*}</code></p> <p>to give the coefficients $b_k,c_k$ for the recursion relation</p> <p>$$\phi_{k+1}(x)=(x-b_k)\phi_k(x)-c_k\phi_{k-1}(x)$$</p> <p>Here, the result</p> <p>$$\int_0^1 (1-t^2)^{\alpha-1/2}t^k \mathrm dt=\frac{\Gamma\left(\frac{1+k}{2}\right)\Gamma\left(\alpha+\frac12\right)}{2\Gamma\left(\alpha+\frac{k}{2}+1\right)}$$</p> <p>is useful.</p> <hr> <p>I might as well throw this in. There is an algorithm due to Chebyshev (1859) for determining recursion coefficients from the moments. I've already talked about the algorithm <a href="http://math.stackexchange.com/questions/13174/34047#34047" rel="nofollow">here</a>, so I shall not repeat myself. Instead, I'll reproduce the <em>Mathematica</em> routine I gave in that answer:</p> <pre><code>chebAlgo[mom_?VectorQ, prec_: MachinePrecision] := Module[{n = Quotient[Length[mom], 2], si = mom, ak, bk, np, sp, s, v}, np = Precision[mom]; If[np === Infinity, np = prec]; ak[1] = mom[[2]]/First[mom]; bk[1] = First[mom]; sp = PadRight[{First[mom]}, 2 n - 1]; Do[ sp[[k - 1]] = si[[k - 1]]; Do[ v = sp[[j]]; sp[[j]] = s = si[[j]]; si[[j]] = si[[j + 1]] - ak[k - 1] s - bk[k - 1] v; , {j, k, 2 n - k + 1}]; ak[k] = si[[k + 1]]/si[[k]] - sp[[k]]/sp[[k - 1]]; bk[k] = si[[k]]/sp[[k - 1]]; , {k, 2, n}]; N[{Table[ak[k], {k, n}], Table[bk[k], {k, n}]}, np] ] </code></pre> <p>Here for instance is how to use <code>chebAlgo[]</code> to generate recursion coefficients for the monic half-range Chebyshev polynomials of the first kind:</p> <pre><code>With[{a = 0}, chebAlgo[Table[Gamma[(k + 1)/2] Gamma[a + 1/2]/Gamma[a + k/2 + 1], {k, 0, 10}]/2, Infinity]] // FullSimplify </code></pre> http://mathoverflow.net/questions/81553/orthogonal-polynomials-functions-on-the-interval-0-1-but-with-same-weight-as-ge/81661#81661 Answer by Ben Adcock for Orthogonal polynomials/functions on the interval [0,1] but with same weight as Gegenbauer polynomials Ben Adcock 2011-11-22T22:28:10Z 2011-11-22T22:28:10Z <p>Actually, quite a lot is known about such polynomials, at least in the asymptotic regime $n \rightarrow \infty$, where $n$ is the index of the $n$th orthogonal polynomial. In the <a href="http://www.sciencedirect.com/science/article/pii/S000187080300344X" rel="nofollow">paper</a> there are asymptotic formulae for not only the polynomials themselves, but also the coefficients of the recurrence relation. On could, in theory, use such formulae to compute recurrence coefficients (for large n), combined with a standard algorithm (such as the one posted by J. M.) for coefficients corresponding to small n.</p> http://mathoverflow.net/questions/81553/orthogonal-polynomials-functions-on-the-interval-0-1-but-with-same-weight-as-ge/117801#117801 Answer by JFvD for Orthogonal polynomials/functions on the interval [0,1] but with same weight as Gegenbauer polynomials JFvD 2013-01-01T19:30:29Z 2013-01-01T19:30:29Z <p>The even (degree) Gegenbauer polynomials $C_{2n}^{(\alpha )}$, $n=0,1,2,\ldots$ form an orthogonal basis for the space of square integrable functions over [0,1] with respect to the weight function $(1-x^2)^{\alpha -1/2}$, and so do the odd (degree) Gegenbauer polynomials $C_{2n+1}^{(\alpha )}$, $n=0,1,2,\ldots$</p>