Devoir
JAAFARI AMINA
disp('%%%%%%%------------------ Cours : Circuit RLC ------------------%%%%%%%');
R = 10; %Résistance en Ohm
L = 0.5e-3; %Inductance en H
C = 10e-9; %Capacité en F
w = [Link]; %Plage de variation de la pulsation
% Paramètres pour les trois facteurs d'amortissement
R1 = 1/(2*pi*sqrt(L*C)*0.05); % Résistance pour Q1 = 0.05
R2 = 1/(2*pi*sqrt(L*C)*0.5); % Résistance pour Q2 = 0.5
R3 = 1/(2*pi*sqrt(L*C)*5); % Résistance pour Q3 = 5
Num = 1; %Numérateur de H(p)
Denom = [L*C,R*C,1]; %Dénominateur de H(p)
Denom1 = [L*C,R1*C,1]; % Dénominateur pour Q1 = 0.05
Denom2 = [L*C,R2*C,1]; % Dénominateur pour Q2 = 0.5
Denom3 = [L*C,R3*C,1]; % Dénominateur pour Q3 = 5
H = polyval(Num,1i*w)./polyval(Denom,1i*w);
H1 = polyval(Num,1i*w)./polyval(Denom1,1i*w);
H2 = polyval(Num,1i*w)./polyval(Denom2,1i*w);
H3 = polyval(Num,1i*w)./polyval(Denom3,1i*w);
Feq_Coupure = 2*pi/sqrt(L*C);
Q1 = 1/(R1*C*Feq_Coupure); % Facteur d'amortissement pour Q1
Q2 = 1/(R2*C*Feq_Coupure); % Facteur d'amortissement pour Q2
Q3 = 1/(R3*C*Feq_Coupure); % Facteur d'amortissement pour Q3
disp(['Feq_Coupure = ',num2str(Feq_Coupure),' Hz']);
disp(['Facteur d''amortissement Q1 = ',num2str(Q1)]);
disp(['Facteur d''amortissement Q2 = ',num2str(Q2)]);
disp(['Facteur d''amortissement Q3 = ',num2str(Q3)]);
figure;
subplot(2,1,1);
semilogx(w, 20*log10(abs(H)), 'm', 'linewidth', 3);
hold on;
semilogx(w, 20*log10(abs(H1)), 'r', 'linewidth', 3);
semilogx(w, 20*log10(abs(H2)), 'g', 'linewidth', 3);
semilogx(w, 20*log10(abs(H3)), 'b', 'linewidth', 3);
axis([0, 1e8, -100, 50]);
grid on;
xlabel('Fréquence (Hz)');
ylabel('Gain (dB)');
title('Réponse en gain');
legend('Q = 0.05', 'Q = 0.5', 'Q = 5');
subplot(2,1,2);
semilogx(w, angle(H)*180/pi, 'm', 'linewidth', 3);
hold on;
semilogx(w, angle(H1)*180/pi, 'r', 'linewidth', 3);
semilogx(w, angle(H2)*180/pi, 'g', 'linewidth', 3);
semilogx(w, angle(H3)*180/pi, 'b', 'linewidth', 3);
axis([0, 1e8, -200, 20]);
grid on;
xlabel('Fréquence (Hz)');
ylabel('Phase (°)');
title('Réponse en phase');
legend( 'Q = 0.05', 'Q = 0.5', 'Q = 5');
