Discussion :

[Takou]

Bonjour,
J'ai un signal composée d'une sinusoïde et d'une rampe comme dans la pièce jointe.
Mon but est de filtrer le signal rampe pour centrer la sinusoïde autour de l'axe x sauf que le signal rampe a une transformée de Fourier assez compliquée et je vois pas quel type de filtre utiliser et comment l'implémenter.
Voilà le bout de code que j'ai commencé à écrire:

Code :

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

1
2
3
4
5
6
7
8
L=length(x)                  % length of signal x
NFFT = 2^nextpow2(L); % Next power of 2 from length of x
X = fft(x,NFFT)/L
figure();
plot(abs(X)); 
X=fftshift(X);
figure();
plot(abs(X));

Ceci donne la FFT dans la deuxième pièce jointe.
Quelqu'un a une idée?
Merci d'avance

[compaq610]

Bonsoir,
Je pense que le signal rampe se comporte comme une sinusoïde a une fréquence "tend vers zéro" donc vous pouvez utilisé un filtre passe haut a une fréquence de coupure inférieure a celle de la sinusoïde.
Voilà le code que j'ai modifié

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
clear all
t=0:0.001:1;
y=10*sin(t.*20*pi);
x=20*t+y;
figure(1)
plot(x)
L=length(x);                  % length of signal x
NFFT = 2^nextpow2(L); % Next power of 2 from length of x
X = fft(x,NFFT)/L;
figure(2)
plot(abs(X)); 
X=fftshift(X);
figure(3)
plot(abs(X));
[b,a]=butter(1,0.01,'high');   %filtre passe haut
sig_rec=filter(b,a,x);             %signal reconstruit 
figure(4)
plot(sig_rec)

[Takou]

Bonjour,

Code :

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

1
2
[b,a]=butter(1,0.01,'high');   %filtre passe haut
sig_rec=filter(b,a,x);             %signal reconstruit

Merci pour ta réponse mais je ne dispose pas de la Signal Processing Toolbox et donc je ne peux pas utiliser la fonction "butter" pour calculer les coefficients du filtre.

[FR119492]

Salut!

Une première remarque: ne penses-tu pas qu'à la longue, ton signal va tendre vers l'infini?

Tout cela mis à part, pourquoi vouloir utiliser une transformée de Fourier? Il y a beaucoup plus simple:

1. Tu cherches tous les maxima et les minima de ta fonction.
2. Tu écris les équations des deux droites qui passent par ces points. Ces équations sont de la forme x=A*t+B1 et x=A*t+B2.
3. L'équation de ta rampe est x=A*t+(B1+B2)/2.

Jean-Marc Blanc

[Takou]

Bonjour,
Merci de votre réponse. J'ai pensé aussi à cette solution mais le problème est que mon signal est bruité, il est donc très difficile de trouver les minimas et les maximas de la fonction.

[Takou]

Re,
Finalement j'ai adopté la méthode la plus facile: calculer l'équation de la droite en utilisant le premier et le dernier point du signal.
Merci quand même pour vos réponses.