1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
#Paramétrisation des espaces de départ et d'arrivée
s1 = np.linspace(-5,5,N, False)
s2 = np.linspace(-5,5,N, False)
S1,S2 = np.meshgrid(s1,s2) #points de discrétisation de l'espace de départ
t1 = np.linspace(-5,5,N, False)
t2 = np.linspace(-5,5,N, False)
T1,T2 = np.meshgrid(t1,t2) #points de discrétisation de l'espace de d'arrivé
#Définition des mesures
def n1(x1,x2):
if np.exp(-(x1**2+x2**2)/2) ==0:
return 10**-5
else:
return np.exp(-(x1**2+x2**2)/2)/2*np.pi
for i in range(N):
for j in range(N):
mu[i][j] = n1(s1[i],s2[j])
mu = mu/np.sum(mu)
def n2(x1,x2):
if np.exp(-(x1**2+x2**2)/2) ==0:
return 10**-5
else:
return np.exp(-(x1**2+x2**2)/4)/4*np.pi
for i in range(N):
for j in range(N):
nu[i][j] = n2(t1[i],t2[j])
print(nu)
nu = nu/np.sum(nu) |
