Discussion :

[lord_of_romns]

bonjour tous le monde

bon voila mon probleme jai un montage electronique sous Simulink et je veux decouper le signal de sortie en harmonique (si il ya un bloc qui permet ca ou si il ya une methode pour le faire) afin de savoir l amplitude de l armonique fondamontale et ca frequence afin de mettre un filtre pour l eliminer .

Merci a l avance

[phryte]

Bonjour.
Je ne pense pas qu'il existe un block pour ton problème.
Pour calculer les harmoniques il faut un nombre d'échantillons de ton signal qui respecte Shannon. Puis cela disponible tu peux utiliser MATLAB Fcn avec le programme suivant :

Code :

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

1
2
3
4
5
6
7
function [Pyy,fr,P,F]=SPyy(s,f,N)
Y = fft(s,N);
Pyy = Y.* conj(Y) / N;
fr = f*(0:N/2)/N;
Pyy=Pyy(1:N/2+1);
[P,nF]=max(Pyy);
F=fr(nF);

Puis connaissant la fondamentale F tu peux utiliser butter pour la supprimer.

[FR119492]

Salut!

jai un montage electronique sous Simulink et je veux decouper le signal de sortie en harmonique

Il me semble y avoir une contradiction dans ton exposé: à la base, Simulink est conçu pour intégrer des équations différentielles avec conditions initiales, donc pour étudier un régime transitoire. En revanche, le découpage en harmonique, donc en série de Fourier s'applique à des régimes permanents périodiques.
Jean-Marc Blanc

[lord_of_romns]

merci phryte pour le coup de main et je veux juste que tu m explique un peu plus comment faire integrer cette equation dans un montage en Simulink (si ce n est pas trop demander et merci )

FR119492-Jean-Marc Blanc il n'ya pas de contradiction dans mon expose car Simulink ne sert pas selement a etudier les regimes transitoire mais aussi les regimes permanent .

[phryte]

Bonjour.
Voici un exemple ( pas très optimisé !) :
Programme à lancer en premier (définit le fichier source et FFT) :

Programme en lancer en second (calcule les éléments de la fondamentale) :

Le programme Matlab :

Code :

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

1
2
3
4
function F=harmo(u)
global H
[i,j]=max(H(:,1,end));
F=[i,j];

Mettre aussi "global H" dans Command Window

[lord_of_romns]

Merci infiniment phryte pour ton aide