Çok ölçekli analiz - Multiple-scale analysis

İçinde matematik ve fizik, çok ölçekli analiz (ayrıca çoklu ölçek yöntemi) aynı şekilde geçerli oluşturmak için kullanılan teknikleri içerir yaklaşımlar çözümlerine tedirginlik sorunları, hem küçük hem de büyük değerler için bağımsız değişkenler. Bu, bağımsız bir değişken için hızlı ölçekli ve yavaş ölçekli değişkenler ekleyerek ve ardından bu değişkenleri hızlı ve yavaş, bağımsızlarmış gibi ele alarak yapılır. Bundan sonraki pertürbasyon probleminin çözüm sürecinde, ortaya çıkan ek özgürlük - yeni bağımsız değişkenler tarafından getirilen - kaldırmak için kullanılır (istenmeyen) laik terimler. İkincisi, yaklaşık çözüme kısıtlamalar koyar. çözülebilirlik koşulları.

1980'lerden itibaren matematik araştırması, koordinat dönüşümlerinin ve değişmez manifoldların çok ölçekli modelleme için daha sağlam bir destek sağladığını önermektedir (örneğin, bkz. merkez manifold ve yavaş manifold ).

Örnek: sönümsüz Duffing denklemi

Diferansiyel denklem ve enerji tasarrufu

Çok ölçekli analiz yöntemine bir örnek olarak, sönümlenmemiş ve zorlanmamış Duffing denklemi:^[1]

${ displaystyle { frac {d ^ {2} y} {dt ^ {2}}} + y + varepsilon y ^ {3} = 0,}$   ${ displaystyle y (0) = 1, qquad { frac {dy} {dt}} (0) = 0,}$

ikinci dereceden olan adi diferansiyel denklem tanımlayan doğrusal olmayan osilatör. Bir çözüm y(t) (pozitif) doğrusal olmayan parametrenin küçük değerleri için aranır 0 <ε ≪ 1. Sönümsüz Duffing denkleminin bir Hamilton sistemi:

${ displaystyle { frac {dp} {dt}} = - { frac { kısmi H} { kısmi q}}, qquad { frac {dq} {dt}} = + { frac { kısmi H} { kısmi p}}, quad { text {with}} quad H = { tfrac {1} {2}} p ^ {2} + { tfrac {1} {2}} q ^ {2} + { tfrac {1} {4}} varepsilon q ^ {4},}$

ile q = y(t) ve p = dy/dt. Sonuç olarak, Hamiltoniyen H(p, q) korunan bir miktardır, sabittir, eşittir H = ½ + ¼ ε verilen için başlangıç ​​koşulları. Bu, her ikisinin de y ve dy/dt sınırlandırılmalıdır:

${ displaystyle sol | y (t) sağ | leq { sqrt {1 + { tfrac {1} {2}} varepsilon}} quad { text {ve}} quad sol | { frac {dy} {dt}} right | leq { sqrt {1 + { tfrac {1} {2}} varepsilon}} qquad { text {tümü için}} t.}$ 

Basit pertürbasyon serisi çözümü

Düzenli tedirginlik serisi yaklaşımı soruna sonucu verir:

${ displaystyle y (t) = cos (t) + varepsilon sol [{ tfrac {1} {32}} cos (3t) - { tfrac {1} {32}} cos (t) - underbrace {{ tfrac {3} {8}} , t , sin (t)} _ { text {seküler}} sağ] + { mathcal {O}} ( varepsilon ^ {2 }).}$

Köşeli parantezler arasındaki son terim dünyevidir: sınırsız büyür |t|. Özellikle, ${ displaystyle t = O ( epsilon ^ {- 1})}$ bu terim Ö(1) ve öncü sıra terimiyle aynı büyüklük sırasına sahiptir. Terimler düzensiz hale geldiğinden, seri artık çözümün asimptotik bir genişlemesi değildir.

Çoklu ölçek yöntemi

Ötesinde geçerli bir çözüm inşa etmek ${ displaystyle t = O ( epsilon ^ {- 1})}$yöntemi çok ölçekli analiz kullanıldı. Yavaş ölçeği tanıtın t₁:

${ displaystyle t_ {1} = varepsilon t ,}$

ve çözümü varsayalım y(t) her ikisine de bağlı bir tedirginlik serisi çözümdür. t ve t₁, şu şekilde ele alınır:

${ displaystyle y (t) = Y_ {0} (t, t_ {1}) + varepsilon Y_ {1} (t, t_ {1}) + cdots.}$

Yani:

${ displaystyle { begin {align} { frac {dy} {dt}} & = left ({ frac { kısmi Y_ {0}} { kısmi t}} + { frac {dt_ {1} } {dt}} { frac { kısmi Y_ {0}} { kısmi t_ {1}}} sağ) + varepsilon left ({ frac { kısmi Y_ {1}} { kısmi t} } + { frac {dt_ {1}} {dt}} { frac { kısmi Y_ {1}} { kısmi t_ {1}}} sağ) + cdots & = { frac { kısmi Y_ {0}} { kısmi t}} + varepsilon left ({ frac { kısmi Y_ {0}} { kısmi t_ {1}}} + { frac { kısmi Y_ {1}} { kısmi t}} sağ) + { mathcal {O}} ( varepsilon ^ {2}), end {hizalı}}}$

kullanma dt₁/dt = ε. Benzer şekilde:

${ displaystyle { frac {d ^ {2} y} {dt ^ {2}}} = { frac { kısmi ^ {2} Y_ {0}} { kısmi t ^ {2}}} + varepsilon left (2 { frac { kısmi ^ {2} Y_ {0}} { kısmi t , kısmi t_ {1}}} + { frac { kısmi ^ {2} Y_ {1}} { kısmi t ^ {2}}} sağ) + { mathcal {O}} ( varepsilon ^ {2}).}$

