Pertanyaan yang diberi tag «newton-method»

3

2
Apakah mungkin untuk menyelesaikan PDE nonlinier tanpa menggunakan iterasi Newton-Raphson?
Saya mencoba memahami beberapa hasil dan akan menghargai beberapa komentar umum tentang penanganan masalah nonlinier. Persamaan Fisher (PDE reaksi-difusi nonlinier), ut=duxx+βu(1−u)=F(u)ut=duxx+βu(1−u)=F(u) u_t = du_{xx} + \beta u (1 - u) = F(u) dalam bentuk diskrit, u′j=Lu+βuj(1−uj)=F(u)uj′=Lu+βuj(1−uj)=F(u) u_j^{\prime} = \boldsymbol{L}\boldsymbol{u} + \beta u_j (1 - u_j) = F(\boldsymbol{u}) di mana LL\boldsymbol{L} …

1
Bisakah Jacobian yang diperkirakan dengan perbedaan hingga menyebabkan ketidakstabilan dalam metode Newton?
Saya telah menerapkan pemecah Euler mundur dalam python 3 (menggunakan numpy). Untuk kenyamanan saya sendiri dan sebagai latihan, saya juga menulis fungsi kecil yang menghitung perkiraan perbedaan hingga dari gradien sehingga saya tidak selalu harus menentukan Jacobian secara analitis (jika itu mungkin!). Menggunakan deskripsi yang disediakan dalam Ascher dan Petzold …

2
Strategi untuk Metode Newton ketika Jacobian pada solusinya adalah tunggal
Saya mencoba memecahkan sistem persamaan berikut untuk variabel dan x 2 (semuanya adalah konstanta):P, x1P,x1P,x_1x2x2x_2 A ( 1 - P)2- k1x1= 0A P2- k2x2= 0( 1 - P) ( r1+ x1)4L1- P( r1+ x2)4L2= 0A(1−P)2−k1x1=0AP2−k2x2=0(1−P)(r1+x1)4L1−P(r1+x2)4L2=0\frac{A(1-P)}{2}-k_1x_1=0 \\ \frac{AP}{2}-k_2x_2=0 \\ \frac{(1-P)(r_1+x_1)^4}{L_1}-\frac{P(r_1+x_2)^4}{L_2}=0 Saya dapat melihat bahwa saya dapat mengubah sistem persamaan ini menjadi …

3
Metode penyelesaian sistem difusi-adveksi non-linear di luar Newton-Raphson?
Saya sedang mengerjakan sebuah proyek di mana saya memiliki dua domain digabungkan adv-diff melalui istilah sumber masing-masing (satu domain menambah massa, yang lain mengurangi massa). Untuk singkatnya, saya memodelkan mereka dalam kondisi stabil. Persamaannya adalah persamaan transport adveksi-difusi standar Anda dengan istilah sumber terlihat seperti ini: ∂c1∂t=0=F1+Q1(c1,c2)∂c2∂t=0=F2+Q2(c1,c2)∂c1∂t=0=F1+Q1(c1,c2)∂c2∂t=0=F2+Q2(c1,c2) \frac{\partial c_1}{\partial t} …
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.