Discussion :

[Lord Bouillon]

Bonjour à tous,

Je cherche à réaliser un phase tree (arbre de phase) comme en suivant le lien ci dessous (tout en bas de la page dans la partie Phase Tree for Continuous Phase Modulation):
http://fr.mathworks.com/help/comm/ug...odulation.html

il est posté le code suivant:

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
% Parameters from the CPM Modulator Baseband block
M_ary_number = 2;
modulation_index = 2/3;
pulse_length = 2;
samples_per_symbol = 8;
 
L = 5;  % Symbols to display
pmat = [];
for ip_sig = 0:(M_ary_number^L)-1
    s = de2bi(ip_sig,L,M_ary_number,'left-msb');
    % Apply the mapping of the input symbol to the CPM
    % symbol 0 -> -(M-1), 1 -> -(M-2), etc.
    s = 2*s'+1-M_ary_number;
    sim('doc_phasetree', .9); % Run model to generate x.
    % Next column of pmat
    pmat(:,ip_sig+1) = unwrap(angle(x(:)));
end;
pmat = pmat/(pi*modulation_index);
t = (0:L*samples_per_symbol-1)'/samples_per_symbol;
plot(t,pmat); figure(gcf); % Plot phase tree.

Le truc est qu'il se sert de simulink sur ce schéma et que je ne dispose que d'octave (matlab).

J'ai essayé de remplacer le bloc de modulation par mon propre code de modulation mais je n'arrive pas à obtenir l'arbre...

En particulier je ne comprends pas ce que fais :
s = de2bi(ip_sig,L,M_ary_number,'left-msb');

Je sais ce que fais la fonction de2bi mais je vois pas le but...

[Lord Bouillon]

en fait avec un unwrap(angle(signal)) ça fait completement l'affaire =)

Par contre je me demandais comment faire du recouvrement temporel ? Je m'explique mon graphe représente la phase de mon signal en fonction du temps. Et je voudrais replier l'axe temporel afin d'avoir tous les changements de phases exprimés sur 1 seul temps symbole. Any idea ?