Discussion :

[Rymuss]

Bonjour à toutes et à tous!

Voilà, comme l'intitulé l'indique, je suis une débutante et j'aimerais savoir comment faire pour programmer la transformée de Hilbert sur MATLAB.
J'ai essayé d'en programmer une d'après ce que j'ai compris de la transformée, mais rien n'y faire.

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
t=0:0.1:10;
f=25;
x=cos(2*pi*f*t);
y=abs(hilbert(x)); % Pour comparer avec ma transformée de Hilbert.
plot(t,x), hold on, plot(t,y,'r')
 
a=zeros(1,101);
 
for t=1:1:10
    for T=0:0.1:10
        if t~=T
            a(t)=a(t)+(cos(2*pi*f*T)/(t-T)); % Calculer l'intégrale qui se trouve dans la transformée de Hilbert.
        end
    end
end
 
A=a.*(1/pi); % Calcul de la transformée de Hilbert.
Abs=abs(A); % Calcul de l'enveloppe.
 
figure, plot(t,Abs)

Pouvez-vous m'aider? Merci d'avance.

[Gooby]

Bonjour,

Je suis également débutant, je ne suis pas sûr de pouvoir t'aider efficacement mais je vois néanmoins que tu établis un vecteur t dans un premier temps, puis que tu t'en sers comme variable de boucle, et enfin tu effectues ton plot avec pour abscisse ton vecteur t (qui n'est plus un vecteur mais une unique valeur!).

Je pense que tu devrais changer t pour un autre nom de variable de boucle dans un premier temps.

ou pour faire plus rapide mais moins propre, recopier

Code :

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

t=0:0.1:10;

juste avant ton plot

[Rymuss]

Merci infiniment pour votre réponse, effectivement j'avais complètement oublié que j'avais deux variables avec le même nom, j'ai donc opté pour la première solution (Renommer le t de la boucle) et maintenant j'ai ENFIN un résultat.

Seulement voilà ... J'ai ceci comme résultat ...

A la place de ce qui est en rouge dans celui là ...

Je ne sais vraiment plus comment faire ...