Computando séries ligeiramente oscilatórias com alta precisão?

13

Suponha que eu tenha a seguinte função interessante: Tem algumas propriedades desagradáveis, como sua derivada não sendo contínua em múltiplos racionais de . Eu suspeito que um formulário fechado não existe.

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

Posso computá-lo calculando somas parciais e usando a extrapolação de Richardson, mas o problema é que é muito lento para calcular a função com um bom número de dígitos decimais (100 seria bom, por exemplo).

Existe um método que possa lidar melhor com essa função?

Aqui está um gráfico de com alguns artefatos:f(πx)

Derivada da função, $ f '(\ pi x) $

Kirill
fonte
1
Talvez você possa usar o fato de que , onde é um polinômio de Chebyshev. Então a soma começa a parecer uma série de polinômios racionais. Então, se você puder transformar a série em um polinômio racional em uma base Chebyshev, isso permitiria uma maneira muito eficiente de resumir. Se você não conhece os polinômios e a base de Chebyshev, as Receitas numéricas em C têm uma boa cartilha, assim como: www2.maths.ox.ac.uk/chebfun/ATAP/ATAPfirst6chapters.pdfT k ( x )cos(kx)=Tk(x)Tk(x)
Jay Lemmon
1
er, que deveria dizer cos(kx)=Tk(cos(x))
Jay Lemmon
@JayLemmon Obrigado por esse link. Vou dar uma olhada e ver se isso ajuda.
Kirill
Estou entrando nessa festa um pouco tarde, mas você tentou usar as aproximações de Padé, ou seja, o algoritmo Algoritmo em vez da extrapolação de Richardson? ε
Pedro
Por analogia com o caso de integrais altamente oscilatórias, acho que você não será capaz de fazer um bom trabalho sem ter conhecimento da separação entre partes oscilatórias e não-oscilatórias. Se você tem essa separação, a resposta da série Fourier fornece uma convergência exponencial fácil.
Geoffrey Irving

Respostas:

7

Se as técnicas analíticas não forem permitidas, mas a estrutura periódica for conhecida, aqui está uma abordagem. Seja ser periódica com período de2π, de modo que g(x)=Σjwjeijx onde wj=1

g(x)=cosx2cosx
2π
g(x)=jWjeEujx
Assim, f ( x )
Wj=12π0 02πg(x)e-Eujxdx
Você pode aproximar as integraiswjdiretamente ou calcular vários valores def(x)e usar um DFT. Em qualquer um dos casos, é possível aplicar a extrapolação de Richardson ao resultado. Como no seu caso,g(x)é analítico em uma vizinhança deR, a série final converge exponencialmente mesmo sem Richardson.
f(x)=k1g(kx)kp=k11kpjWjeEujkx=jWjk1(eEujx)kkp=jWjLip(eEujx)
Wjf(x)g(x)R
Geoffrey Irving
fonte
g(x)=porque(x)/(2-porque(x))
3

x=2πuma/buma,b

f(x)=k1porquekxk2(2-porquekx)=k=1bporquekx2-porquekxn0 01(k+bn)2=k=1bporquekx2-porquekxψ1(k/b)b2
ψ1(z)Valores e derivativos para as séries
Geoffrey Irving
fonte
Obrigado. O problema é que eu selecionei essa função específica como modelo para outra função mais complicada que eu realmente queria avaliar, com recursos semelhantes, mas na verdade não os mesmos. Estou ciente do formulário fechado desta pergunta no MSE . Eu quis dizer isso como uma pergunta sobre somar uma série infinita numericamente sem forma fechada.
Kirill
Talvez minha outra resposta seja melhor então?
Geoffrey Irving