cleaning
This commit is contained in:
346
QAM/DDPLL/qam.c
346
QAM/DDPLL/qam.c
@ -1,346 +0,0 @@
|
||||
#include <math.h>
|
||||
#include <stdint.h>
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <complex.h>
|
||||
#include <string.h>
|
||||
|
||||
#define A 10
|
||||
|
||||
struct qam_system_s {
|
||||
int M; // Nombre de symboles M-QAM
|
||||
int k; // Nombre de bits/symboles
|
||||
double Fs; // Fréquence d'échantillionage
|
||||
double Ts; // Temps d'échantillionage
|
||||
int N; // Nombre d'échantillions
|
||||
double Fc; // Fréquence de la porteuse
|
||||
double complex** constellation; // Tableau de symboles I + j Q
|
||||
};
|
||||
typedef struct qam_system_s qam_system;
|
||||
|
||||
// Initialisation de la constellation (double tableau de taille sqrt(M)),
|
||||
// ToDo : changer à un tableau à 1 dimension pour éviter de calculer sqrt(M)
|
||||
void init_constellation (qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
qam->constellation = (double complex**)malloc(sizeof(double complex*) * sm);
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
qam->constellation[i] = (double complex*)malloc(sizeof(double complex) * sm);
|
||||
}
|
||||
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0); // Pour puissance unitaire
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
double complex ip = -(sm - 1) + 2 * i;
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex qp = -(sm - 1) + 2 * j;
|
||||
qam->constellation[i][j] = (ip + I * qp) / norm_factor;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex iq = symbols[k];
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
int idx = k * qam->N + n;
|
||||
s[idx] = A * iq * cexp(2 * I * M_PI * qam->Fc * ((double)idx / qam->Fs));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM
|
||||
void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat, FILE *fp_constel) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)(k * qam->N + n) / qam->Fs)) / A;
|
||||
}
|
||||
r /= qam->N;
|
||||
|
||||
if (fp_constel) {
|
||||
fprintf(fp_constel, "% .8f % .8f\n", creal(r), cimag(r));
|
||||
fflush(fp_constel);
|
||||
}
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation (lent)
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl = 0;
|
||||
int j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// index du symbole (id) : même mappage que dans bits_to_symbols()
|
||||
int id = i_cl * sm + j_cl;
|
||||
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + b] = (id >> (qam->k - 1 - b)) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// PLL pour corriger le déphasage
|
||||
void pll_qam_symbol(qam_system* qam, double complex* symbols_rx, double complex* r_corr, int nb_symbols, double Kp, double Ki, double alpha, FILE* fp_error) {
|
||||
double phase_est = 0.0;
|
||||
double integrator = 0.0;
|
||||
double filtered_error = 0.0;
|
||||
|
||||
int sm = (int)sqrt(qam->M);
|
||||
int N = qam->N;
|
||||
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r_symbol = 0;
|
||||
for (int n = 0; n < N; n++) {
|
||||
int idx = k * N + n;
|
||||
r_symbol += symbols_rx[idx] * cexp(-2.0 * I * M_PI * qam->Fc * ((double)idx / qam->Fs));
|
||||
}
|
||||
r_symbol /= N;
|
||||
|
||||
r_symbol *= cexp(-I * phase_est);
|
||||
|
||||
double min_d = INFINITY;
|
||||
double complex closest = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r_symbol - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
closest = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
double error = carg(r_symbol * conj(closest));
|
||||
|
||||
filtered_error = (1.0 - alpha) * filtered_error + alpha * error;
|
||||
integrator += Ki * filtered_error;
|
||||
phase_est += Kp * filtered_error + integrator;
|
||||
|
||||
// Écriture de l'erreur PLL dans le fichier
|
||||
if (fp_error) {
|
||||
fprintf(fp_error, "%d % .8f\n", k, 100 * filtered_error);
|
||||
fflush(fp_error);
|
||||
}
|
||||
|
||||
for (int n = 0; n < N; n++) {
|
||||
int idx = k * N + n;
|
||||
r_corr[idx] = symbols_rx[idx] * cexp(-I * phase_est);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
double compare_bits(uint8_t* bits1, uint8_t* bits2, int nb_bits) {
|
||||
int errors = 0;
|
||||
for (int i = 0; i < nb_bits; i++) {
|
||||
if (bits1[i] != bits2[i]) errors++;
|
||||
}
|
||||
return (double)errors / nb_bits;
|
||||
}
|
||||
|
||||
// Minimise le BER (si la pll s'est lockée de maniere déphasée de k*pi/2)
|
||||
void demodulate2(qam_system* qam, double complex* r_corr, int nb_symbols, uint8_t* input_bits, uint8_t* output_bits, FILE* fp_constel) {
|
||||
|
||||
int nb_bits = nb_symbols * qam->k;
|
||||
double best_ber = INFINITY;
|
||||
double best_angle = 0.0;
|
||||
uint8_t* temp_bits = (uint8_t*)malloc(nb_bits * sizeof(uint8_t));
|
||||
|
||||
for (int r = 0; r < 4; r++) {
|
||||
double angle = r * M_PI/2;
|
||||
|
||||
double complex* rotated = (double complex*)malloc(sizeof(double complex) * nb_symbols * qam->N);
|
||||
for (int i = 0; i < nb_symbols * qam->N; i++) {
|
||||
rotated[i] = r_corr[i] * cexp(I * angle);
|
||||
}
|
||||
|
||||
demodulate(qam, rotated, nb_symbols, temp_bits, fp_constel);
|
||||
|
||||
double ber = compare_bits(input_bits, temp_bits, nb_bits);
|
||||
|
||||
if (ber < best_ber) {
|
||||
best_ber = ber;
|
||||
best_angle = angle;
|
||||
}
|
||||
|
||||
free(rotated);
|
||||
}
|
||||
|
||||
for (int i = 0; i < nb_symbols * qam->N; i++) {
|
||||
r_corr[i] *= cexp(I * best_angle);
|
||||
}
|
||||
|
||||
demodulate(qam, r_corr, nb_symbols, output_bits, fp_constel);
|
||||
|
||||
free(temp_bits);
|
||||
}
|
||||
|
||||
|
||||
int main () {
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
//qam.Ts = 0.0003;
|
||||
//qam.N = (int)qam.Fs * qam.Ts;
|
||||
qam.Ts = 0.01;
|
||||
qam.N = (int)(qam.Fs * qam.Ts);
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
//int nb_bits = 1000;
|
||||
//int nb_symbols = nb_bits / qam.k;
|
||||
|
||||
//uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
//for (int i = 0; i < nb_bits; i++) {
|
||||
// input_bits[i] = rand() % 2;
|
||||
//}
|
||||
char* texte = "Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux";
|
||||
int nb_chars = strlen(texte);
|
||||
int nb_bits = nb_chars * 8;
|
||||
int nb_symbols = (nb_bits + qam.k - 1) / qam.k;
|
||||
|
||||
// Conversion du texte en bits
|
||||
uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
for(int b = 0; b < 8; b++){
|
||||
input_bits[i*8 + b] = (texte[i] >> (7-b)) & 1;
|
||||
}
|
||||
}
|
||||
|
||||
// Conversion en symboles
|
||||
double complex* symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
bits_to_symbols(&qam, input_bits, nb_bits, symbols);
|
||||
|
||||
// Modulation
|
||||
int total_samples = qam.N * nb_symbols;
|
||||
double complex* s = malloc(sizeof(double complex) * total_samples);
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
// Ajout du bruit
|
||||
double signal_power = (2.0/3.0)*(qam.M-1); // puissance moyenne
|
||||
double snr_dB = 5; // SNR en dB
|
||||
double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
double sigma = sqrt(signal_power / snr_lin);
|
||||
add_noise(s, total_samples, 10);
|
||||
|
||||
|
||||
FILE *fp_ref = fopen("constellation_ref.dat", "w");
|
||||
int sm = (int)sqrt(qam.M);
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
fprintf(fp_ref, "% .8f % .8f\n", creal(qam.constellation[i][j]), cimag(qam.constellation[i][j]));
|
||||
}
|
||||
}
|
||||
fclose(fp_ref);
|
||||
FILE *fp_constel = fopen("constellation.dat", "w");
|
||||
|
||||
// Ajout de dephasage
|
||||
//double phase_offset = M_PI / 6.0; // 30 degrés
|
||||
//for (int i = 0; i < total_samples; i++) {
|
||||
// s[i] *= cexp(I * phase_offset);
|
||||
//}
|
||||
|
||||
// AJout de decalage de fréquence
|
||||
double freq_offset = 1; // Hz de décalage
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
double t = (double)i / qam.Fs;
|
||||
s[i] *= cexp(I * 2 * M_PI * freq_offset * t);
|
||||
}
|
||||
|
||||
// Ajout de decalage entre les symbole
|
||||
//int offset_samples = (int)(0.3 * qam.N); // décalage de 30% d’un symbole
|
||||
//memmove(s + offset_samples, s, (total_samples - offset_samples) * sizeof(double complex));
|
||||
|
||||
double complex* r_corr = malloc(sizeof(double complex) * total_samples);
|
||||
double Kp = 0.2;
|
||||
double Ki = 0.02;
|
||||
double alpha = 0.1;
|
||||
FILE* fp_error = fopen("pll_error.dat", "w");
|
||||
pll_qam_symbol(&qam, s, r_corr, nb_symbols, Kp, Ki, alpha, fp_error);
|
||||
fclose(fp_error);
|
||||
|
||||
// Démodulation
|
||||
uint8_t* output_bits = (uint8_t*)malloc(nb_bits * sizeof(uint8_t));
|
||||
demodulate2(&qam, r_corr, nb_symbols, input_bits, output_bits, fp_constel);
|
||||
//demodulate(&qam, r_corr, nb_symbols, output_bits, fp_constel);
|
||||
|
||||
fclose(fp_constel);
|
||||
|
||||
// Reconstruction du texte
|
||||
char* texte_recup = malloc(nb_chars + 1);
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
char c = 0;
|
||||
for(int b = 0; b < 8; b++){
|
||||
c |= output_bits[i*8 + b] << (7-b);
|
||||
}
|
||||
texte_recup[i] = c;
|
||||
}
|
||||
texte_recup[nb_chars] = '\0';
|
||||
printf("Texte original : %s\n\n", texte);
|
||||
printf("Texte demodulé : %s\n", texte_recup);
|
||||
|
||||
// Calcul du BER
|
||||
double ber = compare_bits(input_bits, output_bits, nb_bits);
|
||||
printf("Taux d'erreur blind QAM: %.4f\n", ber * 100);
|
||||
|
||||
// Libération mémoire
|
||||
free(input_bits);
|
||||
free(output_bits);
|
||||
free(symbols);
|
||||
free(r_corr);
|
||||
free(s);
|
||||
free_constellation(&qam);
|
||||
|
||||
return 0;
|
||||
}
|
||||
File diff suppressed because it is too large
Load Diff
@ -1,16 +0,0 @@
|
||||
-1.34164079 -1.34164079
|
||||
-1.34164079 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 -1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -1.34164079
|
||||
0.44721360 -0.44721360
|
||||
0.44721360 0.44721360
|
||||
0.44721360 1.34164079
|
||||
1.34164079 -1.34164079
|
||||
1.34164079 -0.44721360
|
||||
1.34164079 0.44721360
|
||||
1.34164079 1.34164079
|
||||
@ -1,32 +0,0 @@
|
||||
import numpy as np
|
||||
import matplotlib.pyplot as plt
|
||||
|
||||
def plot_constellations(ref_file, rx_file, title="Constellation comparison"):
|
||||
# Charger et forcer 2D
|
||||
ref_data = np.atleast_2d(np.loadtxt(ref_file))
|
||||
rx_data = np.atleast_2d(np.loadtxt(rx_file))
|
||||
|
||||
x_ref, y_ref = ref_data[:,0], ref_data[:,1]
|
||||
x_rx, y_rx = rx_data[:,0], rx_data[:,1]
|
||||
|
||||
plt.figure(figsize=(6,6))
|
||||
plt.scatter(x_ref, y_ref, color='blue', s=50, marker='o', label='Référence')
|
||||
plt.scatter(x_rx, y_rx, color='red', s=50, marker='x', label='Reçu')
|
||||
|
||||
# Ajustement automatique des limites
|
||||
all_x = np.concatenate([x_ref, x_rx])
|
||||
all_y = np.concatenate([y_ref, y_rx])
|
||||
margin = 0.1 * max(np.ptp(all_x), np.ptp(all_y))
|
||||
plt.xlim(min(all_x)-margin, max(all_x)+margin)
|
||||
plt.ylim(min(all_y)-margin, max(all_y)+margin)
|
||||
|
||||
plt.xlabel('In-phase (I)')
|
||||
plt.ylabel('Quadrature (Q)')
|
||||
plt.title(title)
|
||||
plt.grid(False)
|
||||
plt.gca().set_aspect('equal', adjustable='box')
|
||||
plt.legend()
|
||||
plt.show()
|
||||
|
||||
plot_constellations("constellation_ref.dat", "constellation.dat", title="Constellation QAM")
|
||||
|
||||
@ -1,68 +0,0 @@
|
||||
import numpy as np
|
||||
import matplotlib.pyplot as plt
|
||||
from watchdog.observers import Observer
|
||||
from watchdog.events import FileSystemEventHandler
|
||||
import time
|
||||
import os
|
||||
|
||||
needs_update = True # flag global
|
||||
|
||||
def plot_constellations(ref_file, rx_file, title="Constellation Comparison", save_path=None):
|
||||
ref_data = np.atleast_2d(np.loadtxt(ref_file))
|
||||
rx_data = np.atleast_2d(np.loadtxt(rx_file))
|
||||
|
||||
x_ref, y_ref = ref_data[:,0], ref_data[:,1]
|
||||
x_rx, y_rx = rx_data[:,0], rx_data[:,1]
|
||||
|
||||
plt.clf() # efface la figure précédente
|
||||
plt.scatter(x_ref, y_ref, color='dodgerblue', s=80, marker='o', edgecolors='k', label='Référence')
|
||||
plt.scatter(x_rx, y_rx, color='tomato', s=80, marker='x', alpha=0.6, label='Reçu')
|
||||
|
||||
all_x = np.concatenate([x_ref, x_rx])
|
||||
all_y = np.concatenate([y_ref, y_rx])
|
||||
margin = 0.15 * max(np.ptp(all_x), np.ptp(all_y))
|
||||
plt.xlim(min(all_x)-margin, max(all_x)+margin)
|
||||
plt.ylim(min(all_y)-margin, max(all_y)+margin)
|
||||
|
||||
plt.grid(True, which='both', linestyle='--', linewidth=0.5, alpha=0.7)
|
||||
plt.axhline(0, color='black', linewidth=1)
|
||||
plt.axvline(0, color='black', linewidth=1)
|
||||
plt.xlabel('In-phase (I)', fontsize=12)
|
||||
plt.ylabel('Quadrature (Q)', fontsize=12)
|
||||
plt.title(title, fontsize=14, fontweight='bold')
|
||||
plt.gca().set_aspect('equal', adjustable='box')
|
||||
plt.legend()
|
||||
if save_path:
|
||||
plt.savefig(save_path, dpi=300, bbox_inches='tight')
|
||||
plt.pause(0.1)
|
||||
|
||||
class FileChangeHandler(FileSystemEventHandler):
|
||||
def on_modified(self, event):
|
||||
global needs_update
|
||||
if event.src_path.endswith("constellation_ref.dat") or event.src_path.endswith("constellation.dat"):
|
||||
print(f"{event.src_path} modifié")
|
||||
needs_update = True # on ne fait que signaler
|
||||
|
||||
if __name__ == "__main__":
|
||||
ref_file = "constellation_ref.dat"
|
||||
rx_file = "constellation.dat"
|
||||
|
||||
plt.ion()
|
||||
plt.figure(figsize=(7,7))
|
||||
plot_constellations(ref_file, rx_file)
|
||||
|
||||
event_handler = FileChangeHandler()
|
||||
observer = Observer()
|
||||
observer.schedule(event_handler, path=os.path.dirname(os.path.abspath(ref_file)) or '.', recursive=False)
|
||||
observer.start()
|
||||
|
||||
try:
|
||||
while True:
|
||||
if needs_update:
|
||||
plot_constellations(ref_file, rx_file)
|
||||
needs_update = False
|
||||
time.sleep(0.2) # boucle principale
|
||||
except KeyboardInterrupt:
|
||||
observer.stop()
|
||||
observer.join()
|
||||
|
||||
File diff suppressed because it is too large
Load Diff
@ -1,16 +0,0 @@
|
||||
-3.00000000 -3.00000000
|
||||
-3.00000000 -1.00000000
|
||||
-3.00000000 1.00000000
|
||||
-3.00000000 3.00000000
|
||||
-1.00000000 -3.00000000
|
||||
-1.00000000 -1.00000000
|
||||
-1.00000000 1.00000000
|
||||
-1.00000000 3.00000000
|
||||
1.00000000 -3.00000000
|
||||
1.00000000 -1.00000000
|
||||
1.00000000 1.00000000
|
||||
1.00000000 3.00000000
|
||||
3.00000000 -3.00000000
|
||||
3.00000000 -1.00000000
|
||||
3.00000000 1.00000000
|
||||
3.00000000 3.00000000
|
||||
@ -1,16 +0,0 @@
|
||||
set terminal qt size 600,600 title "Diagramme de constellation"
|
||||
set xlabel "I (In-phase)"
|
||||
set ylabel "Q (Quadrature)"
|
||||
set grid
|
||||
set key off
|
||||
set pointsize 1.5
|
||||
set xrange [-5:5]
|
||||
set yrange [-5:5]
|
||||
|
||||
while (1) {
|
||||
plot \
|
||||
'constellation_ref.dat' using 1:2 with points pt 7 ps 2.5 lc rgb "red", \
|
||||
'constellation.dat' using 1:2 with points pt 7 ps 1 lc rgb "blue"
|
||||
pause 0.15
|
||||
}
|
||||
|
||||
@ -1,226 +0,0 @@
|
||||
#include <math.h>
|
||||
#include <stdint.h>
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <complex.h>
|
||||
#include <string.h>
|
||||
|
||||
#define A 1
|
||||
|
||||
struct qam_system_s {
|
||||
int M; // Nombre de symboles M-QAM
|
||||
int k; // Nombre de bits/symboles
|
||||
double Fs; // Fréquence d'échantillionage
|
||||
double Ts; // Temps d'échantillionage
|
||||
int N; // Nombre d'échantillions
|
||||
double Fc; // Fréquence de la porteuse
|
||||
double complex** constellation; // Tableau de symboles I + j Q
|
||||
};
|
||||
typedef struct qam_system_s qam_system;
|
||||
|
||||
// Initialisation de la constellation (double tableau de taille sqrt(M)),
|
||||
// ToDo : changer à un tableau à 1 dimension pour éviter de calculer sqrt(M)
|
||||
void init_constellation (qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
qam->constellation = (double complex**)malloc(sizeof(double complex*) * sm);
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
qam->constellation[i] = (double complex*)malloc(sizeof(double complex) * sm);
|
||||
}
|
||||
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0); // Pour puissance unitaire
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
double complex ip = -(sm - 1) + 2 * i;
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex qp = -(sm - 1) + 2 * j;
|
||||
qam->constellation[i][j] = (ip + I * qp);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex iq = symbols[k];
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
s[k * qam->N + n] = iq * cexp(2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM
|
||||
void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat, FILE *fp_constel) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
r /= qam->N;
|
||||
|
||||
if (fp_constel) {
|
||||
fprintf(fp_constel, "% .8f % .8f\n", creal(r), cimag(r));
|
||||
fflush(fp_constel);
|
||||
}
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation (lent)
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// index du symbole (id) : même mappage que dans bits_to_symbols()
|
||||
int id = i_cl * sm + j_cl;
|
||||
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + b] = (id >> (qam->k - 1 - b)) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
double compare_bits(uint8_t* bits1, uint8_t* bits2, int nb_bits) {
|
||||
int errors = 0;
|
||||
for (int i = 0; i < nb_bits; i++) {
|
||||
if (bits1[i] != bits2[i]) errors++;
|
||||
}
|
||||
return (double)errors / nb_bits;
|
||||
}
|
||||
|
||||
int main () {
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
qam.Ts = 0.0003;
|
||||
qam.N = (int)qam.Fs * qam.Ts;
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
//int nb_bits = 1000;
|
||||
//int nb_symbols = nb_bits / qam.k;
|
||||
|
||||
//uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
//for (int i = 0; i < nb_bits; i++) {
|
||||
// input_bits[i] = rand() % 2;
|
||||
//}
|
||||
char* texte = "Vif juge, trempez ce blond whisky aqueux";
|
||||
int nb_chars = strlen(texte);
|
||||
int nb_bits = nb_chars * 8;
|
||||
int nb_symbols = (nb_bits + qam.k - 1) / qam.k;
|
||||
|
||||
// Conversion du texte en bits
|
||||
uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
for(int b = 0; b < 8; b++){
|
||||
input_bits[i*8 + b] = (texte[i] >> (7-b)) & 1;
|
||||
}
|
||||
}
|
||||
|
||||
// Conversion en symboles
|
||||
double complex* symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
bits_to_symbols(&qam, input_bits, nb_bits, symbols);
|
||||
|
||||
// Modulation
|
||||
int total_samples = qam.N * nb_symbols;
|
||||
double complex* s = malloc(sizeof(double complex) * total_samples);
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
// Ajout du bruit
|
||||
double signal_power = (2.0/3.0)*(qam.M-1); // puissance moyenne
|
||||
double snr_dB = 5; // SNR en dB
|
||||
double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
double sigma = sqrt(signal_power / snr_lin);
|
||||
add_noise(s, total_samples, sigma);
|
||||
|
||||
|
||||
FILE *fp_ref = fopen("constellation_ref.dat", "w");
|
||||
int sm = (int)sqrt(qam.M);
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
fprintf(fp_ref, "% .8f % .8f\n", creal(qam.constellation[i][j]), cimag(qam.constellation[i][j]));
|
||||
}
|
||||
}
|
||||
fclose(fp_ref);
|
||||
FILE *fp_constel = fopen("constellation.dat", "w");
|
||||
|
||||
// Démodulation
|
||||
uint8_t* output_bits = (uint8_t*)malloc(nb_bits * sizeof(uint8_t));
|
||||
demodulate(&qam, s, nb_symbols, output_bits, fp_constel);
|
||||
|
||||
fclose(fp_constel);
|
||||
|
||||
// Reconstruction du texte
|
||||
char* texte_recup = malloc(nb_chars + 1);
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
char c = 0;
|
||||
for(int b = 0; b < 8; b++){
|
||||
c |= output_bits[i*8 + b] << (7-b);
|
||||
}
|
||||
texte_recup[i] = c;
|
||||
}
|
||||
texte_recup[nb_chars] = '\0';
|
||||
printf("Texte original : %s\n", texte);
|
||||
printf("Texte demodulé : %s\n", texte_recup);
|
||||
|
||||
// Calcul du BER
|
||||
double ber = compare_bits(input_bits, output_bits, nb_bits);
|
||||
printf("Taux d'erreur blind QAM: %.4f\n", ber * 100);
|
||||
|
||||
// Libération mémoire
|
||||
free(input_bits);
|
||||
free(output_bits);
|
||||
free(symbols);
|
||||
free(s);
|
||||
free_constellation(&qam);
|
||||
|
||||
return 0;
|
||||
}
|
||||
@ -1,880 +0,0 @@
|
||||
-0.31202845 -0.36308093
|
||||
-0.21403073 0.38857646
|
||||
-0.38530118 -0.24484730
|
||||
-0.44021783 0.07284252
|
||||
-0.34501041 0.34872590
|
||||
0.35857538 -0.40045539
|
||||
-0.25176904 0.40377600
|
||||
-0.10377581 0.43578972
|
||||
-1.01896226 0.14271390
|
||||
-1.15718524 -1.07104846
|
||||
-0.40883880 0.30174936
|
||||
0.26396095 0.48082155
|
||||
-0.31654178 1.26105861
|
||||
-0.33255136 -0.17032596
|
||||
-0.41545553 0.38068746
|
||||
-0.23973431 1.32252511
|
||||
-0.49928116 0.40812210
|
||||
-0.60891702 -0.47180571
|
||||
-0.98044772 0.33254717
|
||||
1.23532760 -1.08284043
|
||||
-1.01739273 0.23195386
|
||||
-0.95909631 -1.27747201
|
||||
-0.18344231 0.97172160
|
||||
-0.36005886 -1.19561767
|
||||
-0.38733224 1.03655789
|
||||
-1.14703959 0.49592443
|
||||
-0.34627477 0.37436581
|
||||
-0.32594570 -0.36576269
|
||||
-0.33047672 0.45166025
|
||||
1.03227135 -0.26600616
|
||||
-0.24368430 1.13277766
|
||||
-1.10409800 -1.13670229
|
||||
-0.27024330 0.45950443
|
||||
-0.43012606 -0.49137862
|
||||
-0.42938643 0.96752310
|
||||
0.29460881 0.46255766
|
||||
-1.12981200 0.36271133
|
||||
-1.06546596 -0.99914698
|
||||
-0.32964627 0.54802267
|
||||
-0.97768436 1.05949598
|
||||
-0.40125588 0.35060801
|
||||
-0.14555788 -0.29516627
|
||||
-1.16824437 0.33042499
|
||||
-1.38177121 -1.02828516
|
||||
-0.36030711 0.53502663
|
||||
-0.99270774 0.38492966
|
||||
-0.34180208 0.46901106
|
||||
1.02021653 -1.03189533
|
||||
-0.38784460 0.27839052
|
||||
0.98525087 1.10441090
|
||||
-0.37064093 0.37855671
|
||||
0.92249275 0.38692346
|
||||
-0.37782815 0.26512776
|
||||
-0.38059680 -1.10559455
|
||||
-1.08779693 0.25242918
|
||||
-1.31058915 -0.99926776
|
||||
-0.32423909 1.22001515
|
||||
-0.25335804 1.00836191
|
||||
-0.50573269 0.31900345
|
||||
0.42698256 -0.94413229
|
||||
-0.39166269 0.41072366
|
||||
0.25755208 -0.39660996
|
||||
-0.33022334 1.05502122
|
||||
-1.12361886 0.95444663
|
||||
-0.37428566 0.38306311
|
||||
0.21248197 1.04776058
|
||||
-0.47434424 0.99247118
|
||||
0.43507446 -0.36697818
|
||||
-0.95021817 0.15733583
|
||||
-0.99554159 -1.14130590
|
||||
-0.25435034 0.50150512
|
||||
-1.13525219 -0.27992447
|
||||
-0.54906403 0.94151150
|
||||
-1.17729844 -0.37677613
|
||||
-0.29075006 1.14107781
|
||||
-0.55423390 -0.31265259
|
||||
-0.32086652 0.46332639
|
||||
-0.27004284 -0.25876856
|
||||
-0.56225634 0.99336544
|
||||
-0.32896397 -0.32120490
|
||||
-0.42109267 1.09730148
|
||||
0.16066972 -1.09133678
|
||||
-0.46397415 -0.26636357
|
||||
-0.48496755 0.28388862
|
||||
-0.36968179 0.36567458
|
||||
0.26521387 -0.29984385
|
||||
-0.33480946 0.46485409
|
||||
-0.36463789 0.26204947
|
||||
-1.18759281 0.26015798
|
||||
-0.96596592 -0.93117220
|
||||
-0.33784886 0.39164206
|
||||
0.43769194 0.69739352
|
||||
-0.43032293 0.84413229
|
||||
-0.23774796 -0.28911962
|
||||
-0.33130118 0.41821511
|
||||
-0.50013585 0.97965454
|
||||
-0.34689709 0.35409852
|
||||
-0.23010782 -0.12586032
|
||||
-0.97858139 0.38056045
|
||||
1.08280742 -1.25403467
|
||||
-1.08716964 0.37903094
|
||||
-1.06816096 -1.26485276
|
||||
-0.29749093 1.08244304
|
||||
-0.36919521 -1.24051895
|
||||
-0.41806929 1.02179014
|
||||
-0.87352398 0.31284450
|
||||
-0.32826699 0.34535817
|
||||
-0.52300831 -0.18727837
|
||||
-0.25032100 0.44847978
|
||||
0.91304423 -0.42961637
|
||||
-0.13446642 1.15022408
|
||||
-1.15680622 -1.21068591
|
||||
-0.36015459 0.38355671
|
||||
-0.41281658 -0.28563836
|
||||
-0.28095750 1.17150443
|
||||
0.44841940 0.22784693
|
||||
-1.20991785 0.43888552
|
||||
-1.26949348 -1.09066654
|
||||
-0.49699144 0.41646383
|
||||
-1.10251223 1.15954894
|
||||
-0.35344943 0.42981020
|
||||
-0.35460879 -0.40159958
|
||||
-0.96881406 0.42432504
|
||||
-1.12683149 -0.94337597
|
||||
-0.34090486 0.45009494
|
||||
-1.30640116 0.33027100
|
||||
-0.21513225 0.30060400
|
||||
1.23753084 -1.09859072
|
||||
-0.17778666 0.45746800
|
||||
0.98283097 1.00470073
|
||||
-0.44060208 0.46320487
|
||||
1.21931561 0.58300444
|
||||
-0.34443300 0.35025224
|
||||
-0.20455252 -1.20751826
|
||||
-1.05565774 0.30209405
|
||||
-0.90794063 -0.99564640
|
||||
-0.37062355 1.18006149
|
||||
-0.36830585 1.05842914
|
||||
-0.28829158 0.37278172
|
||||
0.51701024 -1.17362256
|
||||
-0.28685147 0.17704574
|
||||
0.43894598 -0.41798899
|
||||
-0.61389745 1.18724310
|
||||
-1.14959160 1.08955659
|
||||
-0.43773208 0.40961614
|
||||
0.00657211 1.11721014
|
||||
-0.45242309 1.04238246
|
||||
0.28752062 -0.13222957
|
||||
-1.29550622 0.31509320
|
||||
-1.06884344 -1.03394430
|
||||
-0.42241401 0.35919199
|
||||
-1.08379928 -0.57802668
|
||||
-0.41458484 0.83139690
|
||||
-1.05365813 -0.42116548
|
||||
-0.41745035 1.19305593
|
||||
-0.44639650 -0.27576769
|
||||
-0.44621577 0.17412965
|
||||
-0.23004430 -0.32898997
|
||||
-0.36165151 1.20449839
|
||||
-0.55117142 -0.48984181
|
||||
-0.16798699 1.04071241
|
||||
0.31011770 -1.36217467
|
||||
-0.46631118 -0.32516904
|
||||
-0.15805322 0.44809717
|
||||
-0.35860554 0.16913916
|
||||
0.41178134 -0.34124789
|
||||
-0.30074668 0.40714965
|
||||
-0.38017496 0.39131278
|
||||
-0.96741058 0.67253560
|
||||
-1.24569648 -0.72436186
|
||||
-0.27101990 0.38548867
|
||||
0.28420613 0.25890492
|
||||
-0.42278819 1.11012306
|
||||
-0.28250890 -0.32853298
|
||||
-0.28400818 0.27501319
|
||||
-0.17076321 1.11189064
|
||||
-0.45023807 0.37402836
|
||||
-0.22213671 -0.41960172
|
||||
-0.94141884 0.30789514
|
||||
1.05357821 -1.01149325
|
||||
-1.05482103 0.52764631
|
||||
-0.99859625 -0.87049613
|
||||
-0.28261505 1.09640290
|
||||
-0.49015364 -1.08402523
|
||||
-0.14498054 1.09090615
|
||||
-0.94911979 0.42267368
|
||||
-0.34357322 0.44569840
|
||||
-0.34921322 -0.35138845
|
||||
-0.31748827 0.31627969
|
||||
1.11463777 -0.28394406
|
||||
-0.46326050 0.89279233
|
||||
-1.06477639 -1.08022949
|
||||
-0.42836832 0.39212624
|
||||
-0.16028271 -0.37561196
|
||||
-0.45372092 1.09581520
|
||||
-0.04101846 0.41543315
|
||||
-1.24776237 0.50372052
|
||||
-1.01046209 -1.20203071
|
||||
-0.34264731 0.39036274
|
||||
-1.19551204 1.04588443
|
||||
-0.35012639 0.43429116
|
||||
-0.36476204 -0.24929817
|
||||
-1.23662846 0.23121207
|
||||
-1.08869701 -1.09800513
|
||||
-0.35657338 0.32142876
|
||||
-1.16706471 0.40350355
|
||||
-0.47757543 0.33968056
|
||||
1.12643846 -1.03342793
|
||||
-0.40429544 0.51235936
|
||||
0.92993726 1.04386743
|
||||
-0.27921326 0.12895368
|
||||
1.07057147 0.41675044
|
||||
-0.45781543 0.34530738
|
||||
-0.15550368 -1.15290701
|
||||
-1.24474912 0.22373701
|
||||
-1.01636406 -1.15984020
|
||||
-0.45241277 0.94067451
|
||||
-0.39666102 1.01551504
|
||||
-0.30415407 0.52966388
|
||||
0.42242219 -1.04720736
|
||||
-0.48521314 0.34940459
|
||||
0.23494422 -0.47659414
|
||||
-0.39947228 1.05426259
|
||||
-1.18041855 1.15411420
|
||||
-0.24193544 0.39310208
|
||||
0.33862100 1.02349548
|
||||
-0.54773004 1.23503633
|
||||
0.10818630 -0.25533904
|
||||
-0.96224435 0.40891113
|
||||
-0.97643294 -1.01562261
|
||||
-0.42844287 0.45214473
|
||||
-1.00810432 -0.28926920
|
||||
-0.40737767 0.96829636
|
||||
-0.89918110 -0.45258226
|
||||
-0.33496466 1.06184687
|
||||
-0.38077794 -0.20763689
|
||||
-0.40201285 0.36673030
|
||||
-0.45952525 -0.25742302
|
||||
-0.25045162 1.07737031
|
||||
-0.50752572 -0.27866299
|
||||
-0.23083068 1.29397369
|
||||
0.39199302 -1.20171519
|
||||
-0.38433807 -0.37865898
|
||||
-0.55990357 0.41664491
|
||||
-0.29794709 0.51366410
|
||||
0.37475812 -0.43001268
|
||||
-0.27729731 0.29944694
|
||||
-0.47488383 0.35239836
|
||||
-0.96654000 0.69215033
|
||||
-1.11110431 -0.82416572
|
||||
-0.20368948 0.49057982
|
||||
0.52308494 0.33483727
|
||||
-0.24850933 1.07443686
|
||||
-0.34791805 -0.32007730
|
||||
-0.25853676 0.51115807
|
||||
-0.17941734 1.20519231
|
||||
-0.11342097 0.25102777
|
||||
-0.27670160 -0.19012193
|
||||
-1.03994582 0.54908843
|
||||
0.86407132 -1.35562044
|
||||
-0.94185647 0.43008755
|
||||
-1.22186569 -0.88452692
|
||||
-0.32106283 1.09166980
|
||||
-0.53694980 -0.85871688
|
||||
-0.24383048 1.07474394
|
||||
-1.05491406 0.53824228
|
||||
-0.27060332 0.34516662
|
||||
-0.32232377 -0.58146394
|
||||
-0.30486319 0.32528566
|
||||
0.90500008 -0.34391269
|
||||
-0.45103429 1.14989139
|
||||
-0.99391132 -1.09294192
|
||||
-0.38275270 0.33166196
|
||||
-0.35915102 -0.54143680
|
||||
-0.32318896 1.00574008
|
||||
0.40764212 0.45276816
|
||||
-1.11114012 0.15238414
|
||||
-0.92815560 -1.22164919
|
||||
-0.43174187 0.15685803
|
||||
-1.14341183 0.99562192
|
||||
-0.41415278 0.22841552
|
||||
-0.35771155 -0.46436053
|
||||
-1.16822694 0.23892475
|
||||
-0.67346115 -1.18765414
|
||||
-0.42589250 0.14605392
|
||||
-1.09038520 0.17955477
|
||||
-0.49389893 0.24608462
|
||||
1.22544028 -0.75196104
|
||||
-0.41209167 0.21904591
|
||||
0.72389705 1.31839934
|
||||
-0.49661047 0.29730266
|
||||
1.22685896 0.67502585
|
||||
-0.35135181 0.25325888
|
||||
-0.00620830 -1.14363823
|
||||
-1.05333122 0.10254708
|
||||
-0.69867409 -1.13965865
|
||||
-0.53128392 1.07896460
|
||||
-0.45948404 1.13887059
|
||||
-0.29932935 0.38962638
|
||||
0.31620505 -1.34874423
|
||||
-0.35708404 0.43969553
|
||||
0.54695336 -0.41951872
|
||||
-0.36976770 0.99489550
|
||||
-1.00577935 1.04650656
|
||||
-0.31839248 0.32966467
|
||||
0.27710048 1.14046577
|
||||
-0.34954991 0.89579091
|
||||
0.46003161 -0.38668754
|
||||
-1.09483649 0.45558599
|
||||
-1.32749255 -1.08655702
|
||||
-0.41384999 0.53388335
|
||||
-0.92586659 -0.23732743
|
||||
-0.11461978 1.22614317
|
||||
-1.13476543 -0.15516704
|
||||
-0.03625266 0.88792143
|
||||
-0.36426959 -0.31917900
|
||||
-0.17525306 0.13101065
|
||||
-0.29878850 -0.40565286
|
||||
0.06367635 1.21715588
|
||||
-0.30032239 -0.25475547
|
||||
0.09898865 1.02720529
|
||||
-0.06611431 -1.15194721
|
||||
-0.39215060 -0.30150449
|
||||
-0.38680782 0.31314671
|
||||
-0.06247327 0.36132059
|
||||
0.20109728 -0.64570689
|
||||
-0.23945673 0.64085696
|
||||
-0.13446949 0.70075405
|
||||
-0.75243730 0.91726734
|
||||
-1.45028028 -0.43988942
|
||||
-0.07926376 0.54618428
|
||||
0.35151647 0.12328557
|
||||
-0.00043103 1.20547449
|
||||
-0.41878158 -0.11462993
|
||||
-0.13188202 0.60069882
|
||||
0.33204624 1.19833977
|
||||
-0.15340704 0.41633354
|
||||
-0.41617623 -0.17790618
|
||||
-0.56285739 0.88577465
|
||||
0.39265124 -1.47117966
|
||||
-0.57457171 0.94599611
|
||||
-1.57813091 -0.27615434
|
||||
0.25275661 1.03110215
|
||||
-0.84474420 -0.90311830
|
||||
0.25303978 1.04767460
|
||||
-0.63961834 1.02757341
|
||||
-0.15838839 0.44901836
|
||||
-0.48000317 -0.10596114
|
||||
-0.25937394 0.42590391
|
||||
0.84493295 -0.70805225
|
||||
0.24423362 1.01963863
|
||||
-1.48665793 -0.14370221
|
||||
0.07376786 0.42983815
|
||||
-0.59202645 -0.01731182
|
||||
0.18705556 1.24782046
|
||||
0.56661178 -0.02647487
|
||||
-0.60722140 0.86222403
|
||||
-1.48738136 -0.43921204
|
||||
-0.15651173 0.61619336
|
||||
-0.19885992 1.58694503
|
||||
0.07385597 0.42852972
|
||||
-0.41682070 -0.14784924
|
||||
-0.77195688 0.87679915
|
||||
-1.40503884 -0.15440256
|
||||
-0.14103875 0.68667030
|
||||
-0.76379382 0.85984795
|
||||
-0.12619804 0.52072503
|
||||
0.12570952 -1.63658614
|
||||
0.09358821 0.46519043
|
||||
1.51702274 0.10571578
|
||||
-0.15718465 0.55544646
|
||||
1.12819516 -0.37212082
|
||||
-0.01459698 0.53777689
|
||||
-1.07351142 -0.64113878
|
||||
-0.57369112 0.88411436
|
||||
-1.46870106 -0.19352991
|
||||
0.35927911 1.03440056
|
||||
0.30424108 1.10045036
|
||||
0.02981222 0.51366749
|
||||
-0.16684004 -1.23339544
|
||||
0.05097045 0.26090896
|
||||
0.09668531 -0.54237611
|
||||
0.38772762 1.03212178
|
||||
-0.30759050 1.47875603
|
||||
-0.02502066 0.42910717
|
||||
0.94844497 0.80507692
|
||||
0.28830246 1.24135339
|
||||
0.14865732 -0.40934806
|
||||
-0.64849870 1.02545190
|
||||
-1.64220884 -0.23073861
|
||||
-0.04115289 0.50218823
|
||||
-1.06285152 0.11311046
|
||||
0.16416052 1.04673165
|
||||
-1.11102027 0.08544000
|
||||
0.30091542 1.08688638
|
||||
-0.54551043 -0.18095572
|
||||
-0.23788564 0.53265900
|
||||
-0.49872342 -0.26678921
|
||||
0.22270801 0.99856952
|
||||
-0.60363472 -0.05824273
|
||||
0.13204250 1.13170488
|
||||
-0.26664132 -1.21950767
|
||||
-0.56062165 -0.26617337
|
||||
0.00727415 0.21197638
|
||||
-0.10690293 0.56904332
|
||||
-0.12173691 -0.55223481
|
||||
-0.13609125 0.38164037
|
||||
-0.22226458 0.71446498
|
||||
-0.74245911 0.85892687
|
||||
-1.55972244 -0.41484829
|
||||
-0.20334295 0.61970357
|
||||
0.31700129 0.07079464
|
||||
0.13449857 1.02302289
|
||||
-0.54238275 -0.31894746
|
||||
-0.28133418 0.52067079
|
||||
0.28060094 1.12630923
|
||||
-0.12162546 0.36634303
|
||||
-0.34561133 -0.17827835
|
||||
-1.03931702 0.70751582
|
||||
0.59423681 -1.29911937
|
||||
-0.75843257 0.81139880
|
||||
-1.57730446 -0.68556803
|
||||
0.10310968 1.03335362
|
||||
-0.67400306 -0.84700429
|
||||
0.06054034 1.15025742
|
||||
-0.77286167 0.47557829
|
||||
-0.29028847 0.54276291
|
||||
-0.27482760 -0.19772142
|
||||
-0.11994968 0.36408189
|
||||
0.78035295 -0.81713160
|
||||
0.21737616 1.09857289
|
||||
-1.18254850 -0.66539897
|
||||
-0.35210721 0.39493518
|
||||
-0.50975569 -0.29188620
|
||||
0.07661012 1.15470996
|
||||
0.35848649 0.26022851
|
||||
-0.77485063 0.69526899
|
||||
-1.22207178 -0.46175769
|
||||
-0.18692438 0.60206436
|
||||
-0.62741906 1.33494975
|
||||
-0.17855601 0.34009484
|
||||
-0.42621663 -0.13057148
|
||||
-0.83850938 0.63075363
|
||||
-1.47199973 -0.51955942
|
||||
-0.23964524 0.59891281
|
||||
-1.01720429 0.84543467
|
||||
-0.25567097 0.34493172
|
||||
0.65782207 -1.35499684
|
||||
-0.13645103 0.42344955
|
||||
1.42139339 0.70402143
|
||||
-0.14134279 0.58163801
|
||||
1.26429849 -0.05123437
|
||||
-0.15573860 0.66425208
|
||||
-0.84664093 -0.89746670
|
||||
-0.66765526 0.65188491
|
||||
-1.33946979 -0.68336950
|
||||
-0.13547805 1.33391456
|
||||
-0.01311918 1.08381016
|
||||
-0.26885915 0.39044922
|
||||
-0.09629192 -1.12985380
|
||||
-0.18477312 0.41780574
|
||||
0.11788432 -0.36930565
|
||||
-0.11400622 1.14616969
|
||||
-0.71164575 1.14762749
|
||||
-0.19434913 0.51072915
|
||||
0.57606657 0.96776987
|
||||
-0.02115547 1.18705119
|
||||
-0.03177567 -0.52325305
|
||||
-1.03661648 0.39131082
|
||||
-1.19200246 -1.02791732
|
||||
-0.44145789 0.43868764
|
||||
-0.96775095 -0.11064332
|
||||
-0.35379494 1.18148992
|
||||
-1.00698678 -0.48198386
|
||||
-0.26498858 1.12476012
|
||||
-0.26143369 -0.37768083
|
||||
-0.51512946 0.36354685
|
||||
-0.19957862 -0.31824547
|
||||
-0.66653588 0.85918826
|
||||
-0.37520507 -0.35376709
|
||||
-0.44493040 0.93295373
|
||||
0.55064254 -1.23389527
|
||||
-0.28115978 -0.30295983
|
||||
-0.54907227 0.21053437
|
||||
-0.50793530 0.37345227
|
||||
0.37024987 -0.25310419
|
||||
-0.35959453 0.11118461
|
||||
-0.50281000 0.16548987
|
||||
-1.26322433 0.05900470
|
||||
-0.64710033 -1.49218071
|
||||
-0.51790185 0.28396723
|
||||
0.31222570 0.27090413
|
||||
-0.71499001 0.96604950
|
||||
-0.16468815 -0.42049513
|
||||
-0.52287459 0.23433390
|
||||
-0.87932719 0.71816438
|
||||
-0.43828506 0.00462493
|
||||
0.02218342 -0.60768376
|
||||
-1.20981804 -0.13172407
|
||||
1.39808402 -0.43885489
|
||||
-0.97813323 -0.26248931
|
||||
-0.25206935 -1.34281904
|
||||
-0.75325400 0.81961764
|
||||
0.14422066 -1.34191266
|
||||
-0.79582225 0.96136586
|
||||
-1.16155513 -0.26329434
|
||||
-0.86013035 0.06741854
|
||||
-0.00965300 -0.55094971
|
||||
-0.49391380 0.00109182
|
||||
1.11157602 0.40557103
|
||||
-1.11950886 0.76793436
|
||||
-0.11150682 -1.78154507
|
||||
-0.43682487 0.08305273
|
||||
-0.12393961 -0.55946423
|
||||
-0.98192766 0.78343527
|
||||
-0.06091643 0.58229269
|
||||
-1.03669652 -0.42149417
|
||||
-0.21401098 -1.64487810
|
||||
-0.50680768 0.13930195
|
||||
-1.33574727 0.14706836
|
||||
-0.55714862 0.12537793
|
||||
0.13380015 -0.61643711
|
||||
-1.09488278 -0.08747533
|
||||
-0.21066788 -1.44769192
|
||||
-0.50229181 -0.04011964
|
||||
-1.13638894 -0.23384381
|
||||
-0.43836050 0.11383238
|
||||
1.39228950 -0.10799639
|
||||
-0.48802878 0.20264974
|
||||
0.47588619 1.62192837
|
||||
-0.27340939 0.19361598
|
||||
0.68448629 0.80864114
|
||||
-0.47195792 0.05420133
|
||||
0.34477091 -1.07786560
|
||||
-0.92004004 -0.25424321
|
||||
-0.37760272 -1.39721129
|
||||
-0.87787868 0.89426868
|
||||
-0.79267036 0.75520147
|
||||
-0.48065034 0.07817515
|
||||
0.85184258 -0.83563910
|
||||
-0.47239528 0.12627174
|
||||
0.48277009 -0.18180439
|
||||
-0.71997171 0.75293982
|
||||
-1.29153937 0.42384830
|
||||
-0.32613062 0.17402085
|
||||
0.10202947 1.04593394
|
||||
-0.71125374 0.74980584
|
||||
0.50484731 -0.16738092
|
||||
-1.16256423 -0.00470142
|
||||
-0.58952940 -1.43917961
|
||||
-0.71851513 0.20431374
|
||||
-0.82581369 -0.70146519
|
||||
-0.43309512 1.05489748
|
||||
-0.87632697 -0.63470868
|
||||
-0.56174046 0.94078784
|
||||
-0.36343466 -0.74030249
|
||||
-0.47009920 0.32511040
|
||||
-0.28949047 -0.46301691
|
||||
-0.61098971 0.83197266
|
||||
-0.30513604 -0.52105098
|
||||
-0.75061540 0.87475951
|
||||
0.50974994 -1.08883784
|
||||
-0.28154826 -0.52078296
|
||||
-0.34873231 0.28167771
|
||||
-0.41149750 0.26937057
|
||||
0.19007326 -0.31381309
|
||||
-0.39222476 0.35812278
|
||||
-0.40291930 0.24467795
|
||||
-1.10924908 0.34485345
|
||||
-0.98085212 -1.05744586
|
||||
-0.27406176 0.26544570
|
||||
0.35012828 0.35381774
|
||||
-0.40087355 1.19550172
|
||||
-0.44650037 -0.43294378
|
||||
-0.29257517 0.33133304
|
||||
-0.13780562 1.23531423
|
||||
-0.51645492 0.53254651
|
||||
-0.45365730 -0.41872486
|
||||
-1.03611147 0.47877429
|
||||
0.79832435 -1.32896077
|
||||
-1.11551489 0.40800484
|
||||
-1.11975639 -0.79847563
|
||||
-0.20611406 1.00724166
|
||||
-0.69838131 -1.03603068
|
||||
-0.23096181 1.08306487
|
||||
-1.07023457 0.66243075
|
||||
-0.15173137 0.52598993
|
||||
-0.38361285 -0.38311364
|
||||
-0.34378062 0.47000297
|
||||
1.03804857 -0.51402522
|
||||
-0.09490212 1.19969154
|
||||
-1.27634648 -0.70130648
|
||||
-0.42548178 0.45803214
|
||||
-0.33217760 -0.25636176
|
||||
-0.36035026 1.01323624
|
||||
0.58868663 0.21384701
|
||||
-1.12508238 0.78429355
|
||||
-1.32778013 -0.87805779
|
||||
-0.31778980 0.37782146
|
||||
-0.88447401 1.32004993
|
||||
-0.31173481 0.18936826
|
||||
-0.46230529 -0.38937533
|
||||
-0.93875591 0.40586591
|
||||
-1.19880808 -0.88001421
|
||||
-0.33507246 0.41148662
|
||||
-1.02904286 0.47779753
|
||||
-0.30145877 0.41586993
|
||||
0.86851117 -1.10918688
|
||||
-0.50406166 0.13349449
|
||||
1.00458594 0.99686515
|
||||
-0.38307432 0.36008430
|
||||
1.00433321 0.37749771
|
||||
-0.50953245 0.15792285
|
||||
-0.39888021 -1.18812231
|
||||
-1.08444469 0.41794647
|
||||
-1.00859605 -1.36982347
|
||||
-0.39929721 1.00858598
|
||||
-0.57043925 0.94134680
|
||||
-0.45270563 0.30986828
|
||||
0.64740169 -0.87173231
|
||||
-0.34851350 0.43318799
|
||||
0.50373214 -0.46736993
|
||||
-0.49324614 0.90839987
|
||||
-1.16083187 0.85756541
|
||||
-0.37995838 0.30665239
|
||||
0.14942910 1.19252524
|
||||
-0.45750199 1.00124895
|
||||
0.39533489 -0.24618649
|
||||
-1.09194749 0.30059665
|
||||
-0.53064464 -1.29745411
|
||||
-0.32829564 0.19645699
|
||||
-1.00491093 -0.49180643
|
||||
-0.69388058 0.98571285
|
||||
-1.07327052 -0.52013451
|
||||
-0.49550616 0.92068615
|
||||
-0.28149445 -0.44098059
|
||||
-0.42975495 0.38440311
|
||||
-0.34010936 -0.33895185
|
||||
-0.54740262 1.13529512
|
||||
-0.25928445 -0.47033113
|
||||
-0.58739701 1.01800468
|
||||
0.48617337 -1.07609991
|
||||
-0.40499489 -0.57781754
|
||||
-0.40174748 0.55155045
|
||||
-0.32286706 0.26588019
|
||||
0.53629286 -0.42122792
|
||||
-0.39920489 0.51713042
|
||||
-0.45359370 0.40737035
|
||||
-0.84548899 0.27785657
|
||||
-1.00667459 -0.96613366
|
||||
-0.26828484 0.25976331
|
||||
0.25295097 0.37607576
|
||||
-0.46126431 1.03260901
|
||||
-0.40279099 -0.24657659
|
||||
-0.37987717 0.12352906
|
||||
-0.44628748 0.86769292
|
||||
-0.29637438 0.45896986
|
||||
-0.39936855 -0.08745092
|
||||
-0.95340646 0.37842294
|
||||
0.90178683 -1.14345809
|
||||
-0.80646092 0.53265844
|
||||
-1.27088485 -0.89749663
|
||||
-0.24802130 0.91333138
|
||||
-0.36572702 -0.99916742
|
||||
-0.36641645 1.05426660
|
||||
-1.10837220 0.52275107
|
||||
-0.44607592 0.15779324
|
||||
-0.23228534 -0.24406293
|
||||
-0.36384164 0.45002486
|
||||
1.03058155 -0.56498483
|
||||
-0.14187980 1.07303759
|
||||
-1.15376965 -1.13637382
|
||||
-0.36264501 0.52815460
|
||||
-0.39051068 -0.24438521
|
||||
-0.07492065 1.14349287
|
||||
0.51684344 0.08621823
|
||||
-0.96172276 0.52619069
|
||||
-1.04576540 -1.11777830
|
||||
-0.34938487 0.45655008
|
||||
-1.02766240 1.18495194
|
||||
-0.52997885 0.55748157
|
||||
-0.15364100 -0.44508072
|
||||
-1.14329676 0.23788243
|
||||
-1.21779719 -1.11349751
|
||||
-0.38546200 0.42075026
|
||||
-1.17918179 0.52604491
|
||||
-0.28574420 0.52199428
|
||||
1.08330303 -0.98066820
|
||||
-0.24802582 0.27160918
|
||||
1.04853566 1.08102404
|
||||
-0.29313743 0.42485423
|
||||
0.96260308 0.26803572
|
||||
-0.51912360 0.34420975
|
||||
-0.11292000 -1.14400559
|
||||
-1.04739492 0.17540820
|
||||
-1.11441406 -1.22883010
|
||||
-0.52639594 0.84986367
|
||||
-0.52710069 1.01735010
|
||||
-0.30885479 0.27388498
|
||||
0.44574757 -0.74425968
|
||||
-0.46382300 0.44814634
|
||||
0.44904213 -0.30813329
|
||||
-0.63767229 1.21418607
|
||||
-1.22556632 0.86395503
|
||||
-0.34214451 0.30574735
|
||||
0.17259234 0.86496400
|
||||
-0.32480619 1.05612537
|
||||
0.43889334 -0.51933604
|
||||
-1.02674349 0.35609282
|
||||
-0.89393514 -1.21049417
|
||||
-0.36846369 0.37941633
|
||||
-0.98641807 -0.47559413
|
||||
-0.44755679 1.21960269
|
||||
-1.03122128 -0.33682392
|
||||
-0.35318795 1.12940664
|
||||
-0.36331006 -0.32896336
|
||||
-0.28306541 0.38647142
|
||||
-0.21812509 -0.27670337
|
||||
-0.44729084 1.10887448
|
||||
-0.52859975 -0.36404097
|
||||
-0.24010453 0.91605871
|
||||
0.41984611 -1.15062846
|
||||
-0.34248795 -0.24005071
|
||||
-0.33397378 0.35878522
|
||||
-0.39831462 0.50795634
|
||||
0.41702786 -0.62163274
|
||||
-0.33805101 0.45650539
|
||||
-0.27859581 0.40216109
|
||||
-1.05774569 0.53779507
|
||||
-1.42415680 -0.87055713
|
||||
-0.27654375 0.38928806
|
||||
0.25901039 0.35652859
|
||||
-0.33082657 1.15409645
|
||||
-0.40685393 -0.20407365
|
||||
-0.23686438 0.56633035
|
||||
-0.26193902 1.21365439
|
||||
-0.20260600 0.49133397
|
||||
-0.45165163 -0.45751703
|
||||
-1.18346583 0.31810438
|
||||
0.89802934 -0.95013768
|
||||
-1.16309863 0.26604245
|
||||
-1.18386506 -1.10311721
|
||||
-0.16784287 1.05364788
|
||||
-0.49077470 -0.99025124
|
||||
0.00171671 1.15509224
|
||||
-0.86606100 0.28942709
|
||||
-0.30559667 0.49226883
|
||||
-0.39902748 -0.41785924
|
||||
-0.33112115 0.32495728
|
||||
0.94443664 -0.36892219
|
||||
-0.38348413 1.14426380
|
||||
-1.18058683 -1.05089383
|
||||
-0.39125284 0.36801023
|
||||
-0.35680139 -0.41245427
|
||||
-0.25121134 1.15931919
|
||||
0.36056736 0.40619567
|
||||
-1.09116775 0.36957303
|
||||
-1.13642958 -1.11281056
|
||||
-0.34056852 0.38007119
|
||||
-1.21838730 0.91970197
|
||||
-0.17884434 0.42245022
|
||||
-0.27551533 -0.38898298
|
||||
-1.25507508 0.38253654
|
||||
-1.21026313 -0.87005343
|
||||
-0.39228963 0.62309083
|
||||
-1.16448806 0.36472580
|
||||
-0.24315635 0.34112968
|
||||
1.07454905 -0.99510410
|
||||
-0.35416482 0.42511567
|
||||
1.03604010 1.07082031
|
||||
-0.39485070 0.40169023
|
||||
1.02018178 0.32303944
|
||||
-0.24714138 0.27513322
|
||||
-0.25443819 -1.04224213
|
||||
-1.00335985 0.33336367
|
||||
-0.95167241 -1.13384472
|
||||
-0.47178371 1.10380719
|
||||
-0.44347102 0.89865945
|
||||
-0.35612593 0.24575863
|
||||
0.59544992 -0.97082846
|
||||
-0.33079820 0.18842791
|
||||
0.24063116 -0.40954664
|
||||
-0.54597573 1.11706897
|
||||
-1.09356272 0.86551165
|
||||
-0.25644651 0.42451965
|
||||
0.31590044 1.01547745
|
||||
-0.57685164 1.03765963
|
||||
0.26467495 -0.26793133
|
||||
-1.02806397 0.49413336
|
||||
-1.07742309 -1.17025750
|
||||
-0.33545549 0.24268994
|
||||
-1.03818129 -0.41749668
|
||||
-0.37066685 0.90622096
|
||||
-1.13805168 -0.46375121
|
||||
-0.41169235 1.07307811
|
||||
-0.47517936 -0.31180552
|
||||
-0.36701432 0.38670108
|
||||
-0.41929830 -0.47138018
|
||||
-0.34942081 1.04098519
|
||||
-0.47106518 -0.32153553
|
||||
-0.37195463 1.05910723
|
||||
0.14637113 -1.06527293
|
||||
-0.10291353 0.42424225
|
||||
0.26947551 -0.43055636
|
||||
-0.40148014 0.38862774
|
||||
-0.40945108 0.39705641
|
||||
-1.04244368 0.49395299
|
||||
-1.27365240 -0.86786814
|
||||
-0.26467394 0.31488035
|
||||
0.31891504 0.35451162
|
||||
-0.04344638 1.10074605
|
||||
-0.39683701 -0.29984324
|
||||
-0.21146179 0.57934039
|
||||
-0.18430679 1.20083875
|
||||
-0.34647055 0.21374240
|
||||
-0.36800694 -0.20463353
|
||||
-0.96803254 0.36613104
|
||||
1.00439662 -1.14286191
|
||||
-1.08184906 0.38984420
|
||||
-0.86737954 -0.96935370
|
||||
-0.21922013 1.04052779
|
||||
-0.45758990 -1.02426336
|
||||
-0.21924203 1.16949224
|
||||
-0.91730740 0.13145528
|
||||
-0.35857284 0.25584683
|
||||
-0.30528757 -0.35552063
|
||||
-0.45500790 0.29050995
|
||||
1.09478106 -0.29988615
|
||||
-0.36959306 0.95649714
|
||||
-0.91490210 -1.26214271
|
||||
-0.38509185 0.39302859
|
||||
-0.35739352 -0.34429725
|
||||
-0.31262994 0.82016967
|
||||
0.38738991 0.37644184
|
||||
-0.98502444 0.18314526
|
||||
-1.05050547 -1.20645447
|
||||
-0.33670658 0.36506942
|
||||
-1.21343766 0.93698166
|
||||
-0.09855738 0.24705256
|
||||
-0.33113277 -0.38521955
|
||||
-0.89739483 0.08841987
|
||||
-0.89465442 -0.94690980
|
||||
-0.42145584 0.28343312
|
||||
-1.06849193 0.47108899
|
||||
-0.39732226 0.44606233
|
||||
1.06930625 -1.11142101
|
||||
-0.29206820 0.39997851
|
||||
0.78073372 0.90899737
|
||||
-0.25421280 0.32608096
|
||||
0.91519669 0.45997640
|
||||
-0.38956838 0.39409618
|
||||
-0.26364644 -1.24426066
|
||||
-1.05191054 0.38146703
|
||||
-1.01747304 -1.08036571
|
||||
-0.38768865 1.09525211
|
||||
-0.38096670 1.04466527
|
||||
-0.36889228 0.35087106
|
||||
0.39210062 -1.23945201
|
||||
-0.24018273 0.20336419
|
||||
0.40875677 -0.57228107
|
||||
-0.41285536 1.16986380
|
||||
-0.98056291 0.94739535
|
||||
-0.26388705 0.51398988
|
||||
0.27895279 1.06031029
|
||||
-0.27752604 1.15191683
|
||||
0.33890427 -0.48476732
|
||||
-1.09037791 0.46464290
|
||||
-1.33508329 -0.81219355
|
||||
-0.39158408 0.44958947
|
||||
-1.24108745 -0.16653263
|
||||
-0.17524396 1.04890332
|
||||
-1.13377063 -0.37194277
|
||||
-0.44164358 1.11606653
|
||||
-0.50379533 -0.21441051
|
||||
-0.42550406 0.40402679
|
||||
-0.41738928 -0.50437718
|
||||
-0.20289748 0.96368079
|
||||
-0.29945932 -0.41636288
|
||||
-0.37797060 0.94457881
|
||||
0.39310239 -1.29272896
|
||||
@ -1,16 +0,0 @@
|
||||
-1.02899151 -1.02899151
|
||||
-1.02899151 -0.34299717
|
||||
-1.02899151 0.34299717
|
||||
-1.02899151 1.02899151
|
||||
-0.34299717 -1.02899151
|
||||
-0.34299717 -0.34299717
|
||||
-0.34299717 0.34299717
|
||||
-0.34299717 1.02899151
|
||||
0.34299717 -1.02899151
|
||||
0.34299717 -0.34299717
|
||||
0.34299717 0.34299717
|
||||
0.34299717 1.02899151
|
||||
1.02899151 -1.02899151
|
||||
1.02899151 -0.34299717
|
||||
1.02899151 0.34299717
|
||||
1.02899151 1.02899151
|
||||
@ -1,361 +0,0 @@
|
||||
#include <math.h>
|
||||
#include <stdint.h>
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <complex.h>
|
||||
#include <string.h>
|
||||
|
||||
#define A 1
|
||||
|
||||
struct qam_system_s {
|
||||
int M; // Nombre de symboles M-QAM
|
||||
int k; // Nombre de bits/symboles
|
||||
double Fs; // Fréquence d'échantillionage
|
||||
double Ts; // Temps d'échantillionage
|
||||
int N; // Nombre d'échantillions
|
||||
double Fc; // Fréquence de la porteuse
|
||||
double complex** constellation; // Tableau de symboles I + j Q
|
||||
};
|
||||
typedef struct qam_system_s qam_system;
|
||||
|
||||
// Initialisation de la constellation (double tableau de taille sqrt(M)),
|
||||
// ToDo : changer à un tableau à 1 dimension pour éviter de calculer sqrt(M)
|
||||
void init_constellation (qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
qam->constellation = (double complex**)malloc(sizeof(double complex*) * sm);
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
qam->constellation[i] = (double complex*)malloc(sizeof(double complex) * sm);
|
||||
}
|
||||
|
||||
//double norm_factor = sqrt((double)(qam->M - 1) / 3.0); // Pour puissance unitaire
|
||||
double norm_factor = sqrt((double)(1.7*(qam->M - 1)/3.0));
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
double complex ip = -(sm - 1) + 2 * i;
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex qp = -(sm - 1) + 2 * j;
|
||||
qam->constellation[i][j] = (ip + I * qp) / norm_factor;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex iq = symbols[k];
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
//s[k * qam->N + n] = A * iq * cexp(2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
int idx = k * qam->N + n;
|
||||
s[idx] = A * iq * cexp(2 * I * M_PI * qam->Fc * ((double)idx / qam->Fs));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM
|
||||
void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat, FILE *fp_constel) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
//r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)(k * qam->N + n) / qam->Fs));
|
||||
|
||||
|
||||
}
|
||||
r /= qam->N;
|
||||
r /= A;
|
||||
|
||||
if (fp_constel) {
|
||||
fprintf(fp_constel, "% .8f % .8f\n", creal(r), cimag(r));
|
||||
fflush(fp_constel);
|
||||
}
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation (lent)
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// index du symbole (id) : même mappage que dans bits_to_symbols()
|
||||
int id = i_cl * sm + j_cl;
|
||||
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + b] = (id >> (qam->k - 1 - b)) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
double compare_bits(uint8_t* bits1, uint8_t* bits2, int nb_bits) {
|
||||
int errors = 0;
|
||||
for (int i = 0; i < nb_bits; i++) {
|
||||
if (bits1[i] != bits2[i]) errors++;
|
||||
}
|
||||
return (double)errors / nb_bits;
|
||||
}
|
||||
|
||||
// ---------- Helpers ----------
|
||||
static inline double complex interp_linear(double complex *x, int N, double t) {
|
||||
// t: position en échantillons (peut être fractionnaire). On suppose 0 <= t <= N-2
|
||||
int i = (int)floor(t);
|
||||
double mu = t - (double)i;
|
||||
if (i < 0) i = 0;
|
||||
if (i >= N-1) return x[N-1];
|
||||
return (1.0 - mu) * x[i] + mu * x[i+1];
|
||||
}
|
||||
|
||||
static int nearest_constellation_index(qam_system* qam, double complex z, int *out_i, int *out_j) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int best_i = 0, best_j = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(z - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d; best_i = i; best_j = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
if (out_i) *out_i = best_i;
|
||||
if (out_j) *out_j = best_j;
|
||||
return best_i * sm + best_j;
|
||||
}
|
||||
|
||||
// --- Démodulation avec M&M + PLL et écriture dans fichier ---
|
||||
void demodulate_real(qam_system* qam, double complex* rx_samples, int total_samples,
|
||||
uint8_t* bits_hat, const char* filename)
|
||||
{
|
||||
double sps = (double) qam->N; // échantillons par symbole
|
||||
double Kp_t = 0.005, Ki_t = 0.001;
|
||||
double Kp_c = 0.01, Ki_c = 0.01; // PLL porteur
|
||||
double timing_integrator = 0.0;
|
||||
double phase_integrator = 0.0, phase_est = 0.0;
|
||||
double t = sps/2.0;
|
||||
|
||||
int max_symbols = (int)(total_samples / sps) + 10;
|
||||
double complex* syms = malloc(sizeof(double complex) * max_symbols);
|
||||
int sym_count = 0;
|
||||
|
||||
double power = 0.0;
|
||||
for(int i = 0; i < total_samples; i++) power += pow(cabs(rx_samples[i]), 2);
|
||||
power /= total_samples;
|
||||
double agc_gain = 1.0;
|
||||
if(power > 0) agc_gain = 1.0 / sqrt(power);
|
||||
|
||||
// Ouverture du fichier pour écrire les symboles reçus
|
||||
FILE* fp = fopen(filename, "w");
|
||||
if(!fp) {
|
||||
perror("Erreur ouverture fichier .dat");
|
||||
free(syms);
|
||||
return;
|
||||
}
|
||||
|
||||
while((int)floor(t) + 1 < total_samples) {
|
||||
double complex yk = interp_linear(rx_samples, total_samples, t) * agc_gain;
|
||||
double complex z = yk * cexp(-I * phase_est);
|
||||
syms[sym_count++] = z;
|
||||
|
||||
// Écriture de chaque symbole reçu dans le fichier
|
||||
fprintf(fp, "% .8f % .8f\n", creal(z), cimag(z));
|
||||
|
||||
if(sym_count >= 3){
|
||||
double complex y_k = syms[sym_count-1];
|
||||
double complex y_k_1 = syms[sym_count-2];
|
||||
double complex y_k_2 = syms[sym_count-3];
|
||||
double e_t = creal((y_k - y_k_2) * conj(y_k_1));
|
||||
timing_integrator += Ki_t * e_t;
|
||||
double timing_adjust = Kp_t * e_t + timing_integrator;
|
||||
t += sps + timing_adjust;
|
||||
} else t += sps;
|
||||
|
||||
// PLL carrier
|
||||
if(sym_count >= 1) {
|
||||
int ii, jj;
|
||||
int idx = nearest_constellation_index(qam, z, &ii, &jj);
|
||||
double complex d = qam->constellation[ii][jj];
|
||||
double e_phase = cimag(z * conj(d));
|
||||
phase_integrator += Ki_c * e_phase;
|
||||
double phase_adj = Kp_c * e_phase + phase_integrator;
|
||||
phase_est += phase_adj;
|
||||
if(phase_est > M_PI) phase_est -= 2*M_PI;
|
||||
if(phase_est < -M_PI) phase_est += 2*M_PI;
|
||||
}
|
||||
}
|
||||
|
||||
// Décision finale et reconstruction bits
|
||||
int nb_symbols = sym_count;
|
||||
int sm = (int)sqrt(qam->M);
|
||||
int nb_bits = nb_symbols * qam->k;
|
||||
for(int k = 0; k < nb_symbols; k++) {
|
||||
double complex z = syms[k];
|
||||
int ii, jj;
|
||||
int id = nearest_constellation_index(qam, z, &ii, &jj);
|
||||
for(int b = 0; b < qam->k; b++) {
|
||||
int bit = (id >> (qam->k - 1 - b)) & 1;
|
||||
if(k * qam->k + b < nb_bits) bits_hat[k * qam->k + b] = (uint8_t)bit;
|
||||
}
|
||||
}
|
||||
|
||||
fclose(fp);
|
||||
free(syms);
|
||||
}
|
||||
|
||||
|
||||
int main () {
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
//qam.Ts = 0.0003;
|
||||
//qam.N = (int)(qam.Fs * qam.Ts);
|
||||
qam.N = 100;
|
||||
qam.Ts = qam.N / qam.Fs;
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
//int nb_bits = 1000;
|
||||
//int nb_symbols = nb_bits / qam.k;
|
||||
|
||||
//uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
//for (int i = 0; i < nb_bits; i++) {
|
||||
// input_bits[i] = rand() % 2;
|
||||
//}
|
||||
char* texte = "VVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxVif juge, trempez ce blond whisky aqueuxif juge, trempez ce blond whisky aqueux";
|
||||
int nb_chars = strlen(texte);
|
||||
int nb_bits = nb_chars * 8;
|
||||
int nb_symbols = (nb_bits + qam.k - 1) / qam.k;
|
||||
|
||||
// Conversion du texte en bits
|
||||
uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
for(int b = 0; b < 8; b++){
|
||||
input_bits[i*8 + b] = (texte[i] >> (7-b)) & 1;
|
||||
}
|
||||
}
|
||||
|
||||
// Conversion en symboles
|
||||
double complex* symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
bits_to_symbols(&qam, input_bits, nb_bits, symbols);
|
||||
|
||||
// Modulation
|
||||
int total_samples = qam.N * nb_symbols;
|
||||
double complex* s = malloc(sizeof(double complex) * total_samples);
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
//double phase_offset = M_PI / 6.0; // 30 degrés
|
||||
//for (int i = 0; i < total_samples; i++) {
|
||||
// s[i] *= cexp(I * phase_offset);
|
||||
//}
|
||||
//double freq_offset = 0; // Hz de décalage
|
||||
//for (int i = 0; i < total_samples; i++) {
|
||||
// double t = (double)i / qam.Fs;
|
||||
// s[i] *= cexp(I * 2 * M_PI * freq_offset * t);
|
||||
//}
|
||||
//int offset_samples = (int)(0.3 * qam.N); // décalage de 30% d’un symbole
|
||||
//memmove(s + offset_samples, s, (total_samples - offset_samples) * sizeof(double complex));
|
||||
|
||||
|
||||
|
||||
// Ajout du bruit
|
||||
//double snr_dB = 5; // SNR en dB
|
||||
//double signal_power = 1.0; // puissance moyenne unitaire après normalisation
|
||||
//double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
//double sigma = sqrt(signal_power / (2.0 * snr_lin)); // /2 car bruit sur I et Q
|
||||
add_noise(s, total_samples, 0.1);
|
||||
|
||||
FILE *fp_ref = fopen("constellation_ref.dat", "w");
|
||||
int sm = (int)sqrt(qam.M);
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
fprintf(fp_ref, "% .8f % .8f\n", creal(qam.constellation[i][j]), cimag(qam.constellation[i][j]));
|
||||
}
|
||||
}
|
||||
fclose(fp_ref);
|
||||
FILE *fp_constel = fopen("constellation.dat", "w");
|
||||
|
||||
// Démodulation
|
||||
uint8_t* output_bits = (uint8_t*)malloc(nb_bits * sizeof(uint8_t));
|
||||
//demodulate(&qam, s, nb_symbols, output_bits, fp_constel);
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
s[i] *= cexp(-2 * I * M_PI * qam.Fc * (i / qam.Fs));
|
||||
}
|
||||
|
||||
demodulate_real(&qam, s, total_samples, output_bits, "constellation.dat");
|
||||
|
||||
|
||||
fclose(fp_constel);
|
||||
|
||||
// Reconstruction du texte
|
||||
char* texte_recup = malloc(nb_chars + 1);
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
char c = 0;
|
||||
for(int b = 0; b < 8; b++){
|
||||
c |= output_bits[i*8 + b] << (7-b);
|
||||
}
|
||||
texte_recup[i] = c;
|
||||
}
|
||||
texte_recup[nb_chars] = '\0';
|
||||
printf("Texte original : %s\n\n", texte);
|
||||
printf("Texte demodulé : %s\n", texte_recup);
|
||||
|
||||
// Calcul du BER
|
||||
double ber = compare_bits(input_bits, output_bits, nb_bits);
|
||||
printf("Taux d'erreur blind QAM: %.4f\n", ber * 100);
|
||||
|
||||
// Libération mémoire
|
||||
free(input_bits);
|
||||
free(output_bits);
|
||||
free(symbols);
|
||||
free(s);
|
||||
free_constellation(&qam);
|
||||
|
||||
return 0;
|
||||
}
|
||||
BIN
QAM/OLD/qam
BIN
QAM/OLD/qam
Binary file not shown.
@ -1,238 +0,0 @@
|
||||
#include <math.h>
|
||||
#include <stdint.h>
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <complex.h>
|
||||
#include "../wav/wav.h"
|
||||
#include "../files/files.h"
|
||||
#include <string.h>
|
||||
|
||||
#define A 1
|
||||
|
||||
struct qam_system_s {
|
||||
int M; // Nombre de symboles M-QAM
|
||||
int k; // Nombre de bits/symboles
|
||||
double Fs; // Fréquence d'échantillionage
|
||||
double Ts; // Temps d'échantillionage
|
||||
int N; // Nombre d'échantillions
|
||||
double Fc; // Fréquence de la porteuse
|
||||
double complex** constellation; // Tableau de symboles I + j Q
|
||||
};
|
||||
typedef struct qam_system_s qam_system;
|
||||
|
||||
// Initialisation de la constellation (double tableau de taille sqrt(M)),
|
||||
// ToDo : changer à un tableau à 1 dimension pour éviter de calculer sqrt(M)
|
||||
void init_constellation (qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
qam->constellation = (double complex**)malloc(sizeof(double complex*) * sm);
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
qam->constellation[i] = (double complex*)malloc(sizeof(double complex) * sm);
|
||||
}
|
||||
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0); // Pour puissance unitaire
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
double complex ip = -(sm - 1) + 2 * i;
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex qp = -(sm - 1) + 2 * j;
|
||||
qam->constellation[i][j] = (ip + I * qp);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex iq = symbols[k];
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
s[k * qam->N + n] = iq * cexp(2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM
|
||||
void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
r /= qam->N;
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation (lent)
|
||||
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// index du symbole (id) : même mappage que dans bits_to_symbols()
|
||||
int id = i_cl * sm + j_cl;
|
||||
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + b] = (id >> (qam->k - 1 - b)) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
// Compare deux tableaux de bits (0/1) et retourne le pourcentage de fiabilité.
|
||||
double compare_bits(const uint8_t *in_bits, size_t nb_bits_in, const uint8_t *out_bits, size_t nb_bits_out, size_t *erreurs) {
|
||||
if (!in_bits || !out_bits) {
|
||||
if (erreurs) *erreurs = (nb_bits_in < nb_bits_out) ? nb_bits_out : nb_bits_in;
|
||||
return 0.0;
|
||||
}
|
||||
|
||||
size_t n_min = (nb_bits_in < nb_bits_out) ? nb_bits_in : nb_bits_out;
|
||||
size_t n_max = (nb_bits_in > nb_bits_out) ? nb_bits_in : nb_bits_out;
|
||||
|
||||
size_t err = 0;
|
||||
for (size_t i = 0; i < n_min; ++i) {
|
||||
if ((in_bits[i] & 1) != (out_bits[i] & 1)) err++;
|
||||
}
|
||||
|
||||
if (n_max != n_min) {
|
||||
err += (n_max - n_min);
|
||||
}
|
||||
|
||||
if (erreurs) *erreurs = err;
|
||||
|
||||
double total_compared = (double)n_max;
|
||||
if (total_compared == 0.0) return 0.0;
|
||||
|
||||
double ber = (double)err / total_compared;
|
||||
double reliability_percent = (1.0 - ber) * 100.0;
|
||||
return reliability_percent;
|
||||
}
|
||||
|
||||
|
||||
int main (int argc, char *argv[]) {
|
||||
if (argc < 2) {
|
||||
fprintf(stderr, "Utilisation: %s <fichier_entree>\n", argv[0]);
|
||||
return 1;
|
||||
}
|
||||
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
qam.Ts = 0.0003;
|
||||
qam.N = (int)qam.Fs * qam.Ts;
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
printf("Lecture du fichier...\n");
|
||||
// Lecture du fichier et conversion en bits
|
||||
const char *input_filename = argv[1];
|
||||
bit_array input_bits = file_to_bits(input_filename);
|
||||
size_t nb_symbols = input_bits.nb_bits / qam.k;
|
||||
|
||||
printf("Mise en forme des symboles...\n");
|
||||
// Mise en forme des symboles
|
||||
double complex *symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
bits_to_symbols(&qam, input_bits.bits, input_bits.nb_bits, symbols);
|
||||
|
||||
printf("Modulation...\n");
|
||||
// Modulation QAM
|
||||
int total_samples = qam.N * nb_symbols;
|
||||
double complex* s = (double complex*)malloc(sizeof(double complex) * total_samples);
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
// Ajout du bruit
|
||||
double signal_power = (2.0/3.0)*(qam.M-1); // puissance moyenne avant échelle
|
||||
double snr_dB = 10; // Signal to noise ratio
|
||||
double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
double sigma = sqrt(signal_power / snr_lin);
|
||||
printf("Ajout du bruit... \n puissance du signal : %f\n SNR db : %f\n sigma : %f\n", signal_power, snr_dB, sigma);
|
||||
add_noise(s, total_samples, 0);
|
||||
|
||||
printf("Demodulation...\n");
|
||||
|
||||
// Demodulation QAM
|
||||
bit_array output_bits;
|
||||
output_bits.nb_bits = input_bits.nb_bits;
|
||||
output_bits.bits = (uint8_t*)malloc(output_bits.nb_bits * sizeof(uint8_t));
|
||||
demodulate(&qam, s, nb_symbols, output_bits.bits);
|
||||
|
||||
printf("Ecriture...\n");
|
||||
// Ecriture du fichier de Demodulation
|
||||
char *output_filename = make_output_filename(input_filename);
|
||||
bits_to_file(output_filename, &output_bits);
|
||||
|
||||
// Affichage du signal dans un .wav
|
||||
/*
|
||||
double* si = (double*)malloc(sizeof(double) * total_samples);
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
si[i] = cimag(s[i]);
|
||||
}
|
||||
write_wav("output.wav", si, total_samples);
|
||||
*/
|
||||
|
||||
size_t erreurs = 0;
|
||||
double fiabilite = compare_bits( input_bits.bits, input_bits.nb_bits, output_bits.bits, output_bits.nb_bits, &erreurs);
|
||||
|
||||
printf("Comparaison :\n");
|
||||
printf(" Bits d'entrée : %zu\n", input_bits.nb_bits);
|
||||
printf(" Bits de sortie: %zu\n", output_bits.nb_bits);
|
||||
printf(" Erreurs : %zu\n", erreurs);
|
||||
printf(" Fiabilité : %.4f %%\n", fiabilite);
|
||||
|
||||
// Libération mémoire
|
||||
free_bit_array(&input_bits);
|
||||
free_bit_array(&output_bits);
|
||||
free(symbols);
|
||||
free(s);
|
||||
free_constellation(&qam);
|
||||
free(output_filename);
|
||||
|
||||
return 0;
|
||||
}
|
||||
@ -1,237 +0,0 @@
|
||||
#include <math.h>
|
||||
#include <stdint.h>
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <complex.h>
|
||||
#include "../wav/wav.h"
|
||||
#include "../files/files.h"
|
||||
#include "../plot/plot_constellation.h"
|
||||
#include <string.h>
|
||||
#include <SDL2/SDL.h>
|
||||
|
||||
#define A 1
|
||||
|
||||
struct qam_system_s {
|
||||
int M; // Nombre de symboles M-QAM
|
||||
int k; // Nombre de bits/symboles
|
||||
double Fs; // Fréquence d'échantillionage
|
||||
double Ts; // Temps d'échantillionage
|
||||
int N; // Nombre d'échantillions
|
||||
double Fc; // Fréquence de la porteuse
|
||||
double complex** constellation; // Tableau de symboles I + j Q
|
||||
};
|
||||
typedef struct qam_system_s qam_system;
|
||||
|
||||
// Initialisation de la constellation (double tableau de taille sqrt(M)),
|
||||
// ToDo : changer à un tableau à 1 dimension pour éviter de calculer sqrt(M)
|
||||
void init_constellation (qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
qam->constellation = (double complex**)malloc(sizeof(double complex*) * sm);
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
qam->constellation[i] = (double complex*)malloc(sizeof(double complex) * sm);
|
||||
}
|
||||
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0); // Pour puissance unitaire
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
double complex ip = -(sm - 1) + 2 * i;
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex qp = -(sm - 1) + 2 * j;
|
||||
qam->constellation[i][j] = (ip + I * qp);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex iq = symbols[k];
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
s[k * qam->N + n] = iq * cexp(2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM
|
||||
void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
r /= qam->N;
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation
|
||||
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
/*
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0);
|
||||
double Ir = creal(r) * norm_factor / A;
|
||||
double Qr = cimag(r) * norm_factor / A;
|
||||
|
||||
int i = (int)round((Ir + (sm - 1)) / 2.0);
|
||||
int j = (int)round((Qr + (sm - 1)) / 2.0);
|
||||
|
||||
i = (i < 0) ? 0 : ((i >= sm) ? sm - 1 : i);
|
||||
j = (j < 0) ? 0 : ((j >= sm) ? sm - 1 : j);
|
||||
|
||||
int id = i * sm + j;
|
||||
*/
|
||||
|
||||
int id = i_cl * sm + j_cl;
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + (qam->k - 1 - b)] = (id >> b) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
double complex* demodulate_points(qam_system* qam, double complex* s, int nb_symbols) {
|
||||
double complex* points = malloc(sizeof(double complex) * nb_symbols);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
r /= qam->N;
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0);
|
||||
double Ir = creal(r);
|
||||
double Qr = cimag(r);
|
||||
|
||||
points[k] = Ir + I * Qr;
|
||||
}
|
||||
return points;
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
int main (int argc, char *argv[]) {
|
||||
if (argc < 2) {
|
||||
fprintf(stderr, "Utilisation: %s <fichier_entree>\n", argv[0]);
|
||||
return 1;
|
||||
}
|
||||
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
qam.Ts = 0.0003;
|
||||
qam.N = (int)qam.Fs * qam.Ts;
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
printf("Lecture du fichier...\n");
|
||||
// Lecture du fichier et conversion en bits
|
||||
const char *input_filename = argv[1];
|
||||
bit_array input_bits = file_to_bits(input_filename);
|
||||
size_t nb_symbols = input_bits.nb_bits / qam.k;
|
||||
|
||||
printf("Mise en forme des symboles...\n");
|
||||
// Mise en forme des symboles
|
||||
double complex *symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
bits_to_symbols(&qam, input_bits.bits, input_bits.nb_bits, symbols);
|
||||
|
||||
printf("Modulation...\n");
|
||||
// Modulation QAM
|
||||
int total_samples = qam.N * nb_symbols;
|
||||
double complex* s = (double complex*)malloc(sizeof(double complex) * total_samples);
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
// Ajout du bruit
|
||||
double signal_power = (2.0/3.0)*(qam.M-1); // puissance moyenne avant échelle
|
||||
double snr_dB = 10; // Signal to noise ratio
|
||||
double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
double sigma = sqrt(signal_power / snr_lin);
|
||||
printf("Ajout du bruit... \n puissance du signal : %f\n SNR db : %f\n sigma : %f\n", signal_power, snr_dB, sigma);
|
||||
add_noise(s, total_samples, sigma);
|
||||
|
||||
printf("Demodulation...\n");
|
||||
|
||||
// Demodulation QAM
|
||||
bit_array output_bits;
|
||||
output_bits.nb_bits = input_bits.nb_bits;
|
||||
output_bits.bits = (uint8_t*)malloc(output_bits.nb_bits);
|
||||
demodulate(&qam, s, nb_symbols, output_bits.bits);
|
||||
|
||||
printf("Ecriture...\n");
|
||||
// Ecriture du fichier de Demodulation
|
||||
char *output_filename = make_output_filename(input_filename);
|
||||
bits_to_file(output_filename, &output_bits);
|
||||
|
||||
|
||||
// Affichage du signal dans un .wav
|
||||
double* si = (double*)malloc(sizeof(double) * total_samples);
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
si[i] = cimag(s[i]);
|
||||
}
|
||||
write_wav("output.wav", si, total_samples);
|
||||
|
||||
// Plot
|
||||
printf("Ploting...");
|
||||
plot_t plot;
|
||||
plot_init(&plot, 1400, 1400, 0.04);
|
||||
plot_draw_constellation(&plot, qam.constellation, qam.M);
|
||||
SDL_Color red = {255, 0, 0, 255};
|
||||
plot_draw_points_animated(&plot, demodulate_points(&qam, s, nb_symbols), nb_symbols, red, 5);
|
||||
|
||||
// Libération mémoire
|
||||
plot_close(&plot);
|
||||
free_bit_array(&input_bits);
|
||||
free_bit_array(&output_bits);
|
||||
free(symbols);
|
||||
free(s);
|
||||
free_constellation(&qam);
|
||||
free(output_filename);
|
||||
|
||||
return 0;
|
||||
}
|
||||
@ -1,324 +0,0 @@
|
||||
#include <math.h>
|
||||
#include <stdint.h>
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <complex.h>
|
||||
#include "../wav/wav.h"
|
||||
#include "../files/files.h"
|
||||
#include <string.h>
|
||||
|
||||
#define A 1
|
||||
|
||||
struct qam_system_s {
|
||||
int M; // Nombre de symboles M-QAM
|
||||
int k; // Nombre de bits/symboles
|
||||
double Fs; // Fréquence d'échantillionage
|
||||
double Ts; // Temps d'échantillionage
|
||||
int N; // Nombre d'échantillions
|
||||
double Fc; // Fréquence de la porteuse
|
||||
double complex** constellation; // Tableau de symboles I + j Q
|
||||
};
|
||||
typedef struct qam_system_s qam_system;
|
||||
|
||||
// Initialisation de la constellation (double tableau de taille sqrt(M)),
|
||||
// ToDo : changer à un tableau à 1 dimension pour éviter de calculer sqrt(M)
|
||||
void init_constellation (qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
qam->constellation = (double complex**)malloc(sizeof(double complex*) * sm);
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
qam->constellation[i] = (double complex*)malloc(sizeof(double complex) * sm);
|
||||
}
|
||||
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0); // Pour puissance unitaire
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
double complex ip = -(sm - 1) + 2 * i;
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex qp = -(sm - 1) + 2 * j;
|
||||
qam->constellation[i][j] = (ip + I * qp);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex iq = symbols[k];
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
s[k * qam->N + n] = iq * cexp(2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM
|
||||
void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
//int i = k * qam->N + n;
|
||||
//r += s[i] * cexp(-2 * I * M_PI * qam->Fc * ((double)i / qam->Fs));
|
||||
}
|
||||
r /= qam->N;
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation
|
||||
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// index du symbole (id) : même mappage que dans bits_to_symbols()
|
||||
int id = i_cl * sm + j_cl;
|
||||
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + b] = (id >> (qam->k - 1 - b)) & 1;
|
||||
}
|
||||
|
||||
/*
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0);
|
||||
double Ir = creal(r) * norm_factor / A;
|
||||
double Qr = cimag(r) * norm_factor / A;
|
||||
|
||||
int i = (int)round((Ir + (sm - 1)) / 2.0);
|
||||
int j = (int)round((Qr + (sm - 1)) / 2.0);
|
||||
|
||||
i = (i < 0) ? 0 : ((i >= sm) ? sm - 1 : i);
|
||||
j = (j < 0) ? 0 : ((j >= sm) ? sm - 1 : j);
|
||||
|
||||
int id = i * sm + j;
|
||||
|
||||
int id = i_cl * sm + j_cl;
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + (qam->k - 1 - b)] = (id >> b) & 1;
|
||||
}
|
||||
*/
|
||||
}
|
||||
}
|
||||
|
||||
double complex* demodulate_points(qam_system* qam, double complex* s, int nb_symbols) {
|
||||
double complex* points = malloc(sizeof(double complex) * nb_symbols);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
int i = k * qam->N + n;
|
||||
r += s[i] * cexp(-2 * I * M_PI * qam->Fc * ((double)i / qam->Fs));
|
||||
}
|
||||
r /= qam->N;
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0);
|
||||
double Ir = creal(r);
|
||||
double Qr = cimag(r);
|
||||
|
||||
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
double complex p = qam->constellation[i_cl][j_cl];
|
||||
points[k] = creal(p) + I * cimag(p);
|
||||
|
||||
//points[k] = (int)Ir + I * (int)Qr;
|
||||
}
|
||||
return points;
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
#include <stddef.h> // pour size_t
|
||||
|
||||
// Compare deux tableaux de bits (0/1) et retourne le pourcentage de fiabilité.
|
||||
double compare_bits_and_get_reliability(const uint8_t *in_bits, size_t nb_bits_in, const uint8_t *out_bits, size_t nb_bits_out, size_t *erreurs) {
|
||||
if (!in_bits || !out_bits) {
|
||||
if (erreurs) *erreurs = (nb_bits_in < nb_bits_out) ? nb_bits_out : nb_bits_in;
|
||||
return 0.0;
|
||||
}
|
||||
|
||||
size_t n_min = (nb_bits_in < nb_bits_out) ? nb_bits_in : nb_bits_out;
|
||||
size_t n_max = (nb_bits_in > nb_bits_out) ? nb_bits_in : nb_bits_out;
|
||||
|
||||
size_t err = 0;
|
||||
for (size_t i = 0; i < n_min; ++i) {
|
||||
if ((in_bits[i] & 1) != (out_bits[i] & 1)) err++;
|
||||
}
|
||||
|
||||
if (n_max != n_min) {
|
||||
err += (n_max - n_min);
|
||||
}
|
||||
|
||||
if (erreurs) *erreurs = err;
|
||||
|
||||
double total_compared = (double)n_max;
|
||||
if (total_compared == 0.0) return 0.0;
|
||||
|
||||
double ber = (double)err / total_compared;
|
||||
double reliability_percent = (1.0 - ber) * 100.0;
|
||||
return reliability_percent;
|
||||
}
|
||||
|
||||
|
||||
void affiche_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex p = qam->constellation[i][j];
|
||||
printf("(%d,%d) ", (int)creal(p), (int)cimag(p));
|
||||
}
|
||||
printf("\n");
|
||||
}
|
||||
}
|
||||
|
||||
void affiche_points(double complex* r, int len, int len_samples) {
|
||||
for (int i = 0; i < len; i += len_samples) {
|
||||
for (int j = 0; j < len_samples; j++) {
|
||||
double complex p = r[i + j];
|
||||
printf("(%f,%f) ", (float)creal(p), (float)cimag(p));
|
||||
}
|
||||
printf("\n");
|
||||
}
|
||||
}
|
||||
|
||||
int main (int argc, char *argv[]) {
|
||||
if (argc < 2) {
|
||||
fprintf(stderr, "Utilisation: %s <fichier_entree>\n", argv[0]);
|
||||
return 1;
|
||||
}
|
||||
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
qam.Ts = 0.0003;
|
||||
qam.N = (int)qam.Fs * qam.Ts;
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
printf("Lecture du fichier...\n");
|
||||
// Lecture du fichier et conversion en bits
|
||||
const char *input_filename = argv[1];
|
||||
bit_array input_bits = file_to_bits(input_filename);
|
||||
size_t nb_symbols = input_bits.nb_bits / qam.k;
|
||||
|
||||
printf("Mise en forme des symboles...\n");
|
||||
// Mise en forme des symboles
|
||||
double complex *symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
bits_to_symbols(&qam, input_bits.bits, input_bits.nb_bits, symbols);
|
||||
|
||||
printf("Modulation...\n");
|
||||
// Modulation QAM
|
||||
int total_samples = qam.N * nb_symbols;
|
||||
double complex* s = (double complex*)malloc(sizeof(double complex) * total_samples);
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
// Ajout du bruit
|
||||
double signal_power = (2.0/3.0)*(qam.M-1); // puissance moyenne avant échelle
|
||||
double snr_dB = 10; // Signal to noise ratio
|
||||
double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
double sigma = sqrt(signal_power / snr_lin);
|
||||
printf("Ajout du bruit... \n puissance du signal : %f\n SNR db : %f\n sigma : %f\n", signal_power, snr_dB, sigma);
|
||||
add_noise(s, total_samples, 5);
|
||||
|
||||
printf("Demodulation...\n");
|
||||
|
||||
// Demodulation QAM
|
||||
bit_array output_bits;
|
||||
output_bits.nb_bits = input_bits.nb_bits;
|
||||
output_bits.bits = (uint8_t*)malloc(output_bits.nb_bits);
|
||||
demodulate(&qam, s, nb_symbols, output_bits.bits);
|
||||
|
||||
printf("Ecriture...\n");
|
||||
// Ecriture du fichier de Demodulation
|
||||
char *output_filename = make_output_filename(input_filename);
|
||||
bits_to_file(output_filename, &output_bits);
|
||||
|
||||
//printf("Constelattion :\n");
|
||||
//affiche_constellation(&qam);
|
||||
//printf("Points de demodulation :\n");
|
||||
//affiche_points(demodulate_points(&qam, s, nb_symbols), nb_symbols, qam.N);
|
||||
|
||||
// Affichage du signal dans un .wav
|
||||
/*
|
||||
double* si = (double*)malloc(sizeof(double) * total_samples);
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
si[i] = cimag(s[i]);
|
||||
}
|
||||
write_wav("output.wav", si, total_samples);
|
||||
*/
|
||||
|
||||
|
||||
size_t erreurs = 0;
|
||||
double fiabilite = compare_bits_and_get_reliability( input_bits.bits, input_bits.nb_bits, output_bits.bits, output_bits.nb_bits, &erreurs);
|
||||
|
||||
printf("Résultat de la comparaison :\n");
|
||||
printf(" Bits d'entrée : %zu\n", input_bits.nb_bits);
|
||||
printf(" Bits de sortie: %zu\n", output_bits.nb_bits);
|
||||
printf(" Erreurs : %zu\n", erreurs);
|
||||
printf(" Fiabilité : %.4f %%\n", fiabilite);
|
||||
|
||||
// Libération mémoire
|
||||
free_bit_array(&input_bits);
|
||||
free_bit_array(&output_bits);
|
||||
free(symbols);
|
||||
free(s);
|
||||
free_constellation(&qam);
|
||||
free(output_filename);
|
||||
|
||||
return 0;
|
||||
}
|
||||
@ -1,184 +0,0 @@
|
||||
#include <math.h>
|
||||
#include <stdint.h>
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <complex.h>
|
||||
#include <string.h>
|
||||
|
||||
#define A 1
|
||||
|
||||
struct qam_system_s {
|
||||
int M; // Nombre de symboles M-QAM
|
||||
int k; // Nombre de bits/symboles
|
||||
double Fs; // Fréquence d'échantillionage
|
||||
double Ts; // Temps d'échantillionage
|
||||
int N; // Nombre d'échantillions
|
||||
double Fc; // Fréquence de la porteuse
|
||||
double complex** constellation; // Tableau de symboles I + j Q
|
||||
};
|
||||
typedef struct qam_system_s qam_system;
|
||||
|
||||
// Initialisation de la constellation (double tableau de taille sqrt(M)),
|
||||
// ToDo : changer à un tableau à 1 dimension pour éviter de calculer sqrt(M)
|
||||
void init_constellation (qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
qam->constellation = (double complex**)malloc(sizeof(double complex*) * sm);
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
qam->constellation[i] = (double complex*)malloc(sizeof(double complex) * sm);
|
||||
}
|
||||
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0); // Pour puissance unitaire
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
double complex ip = -(sm - 1) + 2 * i;
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex qp = -(sm - 1) + 2 * j;
|
||||
qam->constellation[i][j] = (ip + I * qp);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex iq = symbols[k];
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
s[k * qam->N + n] = iq * cexp(2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM
|
||||
void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
r /= qam->N;
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation (lent)
|
||||
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// index du symbole (id) : même mappage que dans bits_to_symbols()
|
||||
int id = i_cl * sm + j_cl;
|
||||
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + b] = (id >> (qam->k - 1 - b)) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
double compare_bits(uint8_t* bits1, uint8_t* bits2, int nb_bits) {
|
||||
int errors = 0;
|
||||
for (int i = 0; i < nb_bits; i++) {
|
||||
if (bits1[i] != bits2[i]) errors++;
|
||||
}
|
||||
return (double)errors / nb_bits;
|
||||
}
|
||||
|
||||
int main () {
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
qam.Ts = 0.0003;
|
||||
qam.N = (int)qam.Fs * qam.Ts;
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
int nb_bits = 1000;
|
||||
int nb_symbols = nb_bits / qam.k;
|
||||
|
||||
uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
for (int i = 0; i < nb_bits; i++) {
|
||||
input_bits[i] = rand() % 2;
|
||||
}
|
||||
|
||||
// Conversion en symboles
|
||||
double complex* symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
bits_to_symbols(&qam, input_bits, nb_bits, symbols);
|
||||
|
||||
// Modulation
|
||||
int total_samples = qam.N * nb_symbols;
|
||||
double complex* s = malloc(sizeof(double complex) * total_samples);
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
// Ajout du bruit
|
||||
double signal_power = (2.0/3.0)*(qam.M-1); // puissance moyenne
|
||||
double snr_dB = 10; // SNR en dB
|
||||
double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
double sigma = sqrt(signal_power / snr_lin);
|
||||
add_noise(s, total_samples, sigma);
|
||||
|
||||
// Démodulation
|
||||
uint8_t* output_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
demodulate(&qam, s, nb_symbols, output_bits);
|
||||
|
||||
// Calcul du taux d'erreur
|
||||
double ber = compare_bits(input_bits, output_bits, nb_bits);
|
||||
printf("Taux d'erreur : %.4f\n", ber);
|
||||
|
||||
// Libération mémoire
|
||||
free(input_bits);
|
||||
free(output_bits);
|
||||
free(symbols);
|
||||
free(s);
|
||||
free_constellation(&qam);
|
||||
|
||||
return 0;
|
||||
}
|
||||
@ -1,328 +0,0 @@
|
||||
#include <math.h>
|
||||
#include <stdint.h>
|
||||
#include <stdio.h>
|
||||
#include <stdlib.h>
|
||||
#include <complex.h>
|
||||
#include "../files/files.h"
|
||||
#include <string.h>
|
||||
|
||||
#define A 1
|
||||
|
||||
struct qam_system_s {
|
||||
int M; // Nombre de symboles M-QAM
|
||||
int k; // Nombre de bits/symboles
|
||||
double Fs; // Fréquence d'échantillionage
|
||||
double Ts; // Temps d'échantillionage
|
||||
int N; // Nombre d'échantillions
|
||||
double Fc; // Fréquence de la porteuse
|
||||
double complex** constellation; // Tableau de symboles I + j Q
|
||||
};
|
||||
typedef struct qam_system_s qam_system;
|
||||
|
||||
// Initialisation de la constellation (double tableau de taille sqrt(M)),
|
||||
// ToDo : changer à un tableau à 1 dimension pour éviter de calculer sqrt(M)
|
||||
void init_constellation (qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
qam->constellation = (double complex**)malloc(sizeof(double complex*) * sm);
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
qam->constellation[i] = (double complex*)malloc(sizeof(double complex) * sm);
|
||||
}
|
||||
|
||||
double norm_factor = sqrt((double)(qam->M - 1) / 3.0); // Pour puissance unitaire
|
||||
|
||||
for (int i = 0; i < sm; i++) {
|
||||
double complex ip = -(sm - 1) + 2 * i;
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double complex qp = -(sm - 1) + 2 * j;
|
||||
qam->constellation[i][j] = (ip + I * qp);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex iq = symbols[k];
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
s[k * qam->N + n] = iq * cexp(2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM
|
||||
void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n] * cexp(-2 * I * M_PI * qam->Fc * ((double)n / qam->Fs));
|
||||
}
|
||||
r /= qam->N;
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation (lent)
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// index du symbole (id) : même mappage que dans bits_to_symbols()
|
||||
int id = i_cl * sm + j_cl;
|
||||
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + b] = (id >> (qam->k - 1 - b)) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Demodulation QAM avec porteuse estimé
|
||||
void demodulate2(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bits_hat) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < qam->N; n++) {
|
||||
r += s[k * qam->N + n];
|
||||
}
|
||||
r /= qam->N;
|
||||
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation (lent)
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
if (d < min_d) {
|
||||
min_d = d;
|
||||
i_cl = i;
|
||||
j_cl = j;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// index du symbole (id) : même mappage que dans bits_to_symbols()
|
||||
int id = i_cl * sm + j_cl;
|
||||
|
||||
for (int b = 0; b < qam->k; b++) {
|
||||
bits_hat[k * qam->k + b] = (id >> (qam->k - 1 - b)) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
/*
|
||||
// Pour j = 1 dans la def de wiki sur l'autocorrelation
|
||||
double fc_autocorrelation_1_rad(double complex* s, int N) {
|
||||
double complex r = 0;
|
||||
for (int n = 0; n < N - 1; n++) {
|
||||
r += s[n] * conj(s[n + 1]);
|
||||
}
|
||||
r /= (N - 1);
|
||||
double om = -carg(r); // rad / sample
|
||||
return om;
|
||||
}
|
||||
|
||||
// Autocorrelation pour trouver la Fc (fréquence de la porteuse)
|
||||
double fc_autocorrelation_multilag_rad(double complex* s, int N, int max_lag) {
|
||||
double sum_phase = 0;
|
||||
for (int k = 1; k <= max_lag; k++) {
|
||||
double complex rk = 0;
|
||||
for (int n = 0; n < N - k; n++) {
|
||||
rk += s[n] * conj(s[n + k]);
|
||||
}
|
||||
rk /= (double)(N - k);
|
||||
sum_phase += -carg(rk) / k;
|
||||
}
|
||||
return sum_phase / max_lag; // rad/sample
|
||||
}
|
||||
|
||||
// Correction du signal pour enlever l'offset de la porteuse
|
||||
void signal_correction(double complex* s, int total_samples, double om_hat) {
|
||||
for (int n = 0; n < total_samples; n++) {
|
||||
s[n] *= cexp(-I * om_hat * n);
|
||||
}
|
||||
}
|
||||
*/
|
||||
|
||||
void blind_carrier(double complex* s, int N, int M) {
|
||||
int k = (int)sqrt(M); // puissance à élever
|
||||
double complex sum = 0;
|
||||
for(int n = 0; n < N; n++)
|
||||
sum += cpow(s[n], k);
|
||||
|
||||
double phi_hat = carg(sum) / k;
|
||||
|
||||
for(int n = 0; n < N; n++)
|
||||
s[n] *= cexp(-I * phi_hat);
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
// Compare deux tableaux de bits (0/1) et retourne le pourcentage de fiabilité.
|
||||
double taux_erreur_bits(const uint8_t *in_bits, size_t nb_bits_in,
|
||||
const uint8_t *out_bits, size_t nb_bits_out) {
|
||||
if (!in_bits || !out_bits)
|
||||
return 1.0;
|
||||
|
||||
size_t n_min = nb_bits_in < nb_bits_out ? nb_bits_in : nb_bits_out;
|
||||
size_t n_max = nb_bits_in > nb_bits_out ? nb_bits_in : nb_bits_out;
|
||||
size_t err = n_max - n_min;
|
||||
|
||||
for (size_t i = 0; i < n_min; i++)
|
||||
err += ((in_bits[i] ^ out_bits[i]) & 1);
|
||||
|
||||
return n_max ? (double)err / n_max : 0.0;
|
||||
}
|
||||
|
||||
|
||||
int main (int argc, char *argv[]) {
|
||||
if (argc < 2) {
|
||||
fprintf(stderr, "Utilisation: %s <fichier_entree>\n", argv[0]);
|
||||
return 1;
|
||||
}
|
||||
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
qam.Ts = 0.3;
|
||||
qam.N = (int)qam.Fs * qam.Ts;
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
printf("Lecture du fichier...\n");
|
||||
// Lecture du fichier et conversion en bits
|
||||
const char *input_filename = argv[1];
|
||||
bit_array input_bits = file_to_bits(input_filename);
|
||||
size_t nb_symbols = input_bits.nb_bits / qam.k;
|
||||
|
||||
printf("Mise en forme des symboles...\n");
|
||||
// Mise en forme des symboles
|
||||
double complex *symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
bits_to_symbols(&qam, input_bits.bits, input_bits.nb_bits, symbols);
|
||||
|
||||
printf("Modulation...\n");
|
||||
// Modulation QAM
|
||||
int total_samples = qam.N * nb_symbols;
|
||||
double complex* s = (double complex*)malloc(sizeof(double complex) * total_samples);
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
// Ajout du bruit
|
||||
double signal_power = (2.0/3.0)*(qam.M-1); // puissance moyenne avant échelle
|
||||
double snr_dB = 10; // Signal to noise ratio
|
||||
double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
double sigma = sqrt(signal_power / snr_lin);
|
||||
printf("Ajout du bruit... \n puissance du signal : %f\n SNR db : %f\n sigma : %f\n", signal_power, snr_dB, sigma);
|
||||
add_noise(s, total_samples, 0);
|
||||
|
||||
|
||||
// Demodulation QAM
|
||||
//printf("Demodulation...\n");
|
||||
//bit_array output_bits;
|
||||
//output_bits.nb_bits = input_bits.nb_bits;
|
||||
//output_bits.bits = (uint8_t*)malloc(output_bits.nb_bits * sizeof(uint8_t));
|
||||
//demodulate(&qam, s, nb_symbols, output_bits.bits);
|
||||
|
||||
//printf("Ecriture...\n");
|
||||
// Ecriture du fichier de Demodulation
|
||||
//char *output_filename = make_output_filename(input_filename);
|
||||
//bits_to_file(output_filename, &output_bits);
|
||||
|
||||
//double erreurs = taux_erreur_bits(input_bits.bits, input_bits.nb_bits, output_bits.bits, output_bits.nb_bits);
|
||||
|
||||
//printf("Comparaison :\n");
|
||||
//printf(" Erreurs : %f\n", erreurs * 100);
|
||||
|
||||
/*
|
||||
double om_hat = fc_autocorrelation_1_rad(s, total_samples);
|
||||
printf("Estimation de fc 1: %f\n", (om_hat * qam.Fs) / (2 * M_PI));
|
||||
|
||||
om_hat = fc_autocorrelation_multilag_rad(s, total_samples, 10);
|
||||
printf("Estimation de fc multilag: %f\n", (om_hat * qam.Fs) / (2 * M_PI));
|
||||
|
||||
// Correction du signal avec Fc (en rad/sample) trouvé
|
||||
signal_correction(s, total_samples, om_hat);
|
||||
|
||||
printf("Demodulation...\n");
|
||||
demodulate2(&qam, s, nb_symbols, output_bits.bits);
|
||||
*/
|
||||
|
||||
// Avant la demodulation
|
||||
printf("Test de blind carrier correction...\n");
|
||||
|
||||
// Paramètres CMA simples
|
||||
double mu = 0.001; // pas d'adaptation
|
||||
int num_iter = 100; // nombre d'itérations
|
||||
|
||||
// Appel de la fonction blind carrier correction
|
||||
blind_carrier(s, total_samples, qam.M);
|
||||
|
||||
// Ensuite, utilise ton démodulateur actuel
|
||||
bit_array output_bits;
|
||||
output_bits.nb_bits = input_bits.nb_bits;
|
||||
output_bits.bits = (uint8_t*)malloc(output_bits.nb_bits * sizeof(uint8_t));
|
||||
|
||||
demodulate2(&qam, s, nb_symbols, output_bits.bits);
|
||||
|
||||
// Comparaison avec bits originaux
|
||||
double erreurs = taux_erreur_bits(input_bits.bits, input_bits.nb_bits, output_bits.bits, output_bits.nb_bits);
|
||||
printf("Taux d'erreur après blind carrier correction: %f %%\n", erreurs * 100);
|
||||
|
||||
// Ecriture du fichier de Demodulation
|
||||
printf("Ecriture...\n");
|
||||
char *output_filename = make_output_filename(input_filename);
|
||||
bits_to_file(output_filename, &output_bits);
|
||||
|
||||
// Libération mémoire
|
||||
free_bit_array(&input_bits);
|
||||
free_bit_array(&output_bits);
|
||||
free(symbols);
|
||||
free(s);
|
||||
free_constellation(&qam);
|
||||
free(output_filename);
|
||||
|
||||
return 0;
|
||||
}
|
||||
@ -1,416 +0,0 @@
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
1.34164079 0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -1.34164079
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 0.44721360
|
||||
0.44721360 1.34164079
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 -1.34164079
|
||||
-0.44721360 0.44721360
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-1.34164079 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -1.34164079
|
||||
@ -1,16 +0,0 @@
|
||||
-1.34164079 -1.34164079
|
||||
-1.34164079 -0.44721360
|
||||
-1.34164079 0.44721360
|
||||
-1.34164079 1.34164079
|
||||
-0.44721360 -1.34164079
|
||||
-0.44721360 -0.44721360
|
||||
-0.44721360 0.44721360
|
||||
-0.44721360 1.34164079
|
||||
0.44721360 -1.34164079
|
||||
0.44721360 -0.44721360
|
||||
0.44721360 0.44721360
|
||||
0.44721360 1.34164079
|
||||
1.34164079 -1.34164079
|
||||
1.34164079 -0.44721360
|
||||
1.34164079 0.44721360
|
||||
1.34164079 1.34164079
|
||||
1340
QAM/pll_error.dat
1340
QAM/pll_error.dat
File diff suppressed because it is too large
Load Diff
@ -24,7 +24,7 @@ def load_error(filename):
|
||||
|
||||
# -------------------- Application --------------------
|
||||
app = QtWidgets.QApplication(sys.argv)
|
||||
win = pg.GraphicsLayoutWidget(show=True, title="Constellation et PLL Error")
|
||||
win = pg.GraphicsLayoutWidget(show=True, title="Constellation")
|
||||
win.resize(900, 900)
|
||||
win.setBackground('#0a0a0a') # fond noir profond
|
||||
|
||||
@ -39,38 +39,6 @@ void init_constellation (qam_system* qam) {
|
||||
}
|
||||
}
|
||||
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
// Modulation QAM
|
||||
void modulate (qam_system* qam, double complex* symbols, int nb_symbols, double complex* s) {
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
@ -99,7 +67,8 @@ void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bit
|
||||
// Distance euclidien de Ir et Qr pour avoir le point le plus proche de la constellation (lent)
|
||||
int sm = (int)sqrt(qam->M);
|
||||
double min_d = INFINITY;
|
||||
int i_cl, j_cl = 0;
|
||||
int i_cl = 0;
|
||||
int j_cl = 0;
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
double d = cabs(r - qam->constellation[i][j]);
|
||||
@ -120,12 +89,40 @@ void demodulate(qam_system* qam, double complex* s, int nb_symbols, uint8_t* bit
|
||||
}
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
// Calcul du bruit gaussien pour un sigma donné
|
||||
// Formule de Box-Muller
|
||||
double gaussian_noise (double sigma) {
|
||||
double u1 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
double u2 = (rand() + 1) / ((double)RAND_MAX + 2);
|
||||
return sigma * sqrt(-2 * log(u1)) * cos(2 * M_PI * u2);
|
||||
}
|
||||
|
||||
// Ajout du bruit
|
||||
void add_noise (double complex* s, int len, double sigma) {
|
||||
//double signal_power = (2.0/3.0)*(qam.M-1);
|
||||
//double snr_dB = 5; // SNR en dB
|
||||
//double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
//double sigma = sqrt(signal_power / snr_lin);
|
||||
for (int i = 0; i < len; i++) {
|
||||
double nr = gaussian_noise(sigma);
|
||||
double ni = gaussian_noise(sigma);
|
||||
s[i] += nr + I * ni;
|
||||
}
|
||||
}
|
||||
|
||||
// Changer le tableau de bits en boolen ou alors la represenation binaire et shifter pour extraire les bits (pas bien si M plus grand)
|
||||
void bits_to_symbols (qam_system* qam, uint8_t* bits, int nb_bits, double complex* symbols) {
|
||||
int nb_symbols = nb_bits / qam->k;
|
||||
int sm = sqrt(qam->M);
|
||||
for (int k = 0; k < nb_symbols; k++) {
|
||||
int id = 0;
|
||||
for (int b = 0 ; b < qam->k; b++) {
|
||||
id = id * 2 + bits[k * qam->k + b];
|
||||
}
|
||||
int i = id / sm;
|
||||
int j = id % sm;
|
||||
symbols[k] = qam->constellation[i][j];
|
||||
}
|
||||
}
|
||||
|
||||
double compare_bits(uint8_t* bits1, uint8_t* bits2, int nb_bits) {
|
||||
@ -136,35 +133,72 @@ double compare_bits(uint8_t* bits1, uint8_t* bits2, int nb_bits) {
|
||||
return (double)errors / nb_bits;
|
||||
}
|
||||
|
||||
int main () {
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
qam.Ts = 0.0003;
|
||||
qam.N = (int)qam.Fs * qam.Ts;
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
void add_dephasage(double complex* s, double phi_offset, int total_samples) {
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
s[i] *= cexp(I * phi_offset);
|
||||
}
|
||||
}
|
||||
|
||||
//int nb_bits = 1000;
|
||||
//int nb_symbols = nb_bits / qam.k;
|
||||
void add_freq(qam_system* qam, double complex* s, double freq_offset, int total_samples) {
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
double t = (double)i / qam->Fs;
|
||||
s[i] *= cexp(I * 2 * M_PI * freq_offset * t);
|
||||
}
|
||||
}
|
||||
|
||||
//uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
//for (int i = 0; i < nb_bits; i++) {
|
||||
// input_bits[i] = rand() % 2;
|
||||
//}
|
||||
char* texte = "Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, ";
|
||||
int nb_chars = strlen(texte);
|
||||
int nb_bits = nb_chars * 8;
|
||||
int nb_symbols = (nb_bits + qam.k - 1) / qam.k;
|
||||
void fill_constellation_data(qam_system* qam, FILE *fp_ref) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
fprintf(fp_ref, "% .8f % .8f\n", creal(qam->constellation[i][j]), cimag(qam->constellation[i][j]));
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Conversion du texte en bits
|
||||
uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
void reconstruction_text(int nb_chars, uint8_t* output_bits, char* texte_recup) {
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
char c = 0;
|
||||
for(int b = 0; b < 8; b++){
|
||||
c |= output_bits[i*8 + b] << (7-b);
|
||||
}
|
||||
texte_recup[i] = c;
|
||||
}
|
||||
texte_recup[nb_chars] = '\0';
|
||||
}
|
||||
|
||||
void conversion_text_to_bits(int nb_chars, int nb_bits, char* texte, uint8_t* input_bits) {
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
for(int b = 0; b < 8; b++){
|
||||
input_bits[i*8 + b] = (texte[i] >> (7-b)) & 1;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Libération de la mémoire
|
||||
void free_constellation(qam_system* qam) {
|
||||
int sm = (int)sqrt(qam->M);
|
||||
for (int i = 0; i < sm; i++)
|
||||
free(qam->constellation[i]);
|
||||
free(qam->constellation);
|
||||
}
|
||||
|
||||
int main () {
|
||||
qam_system qam;
|
||||
qam.M = 16;
|
||||
qam.k = (int)log2((double)(qam.M));
|
||||
qam.Fs = 44100;
|
||||
qam.Ts = 0.01;
|
||||
qam.N = (int)(qam.Fs * qam.Ts);
|
||||
qam.Fc = 2000;
|
||||
init_constellation(&qam);
|
||||
|
||||
// Conversion du texte en bits
|
||||
char* texte = "Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux, Vif juge, trempez ce blond whisky aqueux";
|
||||
int nb_chars = strlen(texte);
|
||||
int nb_bits = nb_chars * 8;
|
||||
int nb_symbols = (nb_bits + qam.k - 1) / qam.k;
|
||||
uint8_t* input_bits = malloc(nb_bits * sizeof(uint8_t));
|
||||
conversion_text_to_bits(nb_chars, nb_bits, texte, input_bits);
|
||||
|
||||
// Conversion en symboles
|
||||
double complex* symbols = malloc(sizeof(double complex) * nb_symbols);
|
||||
@ -176,53 +210,29 @@ int main () {
|
||||
modulate(&qam, symbols, nb_symbols, s);
|
||||
|
||||
// Ajout du bruit
|
||||
double signal_power = (2.0/3.0)*(qam.M-1); // puissance moyenne
|
||||
double snr_dB = 5; // SNR en dB
|
||||
double snr_lin = pow(10.0, snr_dB / 10.0);
|
||||
double sigma = sqrt(signal_power / snr_lin);
|
||||
add_noise(s, total_samples, sigma);
|
||||
|
||||
add_noise(s, total_samples, 0);
|
||||
|
||||
FILE *fp_ref = fopen("constellation_ref.dat", "w");
|
||||
int sm = (int)sqrt(qam.M);
|
||||
for (int i = 0; i < sm; i++) {
|
||||
for (int j = 0; j < sm; j++) {
|
||||
fprintf(fp_ref, "% .8f % .8f\n", creal(qam.constellation[i][j]), cimag(qam.constellation[i][j]));
|
||||
}
|
||||
}
|
||||
fill_constellation_data(&qam, fp_ref);
|
||||
fclose(fp_ref);
|
||||
FILE *fp_constel = fopen("constellation.dat", "w");
|
||||
|
||||
double phase_offset = M_PI / 6.0; // 30 degrés
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
s[i] *= cexp(I * phase_offset);
|
||||
}
|
||||
double freq_offset = 0; // Hz de décalage
|
||||
for (int i = 0; i < total_samples; i++) {
|
||||
double t = (double)i / qam.Fs;
|
||||
s[i] *= cexp(I * 2 * M_PI * freq_offset * t);
|
||||
}
|
||||
//int offset_samples = (int)(0.3 * qam.N); // décalage de 30% d’un symbole
|
||||
//memmove(s + offset_samples, s, (total_samples - offset_samples) * sizeof(double complex));
|
||||
|
||||
// Ajout de dephasage
|
||||
//add_dephasage(s, M_PI / 6.0, total_samples);
|
||||
|
||||
// AJout de decalage de fréquence
|
||||
//add_freq(&qam, s, 1, total_samples);
|
||||
|
||||
// Démodulation
|
||||
FILE *fp_constel = fopen("constellation.dat", "w");
|
||||
uint8_t* output_bits = (uint8_t*)malloc(nb_bits * sizeof(uint8_t));
|
||||
demodulate(&qam, s, nb_symbols, output_bits, fp_constel);
|
||||
|
||||
fclose(fp_constel);
|
||||
|
||||
// Reconstruction du texte
|
||||
char* texte_recup = malloc(nb_chars + 1);
|
||||
for(int i = 0; i < nb_chars; i++){
|
||||
char c = 0;
|
||||
for(int b = 0; b < 8; b++){
|
||||
c |= output_bits[i*8 + b] << (7-b);
|
||||
}
|
||||
texte_recup[i] = c;
|
||||
}
|
||||
texte_recup[nb_chars] = '\0';
|
||||
printf("Texte original : %s\n", texte);
|
||||
reconstruction_text(nb_chars, output_bits, texte_recup);
|
||||
|
||||
printf("Texte original : %s\n\n", texte);
|
||||
printf("Texte demodulé : %s\n", texte_recup);
|
||||
|
||||
// Calcul du BER
|
||||
Reference in New Issue
Block a user