Noah main sim maintenence

.gitignore

@@ -213,3 +213,4 @@ Projects/.DS_Store
 .fake
 *.joblib
 FullVehicleSim/simulation_output.parquet
+.vscode/settings.json

Data/FSLib/IntegralsAndDerivatives.py

@@ -7,6 +7,7 @@
 from scipy.integrate import RK45
 # from scipy.integrate import solve_ivp
 import matplotlib.pyplot as plt
+from scipy.integrate import cumulative_trapezoid
 
 
 # lat = "VDM_GPS_Latitude"
@@ -67,20 +68,20 @@ def mag (a, b, c):
     # print(f"mag = {np.sqrt(a**2 + b**2 + c**2)}")
     return np.sqrt(a**2 + b**2 + c**2)
 
-def in_place_integrate (derivative, dt=0.01): #Rimann sum that returns an array of the intergral at that point. Assumes dt of 0.01
-    if len(derivative.shape) == 1:
-        width = 1
-    else:
-        width = derivative.shape[1]
-    container = derivative[0]
-    out = np.zeros((derivative.shape[0], width))
-    out[0] = container*dt
-    for i in range(1, derivative.shape[0]):
-        container += dt*derivative[i]
-        # print(f"added {0.01*derivative[i]} to {container}")
-        out[i] = container
-    # print(f"out of integral is {out}")
-    return out
+def in_place_integrate (derivative, t, initial=0): #Rimann sum that returns an array of the intergral at that point. Assumes dt of 0.01
+    # if len(derivative.shape) == 1:
+    #     width = 1
+    # else:
+    #     width = derivative.shape[1]
+    # container = derivative[0]
+    # out = np.zeros((derivative.shape[0], width))
+    # out[0] = container*dt
+    # for i in range(1, derivative.shape[0]):
+    #     container += dt*derivative[i]
+    #     # print(f"added {0.01*derivative[i]} to {container}")
+    #     out[i] = container
+    # # print(f"out of integral is {out}")
+    return cumulative_trapezoid(derivative, t, initial)
 
 def integrate_with_tCol (col, timeCol):
     if col.len() != timeCol.len():

Data/README_TRACTIVE_THERMAL_MODELING.md

@@ -0,0 +1,32 @@
+# Run Thermal Modeling
+Todo
+
+### Run Tests
+
+```bash
+python -m unittest discover -s test
+```
+
+#### Plot the curve
+```bash
+python ./TractiveBatteryThermalModelViewer.py -h
+usage: TractiveBatteryThermalModelViewer.py [-h] --path_parquet PATH_PARQUET [--column-name COLUMN_NAME] [--t-end T_END] [--initial-temp INITIAL_TEMP]
+
+View tractive battery thermal model as a function of current draw.
+
+options:
+  -h, --help            show this help message and exit
+  --path_parquet PATH_PARQUET
+                        Path to the Parquet file containing current data.
+  --column-name COLUMN_NAME
+                        Name of the current column in the Parquet file (default: SME_TEMP_BusCurrent).
+  --t-end T_END         End time in seconds for the simulation (default: 60).
+  --initial-temp INITIAL_TEMP
+                        Initial temperature in °C (default: 22).
+```
+
+Example Command :
+
+```
+python TractiveBatteryThermalModelViewer.py --path_parquet <file_name>
+```

Data/SensorFrequencyAnalysis/apps.py

