Was sind diese Schwingungen?

Ich habe eine numerisch definierte Funktion die zwischen einem Gaußschen und einem Lorentzschen liegt. Es zerfällt viel langsamer als ein Gaußscher, aber immer noch schneller als eine einfache inverse Kraft. $g(x)$

Ich muss seine Fourier-Transformation für großes berechnen . Da Funktionsaufrufe von rechenintensiv sind, definiere ich eine Interpolation von - nenne es - für einen großen Bereich von , , und benutze das für mein Integral. $f(t)\equiv \mathcal{F}[g(x)](t)$ $t$ $g(x)$ $g(x)$ $g_{\text{int}}(x)$ $x$ $-40<x<40$

f (t) = \int_{- - \infty}^{\infty} \cos (t x) G (x) d x \underset{\approx}{⟶} \int_{- - L.}^{L.} \cos (t x) G_{int} (x) d x

$f(t)=\int_{-\infty}^{\infty}\cos (tx)g(x)\,dx\,\,\underset{\approx}{\longrightarrow}\,\,\,\int_{-L}^{L}\cos(tx)g_{\text{int}}(x)\,dx$

Wenn ich jedoch eine Annäherung an die Fourier-Transformation berechne, erhalte ich einige merkwürdige Schwingungen, die ich zunächst nicht erwarten würde.

Wie ich im obigen Bild angegeben habe, haben die Schwingungen eine "Periode" von ungefähr 15,7. Meine erste Vermutung wäre, dass dies ein Artefakt der alternierenden Natur der Aufhebung des Integrals sein könnte, aber das würde die beobachtete "Periode" von 15,7 nicht erklären.

{T.}_{vermuten} = \frac{2 π}{L.} \approx 0,157 \dots

$T_{\text{guess}}=\frac{2\pi}{L}\approx 0.157\ldots$

Das ist genau ein Faktor von 100, der sich von dem unterscheidet, was ich beobachte (ja, ich habe überprüft, ob ich meine Integrale und horizontalen Achsen richtig definiert habe). Wie konnte das sein?

Edit # 1: Interpolationsdetails

Ich interpoliere mit Mathematicas integriertem InterpolationElement, das zwischen aufeinanderfolgenden Punkten mit einer kubischen Kurve interpoliert (also wird an jedem Punkt bis zur Ableitung 2 definiert). Ich interpoliere spezifisch die Funktion über den Bereich in Schritten von . $^{\text{nd}}$ $g(x)$ $-40<x<40$ $dx=40/100=0.4$

In der Tat, jetzt, wo ich das schreibe, ist mir klar, dass es sehr gut ein Artefakt meiner endlichen Stichprobe sein könnte, weil:

{T.}_{Vermutung # 2} = \frac{2 π}{d x} = \frac{2 π}{0,4} = 15.7 \dots

$T_{\text{guess #2}}=\frac{2\pi}{dx}=\frac{2\pi}{0.4}=15.7\ldots$

Ich würde mich über weitere Hilfe freuen, insbesondere über einen guten Weg, um dieses Problem zu lösen.

Edit # 2: Die Funktion $g(x)$

h[x_?NumericQ, En_?NumericQ, pz_?NumericQ] := 
 1./(En^2 + pz^2 + 0.24^2)*
  NIntegrate[((Sqrt[
      0.316/(1. + 
         1.2*((k4 + 0.5*En)^2 + kp + (x*pz)^2))^1.*0.316/(1. + 
         1.2*((k4 - 0.5*En)^2 + kp + ((1. - x)*pz)^2))^1.])*((1. - 
         x)*0.316/(1. + 1.2*((k4 + 0.5*En)^2 + kp + (x*pz)^2))^1. + 
      x*0.316/(1. + 
         1.2*((k4 - 0.5*En)^2 + kp + ((1. - x)*pz)^2))^1.))/(((k4 + 
        0.5*En)^2 + 
      kp + (x*pz)^2 + (0.316/(1. + 
         1.2*((k4 + 0.5*En)^2 + kp + (x*pz)^2))^1.)^2)*((k4 - 
        0.5*En)^2 + 
      kp + ((1. - x)*
        pz)^2 + (0.316/(1. + 
         1.2*((k4 - 0.5*En)^2 + 
            kp + ((1. - x)*
              pz)^2))^1.)^2)), {k4, -\[Infinity], \[Infinity]}, {kp, 
    0, \[Infinity]}, Method -> "LocalAdaptive", 
   MaxRecursion -> 
    100]; (*LocalAdaptive seems to work slightly faster *)

g[x_]:=h[0.5,x,2.]; (*this is the function*)

interpolation integration mathematica oscillations Arturo don Juan
quelle

Welche Art von Interpolation verwenden Sie? Es ist eigentlich nicht die einfachste Funktion, auf einer so großen Domäne genau zu interpolieren.

Anton Menshov

@AntonMenshov Siehe meine letzte Bearbeitung.

Arturo Don Juan

Ich möchte darauf hinweisen, dass jede numerische Integrationsmethode im Wesentlichen einen Interpolanten Ihrer Funktion erstellt und diesen genau integriert. Wenn Sie also die Funktion zuerst selbst interpolieren, ersetzen Sie Integratetatsächlich die adaptive Konstruktion von Interpolanten (und all den Gedanken, der in ihre Funktion eingeflossen ist) Design) mit Ihrem eigenen, was wahrscheinlich eine schlechte Idee ist.

Kirill

g

$g$

\log g (x)

$\log g(x)$

[- L, L]

$[-L,L]$ Module[{n=32,L=10},Exp[InterpolatingPolynomial[Table[{x,Log[g[x]]},{x,-L Cos[(π N@Range[1,2n-1,2])/(2n)]}],x]]]

Ich denke auch, dass es hilfreich wäre, den Integranden ( gist.github.com/ikirill/931dbf6fb3581f7d66353b74c4743cd3 ) zu kompilieren , unabhängig von all den schwierigeren Problemen, die es für mich um den Faktor ~ 8 beschleunigt.

Kirill

Was sind diese Schwingungen?

Edit # 1: Interpolationsdetails

Edit # 2: Die Funktion G( x )G(x)g(x)

Antworten:

Edit # 2: Die Funktion $g(x)$