Resume TP Tns

TP1 : Initiation à MATLAB et Génération des signaux :
Objectif
.Générer les différents types de signaux exponentiels, sinusoïdales et d’autres signaux usuels, les
visualiser en forme continue et discrète.
x=
. the command { k = find(X) } returns the indices of the -1.2000 0 3.1000 6.2000 -3.3000 -2.1000
array X >> inds = find (x < 0)
inds =
156
>> y = x(inds)
3. Tableaux et Matrices y=
-1.2000 -3.3000 -2.1000
. On utilise les crochets { M=[ et ;et ] } pour définir le %autre méth%
début et la fin d’une matrice >>y=[x > 0]
y=
. Pour accéder ou modif un élém de la matr { >>M(l,c)=0 } 0 0 1 1 0 0
. le transposer d’une matrice { M' }
. fonctions zeros(l,c) et ones(l,c) , matrices identité en utilisant la fonction eye(f).
. X=valeur début : pas : valeur fin, peut être remplacée par la fonction
. linespace(val du début, val de fin,nb de val entre les deux)
. square, génère un signal carrée de periode 2π
. Sawtooth, génère un signal triangulaire de periode 2π

. stem(x,y) signal-Discret o
. axis([xmin xmax ymin ymax]).
Input , demande à l’util de saisir une val
. Impulsion :
I(t) = { 0 si t ≠0 ; 1 si t=0
. Echelon : E(t) = { 0 si t <0 ; 1 si t>=0
TP2 : La transformée en Z
Objectifs : Calcul F Trsf H, Calcul et repr ls pôles/zéros de H, Trac reponse en fréq de H
1- commandes utiles
. lookfor sert à rechercher une fonction, ex: t lookfor z transform

. doc ouvrir un browser, et avoir plus de détaille
. disp. Pour afficher un message
. syms sert à créer des variables symboliques (>> syms n z)
2- Transformée en Z
Pour calcul la transformée en Z d’un signal x(n) on utilise ztrans. ( >> X1(z)=ztrans(x1) )
3- Produit de convolution
. conv, La convolution se transforme à une multiplication dans le domaine en Z. ( >> x3=conv(x1,x2)

)
4- Fonction du transfert : H=Y /X
. tf(Y,X,fréqE) fréq reçoit -1 si la fréquence n’est pas connue, Pour déclarer une fonction de
transfert
>> n=[1 0];
>> d=[1 -0.9];
>> tf(n,d,-1)
5- Les pôles et les zéros
On peut tracer les pôles et les zéros sur un graphique zplane(a,b) et pzmap (H).
TP3: Analyse spectrale

Objectifs : calcul transf de Fourier, Rep un signal temporel dans le domaine fréquentiel, Calcul
transf de Fourier inverse.
Commandes utiles
. Diag retourne la diagonale d’une matrice passée en argument à la commande
. conj(z) retourne le conjuguer de z, angle(z) retourne l’argument de z
. real(z) la partie réel , imag(z) la partie imaginaire
2- Domaine fréquentiel : montrera quelle proportion du signal appartient à telle ou telle bande de
fréquence
La représentation d’un signal temporal dans le domaine fréquentiel s’appelle le spectre fréquentiel.
Ce spectre peut être généré par la transformée de Fourier du signal.
3- Transformée de Fourier :
fft(X), calcul de la transformé de Fourier, il reçoit le signal xn, retourne vect de TF
bruit=2*randn(size(t));
plot(Fe*t(1:100),X(1:100)) % pr trac 100 pre
%nn centr
Y = fft(X); % un vect de N pts content coeff Fourier
P1 = (abs(Y))/L; % estimation du spectre
P2=fftshift(P1); % pour centrer le spectre
f = (Fe/L)*(0:(L/2)); % on va prendre que la partie positif du spectre
plot(f,P2(L/2:end))
4- Transformée de Fourier inverse
. ifft, Pr conversion d’un spectre ds l dom fréq à un sig ds dom temporel
il faut utiliser ifftshift avant de procéder à transf inverse afin de rendre au spectre sa nature initiale.
TP4: Filtrage numérique

Objectifs :Calcul les para d’un filtr, Créer filtr num, Filtrer un sig avec filtr num.
Passe bas : Un filtre qui laisse passer les basses fréq et qui atténue les fréq sup à la fréq de coupure
sont rejetées.
Passe haut : Un filtre qui laisse passer les hautes fréq et qui atténue les fréq inf à la fréqde coupure
sont Éliminées
Passe bande : Un filtre qui ne laisse passer qu’une bande ou un intervalle de fréq compris entre une
fréq de coupure basse et une fréq de coupure haute.
Coupe bande : est un filtre empêchant le passage d'une partie des fréq
. Y = filter (b, a, x) filtrer un signal discret avec un filtre num
. fvtool, sert à visualiser et tracer la réponse
. filtfilt m traitt que filter mais en éliminant la phase(=0)
. freqs donne la réponse fréq d’un filtre analog
. freqz donne la réponse fréq d’un filtre num

Filtre à réponse finie(RIF) :
. fir1, fonctions qui servent à créer FIR. [b,a] = fir1(n,fn,'low');

freqz(b,a);
. freqz, a réponse fréq d’un filtre num
. n : ordre de filtr ; fn : fréq de coupure normalisée (fn=2*f/fe)
‘ low’ : pass_bas, ‘high’ , ‘band pass’, ‘stop’
Filtre à réponse infinie
filtre de Butterworth Filtres de Tchebychev

% Type2
[b1,a1] = butter(n,fn,'high'); %Type1 [b3,a3]=cheby2(n,15,fn,'high');
[h1,w1]=freqz(b1,a1); [b2,a2]=cheby1(n,2,fn,'high'); [h3,w3]=freqz(b3,a3);
[h2,w2]=freqz(b2,a2);
plot(w1*fe/(2*pi),abs(h1))

Resume TP Tns

Transféré par

Droits d'auteur :

Formats disponibles

Resume TP Tns

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Resume TP Tns

Transféré par

Droits d'auteur :

Formats disponibles

TP1 : Initiation à MATLAB et Génération des signaux :

. le transposer d’une matrice { M' }

. fonctions zeros(l,c) et ones(l,c) , matrices identité en utilisant la fonction eye(f).

. linespace(val du début, val de fin,nb de val entre les deux)

. square, génère un signal carrée de periode 2π

. Sawtooth, génère un signal triangulaire de periode 2π

. axis([xmin xmax ymin ymax]).

Input , demande à l’util de saisir une val

. Echelon : E(t) = { 0 si t <0 ; 1 si t>=0

. lookfor sert à rechercher une fonction, ex: t lookfor z transform

. disp. Pour afficher un message

. syms sert à créer des variables symboliques (>> syms n z)

. conv, La convolution se transforme à une multiplication dans le domaine en Z. ( >> x3=conv(x1,x2)

4- Fonction du transfert : H=Y /X

>> n=[1 0];

>> d=[1 -0.9];

5- Les pôles et les zéros

TP3: Analyse spectrale

. Diag retourne la diagonale d’une matrice passée en argument à la commande

. conj(z) retourne le conjuguer de z, angle(z) retourne l’argument de z

. real(z) la partie réel , imag(z) la partie imaginaire

plot(Fe*t(1:100),X(1:100)) % pr trac 100 pre

Y = fft(X); % un vect de N pts content coeff Fourier

P1 = (abs(Y))/L; % estimation du spectre

P2=fftshift(P1); % pour centrer le spectre

f = (Fe/L)*(0:(L/2)); % on va prendre que la partie positif du spectre

4- Transformée de Fourier inverse

. ifft, Pr conversion d’un spectre ds l dom fréq à un sig ds dom temporel

TP4: Filtrage numérique

. Y = filter (b, a, x) filtrer un signal discret avec un filtre num

. fvtool, sert à visualiser et tracer la réponse

. filtfilt m traitt que filter mais en éliminant la phase(=0)

. freqs donne la réponse fréq d’un filtre analog

. freqz donne la réponse fréq d’un filtre num

. fir1, fonctions qui servent à créer FIR. [b,a] = fir1(n,fn,'low');

. n : ordre de filtr ; fn : fréq de coupure normalisée (fn=2*f/fe)

‘ low’ : pass_bas, ‘high’ , ‘band pass’, ‘stop’

Filtre à réponse infinie

filtre de Butterworth Filtres de Tchebychev

Vous aimerez peut-être aussi