simu_telecom_numerique/generate_fig_recepteur.m at master · PPaques/simu_telecom_numerique · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
%% generate_fig_emetteur : genere toutes les figures du recepteur
% Script de simulation d'une chaine complete de telecommunications
% Par Nicolas Segui et Pierre Paques
% 2012-2013

%% filtrage
% Affichage de la reponse en frequence des filtres %%OK!!
figure(30)
plot(f_freqs,(abs(filtres_coeff)));
title(['Reponse en frequence des filtres generes']);
axis([0 (2*N/T_b) 0 1]); grid;
xlabel('Frequence (Hz)');
ylabel('Amplitude du spectre') ;

% affichage des filtres en temporel
for i = 1 : N
    figure(30+i)
    plot(filtres_a_temporel(:,i))
    title(['Reponse impulsionnelle du canal ',int2str(i)])
    grid;
    xlabel('Nombre d''echantillons');
    ylabel('Amplitude');
end


    % affichage du filtrage du signal
	figure
    amp = max(abs(fft(canal_final)));
    %plot(f_freqs,2*amp*abs(f_tronc),f_freqs*2,abs(fft(canal_final)))
    plot(f_freqs,amp*abs(f_tronc),f_freqs,abs(fft(canal_final(1:(size(canal_final)-100),:))));
    title(['Filtrage aux differentes porteuses'])
    axis([0 (2*N/T_b) 0 amp*1.5])
    xlabel('Frequence (Hz)')
    ylabel('Amplitude du spectre')
    grid

    % signal filtre aux differentes porteuses
    figure
    amp = max(abs(fft(canal_final)));
    signal_fft_tronq = abs(fft(recepteur_signal_anal_conv));
    signal_fft_tronq = signal_fft_tronq(1:1000,:);
    freq = (Fa/2)*(0:1:(length(recepteur_signal_anal_conv)-1))/length(recepteur_signal_anal_conv);
    plot(freq(1,1:1000),signal_fft_tronq)
    title(['Signal filtre aux differentes porteuses'])
    axis([0 (2*(N-1)/T_b+b_pass/(2*pi)) 0 amp])
    xlabel('Frequence (Hz)')
    ylabel('Amplitude du spectre')
    grid


%% synchronisation
% Filtre de la sequence pilote genere dans la synchro
figure(40)
plot(recepteur_sync_filtre_n);
xlabel('Echantillons');
title('Filtres genere pour la sequence pilote (Temporel)');
legend('canal 1','canal 2', 'canal 3');

% filtre pilote en frequentiel
figure(31);
plot(fftshift(abs(fft(recepteur_sync_filtre_n))));
xlabel('Echantillons');
title('Filtres genere pour la sequence pilote (frequentiel)');
legend('canal 3','canal 2', 'canal 1');


% affichage de la correlation entre le signal et la sequence pilote et
% prise du poinr maximal
for i = 1:N
    amp = max(recepteur_xcorr_x(:,i))*1.20;
    figure(41+i)
    plot(recepteur_xcorr_y(:,i)*T_n,recepteur_xcorr_x(:,i),'',recepteur_xcorr_y(recepteur_xcorr_max_indice(1,i),i)*T_n,recepteur_xcorr_x(recepteur_xcorr_max_indice(1,i),i),'c*',recepteur_xcorr_y(recepteur_xcorr_max_indice(1,i),i),recepteur_xcorr_x(recepteur_xcorr_max_indice(1,i),i),'co')
    title(['Correlation du signal ',int2str(i),' avec la sequence pilote avec un Delai de ',num2str(recepteur_retards(1,i)*1000*T_n),'ms'])
    axis([recepteur_xcorr_y(1,i)*T_n recepteur_xcorr_y(end,i)*T_n -amp amp])
    xlabel('Temps (s)');
    ylabel('Amplitude du signal');
    grid
end


% affichage de la synchronisation de la sequence pilote sur le signal
amp = max(abs(signal_FA(:,1)))*1.2;
amp2 = max(abs(recepteur_sync_seq_pilote(:,1)));
figure(50)
plot(T_n_FA,signal_FA(:,1),'',recepteur_sync_filtre_n_t+recepteur_retards(1,1)*T_n,amp*recepteur_sync_seq_pilote(:,1)/amp2,'');
axis([recepteur_retards(1,1)*T_n (length(recepteur_sync_seq_pilote)+recepteur_retards(1,1))*T_n -amp amp]);
title(['Synchronisation de la sequence pilote sur le signal R tot 0(t) filtre']);
xlabel('Temps (s)');
ylabel('Amplitude du signal') ;
legend('Signal recu','Sequence pilote');
grid

% affichage des instants d'echantillonnage
for i = 1:N
    amp = max(abs(signal_FA(:,i)))*1.20;
    figure (50+i)
    plot(T_n_FA,signal_FA(:,i),'',T_n_FA(indices_delai(:,i)),signal_FA(indices_delai(:,i),i),'*');
    title(['Estimation des symboles sur le signal R tot ',int2str(i),'(t)']);
    xlabel('Temps (s)'); ylabel('Amplitude du signal') ;
    legend('Signal recu','Symboles estimes');
    grid
end