7 svar
139 visningar
AlexMu behöver inte mer hjälp
AlexMu 1028
Postad: 22 dec 2025 03:13 Redigerad: 22 dec 2025 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 2766
Postad: 22 dec 2025 08:55 Redigerad: 22 dec 2025 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 2766
Postad: 22 dec 2025 09:15 Redigerad: 22 dec 2025 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 2766
Postad: 22 dec 2025 09:38 Redigerad: 22 dec 2025 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.

AlexMu 1028
Postad: 22 dec 2025 23:31 Redigerad: 22 dec 2025 23:35

Intressanta ideer! Ja, Vietes formler är inte särskilt snygga i ett mer generellt sammanhang. Jag har ingen aning om det finns någon snygg formel, vilket i sig är rätt intressant. Ofta brukar ju rötterna vara det mest önskvärda hos ett polynom. Men här är det mycket svårare att finna derivatorna med hjälp av rötterna än om vi hade alla koefficienter (snyggt) givna. Den slutna formen för integralen jag håller på med är en rätt stor formel:

In=2(n-1)!k=1nnn-k=1n-1Pk(2)1-1k-Pk(2)0π2k2+1-n1-2-1-2ζ2+1\displaystyle I_n = \frac{2}{(n-1)!}\sum_{k=1}^{n}\binom{n}{n-k}\sum_{\ell = 1}^{n} \frac{\left(-1\right)^{\ell}\left(P_k^{(2\ell)}\left(1\right)\left(-1\right)^k-P_k^{(2\ell)}\left(0\right)\right)}{\pi^{2\ell}k^{2\ell + 1-n}}\left(1-2^{-1-2\ell}\right)\zeta\left(2\ell + 1\right)


Det kanske skulle vara lite komiskt att använda Vietes, skulle lägga till några fler summor och produkter!

Jag ska leka runt lite mer med rekursionen för RmR_m, jag tyckte det kändes halvlovande. Det är i alla fall lite enklare uttryck att jobba med. 

SeriousCephalopod 2766
Postad: 23 dec 2025 07:01 Redigerad: 23 dec 2025 07:04

Jag har mer hopp om att kollapsa själva summan eller åtminstine täljarsifferenserna i din integralformel till något enklare uttryck än att erhålla alla individuella derivatavärden

Summan av alla 2n-derivator hos ett polynom är ju effektivt summan av alla jämna koefficienter och summan av alla jämna koefficienter till ett polynom (plus fakulteter) med endast jänna faktorer är bokstavligen p(0)

Man kan även lätt jämnifiera (symmetrisera) ett polynom genom att konstruera p(x) + p(-x) eftersom de udda termerna tar ut varandra.

Detta känns i alla fall som ideer men kan testa att utforska. 

AlexMu 1028
Postad: 25 dec 2025 01:39 Redigerad: 25 dec 2025 03:45

Lite progress med det ursprungliga jag tog upp: 

Det finns en formel som heter Faà di Bruno's formel, som är en generaliserad kedjeregel för den nn:te derivatan. Enligt mina definitioner tidigare har vi att (ignorerar index på PP här)

ddxlnPx=Q0x\displaystyle \frac{d}{dx}\ln P\left(x\right) = Q_0\left(x\right) och dmdxmlnPx=Qmx\displaystyle \frac{d^m}{dx^m} \ln P\left(x\right) = Q_m\left(x\right)

Faà di Bruno's formel (speciallfall längre ned på wikipediasidan) ger

dmdxmPx=dmdxmelnPx=PxBmQ0x,Q1x,,Qm-1x\displaystyle \frac{d^m}{dx^m}P\left(x\right) = \frac{d^m}{dx^m}e^{\ln P\left(x\right)} = P\left(x\right)B_m\left(Q_0\left(x\right), Q_1\left(x\right),\dots, Q_{m-1}\left(x\right)\right)

Där Bm(x1,,xm)B_m(x_1, \dots, x_m) är det mm:te "complete exponential bell polynomial". Detta ger oss då att den "slutna" formeln för RmR_m är BmQ0x,Q1x,,Qm-1xB_m\left(Q_0\left(x\right), Q_1\left(x\right),\dots, Q_{m-1}\left(x\right)\right).

