def find_symbolic_bodes(self, save=1): self.Usys = self.U6*(self.U5*(self.U4*(self.U3*\ (self.U2*(self.U1*self.U0))))) z_b = sympy_TMM.find_base_vector(self.Usys) z_enc = self.U2*(self.U1*(self.U0*z_b)) enc_gain = 180.0/pi*1024.0/360.0 theta = z_enc[1]*enc_gain z_a = self.U5*(self.U4*(self.U3*z_enc)) atip = s**2*z_a[0]*a_gain self.sym_bodes = [theta, atip] self.Usys = self.Usys if save: self.cse_to_file() return self.sym_bodes
def find_symbolic_bodes(self, save=1): U0 = AVS.Get_Aug_Mat(s) U1 = TSD.Get_Aug_Mat(s) U2 = beam.Get_Aug_Mat(s) Usys = U2*(U1*U0) z_b = sympy_TMM.find_base_vector(Usys) z_enc = U1*(U0*z_b) enc_gain = 180.0/pi*1024.0/360.0 theta = z_enc[1]*enc_gain z_a = U2*z_enc atip = s**2*z_a[0]*a_gain self.sym_bodes = [theta, atip] self.Usys = Usys if save: self.cse_to_file() return self.sym_bodes
def find_symbolic_bodes(self, save=1): ## self.list = [self.avs, self.base_mass, \ ## self.clamp_spring, \ ## self.beam, self.accel, self.beam2] U0 = AVS_ThetaFB.Get_Aug_Mat(s) U1 = Base_Mass.Get_Aug_Mat(s) U2 = TSD_clamp.Get_Aug_Mat(s) U3 = beam1.Get_Aug_Mat(s) U4 = Accel_Mass.Get_Aug_Mat(s) U5 = beam2.Get_Aug_Mat(s) Usys = U5*(U4*(U3*(U2*(U1*U0)))) z_b = sympy_TMM.find_base_vector(Usys) z_enc = U1*(U0*z_b) enc_gain = 180.0/pi*1024.0/360.0 theta = z_enc[1]*enc_gain z_a = U4*(U3*(U2*z_enc)) atip = s**2*z_a[0]*a_gain self.sym_bodes = [theta, atip] self.Usys = Usys if save: self.cse_to_file() return self.sym_bodes