Estoy leyendo un artículo sobre un modelo de osciladores (modelo Kuramoto), y me encontré con la siguiente ecuación de Fokker-Planck.
r(t) se define así: reiω=∫∫eiθρ(t,θ.ω)g(ω)dωdθ rho es una densidad (que tenía la restricción de integrarse a la unidad) y ahora está perturbada de esta manera: ρ(ω,t,θ)=ϵη(ω,t,θ)+1/2π and ϵ<<1
ahora la ecuación de Focker Planck dice : ε∂η∂t=εD∂2η∂θ2−∂∂θ[(12π+ϵη)v(ω,t,θ)]
...luego dice "Consideremos esta ecuación en el orden más bajo de η . Para encontrar el valor de O( ϵ ) en el término entre corchetes, observamos que r(t) es O( ϵ ). Más concretamente, encontramos : r(t)=ϵr1(t)+O(ϵ2)where:r1eiψ=∫∫eiθη(t,θ.ω)g(ω)dωdθ
Así que Ok, reescriben la primera ecuación enchufando la densidad modificada (eqn.2). Esto es lo que obtengo:
r(t)eiψ=∫∫eiθρ(θ,t,ω)g(ω)dωdθ=∫∫eiθ[12π+ϵη(θ,t,ω)]g(ω)dωdθ=ϵ∫∫eiθη(θ,t,ω)g(ω)dωdθ+12π∫∫eiθg(ω)dωdθ
⇔r(t)=ϵe−iψ∫∫eiθη(θ,t,ω)g(ω)dωdθ+e−iψ2π∫+∞−∞∫2π0g(ω)eiθdωdθ=ϵr1(t)+e−iψ2π∫+∞−∞g(ω)dω∫2π0eiθdθ ...nada es de orden épsilon al cuadrado aquí? Además no voy a apostar demasiado pero me parece que la segunda integral se desvanece por la [eiθ]2π0 ...
Gracias.