import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
df = pd.read_csv('rawData.csv')
df
Vps = df.Vps
Vpr = df.Vpr
R = 1000
I = (Vps - Vpr)/R
m, b = np.polyfit(I, Vpr, 1)
Vth = m*I + b
plt.title("Tinkercad data")
plt.xlabel("Current (A)")
plt.ylabel("Power supply voltage (V)")
plt.plot(I, Vpr, 'b.', label='data')
plt.plot(I, Vth, 'r-', label='fit')
plt.legend()
df2 = pd.read_csv('P0-LTSpice.txt',delimiter='\t')
df2
df2.columns
Vps = df2['V(n001)']
Vpr = df2['V(n002)']
Vps, Vpr
Vps = df2['V(n001)']
Vpr = df2['V(n002)']
R = 1000
I = (Vps - Vpr)/R
m, b = np.polyfit(I, Vpr, 1)
Vth = m*I + b
plt.title("LTSpice data")
plt.xlabel("Current (A)")
plt.ylabel("Power supply voltage (V)")
plt.plot(I, Vpr, 'b.', label='data')
plt.plot(I, Vth, 'r-', label='fit')
plt.legend()