// *************************************************
// Tracé expérimental diagramme de Bode
// *************************************************
//
// *************************************************
// Saisie des valeurs relevées
// *************************************************
w=[2 5 10 20 30 40 50 60 70 80]; // Pulsation en rad/s
GdB=[13 12 12 11 11 9.4 8.1 6.4 5.8 4.3]; // Gain en dB
phi=-[5 11 17 28 48 60 63 96 104 115]; // Phase en degré


//
// Affichage des courbes
scf(0);clf(0);
subplot(211);
plot2d(w,GdB)
a=gca()
a.log_flags=('lnn');
a.grid=[1,1]
a.sub_ticks=[1,4];
a.x_location="origin";
a.background=color(240,255,240)
a.x_label.text="$\omega\ (rad/s)$"
a.x_label.font_size=5
a.y_label.text="$Gain (dB)$"
a.y_label.font_size=5
a.children.children.thickness = 3
subplot(212);plot2d(w,phi)
a=gca()
a.log_flags=('lnn');
a.grid=[1,1]
a.sub_ticks=[1,4];
a.x_location="origin";
a.background=color(240,255,240)
a.y_label.text="$Phase (°)$"
a.y_label.font_size=5
a.children.children.thickness = 3