Come uso il metodo di NEWTON RAPHSON per un sistema di eq. non lineari?
voglio rimuovere {ξ}={ξ₁,ξ₂,...,ξm} tali che fi(ξ₁,ξ₂,...) = Yi
Applico Taylor prendendo una n-upla iniziale {X₁₀,X₂₀,...Xm₀} a piacere (di solito metto E=Em e Θ=Θ 'come valori iniziali')
lo trascuro
Yi=fi(X₁₀,X₂₀,...,Xm₀) + ∂fi/∂x₁ ΔX1*+...+∂fi/∂xm ΔXm*
In forma matriciale:
{Δy1
- Δy1
- Δy2
- Δym
{ΔΘ
{ΔE
- = matice dei termini noti
- Ϟ = matrice dei coefficienti del sistema linearizzato
Dunque trovo ΔXm+
{ fi(X₁,X₂,...) - Yi } < Ѓ
in generale | fi(X1ⁿいm) - Yi| < Ѓ
Tema: uso il metodo di Newton Raphson per un sistema di eq. non lineari?
y1 = f1(x1, x2, ..., xm)y2 = f2(x1, x2, ..., xm)...ym = fm(x1, x2, ..., xm)
Voglio trovare {Ξ} = {Ξ1, Ξ2, ..., Ξm} tali che fi(Ξ1, Ξ2, ...) = yi
Applico Taylor prendendo una n-upla iniziale {x10, x20, ..., xm0} a piacere (di fatto metto E = Em e Θ = 0 come valori iniziali)
fi(x10 + Δx1, x20 + Δx2, ..., xm0 + Δxm) = yi
fi(x10 + Δx1, ..., xm0 + Δxm) = fi(x10, x20, ..., xmo) + ∂fi/∂x1 Δx1* + ... + Ri(xio, xm0)
(lo trascuro)... + ∂fi/∂xk Δxk* +...
∂fi/∂xm Δxm*
yi - fi(x10, x20, ..., xmo) = Σj=1n ∂f/∂xj Δxj*
(valutate con la n-upla iniziale)
In forma matriciale:
| y1 - f10 | | ∂f1/∂x1 | ... | ∂f1/∂xm | | Δx1 || y2 - f20 | = | ∂f2/∂x1 | ... | ∂f2/∂xm | | Δx2 || ym - fm0 | | ∂fm/∂x1 | ... | ∂fm/∂xm | | Δxm |
n × 1 m × m
m vettore dei MATRICE DEI COEFFICIENTI
iincrementi DEL SISTEMA LINEARIZZATO
(Jacobiano)
Dunque trovo Δxi+1 → xi+1 = xi0 + Δxi x24 = x20 - Δx2
| fi(x1, x2, ...) - yi | < ε → SI ← le soluzioni
NO
↓ ricomincio con in generale | fi(xm) - yi | < ε
le iterazioni il piccolo
Nel nostro caso all'iterazione i-esima:
| P0 - P0 | | ∂P/∂Θ | | | [ΔΘ]i || Pm+1,n | = | | [ΔE]i |
| Qm+1,n | | |
[ΠΩ,n],Θ,n | ∂Q/∂E | | |
ESEMPIO:
Scrivo dunque le quattro equazioni:
P1 = 3E2E1y21cos(θ2-θ1-χ21) + 3E2y22cos(χ22) + 3E2E3y23cos(θ2-θ3-χ23)
P3 = 3E3E2y32cos(θ3-θ2-χ32) + 3E3y33cos(χ33)
Q1 = 3E2E1y21sin(θ2-θ1-χ21) + 3E2y22sin(-χ22) + 3E2E3y23sin(θ2-θ3-χ23)
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
-
Appunti di Elementi di impianti e centrali elettriche (parte 3)
-
Appunti di Elementi di impianti e centrali elettriche (parte 1)
-
Elementi di impianti e centrali elettriche
-
Esercizi degli Esami passati di Elementi di impianti e centrali elettriche