Discussion :

[bucher0n]

Bonjour a tous,

Je travaille actuellement dans le cadre de mon école d'ingénieurs sur un gros projet annuel qui se finit dans très peu de temps.
Le but de mon projet est de simuler une chaîne de transmission sous Matlab (sans utiliser Simulink je précise).
Les elements de la chaîne sont les suivants :
- signal sinusoïdal de frequence 200kHz en entrée
- échantillonnage, quantification
- code convolutionnel
- codage en ligne : Manchester
- modulation FSK
- filtre d'émission de Nyquist (en cosinus surélevé)
- canal (bruit)
et partie réceptrice, qui permet de reconstituer le signal de départ.
J'ai pour l'instant tout programmé sauf le décodage Manchester, dans la partie receptrice.

Est-ce que quelqu'un aurait une idee pour implementer ce code, sachant que j'ai en entrée un signal binaire prenant les valeurs -1 et 1 ??

Merci d'avance

[bucher0n]

Je rappelle que le decodage de Manchester doit renvoyer un 1 sur une durée T lorsqu'un -1 (pendant T/2) puis un 1 (pendant T/2) sont en entrée, et 0 lorsqu'un 1 (pendant T/2) puis un -1 (pendant T/2) sont en entrée.
J'ai deja essayé un algorithme simple, mais celui-ci ne porte pas ses fruits :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
for j=0:1:T/2,
        if fskd(j) == 0
            res(j) = 1;
        else res(j)=1;
        end
end
for k=(T/2)+1:1:T,
        if fsk(k)==1
            res(k) = 1;
        else res(k) = 0;
        end
end

[Jerome Briot]

Je n'y connais rien mais ça doit donner quelque chose comme ceci, non ?

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
%Data
Data = [1 0 1 0 0 1 1 1 0 0 1];
Data = [Data;Data];
 
% Clock
Clock = repmat([1;0],numel(Data)/2,1);
 
% Manchester code
MC = (Data(:)&Clock(:)) | (~Data(:)&~Clock(:));
 
% Data from Manchester code
DM = xor(Clock(:),~MC (:));
 
% Period
T = 5;
t = 0:T/4:T/4*numel(MC )-1;
 
figure
 
subplot(411)
stairs(t(:),Clock(:),'r-')
title('Clock')
set(gca,'ylim',[-.5 1.5],'xlim',[0 max(t)], ...
    'xtick',t(1:2:end),'xticklabel','',...
    'xgrid','on')
 
subplot(412)
stairs(t(:),Data(:),'r-')
title('Data')
set(gca,'ylim',[-.5 1.5],'xlim',[0 max(t)], ...
    'xtick',t(1:2:end),'xticklabel','',...
    'xgrid','on')
 
subplot(413)
stairs(t(:),MC (:),'r-')
title('Manchester code')
set(gca,'ylim',[-.5 1.5],'xlim',[0 max(t)], ...
    'xtick',t(1:2:end),'xticklabel','',...
    'xgrid','on')
 
subplot(414)
stairs(t(:),DM(:),'r-')
title('Data from Manchester code')
set(gca,'ylim',[-.5 1.5],'xlim',[0 max(t)], ...
    'xtick',t(1:2:end),'xticklabel','',...
    'xgrid','on')

(Source : wikipedia )

[bucher0n]

Merci pour cette réponse si rapide.
Oui en effet c'est exactement ça le principe.
J'ai essayé d'integrer ce programme dans le mien pour le decodage comme ceci :

Code :

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

1
2
3
4
5
6
7
8
T=10;
t2 = 0:T/4:T/4*numel(fskd )-1;
h = repmat([1;0],numel(fskd)/2,1);
ndims(h)
ndims(fskd)
length(h)
length(fskd)
res = xor(h(:),~fskd(:));

Sachant que :
*la durée d'un bit est de 5
*mon signal a decoder est fskd et equivaut au MC de ton programme
*res correspond au DM de ton programme

Probleme : lorsque je compile, matlab me retourne une erreur : "matrix dimensions must agree" pour la ligne du xor.
C'est pour ça que j'ai mis ndims(h) et ndims(fskd) pour voir les dimensions, mais il me dit que ces 2 matrices ont une dimension de 2. Les longueurs (length) sont 124 pour h et 125 pour fskd. Comment faire ?

[Jerome Briot]

Les longueurs (length) sont 124 pour h et 125 pour fskd. Comment faire ?

Enleve une valeur à fskd, non ?

[bucher0n]

oui , mais comment faire ?

[Jerome Briot]

=> supprimer un élément d'une matrice

[bucher0n]

merci beaucoup, je vais essayer ça :

Code :

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

fskd(125) = []

[bucher0n]

J'ai donc essayé

Code :

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

fskd(125) = []

, ca m'enleve l'erreur mais ca me donne pas le bon signal. J'ai bien un signal carré mais les valeurs ne correspondent pas.
En pieces jointes l'allure de mon signal a decoder (demod) et le signal que j'obtiens (decod)

Voici le code qui m'a généré ceci :

Code :

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

1
2
3
4
5
6
7
T=10;
t2 = 0:T/4:T/4*numel(fskd )-1;
h = repmat([1;0],numel(fskd)/2,1);
fskd(125)=[];
res = xor(h(:),~fskd(:));
 
stairs(t2(:),res (:),'r-')

Merci

[Jerome Briot]

Tu pourrais attacher un exemple de signal sous forme d'un fichier mat compressé en zip ?

[bucher0n]

J'ai attaché mon programme. Le fichier principal est main.m, a partir duquel on lance la simulation.
Il faudrait juste que tu jettes un coup d'oeil a la partie decodage manchester, qui recupere le resultat de la partie demodulation fsk. (fskd)