Echantionnage Filtres-Analogiques

Université des Sciences et de la Technologie Houari Boumediene (USTHB) REALISATION
Faculté d'Electronique et d'Informatique (FEI) Manip1 :

Département de Télécommunication Taper le script suivant :
Section: Licence Telecom, 3ème année, Semestre (5)- Section A %clear all;close all,clc
TP2 en Traitement du signal
Ts=100; % Période du sinus
Duree=200; % Durée d'observation signal 2 périodes
dt=50; % Pas ou période d'échantillonnage temporel du signal
Objectif:
N=Duree/dt; % Nombre total d'échantillons
n=0:N; % "Vectorisation du temps" (échantillonnage)
1. Echantillonnage des signaux : Le but de la 1ère partie est d'étudier l’opération
d'échantillonnage indispensable avant tout traitement numérique d'un signal. t=n*dt; % Définition du temps
2. Analyse et synthèse des filtres analogiques : De synthétiser un filtre y=sin(2*pi*t/Ts); % Génération du sinus
figure();plot(t,y,'*');
analogique sous Matlab
title('le siganl échantillonné');
xlabel('t ');ylabel('y');
Partie I. Echantillonnage des signaux grid on;
%hold on;
Le signal analogique x(t), continu dans le temps, est représenté par un ensemble de
valeur discrètes xe(t) = x(nTe) avec n : entier et Te : période d’échantillonnage. 2. Exécuter le même programme pour : dt=1, dt=10, dt=50, dt=100.
3. Que remarquez-vous ?
Manip2 :
On considère le signal sinusoïdal x = cos(2*pi*50*t) généré par le code suivant :
clear all; close all

t= 0:1/1000:.04; % Durée d'observation signal
x = cos(2*pi*50*t); % approx. to continuous-time
figure(1);plot(t,x);
xlabel('Time (sec)'); ylabel('Amplitude');
Préparation théorique
1. Un signal est constitué de 3 sinusoïdes de fréquences et d'amplitudes différentes: 1. Quelle est la fréquence du signal f?
( ) = cos(2 ) + 2 cos(2 ) + cos(2 ) 2. Quelle est la fréquence d’échantillonnage de ce signal?
avec = 0.5 , = 1.5 , et = 1.6 3. Le théorème de Shannon d’échantillonnage est-il respecté ? Quelle conclusion en
 Représenter le spectre d'amplitude du signal échantillonné pour les 2 tirez-vous,
fréquences d'échantillonnage suivantes: a) = 2 et b) = 4 4. Déterminer le nombre total d'échantillons (N)
 b. Toujours pour les deux cas de fréquences d'échantillonnage, représenter le 5. ajouter les lignes suivantes et commenter chaque ligne d'instruction
spectre d'amplitude du signal continu obtenu pour reconstruction idéal à
partir du signal échantillonné. n=0:(N-1); %Vectorisation du temps" (échantillonnage)
2. Donner et expliquer le théorème d'échantillonnage. Vous illustrerez vos propos en tn=n*Te; % Définition du temps discret
traitant l'échantillonnage pour Fe=150 Hz et Fe=500 Hz d'un signal x(t) dont le xn= cos(2*pi*50*tn); % signal echantillonnée (discret-time)
spectre est donné par la figure ci-dessous: figure(2);stem(tn,xn); % visualisation du signal x(tn)
xlabel('Time (sec)'); ylabel('Amplitude');
figure(3);stem(n,xn); % visualisation du signal x[n]
xlabel('Time Index (n)')
ylabel('x[n]')
6. Exécuter le programme,
Boularbah. SOUISSI Traitement du Signal : L3-SECTION B email:souissib@hotmail.com Boularbah. SOUISSI Traitement du Signal : L3-SECTION B email:souissib@hotmail.com
7. Comparer entre les trois figures (1, 2, 3). 3. Quelle sont la différence entre prog 1 et prog 2 (Programme et résultats)?
Commenter vos résultats 4. Dans les 02 figures, mesurer le module pour la fréquence de coupure, que concluez-
vous?
5. Le Signal d’entrée Ve(t) est signal sinusoïdal d’amplitude A = 5V et de la pulsation
Partie II. Analyse et synthèse des filtres analogiques wm =2πfm= 2 rad/s,
 Dans les 2 figures, déterminer H(wm), et le déphasage (wm)-
Instructions: freqs, loglog, semilogx, buttord, butter, buttord,cheby1, cheb1ord  Calculer manuellement l'amplitude de sortie,
 Tracer sous Matlab en fonction du temps sur le même graphique le signal
Exo1 (Voir Exo1_TD2): d'entrée Ve(t) et le signal de sortie Vs(t)
Visualisons sous Matlab la réponse en amplitude, la réponse en phase d’un filtre dont on 6. Même question (4) pour wm=0.05 rad/s,
connaît la fonction de transfert opérationnelle : ( )= , et ( )= avec s=jw 7. Que concluez-vous?
8. Faire les changements nécessaires au programme ci-dessous pour afficher le
1. Préparation théorique: module et la phase de la fonction de transfert suivante : ( ) = ,
Quelle conclusion en tirez-vous?
1. Quel est le type et l'ordre de ce filtre?
2. Déterminer la pulsation de coupure wc=2πfc (où fc est la fréquence de coupure) Si on applique le même signale, Déterminer l'amplitude, le temps de déphasage
3. Donner le module (H et Hdb=20log(H)) et la phase entre le signal d'entrée et de sortie.
4. Tracer les deux diagrammes de Bode (Amplitude et phase)
REALISATION
Manip3 :
2. Ecrire le programme suivant sous Matlab
%freqs([1],[1 1]); %montre la réponse de 1/(s+1) entre 10e-2 et 10e1 rad/s
. .
%** ( ) = .
%**a = [1 0.4 1];
%**b = [0.2 0.3 1];
%**w = logspace(-1,1);
%**freqs(b,a,w)
clear all; close all

%**********************prog 1**************************
b=[1];a=[1 1];
freqs(b,a);
%**********************prog 2**************************
[H,w]=freqs(b,a);
%[H,w]=freqs([1],[1 1]);
figure();
subplot(2,1,1)
semilogx(w,20*log10(abs(H)));
xlabel('Frequency (radians)'); ylabel('Attenuation (dB)'); grid;
subplot(2,1,2)
%semilogx (w, angle(H)*180/pi); %calcul de phase en degree
semilogx (w, angle(H)); %calcul de phase en radians
xlabel('Frequency (radians)'); ylabel('Phase in degrees'); grid;
1. Commenter chaque instruction

2. Exécuter le programme et interpréter les résultats des deux figures
Boularbah. SOUISSI Traitement du Signal : L3-SECTION B email:souissib@hotmail.com Boularbah. SOUISSI Traitement du Signal : L3-SECTION B email:souissib@hotmail.com

Echantionnage Filtres-Analogiques

Transféré par

Droits d'auteur :

Formats disponibles

Echantionnage Filtres-Analogiques

Transféré par

Informations du document

Description originale:

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

Echantionnage Filtres-Analogiques

Transféré par

Droits d'auteur :

Formats disponibles

Université des Sciences et de la Technologie Houari Boumediene (USTHB) REALISATION

Faculté d'Electronique et d'Informatique (FEI) Manip1 :

clear all; close all

%freqs([1],[1 1]); %montre la réponse de 1/(s+1) entre 10e-2 et 10e1 rad/s

clear all; close all

1. Commenter chaque instruction

Vous aimerez peut-être aussi