zoukankan      html  css  js  c++  java
  • Clausen Functions (and related series, functions, integrals)

    Since the Clausen functions are intimately related to a number of other important special functions, such as Inverse Tangent Integrals, Polylogarithms, Polygamma Functions, Zeta Functions, and more besides - many of which are at the forefront of modern mathematical research.
    Define the Clausen Functions (displaystyle ext{Cl}_m( heta)) and (displaystyle ext{Sl}_m( heta)) as follows:

    [ ext{Cl}_m( heta) = egin{cases} displaystylesum_{k=1}^{infty}frac{sin k heta}{k^m} & ext{if }m ext{ is even} \ displaystylesum_{k=1}^{infty}frac{cos k heta}{k^m} & ext{if }m ext{ is odd} end{cases}~,~ ext{Sl}_m( heta) = egin{cases} displaystylesum_{k=1}^{infty}frac{cos k heta}{k^m} & ext{if }m ext{ is even} \ displaystylesum_{k=1}^{infty}frac{sin k heta}{k^m} & ext{if }m ext{ is odd} end{cases} ]

    As one would therefore expect, basic trigonometry can be used to develop all sorts of Clausen function properties. For example, consider the difference

    [egin{align*} & ext{Cl}_{2m}( heta)- ext{Cl}_{2m}(pi- heta)=sum_{k=1}^{infty}frac{sin k heta}{k^{2m}}-sum_{k=1}^{infty}frac{sin k(pi- heta)}{k^{2m}}\ &=sum_{k=1}^{infty}frac{sin k heta}{k^{2m}}-sum_{k=1}^{infty}frac{(sinpi kcos k heta-cospi ksin k heta)}{k^{2m}}\ &=sum_{k=1}^{infty}frac{sin k heta}{k^{2m}}+sum_{k=1}^{infty}frac{(cospi ksin k heta)}{k^{2m}}=sum_{k=1}^{infty}frac{sin k heta}{k^{2m}}+sum_{k=1}^{infty}(-1)^kfrac{sin k heta}{k^{2m}}\ &=2\,left(frac{sin 2 heta}{2^{2m}}+frac{sin 4 heta}{4^{2m}}+frac{sin 6 heta}{6^{2m}}+cdots\, ight)\ &=displaystyle frac{2}{2^{2m}}\,left(sin 2 heta+frac{sin 4 heta}{2^{2m}}+frac{sin 6 heta}{3^{2m}}+cdots\, ight)=frac{1}{2^{2m-1}} ext{Cl}_{2m}(2 heta) end{align*}]

    We now have the duplication formula for the CL-type Clausen function of even order:
    (Largemathbf{color{Purple}{Result ~1:}})

    [ ext{Cl}_{2m}(2 heta)=2^{2m-1}Big[ ext{Cl}_{2m}( heta)- ext{Cl}_{2m}(pi- heta)Big] ]


    Next, let's take the above duplication formula, replace ( heta) with the variable (x), and integrate both sides:

    [int_0^{varphi} ext{Cl}_{2m}(2x)\,mathrm{d}x=2^{2m-1}left[int_0^{varphi} ext{Cl}_{2m}(x)\,mathrm{d}x-int_0^{varphi} ext{Cl}_{2m}(pi-x)\,mathrm{d}x ight] ]

    The L.H.S. is equal to

    [egin{align*} &sum_{k=1}^{infty}frac{1}{k^{2m}}\,int_0^{varphi}sin 2kx\,mathrm{d}x=-frac{1}{2}\,sum_{k=1}^{infty}frac{1}{k^{2m+1}}Big[cos 2kxBig]_0^{varphi}\ &=-frac{1}{2}\,sum_{k=1}^{infty}frac{cos 2kvarphi}{k^{2m+1}}+frac{1}{2}\,sum_{k=1}^{infty}frac{1}{k^{2m+1}}=frac{1}{2}\,Big[zeta(2m+1)- ext{Cl}_{2m+1}(2varphi)Big] end{align*}]

    Whereas the difference of the two integrals on the R.H.S. is

    [egin{align*} &2^{2m-1}\,sum_{k=1}^{infty}frac{1}{k^{2m}}int_0^{varphi}left[sin kx-sin k(pi-x) ight]\,mathrm{d}x\ =&2^{2m-1}\,sum_{k=1}^{infty}frac{1}{k^{2m+1}}Big[-cos kx+cos k(pi-x)Big]_0^{varphi}\ =&2^{2m-1}\,sum_{k=1}^{infty}frac{[cos k(pi-varphi)-cos kvarphi]}{k^{2m+1}}=2^{2m-1}[ ext{Cl}_{2m+1}(pi-varphi)- ext{Cl}_{2m+1}(varphi)] end{align*}]

    We now have the duplication formula for a CL-type Clausen function of odd order:
    (Largemathbf{color{Purple}{Result ~2:}})

    [ ext{Cl}_{2m+1}(2 heta)=zeta(2m+1)+2^{2m}Big[ ext{Cl}_{2m+1}( heta)- ext{Cl}_{2m+1}(pi- heta)Big] ]


    Let's return to the first of the two duplication formulae, and derive a few particular values and correlations:

    [ ext{Cl}_{2m}(2 heta)=2^{2m-1}Big[ ext{Cl}_{2m}( heta)- ext{Cl}_{2m}(pi- heta)Big] ]

    Let ( heta =0), to obtain the obvious results ((n=) integer):

    [ ext{Cl}_{2m}(0) = 0~,~ ext{Cl}_{2m}(pi n) = 0 ]

    Let ( heta =pi/3), to obtain

    [ ext{Cl}_{2m}left ( frac{pi }{3} ight )=frac{(1+2^{2m-1})}{2^{2m-1}} ext{Cl}_{2m}left ( frac{2pi }{3} ight ) ]

    Similarly, let ( heta =pi/4), to obtain

    [ ext{Cl}_{2m}left ( frac{pi }{2} ight )=2^{2m-1} ext{Cl}_{2m}left ( frac{pi }{4} ight )-2^{2m-1} ext{Cl}_{2m}left ( frac{3pi }{4} ight ) ]

    The leftmost Clausen function is expressible in terms of the Dirichlet Beta function:

    [eta(x)=sum_{k=0}^{infty}frac{(-1)^k}{(2k+1)^x} ]

    Since

    [ ext{Cl}_{2m}left ( frac{pi }{2} ight )=sum_{k=1}^{infty}frac{sin (pi k/2)}{k^{2m}}=1-frac{1}{3^{2m}}+frac{1}{5^{2m}}-frac{1}{7^{2m}}+\,cdots\, =eta(2m) ]

    So, in summation, thus far we have the following useful values/relations:
    (Largemathbf{color{Purple}{Result ~3:}})

    [ ext{Cl}_{2m}(0) = 0~,~ ext{Cl}_{2m}(pi n) = 0~,~ ext{Cl}_{2m}left ( frac{pi }{2} ight ) = eta(2m) ]

    [ ext{Cl}_{2m}left ( frac{pi }{3} ight )=frac{(1+2^{2m-1})}{2^{2m-1}} ext{Cl}_{2m}left ( frac{2pi }{3} ight )~,~ ext{Cl}_{2m}left ( frac{pi }{4} ight )=frac{eta(2m)}{2^{2m-1}}+ ext{Cl}_{2m}left ( frac{3pi }{4} ight ) ]


    Continuing on with the first of the two duplication formulae, note that we can apply it to itself as follows:

    [egin{align*} & ext{Cl}_{2m}(2 heta)=2^{2m-1}left[ ext{Cl}_{2m}( heta)- ext{Cl}_{2m}(pi- heta) ight]\ Rightarrow & ext{Cl}_{2m}(4 heta)=2^{2m-1}left[ ext{Cl}_{2m}(2 heta)- ext{Cl}_{2m}(pi-2 heta) ight]\ =&2^{2m-1}left[2^{2m-1}left( ext{Cl}_{2m}( heta)- ext{Cl}_{2m}(pi- heta) ight) ight]-2^{2m-1} ext{Cl}_{2m}(2(pi/2- heta))\ =&4^{2m-1}left[ ext{Cl}_{2m}( heta)- ext{Cl}_{2m}(pi- heta) ight]-4^{2m-1}left[ ext{Cl}_{2m}(pi/2- heta)- ext{Cl}_{2m}(pi/2+ heta) ight] end{align*}]

    (Largemathbf{color{Purple}{Result ~4:}})

    [ ext{Cl}_{2m}(4 heta)=4^{2m-1}Bigg[ ext{Cl}_{2m}( heta)- ext{Cl}_{2m}(pi- heta) - ext{Cl}_{2m}left(frac{pi}{2}- heta ight)+ ext{Cl}_{2m}left(frac{pi}{2}+ heta ight)Bigg] ]


    By applying the exact same process to the second of the two duplication formulae, then for CL-type Clausen function of odd index we have:
    (Largemathbf{color{Purple}{Result ~5:}})

    [egin{align*} & ext{Cl}_{2m+1}(4 heta)\ &=zeta(2m+1)+4^{2m}Bigg[ ext{Cl}_{2m+1}( heta)- ext{Cl}_{2m+1}(pi- heta) - ext{Cl}_{2m+1}left(frac{pi}{2}- heta ight)+ ext{Cl}_{2m+1}left(frac{pi}{2}+ heta ight)Bigg] end{align*}]


    In (mathbf{color{Purple}{Result~4}}), set (displaystyle heta=pi/8\,) to obtain:

    [ ext{Cl}_{2m}left(frac{pi}{2} ight)=4^{2m-1}left[ ext{Cl}_{2m}left(frac{pi}{8} ight)- ext{Cl}_{2m}left(frac{7pi}{8} ight) - ext{Cl}_{2m}left(frac{3pi}{8} ight)+ ext{Cl}_{2m}left(frac{5pi}{8} ight) ight] ]

    We've already seen that the leftmost Clausen function is expressible in terms of the Dirichlet Beta function, so we quickly conclude that:
    (Largemathbf{color{Purple}{Result ~6:}})

    [ ext{Cl}_{2m}left(frac{pi}{8} ight) - ext{Cl}_{2m}left(frac{3pi}{8} ight) + ext{Cl}_{2m}left(frac{5pi}{8} ight) - ext{Cl}_{2m}left(frac{7pi}{8} ight)=frac{eta(2m)}{4^{2m-1}} ]

    And in particular:

    [ ext{Cl}_{2}left(frac{pi}{8} ight) - ext{Cl}_{2}left(frac{3pi}{8} ight) + ext{Cl}_{2}left(frac{5pi}{8} ight) - ext{Cl}_{2}left(frac{7pi}{8} ight)=frac{eta(2)}{4}=frac{mathbf{G}}{4} ]


    To express a CL-type Clausen function of rational argument and even order as a sum of polygamma functions, proceed as follows:

    [ ext{Cl}_{2m}left(frac{ppi}{q} ight)=sum_{k=1}^{infty}frac{sin(ppi/q)}{k^{2m}} ]

    Split this into (q) parts

    [egin{align*} &sum_{k=0}^{infty}frac{sinleft[(kq+1)dfrac{ppi}{q} ight]}{(kq+1)^{2m}}+ sum_{k=0}^{infty}frac{sinleft[(kq+2)dfrac{ppi}{q} ight]}{(kq+2)^{2m}}+ sum_{k=0}^{infty}frac{sinleft[(kq+3)dfrac{ppi}{q} ight]}{(kq+3)^{2m}}+ \,cdots\,+\ &sum_{k=0}^{infty}frac{sinleft[(kq+(q-1))dfrac{ppi}{q} ight]}{(kq+(q-1))^{2m}}+ sum_{k=0}^{infty}frac{sinleft[(kq+q)dfrac{ppi}{q} ight]}{(kq+q)^{2m}} end{align*}]

    Convert into a finite sum:

    [egin{align*} &sum_{j=1}^{j=q}sum_{k=0}^{infty}frac{sinleft(kppi+dfrac{jppi}{q} ight)}{(kq+j)^{2m}}=frac{1}{q^{2m}}sum_{j=1}^{j=q}sum_{k=0}^{infty}frac{sinleft(kppi+dfrac{jppi}{q} ight)}{(k+j/q)^{2m}}=frac{1}{q^{2m}}sum_{j=1}^{j=q}sum_{k=0}^{infty}frac{sinleft(kppi+dfrac{jppi}{q} ight)}{(k+j/q)^{2m}}\ &=frac{1}{q^{2m}}sum_{j=1}^{j=q}sum_{k=0}^{infty}frac{cos(kppi)sinleft(dfrac{jppi}{q} ight)}{(k+j/q)^{2m}}=frac{1}{q^{2m}}sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,sum_{k=0}^{infty}frac{(-1)^{kp}}{(k+j/q)^{2m}} end{align*}]

    We now have two distinct cases to consider, since

    [(-1)^{kp} = egin{cases} 1 & ext{if }p ext{ is even} \ (-1)^k & ext{if }p ext{ is odd} end{cases} ]

    (mathbf{color{Teal}{Case ~1. ~when ~oldsymbol{p} ~is ~even:}})
    This is, naturally, the simpler of the two, since the result above reduces to:

    [frac{1}{q^{2m}}sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,sum_{k=0}^{infty}frac{1}{(k+j/q)^{2m}} ]

    We then apply the series definition for the polygamma function

    [psi_n(z)= (-1)^{n+1}\,n!\,sum_{k=0}^{infty}frac{1}{(z+k)^{n+1}} ]

    to obtain

    [frac{1}{q^{2m}}sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,frac{psi_{2m-1}left ( dfrac{j}{q} ight )}{(2m-1)!} ]

    This, in turn, can be expressed in the infinitely more elegant form:
    (Largemathbf{color{Purple}{Result ~7A:}})

    [(2m-1)!\,q^{2m}\, ext{Cl}_{2m}left(frac{ppi}{q} ight)=sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,psi_{2m-1}left(frac{j}{q} ight) ]


    (mathbf{color{Teal}{Case ~2. ~when ~oldsymbol{p} ~is ~odd:}})
    When (p) is odd, we are faced with an alternating series, which must be split in two:

    [egin{align*} &frac{1}{q^{2m}}sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,sum_{k=0}^{infty}frac{(-1)^{k}}{(k+j/q)^{2m}}\ =&frac{1}{q^{2m}}sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,left[ sum_{k=0}^{infty}frac{1}{(2k+j/q)^{2m}}- sum_{k=0}^{infty}frac{1}{(2k+1+j/q)^{2m}} ight]\ =&frac{1}{(2q)^{2m}}sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,left[ sum_{k=0}^{infty}frac{1}{[k+(j/2q)]^{2m}}- sum_{k=0}^{infty}frac{1}{[k+((j+q)/2q)]^{2m}} ight]\ =&frac{1}{(2q)^{2m}}sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,left[ frac{psi_{2m-1}left(dfrac{j}{2q} ight)-psi_{2m-1}left(dfrac{j+q}{2q} ight) }{(2m-1)!} ight] end{align*}]

    From which we obtain the final result:
    (Largemathbf{color{Purple}{Result ~7B:}})

    [displaystyle (2m-1)!\,(2q)^{2m}\, ext{Cl}_{2m}left(frac{ppi}{q} ight)=sum_{j=1}^{j=q}\,sinleft(frac{jppi}{q} ight)\,left[ psi_{2m-1}left(frac{j}{2q} ight)-psi_{2m-1}left(frac{j+q}{2q} ight) ight] ]


    Here's an extremely elementary application of trigonometry to express a couple of trig series in terms of Clausen functions. Start of with a CL-type Clausen function of odd order and double argument:

    [ ext{Cl}_{2m+1}(2 heta)=sum_{k=1}^{infty}frac{cos 2k heta}{k^{2m+1}} ]

    Next, split the sum using the double angle formula for the cosine:

    [ ext{Cl}_{2m+1}(2 heta)=sum_{k=1}^{infty}frac{(cos^2 k heta-sin^2 k heta)}{k^{2m+1}}= sum_{k=1}^{infty}frac{cos^2 k heta}{k^{2m+1}}- sum_{k=1}^{infty}frac{sin^2 k heta}{k^{2m+1}} ]

    Again, using nothing but the most basic level of trigonometry, ie (displaystyle sin^2 x+cos^2 x=1\,), we see that the sum of the last two series must be

    [sum_{k=1}^{infty}frac{cos^2 k heta}{k^{2m+1}}+ sum_{k=1}^{infty}frac{sin^2 k heta}{k^{2m+1}}=sum_{k=1}^{infty}frac{1}{k^{2m+1}}=zeta(2m+1) ]

    Taking the sum/difference of these two results gives:
    (Largemathbf{color{Purple}{Result ~8:}})

    [sum_{k=1}^{infty}frac{sin^2 k heta}{k^{2m+1}}=frac{1}{2}Big[zeta(2m+1)- ext{Cl}_{2m+1}(2 heta)Big] ]

    [sum_{k=1}^{infty}frac{cos^2 k heta}{k^{2m+1}}=frac{1}{2}Big[zeta(2m+1)+ ext{Cl}_{2m+1}(2 heta)Big] ]


    As mentioned right at the start, there are deep connections between the Clausen functions and - amongst others - polylogarithms. Here's a simple example.
    For (displaystyle |z|le 1), the Polylogarithm of order (m) has the series expansion:

    [ ext{Li}_m(z)=sum_{k=1}^{infty}frac{z^k}{k^m} ]

    Setting (displaystyle z=e^{i heta}\,) in the series above gives:

    [egin{align*} ext{Li}_m(e^{i heta})&=sum_{k=1}^{infty}frac{(e^{i heta})^k}{k^m}=sum_{k=1}^{infty}frac{(cos heta+isin heta)^k}{k^m}\ &=sum_{k=1}^{infty}frac{(cos k heta+isin k heta)}{k^m}=sum_{k=1}^{infty}frac{cos k heta}{k^m}+i\,sum_{k=1}^{infty}frac{sin k heta}{k^m} end{align*}]

    So, for a polylogarithm of even order, we have:

    [ ext{Li}_{2m}(e^{i heta})=sum_{k=1}^{infty}frac{cos k heta}{k^{2m}}+i\,sum_{k=1}^{infty}frac{sin k heta}{k^{2m}}= ext{Sl}_{2m}( heta)+i\, ext{Cl}_{2m}( heta) ]

    And similarly, for a polylogarithm of odd order:

    [ ext{Li}_{2m+1}(e^{i heta})=sum_{k=1}^{infty}frac{cos k heta}{k^{2m+1}}+i\,sum_{k=1}^{infty}frac{sin k heta}{k^{2m+1}}= ext{Cl}_{2m+1}( heta)+i\, ext{Sl}_{2m+1}( heta) ]

    (Largemathbf{color{Purple}{Result ~9:}})

    [egin{align*} & ext{Li}_{2m}(e^{i heta})= ext{Sl}_{2m}( heta)+i\, ext{Cl}_{2m}( heta)~,~ ext{Li}_{2m+1}(e^{i heta})= ext{Cl}_{2m+1}( heta)+i\, ext{Sl}_{2m+1}( heta)\ & ext{Li}_{2m}(e^{-i heta})= ext{Sl}_{2m}( heta)-i\, ext{Cl}_{2m}( heta)~,~ ext{Li}_{2m+1}(e^{-i heta})= ext{Cl}_{2m+1}( heta)-i\, ext{Sl}_{2m+1}( heta) end{align*}]


    As later posts in this thread will demonstrate, the (definite integral) moments of Clausen functions have a number of important applications, so it would, therefore, be quite useful to find a closed form for these moments:

    [mathbf{Cl}_{(m, n)}( heta)=int_0^{ heta}x^n\, ext{Cl}_{m}(x)\,mathrm{d}x ]

    In order to find these, we start off with the integral definition of (displaystyle ext{Cl}_2( heta)):

    [ ext{Cl}_2( heta)=-int_0^{ heta}lnBig|2sinfrac{x}{2}Big|\,mathrm{d}x ]

    We could just as easily use the series definition, but it important to observe that, within the range (displaystyle 0 < hetale 2pi\,), the absolute value of the logsine term can be omitted. Hence while within this range, we can write

    [ ext{Cl}_2( heta)=-int_0^{ heta}lnleft(2sinfrac{x}{2} ight)\,mathrm{d}x~,~ ext{Cl}_1( heta)=-lnleft(2sinfrac{ heta}{2} ight) ]

    with no loss of generality. Note that the omission of ( heta =0) was to avoid the divergence of (displaystyle ext{Cl}_1( heta)\,) at this point (although, in the integrals presented below, this singularity is of no consequence).
    The generalized moments of Clausen functions will be dealt with shortly, but for the moment, let's consider the special - and perhaps most important - case of the moments of (displaystyle ext{Cl}_1( heta)\,).

    [mathbf{Cl}_{(1, n)}( heta)=int_0^{ heta}x^n\, ext{Cl}_{1}(x)\,mathrm{d}x= -int_0^{ heta}x^nlnleft(2sinfrac{x}{2} ight)\,mathrm{d}x ]

    Applying the series definition for (displaystyle ext{Cl}_1( heta)) this can be written as

    [mathbf{Cl}_{(1, n)}( heta)=-sum_{k=1}^{infty}frac{1}{k}int_0^{ heta}x^ncos kx\,mathrm{d}x ]

    Although tedious to derive, the closed form of that last trigonometric integral is easily deduced. Define

    [mathcal{I}_{(p)}=int_0^{ heta}x^pcos kx\,mathrm{d}x ]

    Then we note that, for (p > 2):

    [mathcal{I}_{(p)}=frac{x^p}{k}sin kx\,Biggr|_0^{ heta}+frac{px^{p-1}}{k^2}cos kx\,Biggr|_0^{ heta}-frac{p(p-1)}{k^2}mathcal{I}_{(p-2)} ]

    A similarly, for (p > 4):

    [egin{align*} mathcal{I}_{(p)}&=frac{x^p}{k}sin kx\,Biggr|_0^{ heta}+frac{px^{p-1}}{k^2}cos kx\,Biggr|_0^{ heta}- frac{p(p-1)x^{p-2}}{k^3}sin kx\,Biggr|_0^{ heta}\ &~~~-frac{p(p-1)(p-2)x^{p-3}}{k^4}cos kx\,Biggr|_0^{ heta}+frac{p(p-1)(p-2)(p-3)}{k^4}mathcal{I}_{(p-4)} end{align*}]

    Using the iteration process above, and creating two sums - one for sine terms and the other for cosine terms - we get the desired result:

    [mathcal{I}_{(p)}=p!\,left[sum_{j=0}^{lfloor {p/2} floor}(-1)^jfrac{x^{p-2j}sin kx}{k^{2j+1}(p-2j)!}Biggr|_0^{ heta}+ sum_{j=0}^{lfloor {(p-1)/2} floor}(-1)^jfrac{x^{p-2j-1}cos kx}{k^{2j+2}(p-2j-1)!)}Biggr|_0^{ heta} ight] ]

    Inserting this into our moment integral gives:

    [egin{align*} &mathbf{Cl}_{(1, p)}( heta)=int_0^{ heta}x^p\, ext{Cl}_{1}(x)\,mathrm{d}x= -int_0^{ heta}x^plnleft(2sinfrac{x}{2} ight)\,mathrm{d}x=-sum_{k=1}^{infty}frac{1}{k}int_0^{ heta}x^pcos kx\,mathrm{d}x\ &=-p!\,sum_{k=1}^{infty}frac{1}{k}\,left[sum_{j=0}^{lfloor {p/2} floor}(-1)^jfrac{x^{p-2j}sin kx}{k^{2j+1}(p-2j)!}Biggr|_0^{ heta}+ sum_{j=0}^{lfloor {(p-1)/2} floor}(-1)^jfrac{x^{p-2j-1}cos kx}{k^{2j+2}(p-2j-1)!)}Biggr|_0^{ heta} ight] end{align*}]

    For the lower bound - (x=0) - the terms in the (finite!) sine series vanish. The same is true for all the lower bound term in the (finite) cosine series, except for the final term - containing (x_0), which is present only when (p=2m+1) is odd, in which case this final cosine term is (displaystyle (-1)^{lfloor (p-1)/2 floor +1} / k^{p+1}\,). To account for this term, we introduce the function

    [frac{[1+(-1)^{p+1}]}{2} = egin{cases} 0, & ext{if }p ext{ is even} \ 1, & ext{if }p ext{ is odd} end{cases} ]

    Finally, expressing our sum in terms of Clausen functions, we arrive at:
    (Largemathbf{color{Purple}{Result ~10:}})

    [egin{align*} &int_0^{ heta}x^p\, ext{Cl}_{1}(x)mathrm{d}x\ &=-p!\,left[ sum_{j=0}^{lfloor {p/2} floor}(-1)^jfrac{{ heta}^{\,p-2j}}{(p-2j)!} ext{Cl}_{2j+2}( heta)+ sum_{j=0}^{lfloor {(p-1)/2} floor}(-1)^jfrac{{ heta}^{\,p-2j-1}}{(p-2j-1)!)} ext{Cl}_{2j+3}( heta) ight]\ &~~~+p! (-1)^{lfloor (p-1)/2 floor} frac{[1+(-1)^{p+1}]}{2}zeta(p+2) end{align*}]


    A number of simple - but nonetheless important - trigonometric integrals follow immediately from the previous evaluation. For example

    [egin{align*} &int_0^{ heta}x^p\, ext{Cl}_{1}(x)\,mathrm{d}x= -int_0^{ heta}x^plnleft(2sinfrac{x}{2} ight)\,mathrm{d}x\ &=-frac{ heta^{\,p+1}}{p+1}lnleft(2sinfrac{ heta}{2} ight)+frac{1}{2(p+1)}\,int_0^{ heta}x^{p+1}cotfrac{x}{2}\,mathrm{d}x end{align*}]

    Which can be re-written in the more convenient form:
    (Largemathbf{color{Purple}{Result ~11:}})

    [egin{align*} &int_0^{phi}x^{p+1}cot x\,mathrm{d}x\ &=(p+1)! frac{(-1)^{lfloor (p-1)/2 floor}[1+(-1)^{p+1}]}{2^{p+2}}zeta(p+2)+phi^{p+1}ln(2sinphi)\ &~~~-(p+1)!\,sum_{j=0}^{lfloor {p/2} floor}(-1)^jfrac{{phi}^{\,p-2j}}{2^{2j+1}(p-2j)!} ext{Cl}_{2j+2}(phi)\ &~~~-(p+1)!\,sum_{j=0}^{lfloor {(p-1)/2} floor}(-1)^jfrac{{phi}^{\,p-2j-1}}{2^{2j+2}(p-2j-1)!} ext{Cl}_{2j+3}(phi) end{align*}]


    (mathbf{color{DarkOrange}{Logcosine ~moments - ~part ~1:}})
    Following on from the logsine moments above, we find that the equivalent logcosine moments are slightly trickier - which was to be expected - although they are far richer, since the complex parts also yield useful information.
    By analogy, we start off with:

    [egin{align*} &int_0^{ heta}x^mlnleft(2cosfrac{x}{2} ight)\,mathrm{d}x=int_0^{ heta}x^mlnleft(frac{1+e^{-ix}}{e^{-ix/2}} ight)\,mathrm{d}x\ &=int_0^{ heta}x^mln(1+e^{-ix})\,mathrm{d}x+i\,frac{ heta^{m+2}}{2(m+2)} end{align*}]

    We'll ignore that final complex function of ( heta) for now, and continue with the evaluation of the complex logarithmic integral part:

    [egin{align*} &int_0^{ heta}x^mln(1+e^{-ix})\,mathrm{d}x=sum_{k=1}^{infty}frac{(-1)^{k+1}}{k}int_0^{ heta}x^m(cos kx -i\,sin kx)\,mathrm{d}x\ &=sum_{k=1}^{infty}frac{(-1)^{k+1}}{k}int_0^{ heta}x^mcos kx\,mathrm{d}x-i\,sum_{k=1}^{infty}frac{(-1)^{k+1}}{k}int_0^{ heta}x^msin kx\,mathrm{d}x\ &=int_0^{ heta}x^mleft[ ext{Cl}_1(x)-frac{1}{2} ext{Cl}_1(2x) ight]\,mathrm{d}x-i\,int_0^{ heta}x^mleft[ ext{Sl}_1(x)-frac{1}{2} ext{Sl}_1(2x) ight]\,mathrm{d}x\ &sum_{k=1}^{infty}frac{1}{k}int_0^{ heta}x^mleft(cos kx-frac{cos 2kx}{2} ight)\,mathrm{d}x-i\,sum_{k=1}^{infty}frac{1}{k}int_0^{ heta}x^mleft(sin kx-frac{sin 2kx}{2} ight)\,mathrm{d}x end{align*}]

    We already have the closed form for the two leftmost CL-type integrals, so it remains to find the remaining two SL-types. [The real logcosine moments are omitted below, since they are easily deduced from the previous result]. Let

    [mathcal{I}_{(p)}=int_0^{ heta}x^psin kx\,mathrm{d}x=-frac{1}{k}cos kxBiggr|_0^{ heta}+frac{px^{p-1}}{k^2}sin kxBiggr|_0^{ heta}-frac{p(p-1)}{k^2}mathcal{I}_{(p-2)}=\,cdots ]

    As before, a similar iteration process leads to the general result:

    [egin{align*} &int_0^{ heta}x^psin kx\,mathrm{d}x= \ &p!\,left[sum_{j=0}^{lfloor {p/2} floor}(-1)^{j+1}frac{x^{p-2j}}{k^{2j+1}(p-2j)!}cos kx \,Biggr|_0^{ heta}+sum_{j=0}^{lfloor {(p-1)2} floor}(-1)^{j+1}frac{x^{p-2j-1}}{k^{2j+2}(p-2j-1)!}sin kx\,Biggr|_0^{ heta} ight] end{align*}]

    So the imaginary part of (displaystyle int_0^{ heta}x^plnleft(2cosfrac{x}{2} ight)\,mathrm{d}x) yields

    [egin{align*} &p!\,sum_{k=0}^{infty}frac{1}{k}\,sum_{j=0}^{lfloor {p/2} floor}(-1)^{j+1}frac{x^{p-2j}}{k^{2j+1}(p-2j)!}cos kx \,Biggr|_0^{ heta}\ &+p!\,sum_{k=0}^{infty}frac{1}{k}\,sum_{j=0}^{lfloor {(p-1)2} floor}(-1)^{j+1}frac{x^{p-2j-1}}{k^{2j+2}(p-2j-1)!}sin kx\,Biggr|_0^{ heta}\ &+frac{p!}{2}\,sum_{k=0}^{infty}frac{1}{k}\,sum_{j=0}^{lfloor {p/2} floor}(-1)^{j+1}frac{x^{p-2j}}{(2k)^{2j+1}(p-2j)!}cos 2kx \,Biggr|_0^{ heta}\ &-frac{p!}{2}\,sum_{k=0}^{infty}frac{1}{k}\,sum_{j=0}^{lfloor {(p-1)2} floor}(-1)^{j+1}frac{x^{p-2j-1}}{(2k)^{2j+2}(p-2j-1)!}sin 2kx\,Biggr|_0^{ heta}+frac{ heta^{m+2}}{2(m+2)}=0 end{align*}]

    [egin{align*} &p!\,sum_{j=0}^{lfloor {p/2} floor}(-1)^{j+1}frac{x^{p-2j}}{(p-2j)!} ext{Sl}_{2j+2}(x)\,Biggr|_0^{ heta}+p!\,sum_{j=0}^{lfloor {(p-1)2} floor}(-1)^{j+1}frac{x^{p-2j-1}}{(p-2j-1)!} ext{Sl}_{2j+3}(x)\,Biggr|_0^{ heta}\ &-frac{p!}{2}\,sum_{j=0}^{lfloor {p/2} floor}(-1)^{j+1}frac{x^{p-2j}}{2^{2j+1}(p-2j)!} ext{Sl}_{2j+2}(2x)\,Biggr|_0^{ heta}\ &-frac{p!}{2}\,sum_{j=0}^{lfloor {(p-1)2} floor}(-1)^{j+1}frac{x^{p-2j-1}}{2^{2j+2}(p-2j-1)!} ext{Sl}_{2j+3}(2x)\,Biggr|_0^{ heta}+frac{ heta^{m+2}}{2(m+2)}=0 end{align*}]

    As promised before - on other threads, and indeed other forums - I'll start to find closed form expressions for various polygamma functions, at the rational arguments $1/2, 1/3, 2/3, 1/4 (,) 3/4, 1/6$, and (5/6). This might take a while, and be posted in stages.
    To start with, let's consider the following particular Clausen function of (arbitrary) odd order:

    [ ext{Cl}_{2m+1}left(frac{pi}{3} ight)=sum_{k=1}^{infty}frac{cos (pi k/3)}{k^{2m+1}} ]

    We want to split this into six sums, where the first sum contains the first of every six terms, the second contains the second of every six terms, and so on. We also change summation index so our new series start at (k=0), rather than (k=1) above.

    [egin{align*} & ext{Cl}_{2m+1}left(frac{pi}{3} ight)=sum_{k=1}^{infty}frac{cos (pi k/3)}{k^{2m+1}}=\ &sum_{k=0}^{infty}frac{cosdfrac{pi}{3}(6k+1)}{(6k+1)^{2m+1}}+ sum_{k=0}^{infty}frac{cosdfrac{pi}{3}(6k+2)}{(6k+2)^{2m+1}}+ sum_{k=0}^{infty}frac{cosdfrac{pi}{3}(6k+3)}{(6k+3)^{2m+1}}+\ &sum_{k=0}^{infty}frac{cosdfrac{pi}{3}(6k+4)}{(6k+4)^{2m+1}}+ sum_{k=0}^{infty}frac{cosdfrac{pi}{3}(6k+5)}{(6k+5)^{2m+1}}+ sum_{k=0}^{infty}frac{cosdfrac{pi}{3}(6k+6)}{(6k+6)^{2m+1}} end{align*}]

    Simplify the trig term in each series:

    [egin{align*} cos frac{pi}{3}(6k+n)&=cosleft(2pi k+frac{pi n}{3} ight)\ &=cos 2pi kcosfrac{pi n}{3}-sin 2pi ksinfrac{pi n}{3}equiv cosfrac{pi n}{3} end{align*}]

    Our new sextet of series is thus

    [egin{align*} & ext{Cl}_{2m+1}left(frac{pi}{3} ight)=cosleft(frac{pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(6k+1)^{2m+1}}+ cosleft(frac{2pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(6k+2)^{2m+1}}\ &+cosleft(pi ight)\, sum_{k=0}^{infty}frac{1}{(6k+3)^{2m+1}}+ cosleft(frac{4pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(6k+4)^{2m+1}}\ &+cosleft(frac{5pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(6k+5)^{2m+1}}+ cosleft(2pi ight)\, sum_{k=0}^{infty}frac{1}{(6k+6)^{2m+1}} end{align*}]

    Multiply both sides by (6^{2m+1}), and then subtract the third and sixth series on the RHS from the Clausen term on the LHS (using (displaystyle cospi = -1\,) and (displaystyle cos 2pi=1\,) ) to obtain:

    [egin{align*} &6^{2m+1} ext{Cl}_{2m+1}left(frac{pi}{3} ight)+sum_{k=0}^{infty}frac{1}{(k+1/2)^{2m+1}}-sum_{k=0}^{infty}frac{1}{(k+1)^{2m+1}}\ &=cosleft(frac{pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(k+1/6)^{2m+1}}+ cosleft(frac{2pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(k+1/3)^{2m+1}}\ &~~~+cosleft(frac{4pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(k+2/3)^{2m+1}}+cosleft(frac{5pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(k+5/6)^{2m+1}} end{align*}]

    Express the cosine terms on the RHS in real/rational form to make the RHS

    [egin{align*} &frac{1}{2}\, sum_{k=0}^{infty}frac{1}{(k+1/6)^{2m+1}}- frac{1}{2}\, sum_{k=0}^{infty}frac{1}{(k+1/3)^{2m+1}}\ -&frac{1}{2}\, sum_{k=0}^{infty}frac{1}{(k+2/3)^{2m+1}}+frac{1}{2}\, sum_{k=0}^{infty}frac{1}{(k+5/6)^{2m+1}} end{align*}]

    Now use

    [psi_{nge 1}(x)=(-1)^{n+1}n!sum_{k=0}^{infty}frac{1}{(k+x)^{n+1}} ]

    to re-write the RHS as:

    [egin{align*} &frac{1}{2}left(frac{(-1)^{2m}}{(2m)!} ight)Bigg{ psi_{2m}left( frac{1}{6} ight) -psi_{2m}left( frac{1}{3} ight) -psi_{2m}left( frac{2}{3} ight) +psi_{2m}left( frac{5}{6} ight) Bigg}\ =&frac{1}{2\,(2m)!}\, Bigg{ psi_{2m}left( frac{1}{6} ight) -psi_{2m}left( frac{1}{3} ight) -psi_{2m}left( frac{2}{3} ight) +psi_{2m}left( frac{5}{6} ight) Bigg} end{align*}]

    Next, apply the same process to the two series on the LHS (with the Clausen term):

    [egin{align*} &6^{2m+1} ext{Cl}_{2m+1}left(frac{pi}{3} ight)+sum_{k=0}^{infty}frac{1}{(k+1/2)^{2m+1}}-sum_{k=0}^{infty}frac{1}{(k+1)^{2m+1}}\ =&6^{2m+1} ext{Cl}_{2m+1}left(frac{pi}{3} ight)+frac{1}{(2m)!}Bigg{psi_{2m}left(frac{1}{2} ight)-psi_{2m}(1)Bigg} end{align*}]

    Multiplying both sides by (2(2m)!) then gives the identity

    [egin{align*} &2\, (2m)! \,6^{2m+1} ext{Cl}_{2m+1}left(frac{pi}{3} ight)+2psi_{2m}left(frac{1}{2} ight)-2psi_{2m}(1)\ =& psi_{2m}left( frac{1}{6} ight) -psi_{2m}left( frac{1}{3} ight) -psi_{2m}left( frac{2}{3} ight) +psi_{2m}left( frac{5}{6} ight) end{align*}]

    Next, repeat all of the above, but this time in terms of the Clasuen function with argument (2pi /3)

    [egin{align*} & ext{Cl}_{2m+1}left(frac{2pi}{3} ight)=sum_{k=1}^{infty}frac{cos (2pi k/3)}{k^{2m+1}}\ &=cosleft(frac{2pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(6k+1)^{2m+1}}+ cosleft(frac{4pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(6k+2)^{2m+1}}+\ &~~~cosleft(2pi ight)\, sum_{k=0}^{infty}frac{1}{(6k+3)^{2m+1}}+ cosleft(frac{8pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(6k+4)^{2m+1}}+\ &~~~cosleft(frac{10pi}{3} ight)\, sum_{k=0}^{infty}frac{1}{(6k+5)^{2m+1}}+ cosleft(4pi ight)\, sum_{k=0}^{infty}frac{1}{(6k+6)^{2m+1}}\ &=-frac{1}{2}\, sum_{k=0}^{infty}frac{1}{(6k+1)^{2m+1}} -frac{1}{2}\, sum_{k=0}^{infty}frac{1}{(6k+2)^{2m+1}}+sum_{k=0}^{infty}frac{1}{(6k+3)^{2m+1}}\ &~~~- frac{1}{2}\, sum_{k=0}^{infty}frac{1}{(6k+4)^{2m+1}}-frac{1}{2}\, sum_{k=0}^{infty}frac{1}{(6k+5)^{2m+1}}+ sum_{k=0}^{infty}frac{1}{(6k+6)^{2m+1}} end{align*}]

    Continue exactly as before, and you get the second relation

    [egin{align*} &2\, (2m)! \,6^{2m+1} ext{Cl}_{2m+1}left(frac{2pi}{3} ight)-2psi_{2m}left(frac{1}{2} ight)-2psi_{2m}(1)\ &=-psi_{2m}left( frac{1}{6} ight) -psi_{2m}left( frac{1}{3} ight) -psi_{2m}left( frac{2}{3} ight) -psi_{2m}left( frac{5}{6} ight) end{align*}]

    Relative to the arguments (1/3, 2/3, 1/6), and (5/6), the arguments (1) and (1/2) are pretty straightforward, so I'll simply state them now and prove them later.

    [psi_{2m}(1)=-(2m)!zeta(2m+1)~,~psi_{2m}left(frac{1}{2} ight)=-(2m)!\,(2^{2m+1}-1)zeta(2m+1) ]

    Now, if you add the final forms or relation (1) and relation (2) you get:

    [egin{align*} &2\, (2m)! \,6^{2m+1}left[ ext{Cl}_{2m+1}left(frac{pi}{3} ight)+ ext{Cl}_{2m+1}left(frac{2pi}{3} ight) ight]-4psi_{2m}(1)\ =&2\, (2m)! \,6^{2m+1}left[ ext{Cl}_{2m+1}left(frac{pi}{3} ight)+ ext{Cl}_{2m+1}left(frac{2pi}{3} ight) ight]+4\,(2m)!zeta(2m+1)\ =&2Bigg{ psi_{2m}left( frac{1}{3} ight) +psi_{2m}left( frac{2}{3} ight)Bigg} end{align*}]

    Or

    [egin{align*} &(2m)! \,6^{2m+1}left[ ext{Cl}_{2m+1}left(frac{pi}{3} ight)+ ext{Cl}_{2m+1}left(frac{2pi}{3} ight) ight]+2\,(2m)!zeta(2m+1)\ =&psi_{2m}left( frac{1}{3} ight) +psi_{2m}left( frac{2}{3} ight) end{align*}]

    On the other hand, the reflection formula for the polygamma function gives:

    [egin{align*} &psi_{2m}(x)-psi_{2m}(1-x)=pifrac{mathrm{d}^{2m}}{mathrm{d}x^{2m}}cotpi x\ Rightarrow &psi_{2m}left( frac{1}{3} ight) -psi_{2m}left( frac{2}{3} ight)=pifrac{mathrm{d}^{2m}}{mathrm{d}x^{2m}}cotpi x\,Biggr|_{x=1/3} end{align*}]

    So

    [egin{align*} &psi_{2m}left( frac{1}{3} ight)\ =&displaystyle frac{(2m)! \,6^{2m+1}}{2}left[ ext{Cl}_{2m+1}left(frac{pi}{3} ight)+ ext{Cl}_{2m+1}left(frac{2pi}{3} ight) ight]+\,(2m)!zeta(2m+1)+pifrac{mathrm{d}^{2m}}{mathrm{d}x^{2m}}cotpi x\,Biggr|_{x=1/3} end{align*}]

    and

    [egin{align*} &psi_{2m}left( frac{2}{3} ight)\ =&frac{(2m)! \,6^{2m+1}}{2}left[ ext{Cl}_{2m+1}left(frac{pi}{3} ight)+ ext{Cl}_{2m+1}left(frac{2pi}{3} ight) ight]+\,(2m)!zeta(2m+1)-pifrac{mathrm{d}^{2m}}{mathrm{d}x^{2m}}cotpi x\,Biggr|_{x=1/3} end{align*}]

  • 相关阅读:
    Kail命令
    [UIScreen mainScreen].applicationFrame与[UIScreen mainScreen].bounds区别
    Excel使用技巧
    vs2012配置OpenGL
    IP首部校验和的计算
    链接保存
    Android精讲--界面编程1(界面编程与视图的组件)
    Intent对象详解
    安卓四大组件之--Service
    android的事件处理机制
  • 原文地址:https://www.cnblogs.com/Renascence-5/p/5479156.html
Copyright © 2011-2022 走看看