Einzelnen Beitrag anzeigen
  #3  
Alt 30.03.10, 15:14
Benutzerbild von richy
richy richy ist offline
Singularität
 
Registriert seit: 01.05.2007
Ort: karlsruhe
Beitr?ge: 4.170
Standard AW: kleine numerische Tricks

ENTKOPPELN EINES allgemeinen 1-D PDE SYSTEMS
**************************************
Das Sytem sei in primitiver Form gegeben :

dw/dt + A(x)*dw/dx + C = 0

Alle Groessen sind Vektoren. A die Uebertragungsmatrix.

TRANSFORMATION AUF DIAGONALGESTALT
********************************
Zunaechst ist es notwendig A(x) auf Diagonalgestalt zu transformieren :
D=S_1*A*S

Dabei ist S eine Matrix, deren Spalten aus den RECHTS-Eigenvektoren der Matrix A bestehen.
S ist die inverse Matrix zu S_1
S_1 besteht aus den transponierten Linkseigenvektoren der Matrix A

Die Elemente der Diagonalmatrix D bestehen dann aus den Eigenwerten der Matrix A

Anwenden auf das PDE System
***********************
Multiplikation von links mit S_1

S_1*dw/dt + S_1*A*dw/x + S_1*C = 0


Trick fuer gemischte Form :
Nun multipliziert man A von rechts Mit der Einheitsmatrix E und drueckt diese aus ueber E=S*S_1

S_1*dw/dt + S_1*A*E*dw/x + S_1*C = 0
S_1*dw/dt + S_1*A*S*S_1*dw/x + S_1*C = 0

Charakteristische Form :
S_1*dw/dt + D*S_1*dw/x + S_1*C = 0
***********************************
Dies ist bereits die entkoppelte PDE in charakteristischen Variablen S_1*w, die sich in einem weiteren Schritt auch in primitiven Variablen umformulieren laesst.
Die Methode arbeitet nur fuer lineare PDE's

Ge?ndert von richy (31.03.10 um 00:03 Uhr)
Mit Zitat antworten