Discussion :

[CechD]

Bonjour,

J'ai des bases en python mais pas suffisante pour faire ce que je souhaite faire, ce serait super si vous pouviez me donner un coup de main.

1) J'aimerais afficher en rouge les points (i,j) par exemple (0.2,0.6) et en bleu des points (k,l) par exemple (0.1,0.3) et les lier par créer une fonction qui relie le point (i,j) au point (k,l) par un trait.
2) De plus j'aimerais afficher les points bleues et rouges sur le même graphique, pour qu'on y voit quelque chose j'aimerais séparer le nuage de point bleue et rouge. Je vous met une image pour que vous voyez ce que je veux dire :



Je vous remercie sincèrement pour toute aide.

[CechD]

Je peux vous donner mes matrices (mu et nu) si vous le souhaitez.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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)