Skip to content

Latest commit

 

History

History
286 lines (150 loc) · 6.85 KB

README.md

File metadata and controls

286 lines (150 loc) · 6.85 KB
Raw

EMpython - documentation not finished

This python lib can be used freely but it's build for solving problems in course 31035 "Anvendt elektromagnetisme" DTU

31035 Anvendt elektromagnetisme

Literature

Install python 3 and pip

sudo apt install python3 python3-pip

Install numpy

pip3 install numpy


  1. printJSON.py

    Print material table and constant table from Ulaby appendix B and FUNDAMENTAL PHYSICAL CONSTANTS

  • constants()

    Prints constant table

  • meterials()

    Prints material table

  • getMu_0()

    Returns mu zero

  • getEpsilon_0()

    Returns epsilon zero

  • getSpeedOfLight()

    Returns speed of light


  1. prefix.py

    Returns magnitude for easy convertion to meter.

    • E_exa()
    • P_peta()
    • T_tera()
    • G_giga()
    • M_mega()
    • k_kilo()
    • m_milli()
    • mu_micro()
    • n_nano()
    • p_pico()
    • f_femto()
    • a_atto()

  1. functions.py

    Helpful functions

  • Transmision line

    1. Find characteristic impedance Z_0 with R', G', L', C' and omega. Ulaby s. 89

      • General case

        TL_find_Z0_RpGpLpCpOmega(Rp, Gp, Lp, Cp, omega)

    2. Find u_p Phase velocity with omega and beta. Ulaby s. 89

      • General case

        TL_find_Up_OmegaBeta(omega, beta)

    3. Find Propagation constant (gamma = alpha + j*beta) with R', G', L', C' and omega. Ulaby s. 89

      • General case

        TL_find_smallGamma_RpGpLpCp(Rp, Gp, Lp, Cp, omega)

    4. Find beta with omega and epsilon. Ulaby s. 89 and omega. Ulaby s. 89

      • Lossless case

        TL_find_beta_omegaEpsilon(omega, epsilon_r)

    5. Find u_p Phase velocity with epsilon. Ulaby s. 89

      • Lossless case

        TL_find_Up_epsilon(epsilon_r)

    6. Find characteristic impedance Z_0 with L' and C'. Ulaby s. 89

      • Lossless case

        TL_find_Z0_LpCp(Lp, Cp)

    7. Find characteristic impedance Z_0 with epsilon, a (inner radius) and b (outer radius). Ulaby s. 89

      • Lossless coaxial

        TL_find_Z0_epsilonAB(epsilon_r, a, b)

    8. Find R's with frequency, mu of conductor and sigma of conductor. Ulaby s. 76

      • General case

        TL_find_Rs_freqMucSigmac(freq, mu_c, sigma_c)

    9. Find R' with R's, a and b. Ulaby s. 76

      • General coaxial case

        TL_find_Rprime_RsAB(Rs, a, b)

    10. Find L' with mu (insulator), a and b. Ulaby s. 76

      • General coaxial case

        TL_find_Lprime_muAB(mu_r, a, b)

    11. Find G' with sigma (insulator), a and b. Ulaby s. 76

      • General coaxial case

        TL_find_Gprime_sigmaAB(sigma, a, b)

    12. Find C' with epsilon (insulator), a and b. Ulaby s. 76

      • General coaxial case

        TL_find_Cprime_epsilonAB(epsilon_r, a, b)

    13. Find standing wave ratio (SWR) with reflection coefficient. Ulaby s. 94

      • General case

        TL_find_SWR_GAMMA(GAMMA)

    14. Find standing wave ratio (SWR) with with V_min and V_max. Ulaby s. 94

      • General case

        TL_find_SWR_VminVmax(V_min, V_max)

    15. Find reflection coefficient (GAMMA) abs value with only SWR. Ulaby s. 94

      • General case

        TL_find_absGAMMA_SWR(SWR)

    16. Find normalized load impedance (z_L) with GAMMA. Ulaby s. 90

      • General case

        TL_find_zL_GAMMA(GAMMA)

    17. Find characteristic impedance with short- and opencircuit impedance. Ulaby s. 103

      • Lossless case

        TL_find_Z0_ZscZoc(Zsc, Zoc)

    18. Find beta with length of TL and short- and opencircuit impedance. Ulaby s. 103.

      • Lossless case

        TL_find_beta_lenZscZoc(lenght, Zsc, Zoc)

  • Plane-wave propagation

    1. Find betahat vector with E_field and H_field. Ulaby s. 369

      • Plane-wave case !!Only returns beta hat vector!!

        find_betaHatVec_EfieldHfield(E_field, H_field)

    2. Find reflection coefficient (big Gamma) with Z_L and Z_0. Ulaby s. 90

      • Plane-wave case

        find_GAMMA_ZLZ0(Z_L, Z_0)

    3. Find Z_L with reflection coefficient and Z_0. Ulaby s. 90 If GAMMA is on polar form, pol needs to be True and GAMMA is magnitude and angle is angle. deg needs to be True if angle is in degrees.

      • Plane-wave case

        find_ZL_GAMMAZ0(GAMMA, Z_0, pol=False, deg=False, angle=0)

  • Wave reflection and transmission

    1. Find reflection coefficient (GAMMA) with complex intrinsic impandance. Ulaby s. 379

      • Normal incidence

        WRT_find_GAMMA_eta1eta2(eta1, eta2)

    2. Find reflection coefficient (GAMMA) with intrinsic impandance. Ulaby s. 379

      • Normal incidence

        WRT_find_GAMMA_eta1eta2(eta1, eta2)

    3. Find complex intrinsic impandance with mu_r, epsilon_r, sigma and omega. Ulaby s. 379

      • Normal incidence !OBS! not tested

        WRT_find_eta_c_murEprSig(mu_r, epsilon_r, sigma, omega)

    4. Find transmitting angle with n1, n2 (index of refraction) and incident angle. If deg=False radians is used. Ulaby s. 385

      • Snell's law

        WRT_find_Tangle_n1n2Angle(n1, n2, angle)

  • Wave polarization

    1. Find E wave polarization with E0 vector. Vector can also be complex. Lection 12 (9, 10, 11)

      • Plane wave

        POL_find_polarization(E0_vector)

  • Wave propagation

    1. Find intrinsic impandance with mu_r, epsilon_r. Ulaby s. 357

      • Lossless case

        WP_find_eta_murEpr(mu_r, epsilon_r)

  • getPolarization(E_0real, E_0imag)

    Input

    E_0real and E_0imag are real vectors

    Should also work with H_0real and H_0imag as well

    Output

    Returns polarization

    Else returns bool False

  • getRightOrLeftPolE_field(E_0real, E_0imag, betaVec)

    Input

    E_0real, E_0imag and betaVec are real vectors

    Only works with Electric field

    Output

    Returns left or right polarization

    Else returns error message

  • getRightOrLeftPolH_field(H_0real, H_0imag, betaVec)

    Input

    H_0real, H_0imag and betaVec are real vectors

    Only works with Magnetic field

    Output

    Returns left or right polarization

    Else returns error message