# Génération d'une courbe cible de type Linkwitz-Riley
# Pierre Ragnarsson Avril 2020
# inspiré du programme Python de Jean-Marc Plantefeve pour les cyble filtre Duelund

import math
import cmath

fc=1600                       # fréquence caractéristique
inf=16                        # fréquence inférieure
sup=16000                     # fréquence supérieure
N=100                         # nombre de points
offset=65                     # décalage en dB

file=open("LPLR2.frd","w")       # nom du fichier .frd

for i in range(0,N+1):
    f=inf*(sup/inf)**(i/N)
    LP=1/(1+2j*f/fc+(1j*f/fc)**2)          # fonction de transfert
    file.write(str(f))
    file.write("\t")
    file.write(str(20*math.log10(abs(LP))+offset))
    file.write("\t")
    file.write(str(57.29*cmath.phase(LP)))
    file.write("\n")
file.close()

file=open("HPLR2.frd","w")
for i in range(0,N+1):
    f=inf*(sup/inf)**(i/N)
    HP=((1j*f/fc)**2)/(1+2j*f/fc+(1j*f/fc)**2)          # fonction de transfert
    file.write(str(f))
    file.write("\t")
    file.write(str(20*math.log10(abs(HP))+offset))
    file.write("\t")
    file.write(str(57.29*cmath.phase(HP)))
    file.write("\n")
file.close()