fc=1000      //fréquence caractéristique du système 3 voies "Synkron"
a=4          //amortissement du système 3 voies "Synkron"
inf=20       //fréquence basse du balayage fréquentiel
sup=20000    //fréquence haute du balayage fréquentiel
N=120        //nombre de points au balayage (aux environs du 1/12eme d'octave)
offset=90    //offset en dB pour rejoindre la sensibilité des HPs

for x = 1:N                          //boucle pour créer 3 matrices de 120 lignes
    f=inf*(sup/inf)**((x-1)/(N-1))
    low=1/(1+a*%i*f/fc-(f/fc)**2)**2
    mid=(a**2-2)*low*(f/fc)**2                    //définition des fonctions LP, BP et HP
    high=low*(f/fc)**4
    Mlow(x,1)=f
    Mlow(x,2)=offset+20*log10(abs(low))           //calculs des 3 colonnes LP (Hz,dB,°)
    Mlow(x,3)=57.3*atan(imag(low),real(low))
    Mmid(x,1)=f
    Mmid(x,2)=offset+20*log10(abs(mid))           //calculs des 3 colonnes BP (Hz,dB,°)
    Mmid(x,3)=57.3*atan(imag(mid),real(mid))
    Mhigh(x,1)=f
    Mhigh(x,2)=offset+20*log10(abs(high))         //calculs des 3 colonnes HP (Hz,dB,°)
    Mhigh(x,3)=57.3*atan(imag(high),real(high))
end

fprintfMat("low.frd",Mlow)
fprintfMat("mid.frd",Mmid)            //écriture des 3 matrices cibles dans 3 fichiers .frd
fprintfMat("high.frd",Mhigh)

clf()
subplot(2,1,1); xgrid(4)
plot2d("ln",Mlow(:,1),[Mlow(:,2),Mmid(:,2),Mhigh(:,2)],style=[2,15,5],rect=[20,offset-40,20000,offset])    //tracés des 3 gains pour visu
subplot(2,1,2); xgrid(4)
plot2d("ln",Mlow(:,1),[Mlow(:,3),Mmid(:,3),Mhigh(:,3)],style=[2,15,5],rect=[20,-180,20000,180])    //tracés des 3 phases pour visu
