Domande taggate «newton-method»

3


1
Un giacobino approssimativo con differenze finite può causare instabilità nel metodo Newton?
Ho implementato un solutore all'indietro-Euler in Python 3 (usando numpy). Per mia comodità e come esercizio, ho anche scritto una piccola funzione che calcola un'approssimazione di differenza finita del gradiente in modo che non debba sempre determinare analiticamente il giacobino (se è anche possibile!). Usando le descrizioni fornite in Ascher …

2
Strategie per il metodo di Newton quando il giacobino alla soluzione è singolare
Sto cercando di risolvere il seguente sistema di equazioni per le variabili e (tutte le altre sono costanti):P,x1P,x1P,x_1x2x2x_2 A(1−P)2−k1x1=0AP2−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 Vedo che posso trasformare questo sistema di equazioni in una singola equazione di una singola variabile risolvendo le equazioni 1 e 2 rispettivamente per e e sostituendole …

3
Metodi per risolvere sistemi di diffusione dell'avvocatura non lineari oltre Newton-Raphson?
Sto lavorando a un progetto in cui ho due domini accoppiati adv-diff attraverso i rispettivi termini sorgente (un dominio aggiunge massa, l'altro sottrae massa). Per brevità, li sto modellando in stato stazionario. Le equazioni sono le equazioni di trasporto standard di diffusione dell'avvezione con un termine sorgente simile al seguente: …
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.