%Programme Permettant dimensionnement des filtres

%Structure de Rauch - Passe bande

%1er Partie : Calcule des composants

    %definition des fréquences de chacune des cordes (E1 A1 D2 G2 S2 E3):
    fm = [82.41;110;146.83;196;246.94;329.63];
    %fréquence centrale du filtre
    wn = 2*pi.*fm;
    %on pose le facteur de qualité:
    Q=[150;200;300;400;500;600];
    %Gain à la fréquence de coupure
    K = -2;
    %on fixe les condensateurs (on impose C=C1=C)
    C=10^(-6);
     
    %résistances :
    R2 = (2.*Q)./(wn*C)
    R1 = -R2/(2*K)
    R3 = -(K*R1)./(2.*Q.^2+K)

    %Bande passante du filtre 
    BW = 2./(R2*C)



%---------------------------------------------------



 %2ème partie : Calcule de la fonction de transfert & diagramme de Bode 
    
    %Fonction de transfert sous la forme H(p) = (z*p)/(ap²+bp+c)
   
    %On va créer une boucle qui permettra de calculer la fonction de transfert 
    %associés à chacune des résistances que nous avons obtenue dans la 1er partie 
    for i=1:6
        z = C*R2(i)*R3(i)/(R1(i)+R3(i));
        a = (C^2)*(R1(i)*R2(i)*R3(i))/(R1(i)+R3(i));
        b = 2*C*R1(i)*R3(i)/(R1(i)+R3(i));
        c = 1;
        num = [-z 0];
        den= [a b c];
        sys(i) = tf(num,den);   % /!\ tf = transfert function (et non transformé de fourier)
    end
    
    %Trace les 6 diagrammes de Bode 
    bode(sys(1),sys(2),sys(3),sys(4),sys(5),sys(6));
    %Les légendes sinon on ne comprend rien, surtout qu'il y en 6
    legend('Filtre Mi','Filtre La','Filtre Re','Filtre Sol','Filtre Si','Filtre Mi');