Subversion Repositories pentevo

Rev

Blame | Last modification | View Log | Download | RSS feed

  1. #!/usr/bin/env python
  2.  
  3. import sys
  4. from sympy import *
  5.  
  6.  
  7.  
  8.  
  9. def main():
  10.         u,v,w = symbols('u v w')
  11.         r,c = symbols('r c')
  12.         n,m = symbols('n m')
  13.         omega = symbols('omega')
  14.         i1,i2,i3 = symbols('i1 i2 i3')
  15.         omega_c,q,rc = symbols('omega_c q rc')
  16.  
  17.         eq1=Eq(i1,(v-w)/(m*r))
  18.         eq2=Eq(i2,(u-w)/(1/(I*omega*n*c)))
  19.         eq3=Eq(i1+i2,i3)
  20.         eq4=Eq(w,i3*(r+1/(I*omega*c)))
  21.         eq5=Eq(u,i3*(1/(I*omega*c)))
  22.  
  23. #       eq6=Eq(rc,r*c)
  24.         eq7=Eq(omega_c,1/(r*c*sqrt(n*m)))
  25.         eq8=Eq(q,sqrt(n*m)/(m+1))
  26.  
  27.         s = solve([eq1,eq2,eq3,eq4,eq5,eq7,eq8],[w,u,i1,i2,i3,r,c],dict=True,domain=S.Complexes)
  28.  
  29.         if len(s)!=1:
  30.                 sys.stderr.write('several solutions???\n')
  31.                 exit(1)
  32.  
  33.         sol = s[0]
  34.  
  35.         uexpr = sol[u]
  36.         h = uexpr/v
  37.  
  38.         print(h)
  39.  
  40.         breakpoint()
  41.  
  42.  
  43.  
  44.  
  45.  
  46.  
  47. if __name__=="__main__":
  48.         main()
  49.