«numerical-analysis» etiketlenmiş sorular

Sürekli problemlerin yaklaşık ayrık çözümlerini hesaplamak için algoritmaların oluşturulması ve analizi. Kanonik bir örnek, türevlerin fark katsayıları ile yakınlaştırılmasıdır.

17
Python için yüksek kaliteli bir doğrusal olmayan programlama çözücü var mı?
Çözülmesi gereken birkaç dışbükey olmayan küresel optimizasyon problemim var. Şu anda , oldukça etkili olan MATLAB's Optimizasyon Araç Kutusu'nu (özellikle, fmincon()algoritma = ile 'sqp') kullanıyorum . Ancak, kodumun çoğu Python'da ve optimizasyonu Python'da da yapmak isterim. Rekabet edebilecek Python bağlarına sahip bir NLP çözücü var mı fmincon()? O olmalı Doğrusal …

3
Paralel ODE yöntemlerinde sanatın durumu nedir?
Şu anda ODE entegrasyonu için paralel yöntemler arıyorum. Geniş bir yaklaşım yelpazesini tanımlayan pek çok yeni ve eski literatür var, ancak konuyu genel olarak anlatan herhangi bir son anket veya genel bakış makalesi bulamadım. Burrage [1] adlı kitabın kitabı var, ancak neredeyse 20 yaşında ve bu nedenle parareal algoritması gibi …


3


1
Sayısal kodda hata ayıklama, bu salınım hatasının kaynağı ne olabilir?
Sessizce çok fazla fikir edinilebilir, daha önce böyle bir şey görüp görmediğini merak ediyordum. Grafik, adveksiyon-difüzyon denklemi için başlangıç ​​koşulunu (yeşil), daha sonra iterasyon 200 (mavi) ve sonra tekrar iterasyon 400 (kırmızı) 'daki çözeltiyi gösterir. Adveksiyon-difüzyon denkleminin çözeltisi birkaç tekrardan sonra patlar. Péclet numarası ve CFL koşulu karşılanır, , bu …

1
Sonlu hacimler yöntemi kullanılırken sınır koşulları nasıl uygulanmalıdır?
Benim aşağıdaki önceki soruya ben, bu düzgün olmayan sonlu hacim örgü için sınır koşulları kullanılarak çalışıyorum Etki alanının lhs ( bir Robin türü sınır koşulu uygulamak istiyorum , böylece,x=xL)x=xL)x=x_L) σL=(dux+au)∣∣∣x=xLσL=(dux+au)|x=xL \sigma_L = \left( d u_x + a u \right) \bigg|_{x=x_L} burada sınır değeridir; a , d sırasıyla sınır, kabul ve …

3
Bilimsel Programlama Yarışmaları
Sınırlı bir zaman dilimi içinde kendi kodunuz ve problem çözme becerilerinizle zorlu algoritmik problemleri çözdüğünüz "Programlama Yarışmaları" na düzenli olarak katılıyorum. Bunların neye benzeyebileceğine ilişkin referans örnekleri için Google Code Jam veya ACM-ICPC gibi yarışmalar arayın. (Programlama yarışmalarının ne olduğunu biliyorsanız, aşağıdaki paragrafı atlayabilirsiniz) Bu yarışmalarda, bireysel olarak veya takımlar …

1
Beşinci dereceden Runge-Kutta yönteminin istikrar bölgesi hakkında şaşırtıcı açıklama
Gazetede kafa karıştırıcı bir ifadeyle karşılaştım PJ van der Houwen, Kısmi diferansiyel denklemler için Runge-Kutta yöntemlerinin geliştirilmesi, Appl. Num. Matematik. 20: 261, 1996 264. Sayfadaki 8ff satırlarında, van der Houwen şöyle yazıyor: "Taylor polinomları için bu hayali kararlılık aralığı için boş olduğunu ima p = 1 , 2 , 5 …

3
Matris karekök tersinin verimli hesaplanması
İstatistiklerde sık karşılaşılan bir problem, simetrik pozitif belirli bir matrisin ters kökünün hesaplanmasıdır. Bunu hesaplamanın en etkili yolu ne olabilir? Burada (henüz okumadım) bazı literatür ve rastlantısal olarak burada çoğaltacağım bazı rastlantısal R kodlarıyla karşılaştım. # function to compute the inverse square root of a matrix fnMatSqrtInverse = function(mA) { …

2
Doğrusal olmayan PDE'leri Newton-Raphson yinelemesi kullanmadan çözmek mümkün müdür?
Bazı sonuçları anlamaya çalışıyorum ve doğrusal olmayan problemlerle ilgili bazı genel yorumları takdir ediyorum. Fisher denklemi (doğrusal olmayan reaksiyon difüzyon PDE), ut= dux x+ βu ( 1 - u ) = F( u )ut=duxx+βu(1-u)=F(u) u_t = du_{xx} + \beta u (1 - u) = F(u) takdirine bağlı olarak, u'j= L …


1
Python işlevinde felaket iptali nasıl önlenir?
Bir işlevi sayısal olarak uygulamakta sorun yaşıyorum. Büyük giriş değerlerinde sonucun çok büyük bir sayının çok küçük bir sayının çarpımı olduğu gerçeğinden muzdariptir. Felaket iptalinin doğru terim olup olmadığından emin değilim, lütfen öyleyse düzeltin. Bir şeyin ters gittiğine dair kanıtlar: 6'nın daha büyük girdileri için salınımları ve 0.0 değerini atamayı …

2
Sonlu fark yöntemleri için von neumann kararlılık analizine alternatifler
Birleştirilmiş tek boyutlu poroelastisite denklemlerini ( biot'un modeli) çözmeye çalışıyorum : ∂- ( λ + 2 μ ) ∂2u∂x2+ ∂p∂x= 0−(λ+2μ)∂2u∂x2+∂p∂x=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} = 0 etki alanındaΩ=(0,1)ve sınır koşulları: ∂∂t[ γp + ∂u∂x] - κη[ ∂2p∂x2] = q( x , t )∂∂t[γp+∂u∂x]−κη[∂2p∂x2]=q(x,t)\frac{\partial}{\partial t} \left[ …

2
birçok değişkende sayısal entegrasyon
Let ve m ( → x ) : [ 0 , 1 ] N → Cı bu değişkenlerin bir fonksiyonu olsun.x⃗ =(x1,x2,…,xn)∈[0,1]nx→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nf(x⃗ ):[0,1]n→Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} Bu yinelenen integral için özyinelemeli bir plan var mı? ∫[0,1]n∏dxif(x⃗ )∫[0,1]n∏dxif(x→)\int_{[0,1]^n} \prod dx_i \;f(\vec{x}) Eğer ve [ …

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.