import pandas as pd
import matplotlib.pyplot as plt
import numpy as np
df = pd.read_csv("TinkerCAD.csv")
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.ylabel("Voltage")
plt.xlabel("Current(mA)")
plt.plot(I, Vpr, 'b.', label='data')
plt.plot(I, Vth, 'r-', label='fit')
plt.legend()
lt = pd.read_csv("LtSpice.csv")
Vps2 = lt.v1
Vpr2 = lt.V2
R = 1000
I = (Vps2 - Vpr2)/R
m, b = np.polyfit(I, Vpr2, 1)
Vth = m*I + b
plt.title("LTSpice data")
plt.ylabel("Voltage")
plt.xlabel("Current(mA)")
plt.plot(I, Vpr2, 'b.', label='data')
plt.plot(I, Vth, 'r-', label='fit')
plt.legend()
m, b = np.polyfit(I, Vpr2, 1)