Coolt, men kanske inte särskilt lätt att jobba med. Jag ska leka lite med definitionerna av dessa polynom. Det kanske kan, på något under, förenklas till något trevligt. Det jobbiga är däremot att Qi(a)Q_i(a) där a=0,1a=0, 1 ej är definierade. Man måste egentligen ta gränsvärden på grund av att de är rötter till PP

Ska testa dina förslag efter jag lekt runt med detta!


Tillägg: 25 dec 2025 14:47

Rättning:

dmdxmlnP(x)=Qm-1(x)\displaystyle \frac{d^m}{dx^m}\ln P(x) = Q_{m-1}(x)


Tillägg: 26 dec 2025 17:26

Angående att hålla på integralformeln direkt och skriva om summan av derivatorna tycker jag det känns lite svårt. Varje derivata är har ju π\pi till någon exponent från nämnaren och sedan även zetafunktionen som multipliceras ihop. Om man byter summeringsordning kanske det är lättare, då får man inga sådana problem i alla fall, bara rationella tal. Däremot summerar vi då över indexet kk och polynomen är inte suuper länkade med varandra längre. 

AlexMu 1028
Postad: 26 dec 2025 17:18 Redigerad: 26 dec 2025 19:22

Tillbaka till att utveckla polynomet till a0++a2nx2na_0 + \dots + a_{2n}x^{2n}:

Vi har att

Px-n+kk=x-n+kkn-1j=0nj-x\displaystyle P\left(\frac{x-n+k}{k}\right)=\left(\frac{x-n+k}{k}\right)^{n-1}\prod_{j=0}^{n}\left(j-x\right)

Nu är produkten definitionen av pochhammer symbolen, (x)n=x(x+1)(x+n-1)(x)_n = x(x+1)\cdots(x+n-1) detta är trevligt pga att de har en "känd utveckling". 

Pochhamer definitionen ger

Px-n+kk=x-n+kkn-1-xn+1\displaystyle P\left(\frac{x-n+k}k\right) = \left(\frac{x-n+k}{k}\right)^{n-1}\left(-x\right)_{n+1}
Vilket ger att 

wikipediasidan för nk{n\brack k}, stirlingtalen av första slaget får vi (notation för pochhammer är annorlunda där än jag är van vid)

Px-n+kk=x-n+kkn-1j=0n+1-1jn+1jxj\displaystyle P\left(\frac{x-n+k}k\right) = \left(\frac{x-n+k}{k}\right)^{n-1}\sum_{j=0}^{n+1}\left(-1\right)^{j}{{n+1}\brack j}x^{j}

Trevligt, då har vi nästan en parentesutveckling för PP. Jag är nöjd om man kan skriva derivatorna uttryckta i stirlingtal. Det jobbiga är nu faktorn framför summan. 


Tillägg: 26 dec 2025 19:18

Detta fungerade rätt bra!

Man får att, för aa positivt heltal,

P(a)0=a!n-kn-1-akn-a-1j=0n+1m=0n-1-1j+n-1-mn+1jn-1mj+man-kj\displaystyle P^{(a)}\left(0\right) = \frac{a!\left(n-k\right)^{n-1-a}}{k^{n-a-1}}\sum_{j=0}^{n+1}\sum_{m=0}^{n-1}\left(-1\right)^{j+n-1-m}{n+1\brack j}\binom{n-1}{m}\binom{j+m}{a}\left(n-k\right)^{j}

P(a)1=a!n-kn-1kn-a-1naj=0n+1m=0n-1-1j+n-1-mn+1jn-1mj+manj+mn-k-m\displaystyle P^{(a)}\left(1\right) = \frac{a!\left(n-k\right)^{n-1}}{k^{n-a-1}n^{a}}\sum_{j=0}^{n+1}\sum_{m=0}^{n-1}\left(-1\right)^{j+n-1-m}{n+1\brack j}\binom{n-1}{m}\binom{j+m}a n^{j+m}\left(n-k\right)^{-m}

Jag hoppas det finns en del förenkling att göra, kanske när man lagt in detta i integralformeln. Vi får se! Tack.

Svara
Close