@@ -0,0 +1,41 @@
+import scipy
+import numpy as np
+import polars as pl
+import matplotlib.pyplot as plt
+
+df = pl.read_parquet('/Users/aanyajain/Documents/GitHub/fs-data/FS-3/01112026/011026-23.parquet')
+df = df.fill_null(strategy="forward").fill_null(strategy="backward")
+signal = df["ETC_STATUS_PEDAL_TRAVEL"].to_numpy()
+
+fft = scipy.fft.fft(signal)
+freqs = scipy.fft.fftfreq(len(signal), d=0.01)
+
+attempts = [6, 6.5, 7]
+ffts = []
+iffts = []
+for i, attempt in enumerate(attempts):
+    ffts.append(fft.copy())
+    ffts[-1][np.log(np.abs(fft)) < attempt] = 0
+    ifft = scipy.fft.ifft(ffts[-1])
+    iffts.append(ifft)
+
+RMSs = [np.sqrt(np.mean((ifft.real - signal)**2)) for ifft in iffts]
+plt.plot(attempts, RMSs)
+plt.xlabel("Attempt")
+plt.ylabel("RMS")
+plt.title("APPS - RMS vs Attempt")
+plt.show()
+
+plt.plot(signal, label="Real")
+for ifft, attempt in zip(iffts, attempts):
+    plt.plot(ifft.real, label=f"Attempt {attempt}")
+plt.xlabel("Time")
+plt.ylabel("Amplitude")
+plt.legend()
+plt.show()
+
+plt.scatter(freqs, np.log(np.abs(fft)), s=0.25)
+plt.xlabel("Frequency (Hz)")
+plt.ylabel("Magnitude")
+plt.title("APPS - Frequency Spectrum")
+plt.show()

Data/SensorFrequencyAnalysis/battery_current.py

@@ -0,0 +1,41 @@
+import scipy
+import numpy as np
+import polars as pl
+import matplotlib.pyplot as plt
+
+df = pl.read_parquet('/Users/aanyajain/Documents/GitHub/fs-data/FS-3/11222025/11222025_22.parquet')
+df = df.fill_null(strategy="forward").fill_null(strategy="backward")
+signal = df["ACC_POWER_CURRENT"].to_numpy()
+
+fft = scipy.fft.fft(signal)
+freqs = scipy.fft.fftfreq(len(signal), d=0.01)
+
+attempts = [5, 6, 7]
+ffts = []
+iffts = []
+for i, attempt in enumerate(attempts):
+    ffts.append(fft.copy())
+    ffts[-1][np.log(np.abs(fft)) < attempt] = 0
+    ifft = scipy.fft.ifft(ffts[-1])
+    iffts.append(ifft)
+
+RMSs = [np.sqrt(np.mean((ifft.real - signal)**2)) for ifft in iffts]
+plt.plot(attempts, RMSs)
+plt.xlabel("Attempt")
+plt.ylabel("RMS")
+plt.title("Battery Current - RMS vs Attempt")
+plt.show()
+
+plt.plot(signal, label="Real")
+for ifft, attempt in zip(iffts, attempts):
+    plt.plot(ifft.real, label=f"Attempt {attempt}")
+plt.xlabel("Time")
+plt.ylabel("Amplitude")
+plt.legend()
+plt.show()
+
+plt.scatter(freqs, np.log(np.abs(fft)), s=0.25)
+plt.xlabel("Frequency (Hz)")
+plt.ylabel("Magnitude")
+plt.title("Battery Current - Frequency Spectrum")
+plt.show()

Data/SensorFrequencyAnalysis/battery_tray_temp.py

@@ -0,0 +1,41 @@
+import scipy
+import numpy as np
+import polars as pl
+import matplotlib.pyplot as plt
+
+df = pl.read_parquet('/Users/aanyajain/Documents/GitHub/fs-data/FS-3/11222025/11222025_1.parquet')
+df = df.fill_null(strategy="forward").fill_null(strategy="backward")
+signal = df["ACC_TRAY_TEMPS_BUSBAR"].to_numpy()
+
+fft = scipy.fft.fft(signal)
+freqs = scipy.fft.fftfreq(len(signal), d=0.01)
+
+attempts = [4, 5, 6]
+ffts = []
+iffts = []
+for i, attempt in enumerate(attempts):
+    ffts.append(fft.copy())
+    ffts[-1][np.log(np.abs(fft)) < attempt] = 0
+    ifft = scipy.fft.ifft(ffts[-1])
+    iffts.append(ifft)
+
+RMSs = [np.sqrt(np.mean((ifft.real - signal)**2)) for ifft in iffts]
+plt.plot(attempts, RMSs)
+plt.xlabel("Attempt")
+plt.ylabel("RMS")
+plt.title("Battery Tray Temp - RMS vs Attempt")
+plt.show()
+
+plt.plot(signal, label="Real")
+for ifft, attempt in zip(iffts, attempts):
+    plt.plot(ifft.real, label=f"Attempt {attempt}")
+plt.xlabel("Time")
+plt.ylabel("Amplitude")
+plt.legend()
+plt.show()
+
+plt.scatter(freqs, np.log(np.abs(fft)), s=0.25)
+plt.xlabel("Frequency (Hz)")
+plt.ylabel("Magnitude")
+plt.title("Battery Tray Temp - Frequency Spectrum")
+plt.show()

