Discussion :

[Moustapha2013]

Bonjour, je suis nouveau sur matlab et voudrais de l'aide pour mon projet de session. Je veux savoir comment coder sur matlab la probabilité d'erreur en fonction du rapport signal à bruit pour les modulations 64QAM et 16QAM

[Gooby]

Bonjour,

Si tu veux de l'aide, il faudra commencer par nous montrer ce que tu as fait par toi même.

[Moustapha2013]

Bonjour voici mon code
Je veux tracer le TEB en fonction du rapport signal sur bruit. Le problème que j'ai, c'est au niveau des taux d'erreur binaire. Ça devait être une matrice de 35 éléments mais ça sort juste une valeur. La partie 1 fonctionne bien mais c'est la partie2 qui pose problème
Partie1

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
clc
clear all
% Devoir1 MGR880%
% Déclaration des variables %
f=2600000;
W=5000;
s=-104; %Sensibilité du récepteur
fb = 'fbit';
Eb= 'Energie bit';
Ptx=30;
Lctx=2.5;
Lcrx=2.5;
Gtx=15;
Grx=15;
htx=0.020;
hrx=0.015;
t=20;
T=t+273;
Prx=[];
Loss=[];
K=1.38*10^-23;
N=K*T*W;
n= 'fb/W';
nPSK=1/2; 
nQPSK=1;
n16QAM=2;
n32QAM=2.5;
n64QAM=3;
% Expression de la puissance de réception en fonction de d %
d=1:1:35;
Prx(d)= Ptx+Gtx+Grx+Lctx+Lcrx+20*log10(htx)+20*log10(hrx)-40*log10(d)+30;
Loss(d)= Prx(d)-Ptx-Gtx-Grx;
figure;
subplot(411);
plot(d,Prx);
xlabel('Distance (Km)','fontSize',10);
ylabel('Puissance de réception (dBm)','fontSize',10);
title('Puissance de réception en fonction de la distance','fontSize',10);
subplot(412), plot(d,Loss);
xlabel('Distance (Km)','fontSize',10);
ylabel('Perte totale en fonction de la distance (dB)','fontSize',10);
title('Perte totale en fonction de la distance','fontSize',10);
v=Prx(d)-.10*log10(N)
subplot(413);
plot(d,v)
xlabel('Distance (Km)','fontSize',10);
ylabel(' Rapport signal sur le bruit (S/N)en(dB)','fontSize',10);
title('Rapport signal sur le bruit en fonction de la distance','fontSize',10);
F=Prx(d)-s;
subplot(414);
plot(d,F)
xlabel('Distance (Km)','fontSize',10);
ylabel('Marge évanouissement en(dB)','fontSize',10);
title('Marge évanouissement en fonction de la distance','fontSize',10);

Partie2

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
% Énergie de bit sur la densité spectrale de puissance en dB%
%Déclaration des variables respectif du rapport (Eb/N0) en fonction des
%différents types de modulation
teb1=[];
q1=v-.10*log10(nPSK);% (Eb/N0)PSK
q2=v-.10*log10(nQPSK); % (Eb/N0)QPSK
q3=v-.10*log10(n16QAM); % (Eb/N0)16QAM
q4=v-.10*log10(n32QAM); % (Eb/N0)32QAM
q5=v-.10*log10(n64QAM); % (Eb/N0)64QAM
 
 
 
 
%Énergie de bit sur la densité spectrale de puissance en watt
m1=10.^(q1/10);
m2=10.^(q2/10);
m3=10.^(q3/10);
m4=10.^(q4/10);
m5=10.^(q5/10);
 
%Taux d'erreur binaire pour les différents types de modulation 
teb1=1/2*exp(-m1'.^2)/(1.7725*m1);
teb2=1/2*exp(-m2.^2)/(1.7725*m2);
teb3=3/2*exp(-m3/2.^2)/(1.7725*m3/2);
teb4=1.6464*exp(-m4/2.^2)/(1.7725*m4/2);
teb5=7/4*exp(-m4/2.^2)/(1.7725*m4/2);

[Invité]

Bonsoir,

Sur ces dernières lignes, attention de ne pas oublier d'utiliser l'opérateur ./ (membre à membre) et non / (matriciel) :

Code :

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

1
2
3
4
5
teb1=1/2*exp(-m1.^2)./(1.7725*m1);
teb2=1/2*exp(-m2.^2)./(1.7725*m2);
teb3=3/2*exp(-m3/2.^2)./(1.7725*m3/2);
teb4=1.6464*exp(-m4/2.^2)./(1.7725*m4/2);
teb5=7/4*exp(-m4/2.^2)./(1.7725*m4/2);