fa2338/Masterprojekt-Campusnetz
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Jacobi Matrix numerisch. S, zw, und Richtung überprüfen

Browse Source
This commit is contained in:
fa2338
2025-12-15 16:11:38 +01:00
parent 8e999c80cc
commit 68809a1319
9 changed files with 3561 additions and 3514 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
Funktionales_Modell.py
View File

@@ -63,9 +63,9 @@ class FunktionalesModell:
X_sp, Y_sp, Z_sp = dict_punkt_symbole[standpunkt]
X_zp, Y_zp, Z_zp = dict_punkt_symbole[zielpunkt]
B_sp, L_sp = sp.symbols(f"B{standpunkt} L{standpunkt}")
alpha = sp.symbols(f"alpha{standpunkt}_{zielpunkt}")
zw = sp.symbols(f"zw{standpunkt}_{zielpunkt}")
s = sp.symbols(f"s{standpunkt}_{zielpunkt}")
alpha = sp.symbols(f"{beobachtungenID}_alpha{standpunkt}_{zielpunkt}")
zw = sp.symbols(f"{beobachtungenID}_zw{standpunkt}_{zielpunkt}")
s = sp.symbols(f"{beobachtungenID}_s{standpunkt}_{zielpunkt}")
if beobachtungsart == "tachymeter_distanz":
beobachtungsgleichung = sp.sqrt(
@@ -166,7 +166,7 @@ class FunktionalesModell:
if koordinatenart == "naeherung_us":
db_zugriff = Datenbankzugriff(self.pfad_datenbank)
dict_koordinaten_B_L = self.berechnungen.geometrische_breite_laenge(db_zugriff.get_koordinaten(koordinatenart))
liste_beobachtungen = db_zugriff.get_beobachtungen_from_beobachtungenid()
substitutionen = {}
for punktnummer, vektor in dict_koordinaten_B_L.items():
@@ -178,6 +178,16 @@ class FunktionalesModell:
substitutionen[B_sym] = vektor[1]
substitutionen[L_Sym] = vektor[2]
for standpunkt, zielpunkt, beobachtungenID, tachymeter_richtung, tachymeter_zenitwinkel, tachymeter_distanz in liste_beobachtungen:
alpha = sp.symbols(f"{beobachtungenID}_alpha{standpunkt}_{zielpunkt}")
zw = sp.symbols(f"{beobachtungenID}_zw{standpunkt}_{zielpunkt}")
s = sp.symbols(f"{beobachtungenID}_s{standpunkt}_{zielpunkt}")
substitutionen[alpha] = tachymeter_richtung
substitutionen[zw] = tachymeter_zenitwinkel
substitutionen[s] = tachymeter_distanz
A_numerisch = A_symbolisch.xreplace(substitutionen)
Export.matrix_to_csv(r"Zwischenergebnisse\Jacobi_Matrix_Numerisch_Iteration0.csv", liste_unbekannte,