Discussion :

[kmaniche]

Bonsoir,

Donc je cherche à générer un "Signal sinusoïdal + harmonique" avec Matlab.

Pour un Signal sinusoïdal je prends comme exemple :

Code :

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

1
2
3
 
x = linspace(-pi/2, pi/2) ;
y = sin(x) ;

Alors comment faire pour ajouter des harmoniques à ce signal ?

[Jean Dumoncel]

Salut,

pour rajouter des harmoniques, il te suffit d'ajouter des sinusoides à ton signal. Un exemple sur la décomposition d'un signal carré :
http://fr.wikipedia.org/wiki/Signal_carré

[kmaniche]

je prends le signal défini par :



Comment voir sa représentation graphique

J'ai fais :

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
 
 
% Paramètres du signal
    f = 250 ;
 
% Définition du signal
 
    x=0 ;
    for t=-10:0.1:10
        sm = 0 ;
        for k=1:100
           sm = sm + sin((2*k-1)*2*pi*f*t)/(2*k-1) ;
        end
        x(length(x)+1) = sm ;
    end
 
% Affichage
    plot(x) ;

Mais j'ai pas pu avoir un signal avec harmoniques ?

[tubaas]

Code :

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

1
2
3
4
5
6
7
8
x=0 ;
    for t=-10:0.1:10
        sm = 0 ;
        for k=1:100
           sm = sm + sin((2*k-1)*2*pi*f*t)/(2*k-1) ;
        end
        x(length(x)+1) = sm ;
    end

là tu les a les uns à la suite des autres non?

Code :

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

1
2
3
4
5
6
7
8
 
fs=44100;
f=250;
t=(0:1:44100)/fs;
s=zeros(1,length(t));
for k=1:100
s=s+sin((2*k-1)*2*pi*f*t)/(2*k-1) ;
end

Rq: un signal carré n'a que des harmoniques impaires
+ généralement

Code :

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

1
2
 
s=s+Ak*sin(k*2*pi*f*t+phik);

où Ak et phik amplitude et phase pour l'harmonique k