3 svar
24 visningar
AlexMu 948
Postad: Idag 03:13 Redigerad: Idag 03:44

Hitta jämna derivator i x = 1 och x = 0 av ett 'jobbigt' polynom

Hej! Jag håller på med en integral och den bisarra processen har slutligen lett mig till att jag behöver alla jämna derivator i x=0x=0 och x=1x=1 av polynomet 

Pkx=xn-1j=0nj+k-n-xk\displaystyle P_k\left(x\right) = x^{n-1}\prod_{j=0}^{n}\left(j+k-n-xk\right) 

(nn är ett heltal k\geq k) för att få fram värdet på min integral

Det känns mycket besvärligt att hitta den mm:te derivatan av denna funktion med produktregeln etc, så jag har försökt tänka på andra sätt. En ide jag hade var att först ta log-derivatan, vilket ger att

Pk'(x)Pk(x)=n-1x-kj=0n1j+k-n-xkQ0(x)\displaystyle \frac{P_k^'(x)}{P_k(x)}= \underbrace{\frac{n-1}{x}-k\sum_{j=0}^{n}\frac{1}{j+k-n-xk}}_{Q_0(x)}

Q0Q_0 är lätt att ta ett arbiträrt antal derivator av, så det har jag försökt utnyttja. Vi låter den mm:te derivatan av Q0Q_0 vara QmQ_m. Multiplicera båda led med Pk(x)P_k(x) vilket ger 

Pk'x=PkxQ0x\displaystyle P_k^'\left(x\right) = P_k\left(x\right)Q_0\left(x\right) och derivering ger

Pk(2)x=Pk'xQ0x+PkxQ1x=PkxQ0x2+Q1x\displaystyle P_k^{(2)}\left(x\right) = P_k^'\left(x\right)Q_0\left(x\right) + P_k\left(x\right)Q_1\left(x\right) = P_k\left(x\right)\left(Q_0\left(x\right)^2 + Q_1\left(x\right)\right)

Sedan kan vi derivera igen vilket, efter lite algebra, ger

Pk(3)x=PkxQ0(x)3+3Q0xQ1x+Q2x\displaystyle P_k^{(3)}\left(x\right) = P_{k}\left(x\right)\left(Q_{0}(x)^{3}+3Q_{0}\left(x\right)Q_{1}\left(x\right)+Q_{2}\left(x\right)\right)

Generellt får vi att Pk(m)x=PkxRmxP_k^{(m)}\left(x\right) = P_k\left(x\right)R_m\left(x\right) där RmxR_m\left(x\right) uppfyller rekursionen

R1x=Q0xR_1\left(x\right) = Q_0\left(x\right)
Rm+1x=Q0xRmx+Rm'xR_{m+1}\left(x\right) = Q_0\left(x\right)R_{m}\left(x\right) + R_{m}^'\left(x\right)

Här har jag dock fastnat. Om man kan hitta en sluten form för RmR_m uttryckt i QiQ_i skulle vi kunna få ett slutet uttryck för derivatorna av PP, men jag har inga ideer på hur man kan lösa sådana här rekursionsproblem. Har ni några ideer med rekursionen / ett annat sätt att hitta derivatorna av PP

SeriousCephalopod Online 2754
Postad: Idag 08:55 Redigerad: Idag 08:57

Så Pk definieras av andra raden och du vill bestämma Pk'(0) och Pk'(1)?

Är det problemet?

För x= 0 fallet är det väl enklast att serieutveckla

Pk(x) = a + bx + ...

Och strunta i att utvärdera koeffixienterna förutom den andra eftersom

Pk'(0) = b

För x=1 kan vi derfiniera 

Rk(x) = Pk(x + 1)

och göra samma sak.

Annars har jag missförstått frågan.

Edit:

Ah missade att det var alla jämna derivator och inte bara några specifika.

SeriousCephalopod Online 2754
Postad: Idag 09:15 Redigerad: Idag 09:47

Då påbörjar jag ett nytt försök.

Jag vill börja med att bara utforska omskrivningar av polynomet så att det blir lite mer symmetriskt. Jag lägger till kn som subscript för att betona att n också är ett index.

Jag skriver bara ut några omskrivningar gjort till att börja med. Nedan är inte någon full lösning.

Pkn(x)=xn-1j=0n(j+k-n-xk)P_{kn}(x) = x^{n - 1} \prod_{j = 0}^n (j + k - n - xk)

Först och främst så ser det ut som att vi kan göra en del algebraiska omskrivningar i ett försök att exponera fler symmetrier

Pkn(x)=xn-1j=0n(j -n+k(1-x))P_{kn}(x) = x^{n - 1} \prod_{j = 0}^n (j  - n + k(1 - x))

Eftersom j räknar från 0 till n så kommer j - n = -(n - j) bara att räkna upp -n, -(n - 1), -(n - 2), ..., -1, 0. Dvs de negativa heltalen med belopp mindre eller lika med n, bara i omvänd ordning så vi kan räkna från litet tills stortt vilket blir lättare att läsa

Pkn(x)=xn-1j=0n(k(1-x)-j)P_{kn}(x) = x^{n - 1} \prod_{j = 0}^n (k(1 - x) - j)

Vi bryter ut k så att vi får

Pkn(x)=1knxn-1j=0n((1-x)-jk)P_{kn}(x) = \cfrac{1}{k^n} x^{n - 1} \prod_{j = 0}^n ((1 - x) - \cfrac{j}{k})

eller om vi flyttar om så att vi lättare kan läsa av nollställena

Pkn(x)=(-1)n+1knxn-1j=0n(x-1+jk)P_{kn}(x) = \cfrac{(-1)^{n + 1}}{k^n} x^{n - 1} \prod_{j = 0}^n (x - 1 + \cfrac{j}{k})

Hittills har jag bara gjort omskrivningar och så får vi se om jag kommer vidare därifrån men jag har svårt att tänka mig att någon av dessa inte är användbar.

SeriousCephalopod Online 2754
Postad: Idag 09:38 Redigerad: Idag 09:39

Jag börjar bli skeptisk om vi kommer hitta någon trevlig formel här.

Att bestämma derivatorna P(s)(0)P^{(s)}(0) är ekvivalent med att hitta de jämna koefficienterna till polynomets utvecklade form eftersom P(s)(0)=s!asP^{(s)}(0) = s! a_{s} där asa_s där asa_s är den s:e koefficienten i utvecklingen (förutsatt att s är mindre än polynomets grad)

Att finna koefficienterna till ett polynom givet dess nollställen är rent aritmetiskt trivialt men finns inte i allmännhet slutna formler som är tillfredställande. 

Vi kan exempelvis använda Vietes formler (https://en.wikipedia.org/wiki/Vieta%27s_formulas)

där vi borde få rent konceptuellt att

an-k=(-1)kan1i1i2...iknj=1krija_{n - k} = (-1)^{k} a_n \sum_{1 \leq i_1 \leq i_2 \leq ... i_k \leq n } \prod_{j = 1}^k r_{i_j}

Kanske att våra nollställen här samspelar på något spännande sätt osm gör att de tar ut varandra men jag skulle själv göra en numerisk studie innan jag började hoppas.

Mitt nästa mål är att se om jag kan göra någon ekvivalent till variabelseparation för att försöka separera inflytandet hos k och n.

Svara
Close