Для некоторых демонстрационных проблем я пытаюсь смоделировать ток конденсатора.Моя схема действительно простая.Он имеет 2 синусоидальных источника напряжения (v1 (t) = 2 * cos (2 * pi * 15 * t), v2 (t) = 2 * cos (2 * pi * 30 * t) и один конденсатор 1F. Вы можетесм. мою схему и графику моделирования. На верхней панели находятся мои входные напряжения. А на нижней панели графика - падение напряжения (пурпурный) на конденсаторе и ток конденсатора (коричневый). LTspice sim. И этомой код Python (Чтобы найти ограничение тока. Я просто попытался различить падение напряжения через предел.):
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()
Python sim.
И здесьмой вопрос: все 3 графика (v1, v2, v1-v2) из LTspice и python совпадают, за исключением тока конденсатора. В LTspice колеблется между 150A и -150A, но в питоне ток составляет почти 0A. Если вы мне поможетенайти, где я делаю неправильно, я ценю это.