Tracé polaire d'un diagramme de rayonnement (python)

import numpy as np

from scipy.integrate import odeint

import matplotlib.pyplot as plt

from math import *

thetaListe=[]

EthetaListe=[]

ErListe=[]

sigma=5 #Conductivité

I=1 #(en A) : Intensité de l'antenne à fixer

r=0.2 #(en m) : distance entre l'antenne et le requin

d=0.1 #taille de l'antenne (en m)

theta=0

Er=0

Etheta=0

plt.axes(projection = 'polar')

#Compte tenu l'imparité du sin, le signe de Er et de Etheta importe peu car celui-ci dépend de la direction de mesure. On prend donc la valeur absolue du champ. En effet, nous nous intéressons à la valeur de l'amplitude (une grandeur positive) du champ électrique. Nous tracerons donc indifféremment la valeur absolue des grandeurs adéquates. De fait, le tracé en polaire devra être symétrique par rapport aux valeurs que prennent la variable theta, et ce par symétrie du système émetteur qu'est l'antenne.

for pas in range(361*1):

theta=-180+pas/1

Er=2*cos(theta*pi/180)*I*d/(sigma*2*pi*r**3)

Etheta=sin(theta*pi/180)*I*d/(sigma*2*pi*r**3)

thetaListe.append(theta*pi/180)

ErListe.append(abs(Er))

EthetaListe.append(abs(Etheta))

"""ici polaire"""

EListe=[]

EListe = [sqrt(ErListe[k]**2 + EthetaListe[k]**2) for k in range(len(thetaListe))]

#plt.xlabel("$E=sqrt(Er**2+Etheta**2)$")

plt.plot(thetaListe, ErListe,linewidth=2) #Tracé de Er

plt.plot(thetaListe, EthetaListe,linewidth=2) #Tracé de Etheta

plt.plot(thetaListe, EListe,linewidth=2) #Tracé de E

plt.show()

'''Diagramme polaire'''

#for angle in range (len(thetaListe)):

#