Data/SensorFrequencyAnalysis/battery_voltage.py

@@ -0,0 +1,41 @@
+import scipy
+import numpy as np
+import polars as pl
+import matplotlib.pyplot as plt
+
+df = pl.read_parquet('/Users/aanyajain/Documents/GitHub/fs-data/FS-3/01112026/011026-14.parquet')
+df = df.fill_null(strategy="forward").fill_null(strategy="backward")
+signal = df["ACC_POWER_PACK_VOLTAGE"].to_numpy()
+
+fft = scipy.fft.fft(signal)
+freqs = scipy.fft.fftfreq(len(signal), d=0.01)
+
+attempts = [4.5, 5.5, 6.5]
+ffts = []
+iffts = []
+for i, attempt in enumerate(attempts):
+    ffts.append(fft.copy())
+    ffts[-1][np.log(np.abs(fft)) < attempt] = 0
+    ifft = scipy.fft.ifft(ffts[-1])
+    iffts.append(ifft)
+
+RMSs = [np.sqrt(np.mean((ifft.real - signal)**2)) for ifft in iffts]
+plt.plot(attempts, RMSs)
+plt.xlabel("Attempt")
+plt.ylabel("RMS")
+plt.title("Battery Voltage - RMS vs Attempt")
+plt.show()
+
+plt.plot(signal, label="Real")
+for ifft, attempt in zip(iffts, attempts):
+    plt.plot(ifft.real, label=f"Attempt {attempt}")
+plt.xlabel("Time")
+plt.ylabel("Amplitude")
+plt.legend()
+plt.show()
+
+plt.scatter(freqs, np.log(np.abs(fft)), s=0.25)
+plt.xlabel("Frequency (Hz)")
+plt.ylabel("Magnitude")
+plt.title("Battery Voltage - Frequency Spectrum")
+plt.show()

Data/SensorFrequencyAnalysis/brake_temp.py

@@ -0,0 +1,41 @@
+import scipy
+import numpy as np
+import polars as pl
+import matplotlib.pyplot as plt
+
+df = pl.read_parquet('/Users/aanyajain/Documents/GitHub/fs-data/FS-3/11222025/11222025_1.parquet')
+df = df.fill_null(strategy="forward").fill_null(strategy="backward")
+signal = df["IZZE_BRAKETEMP_S1_CH1"].to_numpy()
+
+fft = scipy.fft.fft(signal)
+freqs = scipy.fft.fftfreq(len(signal), d=0.01)
+
+attempts = [5.5, 6.5, 7.5]
+ffts = []
+iffts = []
+for i, attempt in enumerate(attempts):
+    ffts.append(fft.copy())
+    ffts[-1][np.log(np.abs(fft)) < attempt] = 0
+    ifft = scipy.fft.ifft(ffts[-1])
+    iffts.append(ifft)
+
+RMSs = [np.sqrt(np.mean((ifft.real - signal)**2)) for ifft in iffts]
+plt.plot(attempts, RMSs)
+plt.xlabel("Attempt")
+plt.ylabel("RMS")
+plt.title("Brake Temp - RMS vs Attempt")
+plt.show()
+
+plt.plot(signal, label="Real")
+for ifft, attempt in zip(iffts, attempts):
+    plt.plot(ifft.real, label=f"Attempt {attempt}")
+plt.xlabel("Time")
+plt.ylabel("Amplitude")
+plt.legend()
+plt.show()
+
+plt.scatter(freqs, np.log(np.abs(fft)), s=0.25)
+plt.xlabel("Frequency (Hz)")
+plt.ylabel("Magnitude")
+plt.title("Brake Temp - Frequency Spectrum")
+plt.show()

