Scribd Logo
Importer

Bienvenue sur Scribd !

  • 7 vues

    tp3 Ts 2021

    Transféré par

    You Ssef
    Ce document présente deux méthodes d'intégration numérique, la méthode des rectangles et la méthode des trapèzes. Il donne ensuite deux exemples d'application de ces méthodes pour le calcul numérique de grandeurs caractéristiques de signaux.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd

    tp3 Ts 2021

    Transféré par

    You Ssef
    7 vues4 pages
    Ce document présente deux méthodes d'intégration numérique, la méthode des rectangles et la méthode des trapèzes. Il donne ensuite deux exemples d'application de ces méthodes pour le calcul numérique de grandeurs caractéristiques de signaux.

    Titre original

    tp3_ts_2021

    Copyright

    © © All Rights Reserved

    Formats disponibles

    PDF, TXT ou lisez en ligne sur Scribd

    Avez-vous trouvé ce document utile ?

    Ce contenu est-il inapproprié ?

    Ce document présente deux méthodes d'intégration numérique, la méthode des rectangles et la méthode des trapèzes. Il donne ensuite deux exemples d'application de ces méthodes pour le calcul numérique de grandeurs caractéristiques de signaux.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    7 vues4 pages

    tp3 Ts 2021

    Transféré par

    You Ssef
    Ce document présente deux méthodes d'intégration numérique, la méthode des rectangles et la méthode des trapèzes. Il donne ensuite deux exemples d'application de ces méthodes pour le calcul numérique de grandeurs caractéristiques de signaux.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    sur 4

    UNIVERSITE DE BOUIRA – FACULTE DES SCIENCES TP TRAITEMENT DU SIGNAL

    DEPARTEMENT DE GENIE ELECTRIQUE L3 – ELN – S5 – 2021/2022

    TP3 – CALCUL NUMERIQUE DE GRANDEURS CARACTERISTIQUES


    DES SIGNAUX

    OBJECTIFS DU TP

    Utilisation de l’intégration numérique pour le calcul de certaines grandeurs caractéristiques des signaux
    comme l’énergie, la valeur moyenne, les coefficients de Fourier …

    DEVELOPPEMENT THEORIQUE

    INTEGRATION NUMERIQUE

    L’objectif de l’intégration numérique est de déterminer, en utilisant le calcul numérique, une valeur
    𝑏
    approchée d’une intégrale définie 𝐼 = ∫𝑎 𝑓 (𝑥 ) 𝑑𝑥 , lorsque cette dernière ne peut être évaluée
    analytiquement. Ceci est le cas, par exemple lorsque l’expression de f est un peu compliquée ou que f est
    imparfaitement connue (connue uniquement en un ensemble fini de points).

    Ce dernier cas est souvent rencontré en traitement des signaux physiques dont l’expression peut être
    inconnue, on ne dispose que d’un ensemble de mesures du signal.

    Dans ce cas, le recours à l’intégration numérique devient indispensable pour calculer certaines
    caractéristiques du signal (valeur moyenne, énergie, coefficients de Fourier…) ou effectuer certains
    traitements sur les signaux (corrélation, convolution,…).

    Deux méthodes pour l’intégration numérique sont présentées ci-dessous.

    1- Méthode des rectangles

    Principe

    𝑏−𝑎
    1- Décomposer l’intervalle d’intégration [𝑎, 𝑏] en n sous-intervalles de même longueur ℎ = (h est
    𝑛
    appelé pas d’intégration).
    2- Dans chaque sous-intervalle remplacer la fonction f par une constante égale à la valeur de cette
    fonction à l’une des extrémités du sous-intervalle.
    3- Calculer une valeur approchée de l’intégrale qui n’est rien d’autre que la somme des surfaces des n
    rectangles obtenus de même largeur h et de hauteur variable suivant la valeur de la fonction.

    Formule

    La valeur approchée, notée IR, de l’intégrale est donnée par la formule suivante :
    𝑛−1

    𝐼𝑅 = ℎ ∑ 𝑓 (𝑥𝑖 )
    𝑖=0

    𝑏−𝑎
    Avec ℎ = et 𝑥𝑖 = 𝑎 + 𝑖. ℎ, n représente le nombre de sous-intervalles.
    𝑛

    1
    f(x)

    a=xo x1 x2 xn-1 b=xn x


    h
    2- Méthode des trapèzes

    Principe

    Même principe que la méthode des rectangles sauf que la fonction f est remplacée dans chaque sous-
    intervalle par une droite reliant les 2 points de la fonction se trouvant aux extrémités du sous-intervalle
    (les rectangles hachurés deviennent des trapèzes).

    Formule

    La valeur approchée, notée IT, de l’intégrale est donnée par la formule suivante :
    𝑛−1
    𝑓 (𝑏 ) − 𝑓 (𝑎 )
    𝐼𝑇 = ℎ [ + ∑ 𝑓 (𝑥 𝑖 )]
    2
    𝑖=0

    𝑏−𝑎
    Avec ℎ = et 𝑥𝑖 = 𝑎 + 𝑖. ℎ, n représente le nombre de sous-intervalles.
    𝑛

    TRAVAIL DEMANDE

    −𝑡
    1- On considère le signal 𝑥(𝑡) = {𝑒 𝑠𝑖 0 ≤ 𝑡 ≤ 1
    0 𝑠𝑖𝑛𝑜𝑛
    a)- Calculer analytiquement l’énergie de ce signal.
    b)- Ecrire un programme MATLAB permettant de calculer numériquement cette énergie en utilisant les
    méthodes des rectangles et trapèzes pour différentes valeur du pas d’intégration h = 0.2, 0.1, 0.01 s.
    c)- Comparer les résultats obtenus numériquement avec la valeur exacte de l’énergie. Conclure.

    2- On considère le train d’impulsions rectangulaires


    1 𝑠𝑖 0 ≤ 𝑡 < 0.5
    𝑦 (𝑡 ) = { Périodique de période 𝑇𝑜 = 1 𝑠.
    0 𝑠𝑖 0.5 ≤ 𝑡 < 1
    a)- Ecrire un programme permettant de calculer numériquement les coefficients trigonométriques de
    Fourier du signal jusqu’à 𝑛 = 10. Utiliser la méthode des trapèzes avec un pas h = 0.01 s.
    b)- Comparer les résultats obtenus avec les valeurs exactes en sachant que
    1
    𝑎𝑜 = 0.5 ; 𝑎𝑛 = 0 ; 𝑏𝑛 = [1 − 𝑐𝑜𝑠(𝜋𝑛)]
    𝜋𝑛

    2
    REPONSES

    −𝒕
    1- 𝒙(𝒕) = { 𝒆 𝒔𝒊 𝟎 ≤ 𝒕 ≤ 𝟏
    𝟎 𝒔𝒊𝒏𝒐𝒏

    a)- Energie de ce signal


    +∞ 1
    2( −2𝑡
    1 −2𝑡 1 1 − 𝑒 −2
    𝐸𝑥 = ∫ 𝑥 𝑡)𝑑𝑡 = ∫ 𝑒 𝑑𝑡 = − 𝑒 | = ≈ 0.4323
    2 0 2
    −∞ 0

    b)- Programme

    x=inline('exp(-t)');
    a=0;
    b=1;
    h=0.2;
    n=(b-a)/h;
    s=0;
    for i=0:n-1
    s=s+x(a+i*h)^2;
    end
    Ir=h*s;
    It=h*((x(b)-x(a))/2+s);

    c)- Résultats, comparaison et conclusion

    Méthode des rectangles Méthode des trapèzes


    |𝐸𝑥 − 𝐼𝑅 | |𝐸𝑥 − 𝐼𝑇 |
    𝐼𝑅 𝑒𝑟 = % 𝐼𝑇 𝑒𝑟 = %
    𝐸𝑥 𝐸𝑥
    h = 0.2 0.5245 21.33 % 0.4613 6.71 %
    h = 0.1 0.4770 10.34 % 0.4454 3.03 %
    h = 0.01 0.4367 1.02 % 0.4335 0.28 %

    Les résultats obtenus montrent que :


    - La méthode des trapèzes et meilleur que la méthode des rectangles.
    - En diminuant le pas d’intégration on se rapproche de la valeur exacte.
    - Ce calcul peut être aussi effectué lorsque le signal est défini par un ensemble de valeurs (mesures) à des
    instants donnés.

    𝟏 𝒔𝒊 𝟎 ≤ 𝒕 < 0.5
    2- 𝒚 ( 𝒕) = { de période 𝑻𝒐 = 𝟏 𝒔
    𝟎 𝒔𝒊 𝟎. 𝟓 ≤ 𝒕 < 1
    𝑇0 𝑇0
    1 2 2𝜋𝑛
    𝑎𝑜 = ∫ 𝑥 (𝑡)𝑑𝑡 ; 𝑎𝑛 = ∫ 𝑥 (𝑡)𝑐𝑜𝑠 ( 𝑡) 𝑑𝑡
    𝑇0 𝑇0 𝑇𝑜
    0 0
    𝑇0
    2 2𝜋𝑛
    𝑏𝑛 = ∫ 𝑥(𝑡)𝑠𝑖𝑛 ( 𝑡) 𝑑𝑡
    𝑇0 𝑇𝑜
    0

    3
    a)- Programme

    T0=1;
    h=0.01;
    t=0:h:T0-h;
    x=0.5*square(2*pi/T0*t)+0.5;
    a=0;
    b=T0;
    N=(b-a)/h;
    s=0;
    for i=1:N
    s=s+x(i);
    end
    a0=h*((x(N)-x(1))/2+s)/T0;
    for n=1:10
    sa=0;
    sb=0;
    for i=1:N
    sa=sa+x(i)*cos(2*pi*n*t(i)/T0);
    sb=sb+x(i)*sin(2*pi*n*t(i)/T0);
    end
    a(n)=2*h*((x(N)*cos(2*pi*n*t(N)/T0)-x(1)*cos(2*pi*n*t(1)/T0))/2+sa)/T0;
    b(n)=2*h*((x(N)*sin(2*pi*n*t(N)/T0)-x(1)*sin(2*pi*n*t(1)/T0))/2+sb)/T0;
    end
    disp('Coefficient a0');disp(a0)
    disp('Coefficients an');disp(a)
    disp('Coefficients bn');disp(b)

    Résultat d’exécution

    Coefficient a0
    0.4950
    Coefficients an ( n = 1:10 )
    0.0100 -0.0100 0.0100 -0.0100 0.0100 -0.0100 0.0100 -0.0100 0.0100 -0.0100
    Coefficients bn ( n = 1:10 )
    0.6364 -0.0000 0.2116 -0.0000 0.1263 0.0000 0.0895 -0.0000 0.0688 -0.0000

    b)- Comparaison

    Valeurs exactes
    𝑎𝑜 = 0.5
    𝑎𝑛 = 0
    𝑏𝑛 = 0.6366 0 0.2122 0 0.1273 0 0.0909 0 0.0707 0

    Valeurs approchées obtenues pour h = 0.001


    Coefficient a0
    0.4995
    Coefficients an
    0.0010 -0.0010 0.0010 -0.0010 0.0010 -0.0010 0.0010 -0.0010 0.0010 -0.0010
    Coefficients bn
    0.6366 0.0000 0.2122 -0.0000 0.1273 0.0000 0.0909 -0.0000 0.0707 0.0000
    4

    Vous aimerez peut-être aussi