def plot_exp(startfi=1): exp_data.plot_exp(startfi=startfi)
def plot_exp(self, startfi=1, lt='y-'): exp_data.plot_exp(startfi=startfi, lt=lt)
a_th_mode1 = controls.TransferFunction([g1*w1**2,0,0],[1,2*z1*w1,w1**2]) w2=17.5*2*pi z2=0.03 g2=-0.0005 a_th_mode2 = controls.TransferFunction([g2*w2**2,0,0],[1,2*z2*w2,w2**2]) a_th_tf = a_th_mode1 + a_th_mode2 import exp_data f2 = exp_data.f2 exp_data.plot_exp() a_th_tf.FreqResp(f2, fignum=3, clear=False) wth1=2.8*2*pi zth1=0.1 g_th=0.000003 p_th=10.0*2*pi wth2=19.5*2*pi zth2=0.07 num_th = g_th*p_th*wth1**2*wth2**2 * a_th_tf.den p0 = poly1d([1,p_th,0]) p1 = poly1d([1,2*zth1*wth1,wth1**2]) p2 = poly1d([1,2*zth2*wth2,wth2**2]) den_th = p0*p1*p2