Booster la simulation de courant

le droit est un calcul séquentiel

# coding: utf-8

# In[278]:


get_ipython().magic('matplotlib inline')
import matplotlib.pyplot as plt


# In[279]:


import numpy as np
import math


# In[280]:


points = 360
ac_freq = 50 #Hz
ac_t = 1/ac_freq / 2 #Demi-onde: →/2
ac_volt_real = 220 #V
ac_volt_max = ac_volt_real * math.pow(2,0.5)
reactor_H = 40 #mH
dc_volt = 320
target_dc_volt = 320
career = 40 * math.pow(10,-6)
t_delt = ac_t / points
ep_dead_time_count = 3
dead_time = career * ep_dead_time_count


# In[281]:


def calc_time(i_point):
    return ac_t / points * i_point 

def calc_acv_sin(i_time):
    rad = (i_time / ac_t) * math.pi * 1 #Demi-onde:*2→*1
    acv_sin = ac_volt_max * math.sin(rad)
    return acv_sin

def calc_duty(t, target_dc_volt, dc_volt, acv, aca_sum):
    duty = 0.4
    return duty


# In[282]:


#10ms, 40us, duty:100%, 0.01%Cranté
print("résolution:",10*1000 / 40 * (100/0.1))


# In[283]:


#Avec 250000 points
#Créer une onde sinusoïdale ACV
#40 pour chaque 40us%ON60%OFF, 30%ON70%OFF, ....,En quelque sorte avec des points
#Ou plutôt, c'est impossible avec la carte car il s'agit d'un calcul de droits à partir d'un courant séquentiel ...


# In[284]:


points_list = list(range(0,points))
time_list = list(map(calc_time, points_list))
acv_sin_list = list(map(calc_acv_sin, time_list))


# In[285]:


plt.plot(time_list, acv_sin_list)


# In[286]:


career_count = 0
aca_sum = 0.0
duty = 0
duty_on_count = 0
duty_off_count = 0
duty_on_time = 0
duty_off_time = 0
for t in time_list:
    ac_volt = calc_acv_sin(t)
    if t < dead_time: #Le droit est de 0, mais ACV et ACA sont susceptibles d'être calculés, donc si l'imbrication ou le changement est nécessaire
        duty = 0
    if t >= 0.001 and ac_volt < 30: #1ms,Le droit est de 0, mais ACV et ACA sont susceptibles d'être calculés, donc si l'imbrication ou le changement est nécessaire
        duty = 0
        
    if t >= (career * career_count): #Carrière tous les 40 ans_Sentiment de mettre à jour le compte et le devoir. J'en ai parlé pour calculer les droits pour la première fois
        career_count += 1
        duty = calc_duty(t, target_dc_volt, dc_volt, ac_volt, aca_sum)
        duty_on_count = 0
        duty_off_count = 0
    else:
        if t < ((career * (career_count-1)) + (duty * career)): #Si le devoir est sorti
            duty_off_time = t_delt * duty_off_count
            duty_on_count += 1
        else:
            duty_on_time = t_delt * duty_on_count
            duty_off_count += 1
    
        


# In[ ]: