Discussion :

[minette-87]

bonjour tout le monde
voila j'ai un petit exercice a traiter , mais je n'arrive pas :
l'ennoncé de l'éxercice:
soit un signal échantillonné à 1000 HZ composé de deux sinusoîdes d'amplitude A1=1 et A2=2 et de fréquences f1=150 et f2=140hz
1. ecrire un programme Matlab pour calculer le périodogramme du signal.

bon j'ai essayé mais c'est tout ce que j'ai pu faire , voila mon ptit programme:

Code :

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

1
2
3
4
5
6
7
8
9
A1=1;
A2=2;
f1=150;
f2=140;
fe=1000;
t=1:0.001:(1/fe);
S=A1*sin(2*pi*f1*t)+A2*sin(2*pi*f2*t);
figure(1)
periodogram(S)

[Jerome Briot]

Vérifie bien les valeurs contenues dans le vecteur t...

Et tu as oublié un i dans le nom de la fonction PERIODOGRAM

[minette-87]

Vérifie bien les valeurs contenues dans le vecteur t...

Et tu as oublié un i dans le nom de la fonction PERIODOGRAM

merci pour ta réponse , mais mon grand probleme c'est dans la déclaration du vecteur temps ,

[Jerome Briot]

merci pour ta réponse , mais mon grand probleme c'est dans la déclaration du vecteur temps ,

C'est bine pour ça que je te parlais de t dans mon précédent message...

Un vecteur se définit comme ceci :

Code :

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

V = debut:pas:fin ;

Donc dans ton code, on a :

• debut = 1
• pas = 0.001
• fin = 1/fe soit 0.001

Tu vois ce qui cloche maintenant ?

[minette-87]

C'est bine pour ça que je te parlais de t dans mon précédent message...

Un vecteur se définit comme ceci :

Code :

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

V = debut:pas:fin ;

Donc dans ton code, on a :

• debut = 1
• pas = 0.001
• fin = 1/fe soit 0.001

Tu vois ce qui cloche maintenant ?

ahh oui c'est vrai ,défois je suis vraiment bête , j'ai mis 1 au lieu de mettre un 0.
mercii beaucoup
est ce que le programme est juste STP

[minette-87]

euuhh le probléme c'est que je n'arrive pas a visualiser le graphe de l'addition des deux signaux sinusoîdales , qu'on j'appui sur f5 pour executer le programme , j'obtien une figure vierge , j'arrive pas a trouver le probleme:s

[Jerome Briot]

Toujours le même problème :

Vérifie bien les valeurs contenues dans le vecteur t...

[minette-87]

oui j'ai corriger ,

Code :

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

t=0:0.00001:(1/fe)

[Jerome Briot]

Je pense que tu devrais chercher un peu plus ce que représente la fréquence d'échantillonnage d'un signal et comment on relie la fréquence et la période... mais cela n'a aucun rapport avec MATLAB

Je te rappelle que ton signal est échantillonné à une fréquence de 1000Hz