1 #!/usr/bin/python
   2 # Uloha MAGN1
   3 # Vytvoril Radek Podgorny <radek@podgorny.cz>
   4 
   5 from math import *
   6 
   7 Eps0 = 8.854187817 * 10**-12 #epsilon
   8 mi0 = 4*pi * 10**-7
   9 
  10 
  11 # Tady zadej svoje parametry
  12 A = 0.022
  13 C = 0.0428
  14 D = 0.031
  15 F = 0.263
  16 G = 0.236
  17 H = 0.0719
  18 N1 = 96
  19 N2 = 18
  20 m1 = 0.000791
  21 m2 = 0.00214
  22 mir = 23900
  23 Kp = 0.83
  24 HL = 0.0241
  25 Bm = 0.8
  26 f = 138
  27 In1 = 0.00022
  28 In2 = 0.000624
  29 Iss = 0.00491
  30 
  31 
  32 RA = (H-m2)/(mi0*mir * A*HL) + m2/(mi0 * A*HL)
  33 RC = (H-m1)/(mi0*mir * C*HL) + m1/(mi0 * C*HL)
  34 RD = H/(mi0*mir * D*HL)
  35 RF = F/(mi0*mir * A*HL)
  36 RG = G/(mi0*mir * A*HL)
  37 
  38 Rleft = 2*RG+RD
  39 Rright = 2*RF+RA
  40 R = RC + (Rleft*Rright)/(Rleft+Rright)
  41 L1 = N1**2/R
  42 
  43 print("Vlastni indukcnost L1 je " + str(L1))
  44 
  45 
  46 ### KURVADRAT! Celou dobu to asi fungovalo a ja delal tolik vifikundaci! Vracet uz to nebudu...
  47 #Rright = 2*RF + RA
  48 #R = RD + 2*RG + (RC*Rright)/(RC+Rright)
  49 #L2 = N2**2/R
  50 #M = (L1*L2)**0.5
  51 
  52 ### Tohle je vykradeny z nejakyho prikladu ze skript
  53 R1 = RC
  54 R2 = 2*RF + RA
  55 R3 = 2*RG + RD
  56 M = (N1*N2 * R2)/(R1*R2 + R1*R3 + R2*R3)
  57 
  58 print("Vzajemna indukcnost M je " + str(M))
  59 
  60 
  61 ### Tohle je zatim blbe...
  62 #Rright = 2*RF + RA
  63 #R = RD + 2*RG + (RC*Rright)/(RC+Rright)
  64 I = Bm*R/N1
  65 
  66 print("Maximalni I1 je " + str(I))
  67 
  68 #Spravny vysledky pro zbytek:
  69 #2.42049
  70 #0.00127938