Subversion Repositories pentevo

Rev

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

  1. from sympy import *
  2. import math
  3.  
  4. #w1=184054.140601213
  5. #w2=191540.313016980
  6. #q=1.02654691835174
  7.  
  8. w1=30000*2*math.pi
  9. w2=w1
  10. q=1.0
  11.  
  12.  
  13.  
  14. print('w1={}\nw2={}\nq={}\n\n'.format(w1,w2,q))
  15.  
  16. A=1/w1 + 1/(w2*q)
  17. B=1/(q*w1*w2) + 1/(w2*w2)
  18. C=1/(w1*w2*w2)
  19.  
  20. print('A={}\nB={}\nC={}\n\n'.format(A,B,C))
  21.  
  22. r1 = 1.2e3
  23. r2 = 47e3
  24. r3 = 68e3
  25. c3 = 10e-12
  26.  
  27. r4,c1,c2 = symbols('r4 c1 c2',positive=True)
  28.  
  29. eq1 = Eq( (c1*r1*r2 + c3*r1*r3 + c3*r1*r4 + c3*r2*r3 + c3*r2*r4 + c3*r3*r4)/(r1 + r2), A )
  30. eq2 = Eq( c3*(c1*r1*r2*r3 + c1*r1*r2*r4 + c1*r1*r3*r4 + c2*r1*r3*r4 + c2*r2*r3*r4)/(r1 + r2), B )
  31. eq3 = Eq( c1*c2*c3*r1*r2*r3*r4/(r1 + r2), C )
  32.  
  33. s=solve( [eq1,eq2,eq3], [c1,c2,r4], dict=True )
  34.  
  35. print(s)
  36.  
  37.