Discussion :

[d0n32]

Bonjour,
Je souhaite réaliser la modulation/démodulation d'un signal en BPSK.
J'en suis à la modulation pour le moment, et j'ai écrit :

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
clear all;
close all;
 
%Data brutes UL
X = floor((rand(1,1000)+0.5));  
 
%Codage & Modulation BPSK
alpha = zeros(1,length(X));  
for k=1:length(X)
    if X(k) == 0
        alpha(k) = 1;
    elseif  X(k) == 1
        alpha(k) = -1;
    end
end
phi = (pi/2).*alpha;
%Signal temporel modulé
fe=125*10^6;
t = linspace(0,12,length(phi));
 
s = exp(i.*(2*pi*fe.*t + phi));

Mais quand je plotte s, le résultat n'est pas fou... Est-ce que vous voyez ce qui pourrait ne pas aller ?
Merci

[Invité]

Bonjour d0n32,

Ton problème vient de la génération de ton signal à moduler, et plus particulièrement:

Code :

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

t = linspace(0,12,length(phi));

Je présume que fe est ta fréquence d'échantillonnage, donc fait plutôt comme ceci:

Code :

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

t = linspace(0,12,length(phi))/fe;

[d0n32]

Merci !
Le résultat obtenu n'est-il pas étrange ? La théorie prévoyait des sinusoides qui changeaient de signe pour un valeur de bit différente..

[Invité]

Bonjour d0n32,

La théorie prévoyait des sinusoides qui changeaient de signe pour un valeur de bit différente

Oui c'est le cas, mais tu as peut-être du mal à le voir car ton signal modulateur change practiquement chaque valeur.
Mais par exemple tu peux parfaitement le visualiser en fixant alpha comme ceci:

Code :

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

alpha = [ones(1,floor(length(X)/2)) -ones(1,ceil(length(X)/2))];

[d0n32]

Merci beaucoup

[d0n32]

Re moi,Je n'arrive pas à comprendre pourquoi le code suivant me retourne comme résultat Y2= (0 0 0 0 ... 0) au lieu de ce que je souhaite, à savoir les coordonnées de Y répétées chacune un certain nombre de fois :

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
 
fe=125*10^6;    
fi=70*10^6;
br=10*10^6;   
 t_simu = 1*10^-3;
 
X = floor((rand(1,100)+0.5));
Y = pskmod(X,2);
 
 
phase = zeros(1,t_simu*fe);
 
 
increment = floor((2^48)*br/fe+0.5);
 
 
 
Y2 = zeros(1,t_simu*fe);
k=1;
i=1;
 
 
while k<=length(Y)
    for j=1:length(phase)-1
        phase(j+1)=mod(phase(j)+increment,2^48);
 
        if phase(j+1)<=phase(j)
            if k<=length(Y)
                enable=1;
                Y2(i)=Y(k);
                k=k+1;
            end
        else
            if k<=length(Y)
                enable=0;
                Y2(i)=Y(k);
            end
        end
        if i<length(Y2)
            i=i+1;
        end
    end
end
Y2

Je vous remercie, une fois encore