Subversion Repositories pentevo

Rev

Rev 1247 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed

  1. from sympy import *
  2. import math
  3.  
  4. #r1 = 1.2e3
  5. #r2 = 47e3
  6. #r3 = 68e3
  7. #c3 = 10e-12
  8.  
  9. r1,r2,r3,r4,c1,c2,c3 = symbols('r1,r2,r3,r4,c1,c2,c3')
  10.  
  11. A,B,C,D,E,F,G,H,K,L,M = symbols('A,B,C,D,E,F,G,H,K,L,M')
  12.  
  13. eq1 = Eq( c1*E+r4*F, D )
  14. eq2 = Eq( c1*H + c1*r4*K + c2*r4*L, G )
  15. eq3 = Eq( c1*c2*r4, M )
  16.  
  17. s=solve( [eq1,eq2,eq3], [r4,c1,c2], dict=True, manual=True, simplify=False )
  18.  
  19. print(s)
  20. print('')
  21.  
  22. for ss in s:
  23.     print('\n===')
  24.     print('r4={}\n'.format(ss[r4]))
  25.     print('c1={}\n'.format(ss[c1]))
  26.     print('c2={}\n'.format(ss[c2]))
  27.  
  28.  
  29.  
  30. D = A*(r1+r2)-c3*c1*c2
  31. E = r1*r2
  32. F = c3*(r1+r2+r3)
  33. G = B*(r1+r2)/c3
  34. H = r1*r2*r3
  35. K = r1*(r2+r3)
  36. L = r3*(r1+r2)
  37. M = C*(r1+r2)/(c3*r1*r2*r3)
  38.  
  39.  
  40. A  = (c1*r1*r2 + c3*r1*r3 + c3*r1*r4 + c3*r2*r3 + c3*r2*r4 + c3*r3*r4)/(r1 + r2)
  41. B  = c3*(c1*r1*r2*r3 + c1*r1*r2*r4 + c1*r1*r3*r4 + c2*r1*r3*r4 + c2*r2*r3*r4)/(r1 + r2)
  42. C  = c1*c2*c3*r1*r2*r3*r4/(r1 + r2)
  43.  
  44.