Re: accumulate coefficients of a polynomial
- To: mathgroup at smc.vnet.net
- Subject: [mg102831] Re: [mg102813] accumulate coefficients of a polynomial
- From: Leonid Shifrin <lshifr at gmail.com>
- Date: Fri, 28 Aug 2009 05:43:04 -0400 (EDT)
- References: <200908280444.AAA28832@smc.vnet.net>
Hi,
this will give you the general form of coefficients in terms of rules:
In[1] =
Clear[k];
coeffRules = List @@ ComplexExpand[Re[p[Exp[I Pi/k]]]] /. x_a*y : _ : 1 :>
(x :> y)
Out[1] = {a[10]:>1,a[9]:>Cos[\[Pi]/k],a[8]:>Cos[(2 \[Pi])/k],a[7]:>Cos[(3
\[Pi])/k],a[6]:>Cos[(4 \[Pi])/k],a[5]:>Cos[(5 \[Pi])/k],a[4]:>Cos[(6
\[Pi])/k],a[3]:>Cos[(7 \[Pi])/k],a[2]:>Cos[(8 \[Pi])/k],a[1]:>Cos[(9
\[Pi])/k],a[0]:>Cos[(10 \[Pi])/k]}
This will compute the matrix of vectors of coefficients for 1<=k<=5 (5
instead of 20 just to keep
the output size reasonable):
In[2] = Table[Array[a, {10}] /. coeffRules, {k, 1, 5}]
Out[2] =
{{-1,1,-1,1,-1,1,-1,1,-1,1},{0,1,0,-1,0,1,0,-1,0,1},{-1,-(1/2),1/2,1,1/2,-(1/2),-1,-(1/2),1/2,1},{1/Sqrt[2],1,1/Sqrt[2],0,-(1/Sqrt[2]),-1,-(1/Sqrt[2]),0,1/Sqrt[2],1},{1/4
(1+Sqrt[5]),1/4 (-1+Sqrt[5]),1/4 (1-Sqrt[5]),1/4 (-1-Sqrt[5]),-1,1/4
(-1-Sqrt[5]),1/4 (1-Sqrt[5]),1/4 (-1+Sqrt[5]),1/4 (1+Sqrt[5]),1}}
Regards,
Leonid
On Fri, Aug 28, 2009 at 8:44 AM, BHUPALA <bhupala at gmail.com> wrote:
> I have generated a polynomial as
>
> p[z_] = Sum[a[k]*z^(10 - k), {k, 0, 10}]
>
> to give
>
> z^10 a[0] + z^9 a[1] + z^8 a[2] + z^7 a[3] + z^6 a[4] + z^5 a[5] +
> z^4 a[6] + z^3 a[7] + z^2 a[8] + z a[9] + a[10]
>
> I want to substitute z = Exp[I Pi/k] where k varies from 1 to 20 and
> for each k retain the coefficients of the real part as a vector.
>
> I used the following command for a single iteration (e.g k=6)
>
> ComplexExpand[Re[p[Exp[I Pi/6]]]]]
>
> But how to do it in a loop?
>
> Thanks for any help.
>
> Bhupala
>
>
- References:
- accumulate coefficients of a polynomial
- From: BHUPALA <bhupala@gmail.com>
- accumulate coefficients of a polynomial