Graph theoretical BER for 2 users

Comparaison_Modele_Theorie.py

@@ -106,5 +106,20 @@ def compDiffCarres(loop_range=5, it=1000):
         print(diff)
 
 
+def plotTh2users(P, sigma, g):
+    probas = []
+    puissances = np.array(P)
+    for p in puissances:
+        probas.append(NOMATh.theorie(g, sigma, p, 2)[0])
+
+    plt.plot(puissances[:, -1], probas)
+    plt.show()
+
+
 if __name__ == '__main__':
-    compRepartitionPuissance()
+    # compRepartitionPuissance()
+    P1 = np.arange(0.88, 0.999, 0.001)
+    P = [(1 - p, p) for p in P1]
+    sigma = 0.075
+    g = 1
+    plotTh2users(P, sigma, g)