Supprimer 'projet.py'
This commit is contained in:
parent
578c6ce06c
commit
7c7a131ebf
57
projet.py
57
projet.py
@ -1,57 +0,0 @@
|
|||||||
# -*-coding:Latin-1 -*
|
|
||||||
from math import *
|
|
||||||
print("Entrez la fréquence de coupure:")
|
|
||||||
frequence = input()
|
|
||||||
#print(f)
|
|
||||||
print("Entrez l'impédance d'entrée:")
|
|
||||||
impedance=input()
|
|
||||||
print("Entrez l'amplitude max:")
|
|
||||||
amplitude=input()
|
|
||||||
print("Entrez l'ordre du filtre:")
|
|
||||||
ordre=input()
|
|
||||||
#Calcul des paramètres
|
|
||||||
#Calcule de wc
|
|
||||||
wc=2*pi*frequence
|
|
||||||
#Calcul de b
|
|
||||||
b=log(cosh(amplitude/17.37)/sinh(amplitude/17.37))
|
|
||||||
#Calcul de gamma
|
|
||||||
gamma=sinh(b/(2*ordre))
|
|
||||||
#initialisation i
|
|
||||||
i=1
|
|
||||||
ak=[0]*ordre
|
|
||||||
bk=[0]*ordre
|
|
||||||
gk=[0]*ordre
|
|
||||||
#Cacule de ak
|
|
||||||
for i in range(ordre):
|
|
||||||
ak[i-1]=sin((2*(i-1)*pi)/(2*ordre))
|
|
||||||
#Calcul de bk
|
|
||||||
for i in range(ordre):
|
|
||||||
bk[i-1]=pow(gamma,2)+(sin(i*pi)/ordre)
|
|
||||||
print("test")
|
|
||||||
print(bk)
|
|
||||||
#Calcul de gk
|
|
||||||
|
|
||||||
gk[0] = ((2*ak[1])/gamma)
|
|
||||||
print(ak)
|
|
||||||
print(gk[0])
|
|
||||||
print(bk[0])
|
|
||||||
for i in range(ordre):
|
|
||||||
gk[i]=((4*ak[i-1]*ak[i])/(bk[i-1]*gk[i-1]))
|
|
||||||
#Calcul de ck et lk
|
|
||||||
if (ordre%2)==0:
|
|
||||||
for i in range(ordre):
|
|
||||||
lk[i]=(impedance/wc)*gk[i]
|
|
||||||
for i in range(ordre):
|
|
||||||
ck[i]=(1/impedance)*(1/wc)*gk[i-1]
|
|
||||||
else:
|
|
||||||
for i in range(ordre):
|
|
||||||
lk[i]=(impedance/wc)*gk[i-1]
|
|
||||||
for i in range(ordre):
|
|
||||||
ck[i]=(1/impedance)*(1/wc)*gk[i]
|
|
||||||
#Calcul de Rn
|
|
||||||
if ordre==1:
|
|
||||||
rn=1*impedance
|
|
||||||
else:
|
|
||||||
rn=pow((tanh(b/4)),2)
|
|
||||||
|
|
||||||
print(ck)
|
|
Loading…
Reference in New Issue
Block a user