Daha sonra, Duffing denklemi için çok ölçekli tedirginlik serisinin sıfırıncı ve birinci dereceden problemleri şu hale gelir:

${ displaystyle { başla {hizalı} { frac { kısmi ^ {2} Y_ {0}} { kısmi t ^ {2}}} + Y_ {0} & = 0, { frac { kısmi ^ {2} Y_ {1}} { kısmi t ^ {2}}} + Y_ {1} & = - Y_ {0} ^ {3} -2 , { frac { kısmi ^ {2} Y_ {0}} { kısmi t , kısmi t_ {1}}}. Uç {hizalı}}}$

Çözüm

Sıfırıncı sıra probleminin genel çözümü vardır:

${ displaystyle Y_ {0} (t, t_ {1}) = A (t_ {1}) , e ^ {+ it} + A ^ { ast} (t_ {1}) , e ^ {- o},}$

ile Bir(t₁) bir karmaşık değerli genlik sıfırıncı sıradaki çözüme Y₀(t, t₁) ve ben² = −1. Şimdi, birinci dereceden problemde zorlama sağ taraf diferansiyel denklemin

${ displaystyle sol [-3 , A ^ {2} , A ^ { ast} -2 , i , { frac {dA} {dt_ {1}}} sağ] , e ^ {+ it} -A ^ {3} , e ^ {+ 3it} + cc}$

nerede c.c. gösterir karmaşık eşlenik önceki şartların. Oluşumu laik terimler - henüz bilinmeyen - genlik empoze edilerek önlenebilir Bir(t₁) çözülebilirlik koşulu

${ displaystyle -3 , A ^ {2} , A ^ { ast} -2 , i , { frac {dA} {dt_ {1}}} = 0.}$

Çözülebilirlik koşulunun çözümü, aynı zamanda başlangıç ​​koşullarını da karşılar y(0) = 1 ve dy/dt(0) = 0, şudur:

${ displaystyle A = { tfrac {1} {2}} , exp left ({ tfrac {3} {8}} , i , t_ {1} sağ).}$

Sonuç olarak, çoklu ölçek analizinin yaklaşık çözümü şu şekildedir:

${ displaystyle y (t) = çünkü sol [ sol (1 + { tfrac {3} {8}} , varepsilon sağ) t sağ] + { mathcal {O}} ( varepsilon ),}$

kullanma t₁ = εt ve için geçerlidir εt = O (1). Bu doğrusal olmayan ile aynı fikirde Sıklık kullanılarak bulunan değişiklikler Lindstedt-Poincaré yöntemi.

Bu yeni çözüm şu tarihe kadar geçerlidir: ${ displaystyle t = O ( epsilon ^ {- 2})}$. Birden fazla ölçek yöntemini kullanan daha yüksek dereceli çözümler, ek yavaş ölçeklerin kullanılmasını gerektirir, yani: t₂ = ε² t, t₃ = ε³ t, vb. Bununla birlikte, bu, tedirginlik serisi çözümünde dikkatli bir muamele gerektiren olası belirsizliklere yol açar (bkz. Kevorkian ve Cole 1996; Bender ve Orszag 1999 ).^[2]

Genlik / faz değişkenlerine koordinat dönüşümü

Alternatif olarak, modern ses yaklaşımları bu tür modelleri koordinat dönüşümlerini kullanarak türetmektedir.^[3] aşağıda açıklandığı gibi.

Bir çözüm ${ displaystyle y yaklaşık r cos theta}$ yeni koordinatlarda aranıyor ${ displaystyle (r, theta)}$ genlik nerede ${ displaystyle r (t)}$ yavaşça değişir ve aşama ${ displaystyle theta (t)}$ neredeyse sabit bir oranda değişir, yani ${ displaystyle d theta / dt yaklaşık 1.}$ Basit cebir koordinat dönüşümünü bulur^{[kaynak belirtilmeli ]}

${ displaystyle y = r cos theta + { frac {1} {32}} varepsilon r ^ {3} cos 3 theta + { frac {1} {1024}} varepsilon ^ {2} r ^ {5} (- 21 cos 3 theta + cos 5 theta) + { mathcal {O}} ( varepsilon ^ {3})}$

Duffing'in denklemini yarıçapın sabit olduğu çifte dönüştürür ${ displaystyle dr / dt = 0}$ ve aşama şuna göre gelişir:

${ displaystyle { frac {d theta} {dt}} = 1 + { frac {3} {8}} varepsilon r ^ {2} - { frac {15} {256}} varepsilon ^ { 2} r ^ {4} + { mathcal {O}} ( varepsilon ^ {3}).}$

Yani, Duffing'in salınımları sabit büyüklüktedir ${ displaystyle r}$ ama farklı frekanslara sahip ${ displaystyle d theta / dt}$ genliğe bağlı olarak.^[4]

Daha zor örnekler, karmaşık üstelleri içeren zamana bağlı bir koordinat dönüşümü kullanılarak daha iyi işlenir (önceki çoklu zaman ölçeği yaklaşımında da belirtildiği gibi). Bir web servisi, çok çeşitli örnekler için analizi gerçekleştirecektir.^[5]

Ayrıca bakınız

• Eşleştirilmiş asimptotik genişletme yöntemi
• WKB yaklaşımı

Notlar

Referanslar

Dış bağlantılar

• Carson C. Chow (ed.). "Çoklu ölçekli analiz". Scholarpedia.