Data/SensorFrequencyAnalysis/cell_temp.py

@@ -0,0 +1,41 @@
+import scipy
+import numpy as np
+import polars as pl
+import matplotlib.pyplot as plt
+
+df = pl.read_parquet('/Users/aanyajain/Documents/GitHub/fs-data/FS-3/11222025/11222025_1.parquet')
+df = df.fill_null(strategy="forward").fill_null(strategy="backward")
+signal = df["ACC_SEG0_TEMPS_CELL0"].to_numpy()
+
+fft = scipy.fft.fft(signal)
+freqs = scipy.fft.fftfreq(len(signal), d=0.01)
+
+attempts = [4.5, 5.5, 6.5]
+ffts = []
+iffts = []
+for i, attempt in enumerate(attempts):
+    ffts.append(fft.copy())
+    ffts[-1][np.log(np.abs(fft)) < attempt] = 0
+    ifft = scipy.fft.ifft(ffts[-1])
+    iffts.append(ifft)
+
+RMSs = [np.sqrt(np.mean((ifft.real - signal)**2)) for ifft in iffts]
+plt.plot(attempts, RMSs)
+plt.xlabel("Attempt")
+plt.ylabel("RMS")
+plt.title("Cell Temp - RMS vs Attempt")
+plt.show()
+
+plt.plot(signal, label="Real")
+for ifft, attempt in zip(iffts, attempts):
+    plt.plot(ifft.real, label=f"Attempt {attempt}")
+plt.xlabel("Time")
+plt.ylabel("Amplitude")
+plt.legend()
+plt.show()
+
+plt.scatter(freqs, np.log(np.abs(fft)), s=0.25)
+plt.xlabel("Frequency (Hz)")
+plt.ylabel("Magnitude")
+plt.title("Cell Temp - Frequency Spectrum")
+plt.show()

Data/SensorFrequencyAnalysis/cell_voltage.py

@@ -0,0 +1,41 @@
+import scipy
+import numpy as np
+import polars as pl
+import matplotlib.pyplot as plt
+
+df = pl.read_parquet('/Users/aanyajain/Documents/GitHub/fs-data/FS-3/01112026/011026-14.parquet')
+df = df.fill_null(strategy="forward").fill_null(strategy="backward")
+signal = df["ACC_SEG0_VOLTS_CELL0"].to_numpy()
+
+fft = scipy.fft.fft(signal)
+freqs = scipy.fft.fftfreq(len(signal), d=0.01)
+
+attempts = [1.5, 2.5, 3.5]
+ffts = []
+iffts = []
+for i, attempt in enumerate(attempts):
+    ffts.append(fft.copy())
+    ffts[-1][np.log(np.abs(fft)) < attempt] = 0
+    ifft = scipy.fft.ifft(ffts[-1])
+    iffts.append(ifft)
+
+RMSs = [np.sqrt(np.mean((ifft.real - signal)**2)) for ifft in iffts]
+plt.plot(attempts, RMSs)
+plt.xlabel("Attempt")
+plt.ylabel("RMS")
+plt.title("Cell Voltage - RMS vs Attempt")
+plt.show()
+
+plt.plot(signal, label="Real")
+for ifft, attempt in zip(iffts, attempts):
+    plt.plot(ifft.real, label=f"Attempt {attempt}")
+plt.xlabel("Time")
+plt.ylabel("Amplitude")
+plt.legend()
+plt.show()
+
+plt.scatter(freqs, np.log(np.abs(fft)), s=0.25)
+plt.xlabel("Frequency (Hz)")
+plt.ylabel("Magnitude")
+plt.title("Cell Voltage - Frequency Spectrum")
+plt.show()