s = tf('s');
w0 = 2*pi*300;
Q = 0.9;
G_s = (s^2 + w0^2)/(s^2 + (w0/Q)*s +w0^2);
format long
G_z = c2d(G_s,15e-006,'tustin');
G_z_Numerator= G_z.Numerator{1};
G_z_Denominator=G_z.Denominator{1};
display(G_z_Numerator);
display(G_z_Denominator);
figure('name','Bode of G(s)');
bode(G_s);
评论0