Pertanyaan yang diberi tag «numerical-analysis»

Konstruksi dan analisis algoritma untuk menghitung solusi diskrit perkiraan masalah kontinu. Contoh kanonik adalah perkiraan derivatif melalui quotients perbedaan.




1
Metode numerik untuk membalikkan transformasi integral?
Saya mencoba membalikkan transformasi integral berikut secara numerik: F(y) = ∫∞0yexp[ - 12( y2+ x2) ] Saya0( X y) f( x )d xF(y)=∫0∞yexp⁡[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Jadi untuk saya perlu memperkirakan f ( x ) di mana:F( y)F(y)F(y)f( x )f(x)f(x) dan F ( y ) adalah nyata …

3
Bagaimana seharusnya koefisien non-konstan diperlakukan dengan skema upwind urutan pertama volume terbatas?
Dimulai dengan persamaan advection dalam bentuk konservasi. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x di mana adalah kecepatan yang bergantung pada ruang, dan adalah konsentrasi spesies yang dilestarikan.ua(x)a(x)a(x)uuu Mendiskritkan fluks (di mana fluks , didefinisikan pada tepi sel di antara titik-titik mesh) memberi, u t = 1f=a(x)uf=a(x)uf=a(x)uut=1h(fj−12−fj+12)ut=1h(fj−12−fj+12) u_t = \frac{1}{h}\left( f_{j-{\frac{1}{2}}} - f_{j+{\frac{1}{2}}} …

1
Memecahkan sistem persamaan yang sulit secara numerik
Saya memiliki sistem persamaan -linear yang ingin saya pecahkan secara numerik:nnn f = ( f 1 , … , f n )f( x ) = af(x)=a\mathbf{f}(\mathbf{x})=\mathbf{a} f= ( f1, ... , fn)x =( x1, ... , xn)f=(f1,…,fn)x=(x1,…,xn)\mathbf{f}=(f_1,\dots,f_n)\quad\mathbf{x}=(x_1,\dots,x_n) Sistem ini memiliki sejumlah karakteristik yang membuatnya sangat sulit untuk ditangani. Saya mencari …

2
Tentang perkiraan log yang lebih cepat (x)
Saya telah menulis sebuah kode waktu lalu yang berusaha untuk menghitung tanpa menggunakan fungsi perpustakaan. Kemarin, saya meninjau kode lama, dan saya mencoba membuatnya secepat mungkin, (dan benar). Inilah upaya saya sejauh ini:l o g( x )lHaig(x)log(x) const double ee = exp(1); double series_ln_taylor(double n){ /* n = e^a * …

2
Berapa banyak regularisasi yang ditambahkan untuk membuat SVD stabil?
Saya telah menggunakan Intel MKL SVD ( dgesvdmelalui SciPy) dan memperhatikan bahwa hasilnya sangat berbeda ketika saya mengubah presisi di antara float32dan float64ketika matriks saya dikondisikan secara buruk / bukan peringkat penuh. Apakah ada panduan tentang jumlah minimum regularisasi yang harus saya tambahkan untuk membuat hasil tidak sensitif terhadap float32-> …

2
Apa yang dikatakan analisis stabilitas Von Neumann tentang persamaan perbedaan hingga non-linear?
Saya membaca makalah [1] di mana mereka memecahkan persamaan non-linear berikut menggunakan metode beda hingga. Mereka juga menganalisis stabilitas skema menggunakan analisis stabilitas Von Neumann. Namun, seperti yang disadari oleh penulis, ini hanya berlaku untuk PDE linier. Jadi penulis menyiasatinya dengan "pembekuan" istilah non-linear, yaitu mereka mengganti u u x …

1
Berulang “solver” untuk
Saya tidak bisa membayangkan saya orang pertama yang memikirkan masalah berikut, jadi saya akan puas dengan referensi (tapi jawaban yang lengkap dan terperinci selalu dihargai): Katakanlah Anda memiliki kepastian positif simetris . n dianggap sangat besar, jadi menahan Σ dalam memori tidak mungkin. Anda dapat, bagaimanapun, mengevaluasi Σ x , …

1
Algoritma untuk menghitung eksponensial dari matriks Hessenberg
Saya tertarik dalam menghitung solusi sistem lage ODE menggunakan metode krylov seperti pada [1]. Metode semacam itu melibatkan fungsi-fungsi yang berkaitan dengan fungsi eksponensial (yang disebut fungsi ). Ini pada dasarnya terdiri dari menghitung tindakan fungsi matriks dengan membangun ruang bagian Krylov menggunakan iterasi Arnoldi dan memproyeksikan fungsi pada ruang …


Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.