Leicht oszillierende Reihen mit hoher Präzision berechnen?

13

Angenommen, ich habe die folgende interessante Funktion: Es hat einige unangenehme Eigenschaften, wie die Ableitung, die bei rationalen Vielfachen vonπnicht stetig ist. Ich vermute, dass ein geschlossenes Formular nicht existiert.

f(x)=k1coskxk2(2-coskx).
π

Ich kann es berechnen, indem ich Teilsummen berechne und Richardson-Extrapolation verwende, aber das Problem ist, dass es zu langsam ist, die Funktion auf eine gute Anzahl von Dezimalstellen zu berechnen (100 wäre zum Beispiel schön).

Gibt es eine Methode, die diese Funktion besser beherrscht?

Hier ist eine Darstellung von mit einigen Artefakten:f(πx)

Funktionsableitung $ f '(\ pi x) $

Kirill
quelle
1
Vielleicht kann man die Tatsache , dass verwenden , wobei T k ( x ) ein Tschebyscheff - Polynom ist. Dann beginnt die Summation wie eine Reihe rationaler Polynome auszusehen. Wenn Sie dann die Reihe in ein rationales Polynom auf Chebyshev-Basis umwandeln können, ist dies eine sehr effiziente Methode, um sie zusammenzufassen. Wenn Sie mit Chebyshev-Polynomen und -Basis nicht vertraut sind, haben die numerischen Rezepte in C eine gute Grundierung.cÖs(kx)=Tk(x)Tk(x) www2.maths.ox.ac.uk/chebfun/ATAP/ATAPfirst6chapters.pdf
Jay Lemmon
1
äh, das sollte sagen cos(kx)=Tk(cos(x))
Jay Lemmon
@ JayLemmon Vielen Dank für diesen Link. Ich werde nachsehen, ob es hilft.
Kirill
Ich komme etwas spät zu dieser Party, aber haben Sie versucht, Padé-Approximanten zu verwenden, dh den Algorithmus anstelle der Richardson-Extrapolation? ε
Pedro
In Analogie zu hochoszillierenden Integralen glaube ich, dass Sie ohne Kenntnisse über die Trennung zwischen oszillierenden und nichtoszillierenden Teilen keine gute Arbeit leisten können. Wenn Sie eine solche Trennung haben, gibt Ihnen die Fourier-Reihenantwort eine einfache exponentielle Konvergenz.
Geoffrey Irving

Antworten:

7

Wenn analytische Techniken nicht zulässig sind, aber die periodische Struktur bekannt ist, ist hier ein Ansatz. Sei ist periodisch mit der Periode2π, so dass g(x)=jwjeijx wobei wj=1 ist

G(x)=cosx2-cosx
2π
G(x)=jwjeichjx
Somit ist f ( x )
wj=12π02πG(x)e-ichjxdx
Sie können entweder die Integrale approximierenwjdirekt oder ein Bündel berechnenf(x)Werte und eine DFT verwenden. In beiden Fällen können Sie möglicherweise eine Richardson-Extrapolation auf das Ergebnis anwenden. Da in Ihrem Fallg(x)innerhalb einer Nachbarschaft von analytisch ist
f(x)=k1G(kx)kp=k11kpjwjeichjkx=jwjk1(eichjx)kkp=jwjLip(eichjx)
wjf(x)G(x), konvergiert die endgültige Reihe auch ohne Richardson exponentiell.R
Geoffrey Irving
quelle
G(x)=cos(x)/(2-cos(x))
3

x=2πein/bein,b

f(x)=k1coskxk2(2-coskx)=k=1bcoskx2-coskxn01(k+bn)2=k=1bcoskx2-coskxψ1(k/b)b2
ψ1(z)Werte und Ableitungen für die Serie
Geoffrey Irving
quelle
Vielen Dank. Das Problem ist, dass ich diese spezielle Funktion als Modell für eine andere, kompliziertere Funktion ausgewählt habe, die ich eigentlich bewerten wollte. Sie hat ähnliche Eigenschaften, ist aber nicht die gleiche. Mir ist das geschlossene Formular aus dieser Frage zu MSE bekannt . Ich meinte dies als eine Frage zur numerischen Summierung einer unendlichen Reihe ohne geschlossene Form.
Kirill
Vielleicht ist meine andere Antwort dann besser?
Geoffrey Irving