Проверка симуляции LTspice с помощью решателя SciPy ODE
Для некоторых демонстрационных проблем я пытаюсь смоделировать ток конденсатора. Моя схема действительно простая. Он имеет 2 синусоидальных источника напряжения (v1(t) = 2*cos(2*pi*15*t), v2(t) = 2*cos(2*pi*30*t) и один конденсатор 1F. Вы можете см. мою схему и графику моделирования. На верхней панели есть мои входные напряжения. А на нижней панели графика - падение напряжения (пурпурный) на конденсаторе и ток конденсатора (коричневый). LTspice sim. И это мой код на Python (In Для того, чтобы найти ограничение тока, я только что попытался дифференцировать падение напряжения через предел.):
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
def v1(u,t): # input1
return (2*np.cos(2*np.pi*5*t))
def v2(u,t): #input2
return (2*np.cos(2*np.pi*10*t))
def v(u,t): # voltage drop
return (v1(u,t)-v2(u,t))
n = np.zeros(1000)
t = np.linspace(0,1,1000)
v_ = v(u,t)
v1_ = v1(u,t)
v2_ = v2(u,t)
plt.subplot(2,1,1)
plt.plot(t,v1_,'green',t,v2_,'r')
plt.show()
v0 = [0]
sol = odeint(v,v0,t)
plt.subplot(2,2,1)
plt.plot(t,sol,'brown',t,(v1_-v2_),'blue')
plt.show()
И вот мой вопрос: все 3 графика (v1,v2,v1-v2) от LTspice и python совпадают, кроме тока конденсатора. В LTspice ток колеблется между 150А и -150А, но в питоне ток составляет почти 0А. Если вы поможете мне найти, где я делаю неправильно, я ценю это.