Discussion :

[jimmypoker]

Bonjour,

J'ai un signal (.mat en pièce jointe) et j'aimerais savoir comment, à partir de ce signal, je peux
- Déterminer les composantes dominante
- Enlever les composantes les plus faibles
- Faire différentes approximations

J'ai déjà réalisé un script qui le traite en affichant sa fft (module et phase) et tracé la série de fourier correspondante.

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
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
function four = fourier
 
close all;
clc;
 
load signal.mat s
 
N = length(s);
 
% fréquence d'échantillonage 
fs = N/1.5;
t = 1/fs;
 
% Axe temporel
n = 0:t:(N-1)*t;
S = (fft(s))/N;
 
% Axe fréquenciel
F = ((1/(fs))*N)*(0:N-1);
 
% Filtrage du bruit
eps=10e-9;
S(find(abs(S)<eps))=0;
 
% Décomposition de la FFT  
mod=abs(S);
phaz=angle(S);
 
%Représentation temporelle du signal de départ
subplot(222)
stem(n,s,'.-');
grid on;
title('Signal.mat - Temporel');
 
% Représentation du module et de la phase de la fft
subplot(221)
stem(F,mod,'markerfacecolor', 'blue');
title('Signal.mat - Module de la FFT');
subplot(223)
stem(F,phaz,'markerfacecolor', 'blue');
title('Signal.mat - Phase de la FFT');
 
%Les fréquences et amplifications se trouvent sur le graphe du module de la FFT et
%les déphasages des composantes de Fourier sur son graphe de phase.
A=[2,1,0.8,0.4];
omega=[0,4,8,36];
phi=[0,-2.3562,0,-1.5708];
u=size(n);
Y = zeros(1,u);
 
for i = 1:length(A)   
    Y = Y + A(i)*cos(2*pi*omega(i)*n+phi(i)); 
end
 
% Y = A(1:u)*cos(2*pi*omega(1:u)*n+phi(1:u));
 
%Représentation temporelle du Signal de Fourier
subplot(224)
stem(n,Y,'.-');
grid on;
title('Y=2+cos(2*pi*4*n-2.3562)+0.8*cos(2*pi*8*n)+0.4*cos(2*pi*36*n-1.5708)');
end

Merci de votre aide!