Travaux Pratiques Matlab/Simulink page 1

Exercice 1
Tracer la fonction signe en exécutant le programme suivant :

C(t) = sign(t)

fs=250;
dt=1/fs;
t=-5:1/250:5;

T
T=1
c=(sign(t))
plot(t,c,’r’,’linewidth’,3);
grid on; set(gca,’fontsize’,14);

F
title(’le tracé de c(t)’)
xlabel(’t(ms) ’);
ylim([-1.2 1.2]);
ylabel(’ signe(t)’);

A
Exercice 2
a) Tracer le signal échelon en utilisant la fonction heaviside.
b) Tracer le signal échelon à l’aide de la fonction signe.

R

1 t>0
C(t) =
0 t < 0

Exercice 3

D
a) Tracer le signal rectangulaire à l’aide de la fonction rectpuls.
b) Tracer le signal rectangulaire à l’aide des fonctions heaviside.
c) Tracer le signal rectangulaire à l’aide des fonctions signe.

t
 
C(t) = rect
2

Exercice 4
Calculer le produit de convolution y(t) = x(t) ⊗ h(t)

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 2


0.6

 −1 ≤ t ≤ 0.5
x(t) = 0.3 0.5 ≤ t ≤ 3


0 t < −1 and t > 3

e−t t≥0
h(t) = e−t u(t) =
0 t<0

T
th1=linspace(0,10,1001);
h1=exp(-th1);
h=[0 h1];
th=[0 th1];

F
plot(th,h)
tx=[-1 -1 0.5 0.5 3 3];
x=[0 0.6 0.6 0.3 0.3 0];
figure(1);
plot(tx,x,’:’,th,h)
legend(’x(\tau)’,’h(\tau)’)

A
t=-3;
figure(2);
plot(tx,x,’:’,-th+t,h)
legend(’x(\tau)’,’h(t-\tau)’)

R
Exercice 5
Tracer le signal triangulaire en utilisant la commande tripuls

D
Avec : A = 1, τ = 2 et T = 3

Exercice 6
c) Tracer les signaux Dirac δ(t), δ(t + 3) en utilisant la fonction gauspuls.
d) Tracer les deux signaux dans le même graphe.
e) Tracer les mêmes signaux en utilisant la commande Dirac.

Dr. Ghislain PANDRY

Travaux Pratiques Matlab/Simulink page 3

Exercice 7
d) Tracer le signal sinc(t) en utilisant les commandes linspace et sinc.
e) Tracer le signal sinc(t) en utilisant la commande sin.

sin(πt)
sinc(t) =
πt

Exercice 8

T
Soit le signal périodique :

1 0≤t≤1
x(t) = 
−1 1 ≤ t ≤ 2

F
t1=0:0.01:1;
t2=1:0.01:2;
x1=ones(size(t1));

A
x2=-ones(size(t2));
x=[x1 x2];
xp=repmat(x,1,5);
t=linspace(0,10,length(xp));

R
plot(t,xp)

Le traçage de x(t) sur une période :

D
%%%%%%%%%%%%%%% x(t) sur une période
%%%%%%%%%%%%%%%%%%%%%%
syms t
x=heaviside(t)-2*heaviside(t-1);
ezplot(x,[0 2]);

La décomposition en série de Fourier en forme complexe est :

+∞
Cn ejnω0 t
X
x(t) =
n=−∞

T
1 Z 20
Cn = x(t) e−jnω0 t dt
T0 − T20

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 4

n=-2:2;
t0=0;
T=2;
w=2*pi/T;

c=(1/T)*int(x*exp(-j*n*w*t),t,t0,t0+T);

xx=sum(c.*exp(j*n*w*t))

T
ezplot(xx,[0 10])
title(’Approximation with 5 terms’)
%%%%%%%%%%%%%%%%%%%%%%%

n=-5:5;

F
c=(1/T)*int(x*exp(-j*n*w*t),t,t0,t0+T);
xx=sum(c.*exp(j*n*w*t));
ezplot(xx,[0 10])
title(’Approximation with 11 terms’)
%%%%%%%%%%%%%%%%%%%%%%%

A
n=-10:10;
c=(1/T)*int(x*exp(-j*n*w*t),t,t0,t0+T);
xx=sum(c.*exp(j*n*w*t));
ezplot(xx,[0 10])

R
title(’Approximation with 21 terms’)
%%%%%%%%%%%%%%%%%%%%%%%

n=-30:30;
c=(1/T)*int(x*exp(-j*n*w*t),t,t0,t0+T);

D
xx=sum(c.*exp(j*n*w*t));
ezplot(xx,[0 10])
title(’Approximation with 61 terms’)
%%%%%%%%%%%%%%%%%%%%%%%

Décomposer la fonction x(t) en Séries de Fourier trigonométrique (d’Euler)

∞ 
X 
x(t) = ⟨x⟩ + an cos(nω0 t) + bn sin(nω0 t)
n=1

T
1 Z 20
⟨x⟩ = x(t) dt = x0
T0 − T20

Dr. Ghislain PANDRY

Travaux Pratiques Matlab/Simulink page 5

T T
2 Z 20 2 Z 20
an = x(t) cos(nω0 t) dt bn = x(t) sin(nω0 t) dt
T0 − T20 T0 − T20

NB : APPLIQUER AUX SERIES DE FOURIER DES EXERCICES

DU TD

Exercice 9

T
Soit le signal suivant : 
1 −1 ≤ t ≤ 1
x(t) = 
0 ailleurs

F
1. Tracer la fonction x(t)
2. Calculer et tracer la transformée de Fourier de x(t) en utilisant la fonction
prédéfinie du Matlab fourier.

A
3. Calculer la transformée de Fourier par le calcul mathématique

Exercice 10

R

t + 1 −1 ≤ t ≤ 0
x(t) = 
−t + 1 0 ≤ t ≤ 1

D
1. Tracer la fonction x(t)
2. Calculer et tracer la transformée de Fourier de x(t) en utilisant la fonction
prédéfinie du Matlab fourier.
3. Calculer la transformée de Fourier par le calcul mathématique
NB : APPLIQUER A LA TRANSFORMEE DE FOURIER DES
EXERCICES DU TD

Exercice 11
Pour trouver la transformée de Laplace on utilise l’instruction : laplace().
Pour savoir plus pour la commande Laplace, on écrit :

>>help laplace

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 6

Exécuter les instructions suivantes pour calculer les transformées de Laplace des
fonctions suivantes
%Déclaration des variables symboliques
>>syms a s t w x ;
%calculer la transformée de Laplace
>>laplace(10,t,s)
>>laplace(t^5)
>>laplace(exp(a*s))

T
>>laplace(sin(w*x),t)
>>laplace(sin(w*x),s)
>>laplace(cos(x*w),w,t)
>>laplace(t^(3/2),s)

F
>>laplace(dirac(t),t,s)

Transformée inverse de Laplace

Exercice 11

A
%Déclaration des variables symboliques
>>syms s t w x y;
%calculer la transformée inverse de Laplace
>>ilaplace(1,t)

R
>>ilaplace(1/s)
>>ilaplace(1/(s-1))
>>ilaplace(1/(s-1))
>>ilaplace(1/(t^2+1))
>>ilaplace( y/(y^2+w^2),y,x)

D
>>ilaplace(2*s^(-3))

Développement en fraction simple :

Considérons la fonction de transfert
B(s) num b0 sn + b1 sn−1 + · · · + bn
= = n
A(s) den s + a1 sn−1 + · · · + an
Où certains des constantes, a1 et b1 , peuvent être nuls. Dans Matlab les vecteurs
lignes num et den spécifient les coefficients du numérateur et du dénominateur de
la fonction du transfert. C’est-à-dire :
num = [b0 b1 . . . bn ]
den = [a1 . . . an ]

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 7

Ecriture d’une fonction de transfert avec la commande tf(num,den) : lire et com-

prendre la signification de cette commande en utilisant le help
La commande :
[r,p,k]=residue(num,den)
Calcule les résidus (r), les pôles (p) et le terme direct (k) du développement en
fraction partielle du rapport de deux polynômes B(s) et A(s).
Le développement en fraction partielle de B(s)/A(s) est donné par :

T
B(s) r1 r2 rn
= + + ··· + + k(s)
A(s) s − p1 s − p2 s − pn
r1 , r2 , . . . , rn sont les résidus du développement. p1 , p2 , . . . , pn sont les pôles et k(s)
est le terme direct (reste).

F
La commande residue(B,A) calcule les résidus, les pôles et le terme direct
(reste) de l’expansion du quotient B(s)/A(s). La commande s’écrit :
>>[r,p,k]=residue(B,A)
s+1
Exemple : H1 (s) =

A
s2 + 3s + 1
1) 1er méthode :
Ecrire le programme suivant pour afficher la fonction de transfert H1 (s) :
>>num=[1 1] ; den=[1 3 1] ;
>>Hs1=tf(num,den)

R
2) 2eme méthode :
Faire l’instruction suivante :
>>s=tf(’s’); Hs1=(s+1)/(s^2+3*s+1)

D
Exercice 12 :APPLICATION
Considérons la fonction de transfert suivante :
B(s) 2s3 + 5s2 + 3s + 6
H(s) = = 3
A(s) s + 6s2 + 11s + 6
1. Développer cette fonction en fraction simples par calcul manuel, en calculant
les racines du dénominateur.
2. Trouver ce résultat par la commande [r,p,k]=residue(B,A). Comparer
les résultats.
NB : APPLIQUER A LA TRANSFORMEE DE LAPLACE DES
EXERCICES DU TD

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 8

Exercice 13 :Échantillonnage d’un signal sinusoï-

dal
Le signal à échantillonner est défini par :

y(t) = sin(2πf t).

a) Écrire un programme en Matlab qui trace y(t) et y(t) échantillonné avec la

T
fréquence d’échantillonnage fs = 200hz.
Les étapes à suivre :
1. Définition de la marge de variation de 0 à 0.05 avec un pas de fs−1 (200
valeurs de 0 ; avec un pas de 5ms).

F
2. Définition de la fréquence de y(t).
3. Définition de y(t) signal à étudier
4. Tracer y(t) non échantillonné
5. Crée une nouvelle fenêtre de figure

A
6. Tracer y(t) échantillonné
b) Que remarquez-vous, interpréter vos résultats.
c) Remplacer la première ligne de votre script par fs = 1000 puis par fs = 2000.
d) Interpréter vos résultats.

R
Exercice 14 :Somme de deux sinusoïdes à des fré-
quences différentes et échantillonnées à la même

D
fréquence
Les signaux à échantillonner sont définis par :

y1 (t) = sin(2πf1 t) et y2 (t) = sin(2πf2 t)

Leur somme est

y(t) = y1 (t) + y2 (t).

Prenant la fréquence d’échantillonnage fs = 4000hz. Écrire un programme en

Matlab qui trace y(t).
a) Quelle est la fréquence du signal résultant ?
b) Es ce que le théorème de Shannon est respecté ? Expliquer ?

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 9

Exercice 14 :Produit de convolution

Soit le signal 
1, 0≤n≤2
x[n] = 
0, ailleurs
h[n] = [3, 2, 1] avec 0 ≤ n ≤ 2
1. Tracer les étapes du produit de convolution entre les deux signaux x[n] et

T
h[n]
2. pour chaque étape, calculer mathématiquement le produit de convolution
en utilisant la formule suivante :

F
+∞
X
y[n] = x[k] h[n − k]
k=−∞

La première partie du programme qui trace les deux signaux est donnée par :

x=[1 1 1];

A
h=[3 2 1];
k=0:1:length(x)-1;
figure(1);
stem(n,x,’:’,’y’);
xlim([-4 4])

R
figure(2);
stem(k,h);
xlim([-4 4])

D
Exercice 15 :Fast Fourier Transform
Le but de ce TP est l’analyse spectrale des signaux en utilisant la transformée
de Fourier rapide FFT.
Principe de FFT :
On prend les échantillons x(n) et on les divise en 2 parties :

x(2n) → Échantillons de rang pair.

x(2n + 1) → Échantillons de rang impair

Dr. Ghislain PANDRY

Travaux Pratiques Matlab/Simulink page 10

N −1
2π
x(n) e−j N kn
X
X(k) =
n=0
2π
On note : WN = e−j N

N N N
2
−1 2
−1 2
−1
(2n+1)k
x(n) WNnk x(2n) WN2nk
X X X
X(k) = + + x(2n + 1) WN
n=0 n=0 n=0

T
X(k) = X1 (k) + WNk X2 (k)

Le programme suivant présente la transformée de Fourier rapide FFT d’un

F
signal sinusoïdale de fréquence 5hz, x(t) = sin(2πf t)

Fs = 150; % Sampling frequency

t = 0:1/Fs:1; % Time vector of 1 second

A
f = 5; % Create a sine wave of 1Hz
x = sin(2*pi*t*f);

nfft = 1024; % Length of FFT

R
% Take fft, padding with zeros so that
% length(X) is equal to nfft
X = fft(x,nfft);

D
% FFT is symmetric, throw away second half

X = X(1:nfft/2);

% Take the magnitude of fft of x

mx = abs(X);

% Frequency vector
f = (0:nfft/2-1)*Fs/nfft;

% Generate the plot, title and labels.

figure(1);
plot(t,x);

title(’Sine Wave Signal’);

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 11

xlabel(’Time (s)’);
ylabel(’Amplitude ’)

figure(2);
plot(f,mx);
title(’Power Spectrum of a Sine Wave’);
xlabel(’Frequency (Hz)’);
ylabel(’Power’);

T
1. On amortie la sinusoïde par le signal y(t) = e−0.8t
Que remarquez-vous ?
2. Si on ajoute au signal précédent x(t) une sinusoïde de 50hz

F
Que remarquez-vous ?
3. Si on ajoute au signal précédent x(t) deux sinusoïdes de fréquence 50 et
70 hz
Que remarquez-vous ?

A
Exercice 16
Tracer la transformée de Fourier rapide FFT d’un signal carré périodique.

R
Exercice 17
Soit la séquence causale :

f (n) = {4, 8, 16, . . .} = 2n+2

D
a) Trouver sa transformée en Z.
b) utiliser le programme suivant pour vérifier vos résultats.
Code Matlab :

syms n
f=2^(n+2)
ztrans(f)

c) Faire les mêmes étapes pour la séquence causale :

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 12

f (n) = {0, 1, 2, 3, . . .}

Exercice 2 :
Soit la fonction en z :
z+1
F (z) =
z2
+ 0.3z + 0.02
a) développer F (z) en fraction simple.

T
b) En utilisant la fonction préétabli du Matlab residue, vérifier la décompo-
sition en fraction simple de F (z), calculer la transformée en z inverse.
c) utiliser le programme suivant pour vérifier vos résultats.

F
Code Matlab :

syms z
H=(z+1)/(z^2+0.3*z+0.02);
iztrans(H)

A
d) faire le même travail avec la fonction
1
F (z) = 2
z (z − 0.5)

R
Exercice 18
Soit un circuit RC attaqué par un échelon d’amplitude 24V, avec R = 50Ω,
C = 100µF.

D Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 13

Le condensateur est initialement déchargé.

1. Établissez le modèle simulink de ce montage.
2. Visualisez les courbes en fonction du temps, de la tension et du courant
obtenus au niveau du condensateur.
3. Établissez l’étude théorique et interprétez les résultats obtenus.

Concept : Using Symbolic Workflows

T
Symbolic workflows keep calculations in the natural symbolic form instead of
numeric form. This approach helps you understand the properties of your solution
and use exact symbolic values. You substitute numbers in place of symbolic va-

F
riables only when you require a numeric result or you cannot continue symbolically.
For details, see Choose Numeric or Symbolic Arithmetic. Typically, the steps are :
1. Declare equations.
2. Solve equations.
3. Substitute values.

A
4. Plot results.
5. Analyze results.

Workflow : Solve RLC Circuit Using Laplace Trans-

R
form
Declare Equations
You can use the Laplace transform to solve differential equations with initial

D
conditions. For example, you can solve resistance-inductor-capacitor (RLC) cir-
cuits, such as this circuit.
— Resistances in ohm : R1 , R2 , R3
— Currents in ampere : I1 , I2 , I3
— Inductance in henry : L
— Capacitance in farad : C
— AC voltage source in volts : E(t)
— Capacitor charge in coulomb : Q(t)

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 14

F T
Apply Kirchhoff’s voltage and current laws to get the following equations.

A
I1 = I2 + I3
dI1
L + I1 R1 + I2 R2 = 0
dt
Q
E(t) + I2 R2 − − I3 R3 = 0
C

R
dQ
Substitute the relation I3 = (which is the rate of the capacitor being
dt
charged) to the above equations to get the differential equations for the RLC
circuit.
dI1 R2 dQ R1 + R2

D
− =− I1
dt L dt L
dQ 1 Q R2
 
= E(t) − + I1
dt R2 + R3 C R2 + R3
Declare the variables. Because the physical quantities have positive values, set
the corresponding assumptions on the variables. Let E(t) be an alternating voltage
of 1 V.
syms L C I1(t) Q(t) s
R = sym(’R%d’,[1 3]);
assume([t L C R] > 0)
E(t) = 1*sin(t); % AC voltage = 1 V

Dr. Ghislain PANDRY

Travaux Pratiques Matlab/Simulink page 15

Declare the differential equations.

dI1 = diff(I1,t);
dQ = diff(Q,t);
eqn1 = dI1 - (R(2)/L)*dQ == -(R(1)+R(2))/L*I1

∂ R2 ∂ R1 + R2
I1 (t) − Q(t) = −I1 (t)
∂t L ∂t L
eqn2 = dQ == (1/(R(2)+R(3)))*(E-Q/C) + R(2)/(R(2)+R(3))*I1

T
Q(t)
sin(t) −
∂
Q(t) = C + R2 I1 (t)
∂t R2 + R3 R2 + R3

F
Solve Equations
Compute the Laplace transform of eqn1 and eqn2.
eqn1LT = laplace(eqn1,t,s)

RA syms I1_LT Q_LT


R2 Q(0) − s laplace(Q(t), t, s)
s laplace(I1 (t), t, s)−I1 (0)+
L
=−

(R1 + R2 ) laplace(I1 (t), t, s)

eqn1LT = subs(eqn1LT,[laplace(I1,t,s) laplace(Q,t,s)],[

I1_LT Q_LT])

I1,LT s − I1 (0) +

I1_LT Q_LT])

R2 Q(0) − QLT s
L


C
=−
I1,LT (R1 + R2 )
L
eqn2LT = subs(eqn2LT,[laplace(I1,t,s) laplace(Q,t,s)],[
L
The function solve solves only for symbolic variables. Therefore, to use solve,
first substitute laplace(I1(t),t,s) and laplace(Q(t),t,s) with the variables
I1_LT and Q_LT.

D I1,LT R2 QLT −
QLT s − Q(0) = − s2 + 1
R2 + R3 C (R2 + R3 )
Solve the equations for I1_LT and Q_LT.

Dr. Ghislain PANDRY

Travaux Pratiques Matlab/Simulink page 16

eqns = [eqn1LT eqn2LT];

vars = [I1_LT Q_LT];
[I1_LT, Q_LT] = solve(eqns,vars)

LI1 (0) − R2 Q(0) + CR2 s + Ls2 I1 (0) − R2 s2 Q(0) + CLR2 s3 I1 (0) + CLR3 s3 I1 (0) + CLR2 s I
I1,LT = 
(s2 + 1) R1 + R2 + Ls + CLR2 s2 + CLR3 s2 + CR1 R2 s + CR1 R3 s + CR2 R

Compute I1 and Q by computing the inverse Laplace transform of I1_LT and

T
Q_LT. Simplify the result. Suppress the output because it is long.
I1sol = ilaplace(I1_LT,s,t);
Qsol = ilaplace(Q_LT,s,t);
I1sol = simplify(I1sol);

F
Qsol = simplify(Qsol);

Substitute Values

A
Before plotting the result, substitute symbolic variables by the numeric values
of the circuit elements. Let R1 = 4 Ω, R2 = 2 Ω, R3 = 3 Ω, C = 1/4 F, L = 1.6 H.
Assume that the initial current is I1 (0) = 2 A and the initial charge is Q(0) = 2 C.

R
vars = [R L C I1(0) Q(0)];
values = [4 2 3 1.6 1/4 2 2];
I1sol = subs(I1sol,vars,values)

D
I1 (t) =

Q(t) =
8161

8161

Plot Results
−
+
8161

8161
+
+
− t

− t
81
200 cos(t) 405 sin(t) 16122 e 40 cosh

Qsol = subs(Qsol,vars,values)

924 sin(t) 1055 cos(t) 17377 e 40 cosh

81



√

√
1761
40

1761
40

t
8161

t

−

+
8161

Plot the current I1sol and charge Qsol. Show both the transient and steady
√
742529 1761
14195421

√
sinh

1109425 1761
30600897

state behavior by using two different time intervals : 0 ≤ t ≤ 15 and 2 ≤ t ≤ 25.

Dr. Ghislain PANDRY

√

sinh
1761
40

√
t


1761
40

t
Travaux Pratiques Matlab/Simulink page 17

subplot(2,2,1)
fplot(I1sol,[0 15])
title(’Current’)
ylabel(’I1(t)’)
xlabel(’t’)

subplot(2,2,2)
fplot(Qsol,[0 15])
title(’Charge’)
ylabel(’Q(t)’)
xlabel(’t’)

A FT subplot(2,2,3)
fplot(I1sol,[2 25])
title(’Current’)
ylabel(’I1(t)’)
xlabel(’t’)
text(3,-0.1,’Transient’)
text(15,-0.07,’Steady State’)

subplot(2,2,4)
fplot(Qsol,[2 25])
title(’Charge’)

R
ylabel(’Q(t)’)
xlabel(’t’)
text(3,0.35,’Transient’)
text(15,0.22,’Steady State’)

D
Analyze Results
Initially, the current and charge decrease exponentially. However, over the long
term, they are oscillatory. These behaviors are called “transient” and “steady sta-
te”, respectively. With the symbolic result, you can analyze the result’s properties,
which is not possible with numeric results.
Visually inspect I1sol and Qsol. They are a sum of terms. Find the terms by
using children. Then, find the contributions of the terms by plotting them over
[0 15]. The plots show the transient and steady state terms.
I1terms = children(I1sol);
I1terms = [I1terms{:}];
Qterms = children(Qsol);
Qterms = [Qterms{:}];

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 18

figure;
subplot(1,2,1)
fplot(I1terms,[0 15])
ylim([-0.5 2.5])
title(’Current terms’)

subplot(1,2,2)
fplot(Qterms,[0 15])

T
ylim([-0.5 2.5])
title(’Charge terms’)

F
Exercice 19
Soit le circuit suivant, on donne : E = 24V, R = 1.5 KΩ, C = 10nF, L = 10−4 H.

RA
D
Le condensateur est initialement déchargé.
1. Établissez le modèle Simulink du circuit.
2. Visualisez les tensions E, UC , UR selon le temps.
3. Faire l’étude théorique et établissez les diverses équations.

Exercice 1 : LAPLACE
1. Résoudre les équations suivantes :
′ ′
(a) y ” − 3y + 2y = 4t − 6; y(0) = 1; y (0) = 3
′
(b) y ” + 4y = 0; y(0) = 1; y (0) = 2
′
(c) y ” − 4y = 0; y(0) = 1; y (0) = 2

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 19

(d) y ” + ω 2 y = e(t)
(e) y ” + y = δ(t − τ ); τ > 0
′
(f) y + y = δ(t − τ ); τ > 0
′ ′
(g) y ” (t) + 6y (t) + 9y(t) = e−2t U(t) ; y(0) = 0 ; y (0) = 0
′
(h) y ” (t) − y(t) = (3e−2t + t2 + 1)U(t),avec : y(0) = 0 et y (0) = 0
2. Calculer la transformée de Laplace des fonctions suivantes :
(a) g(t) = U(t − 2)et .

T
(b) f (t) = U(t)sin2 (t).
(c) h(t) =| U(t)sin(ωt) |, ω > 0
(d) k(t) = U(t − 1)et

F
(e) l(t) = sin(t)U(t − π2 )
(f) 
1, si |t| ≤ a
m(t) =
0, sinon

A
(g) f1 (t) = (t2 + t − e−3t )U(t)
(h) f2 (t) = (t + 2)U(t) + (t + 3)U(t − 2)
(i) f3 (t) = (t2 + t + 1)e−2t U(t)
∀t ∈ R, f4 (t) = t2 + 0t sin(t − τ )f4 (τ )dτ
R
(j)
(k) ∀t > 0, f5 (t) = √1t

R
(l) ∀t > 0, f6 (t) = tβ−1 , β > 0
(m) ∀t > 0, f7 (t) = ln(t)
sinx sinx
(n) i. Du développement en série entière de , déduire que L( )=
x x
1

D
arctan .
p
ii. En déduire la transformé de Laplace.
Z x
sint
si(t) = dt
0 t
où si est le sinus intégral.
(o) Exprimer au moyen de la fonction échelon unité, la fonction définie par :

1,

 si 0 ≤ x ≤ a

2, si a ≤ x ≤ 2a

f (x) =




3, si 2a ≤ x ≤ 3a
··· , ··· ·········



Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 20

i. Donnez la représentation graphique.

ii. Donnez l’expression de la transformée de Laplace.

(p) Calculer les transformées de Laplace inverse suivantes :
i. L−1 { p+2
3
− p13 }
−2
ii. L−1 { (p+2) 2}

iii. L−1 { (p+3)(p52 +3p+5) }

T
p
iv. L−1 { (p2 +4p+6) }
p
v. L−1 { (p+1)2 }
vi. L−1 { (2p22p+3
+4p+5)
}

F
Exercice 2 : CONVOLUTION
Calculer le produit de convolution de :
1. U(t)eat ∗ U(t)ebt en distinguant les casa ̸= b et a = b.

A
2. U(t)t ∗ U(t)et
3. On considère le produit de convolution entre 2 signaux x(t) et y(t) :
Z +∞
(x ∗ y)(t) = x(τ )y(t − τ )dτ
−∞

R
Par un changement de variable adéquat, montrer que le produit de
convolution est commutatif.
4. (a) Simplifier les intégrales suivantes :
Z +∞ Z +∞
s(t)δ(t)dt; s(t)δ(t + 1)dt

D
−∞ −∞

où s(t) est un signal quelconque, causal puis non causal.

(b) Calculer la valeur numérique des intégrales suivantes :
Z +∞
r(t)δ(t − 1)dt où r(t) est la fonction rampe.
0
5. Montrer que la convolution d’un signal e(t) avec la fonction rectangle définie
par :

1 t − t0
 
h(t) =
rect
T T
(centrée sur t0 , d’amplitude 1/T et de largeur T ), avec t0 = −T /2,
correspond à un filtrage de type moyenneur.

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 21

6. (a) Déterminer la réponse indicielle (réponse à un signal échelon de Heavi-

side) d’un circuit RC dont la réponse impulsionnelle est définie par :

1 t
 
h(t) = exp −
RC RC
avec t ≥ 0 (0 pour t < 0).
(b) Représenter cette réponse impulsionnelle ainsi que la réponse du circuit.
7. Calculer la réponse d’un circuit RC à une rampe de pente 1, à partir de sa

T
réponse impulsionnelle.
8. On applique à l’entrée d’un filtre passe-bas une impulsion d’amplitude 10V
et de durée 0,00001s. Sa réponse observée en sortie est définie par y(t) =
e−3000t .

F
(a) Calculer l’aire définie par l’impulsion d’entrée et l’axe des abscisses.
L’exprimer en fonction de l’impulsion de Dirac sous la forme c · δ(t).
(b) Représenter y(t), en précisant sa pente à l’origine.
(c) En déduire la "vraie" réponse impulsionnelle du système, que l’on notera
h(t).

A
(d) Déterminer l’expression de la réponse de ce système à une entrée échelon
unité.
(e) Même question pour un signal rampe de pente 1 pour t ≥ 0 et nul pour
t < 0.
9. (a) Soient e et h deux séquences de valeurs provenant respectivement de

R
l’échantillonnage d’un signal et de la réponse impulsionnelle d’un sys-
tème :

e = {0, 0, 0, 1, 1, 1, 1, 0, 0, 0} et h = {1, −1}

D
Calculer la séquence s résultant de la convolution numérique e ∗ h.
(b) Interpréter ces résultats du point de vue des plages de fréquences élimi-
nées et conservées.
(c) Quel est le signal qui permettrait de connaître la séquence h ?
(d) Proposer une séquence h permettant de réaliser un moyennage du signal
e. Même question qu’en b).
(e) Démontrer à l’aide de cet exemple que le produit de convolution est bien
commutatif.

Exercice 3 : CORRÉLATION
1. (a) Calculer la fonction d’autocorrélation d’un signal porte défini par :

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 22

" #
t − T /2
x(t) = rect
T
(b) Conclure sur les propriétés de la corrélation utiles pour la mesure de
ressemblance.
(c) Déterminer l’énergie de ce signal à partir de sa fonction d’autocorréla-
tion.
(d) Calculer la fonction d’autocorrélation du signal carré obtenu par répé-

T
tition de la fonction porte à tous les instants kT ′ , avec T ′ = 2T et k
entier.
(e) Montrer quand dans ce dernier cas, la borne inférieure de l’intégrale
peut être différente de la valeur choisie dans la question précédente.

F
2. Soit un signal défini par :

A
x(t) = t [u(t) − u(t − T )]
T
(a) Représenter ce signal.

A
(b) Calculer sa fonction d’autocorrélation, et la représenter.
(c) Déterminer son énergie à partir de sa fonction d’autocorrélation.

Exercice 4 :SÉRIES DE FOURIER

R
1. On définit une fonction f 2π-périodique sur R.

t, si t ∈ ]−π; π[
f (t) = 
0, si t ∈ {−π, π}

D
(a) Former le développement en série de Fourier de la fonction f .
1
(b) Á l’aide de l’inégalité de Parseval, calculer la série +∞
P
n=1 2 .
n
(c) On définit une fonction g 2π-périodique sur R.

π −t

, si t ∈ ]0; 2π[

g(t) =  2
0, si t ∈ {0, 2π}

Déduire de la question a), le développement en série de Fourier de la

fonction g.
2. Soit f : R → R la fonction 2π périodique définie par

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 23

f (x) = |cos(x)|
(a) Calculer les coefficients de Fourier trigonométriques de f .
P∞ (−1)n+1
(b) En déduire la valeur de n=1 4n2 −1

3. Soit f : R → R la fonction 2π périodique telle que

∀x ∈ [−π; π], f (x) = x
(a) Calculer la série de Fourier de f .

T
(b) Étudier la convergence simple ou uniforme de la série de Fourier de f .
(c) Déterminer de ∞
P 1 P∞ 1
k=0 (2k+1)2 et k=0 (2k+1)4
P∞ 1 P∞ 1
(d) En déduire n=1 n2 et n=1 n4

F
4. Soit la fonction f : R → R 2π périodique définie par :
∀x ∈] − π; π],f (x) = exp(x)
(a) Calculer les coefficients de Fourier exponentiels de f .
(−1)n
(b) En déduire la valeur des sommes +∞
P P+∞ 1
n=1 n2 +1 et n=1 n2 +1

A
5. Soient α ∈ R Z et f : R → R la fonction 2π-périodique définie par :
f (x) = cos(αx)
(a) Déterminer les coefficients de Fourier an et bn de f .
P+∞ (−1)n−1 P+∞ 1

R
(b) En déduire les valeurs des sommes n=1 n2 −α2 et n=1 n2 −α2
(c) En déduire enfin la valeur de +∞ 1
P
n=1 n2

6. Soient α ∈ R∗ et f : R → R la fonction 2π périodique définie par :

f (t) = ch(αt) sur ] − π; π[

D
(a) Déterminer les coefficients de Fourier an et bn de f .
(b) En déduire les valeurs des sommes suivantes :
+∞ +∞
X (−1)n X 1
2 2
et 2 2
n=1 n + α n=1 n + α

7. Soit f : R −→ R la fonction régularisée, 2π périodique, impaire, constante

égale à 1 sur ]0; π[.
(a) Calculer ses coefficients de Fourier trigonométriques.
(b) Etudier la convergence simple ou uniforme de la série de Fourier vers f .
(−1)p 1
(c) En déduire +∞ et +∞
P P
p=0 p=0
2p + 1 (2p + 1)2

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 24

P+∞ 1 P+∞ (−1)n−1

(d) Calculer n=1 et n=1
n2 n2
sin(t) + |sin(t)|
8. On considère la fonction f définie sur R par f (t) =
2
(a) Préciser le mode de convergence de la série de Fourier de f .
(b) En déduire +∞
P 1 P+∞ 1
n=1 4n2 −1 et n=1 (4n2 −1)2
2x
9. f : R → R 2π-périodique paire telle que ∀x ∈ [0; π], f (x) = 1 − π
(a) Développer en série de FOURIER de f ;

T
(b) En déduire +∞ , +∞
P n=1 P 1 P+∞ 1
0 (2n+1)2 n=1 n2 et n=1 n4

10. f : R → R 2π-périodique telle que ∀x ∈] − π; π], f (x) = sin( x2 )

(a) Développer en série de FOURIER de f ;

F
(b) En déduire +∞ n 2n+1
P
n=0 (−1) 16n2 +16n+3
11. f : R → R telle que ∀x ∈ R, f (x) = sup(0, sinx)
(a) Développer en série de FOURIER de f ;
(b) En déduire +∞ 1
P
n=1 4p2 −1

12. Développer en série de FOURIER la fonction f : x 7→ x − E(x) − 1/2

A
13. Soit α ∈ C \ Z. Soit f l’application de R → C 2π-périodique telles que
∀x ∈ [−π; π], f (x) = cos(αx).
(a) Développer la fonction f en série de FOURIER.
(b) En déduire que pour tout z ∈ C \ Z sin(πz)π
= z1 + +∞ n 2z
P
n=1 (−1) z 2 −n2 et

R
πcotan(πz) = z22z−n2
14. Soit f : R −→ R, 2π-périodique, paire, telle que, pour tout t ∈ [0; π] :
π π
f (t) = 1 si 0 ≤ t < , f (t) = 0 si t = ,
2 2

D
π
f (t) = −1 si < t ≤ π.
2
(a) Vérifier f ∈ CM2π et calculer les coefficients de Fourier (trigonomé-
triques) de f .
(b) Étudier les convergences de la série de Fourier de f et préciser sa somme.
(c) En déduire les sommes de séries suivantes :
+∞ +∞ +∞
X (1)p X 1 X 1
, 2
, 2
.
p=0 2p + 1 p=0 (2p + 1) n=1 n

15. Soit f : R −→ R, 2π-périodique, impaire, telle que :

π π
f (t) = t si 0 ≤ t < , f (t) = π − t si ≤ t ≤ π.
2 2
Dr. Ghislain PANDRY
 Travaux Pratiques Matlab/Simulink page 25

(a) Vérifier f ∈ CM2π et calculer les coefficients de Fourier (trigonomé-

triques) de f .
(b) Étudier les convergences de la série de Fourier de f et préciser sa somme.
(c) En déduire les sommes de séries suivantes :
+∞ +∞ +∞ +∞
X 1 X 1 X 1 X 1
2
, 2
, 4
, 4
.
p=0 (2p + 1) n=1 n p=0 (2p + 1) n=1 n

T
16. Soit f : R −→ R, t 7−→ | sin t|.
(a) Vérifier f ∈ CMπ et calculer les coefficients de Fourier (trigonomé-
triques) de f .
(b) Étudier les convergences de la série de Fourier de f et préciser sa somme.
(c) En déduire les sommes de séries suivantes :

F
+∞ +∞ +∞
X 1 X (−1)n X 1
, , .
n=1 4n − 1 n=1 4n − 1 − 1)2
2 2 (4n2
n=1

17. Soit f : R −→ R, 2π-périodique, impaire, telle que : ∀t ∈ [0; π], f (t) =

A
t(π − t).
(a) Vérifier f ∈ CM2π et calculer les coefficients de Fourier (trigonomé-
triques) de f .
(b) Étudier les convergences de la série de Fourier de f et préciser sa somme.
(c) En déduire les sommes de séries suivantes :

R
+∞ +∞ +∞
X (−1)p X 1 X 1
3
, 6
, 6
.
p=0 (2p + 1) p=0 (2p + 1) n=1 n

18. Soit f : [−π; π] −→ C continue telle que : ∀n ∈ Z,

D
Z π
f (t)eint dt = 0.
−π

Montrer : f = 0.
19. Montrer qu’il existe une suite réelle (αn )n∈N telle que : ∀t ∈ R,
+∞
αn cos2 (nt),
X
| cos t| =
n=0

et déterminer une telle suite (αn )n∈N .

20. Soit λ ∈]0; +∞[ fixé. On considère l’application f : R −→ R, 2π-périodique,
telle que : ∀t ∈] − π; π], f (t) = ch(λt).

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 26

(a) Vérifier f ∈ CM2π et calculer les coefficients de Fourier (trigonomé-

triques) de f .
(b) Étudier les convergences de la série de Fourier de f et préciser sa somme.
(c) En déduire les sommes de séries suivantes :
+∞ +∞ +∞
X (−1)n X 1 X 1
2 2
, 2 2
, .
n=1 λ + n n=1 λ + n n=1 (λ2 + n 2 )2
2π
21. Soient T ∈]0; +∞[, ω = , f : R −→ C, T -périodique, de classe C 1 , telle

T
T
que : Z 2π
∀n ∈ {−1, 0, 1}, f (t)einωt dt = 0.
0

Montrer :

F
1
||f ||2 ≤ ||f ′ ||2 ,
2
où : !1
1ZT 2
2
||f ||2 = |f (t)| dt ,
T 0

A
et de même pour ||f ′ ||2 .
22. Soit f : R −→ C, 2π-périodique, continue. On suppose : ∀k ∈ Z,
Z 2π
f (t)ei(2k+1)t dt = 0.
0

R
Montrer que f est π-périodique.
23. Soient T > 0, f : R −→ C, T -périodique, de classe C 1 par morceaux,
continue. Montrer :
2
Z T
T2 Z T ′ 2 1 Z T

D
′ 2
|f | ≤ 2 |f | + f .
0 4π 0 T 0

24. Soit f : R −→ C, 2π-périodique, de classe C 2 par morceaux, de classe C 1 .

Montrer : Z 2π Z 2π Z 2π
′ 2 ′′ 2
4 |f | + 2 |f | ≥ 5 |f ′ |2 .
0 0 0

25. Soit (γn )n∈Z une suite (indexée par Z) à termes dans C.
On note, pour tout p ∈ N : Sp : R −→ C, t 7−→ pk=−p γk eikt .
P

On suppose que la suite (Sp )p∈N converge uniformément sur R vers une
application notée f .
Démontrer que f est 2π-périodique, continue, et que : ∀n ∈ Z, cn (f ) =
γn .

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 27

26. Soient (αn )n≥0 , (βn )n≥1 deux suites réelles telles que la suite d’applications
(Sn )n∈N définie par :
n
α0 X
∀t ∈ R, Sn (t) = + (αk cos(kt) + βk sin(kt))
2 k=1

converge uniformément sur R vers une application notée f .

(a) Montrer que f est 2π-périodique et continue sur R.
(b) Établir : (∀n ≥ 0, an (f ) = αn ) et (∀n ≥ 1, bn (f ) = βn ).

T
1
27. Soient z ∈ C tel que |z| < 1, et f : R −→ C, t 7−→ 1+ze it . Vérifier f ∈ CM2π

et calculer les coefficients de Fourier exponentiels de f .

1
28. Soit a ∈]0; +∞[. On note f : R −→ R, t 7−→ cha+cos t
.
(a) Vérifier f ∈ CM2π et déterminer les coefficients de Fourier (trigonomé-

F
triques) de f .
(b) En déduire, pour tout n ∈ N :
Z π
cos(nt) π(−1)n e−na
dt =
0 cha + cos t sha

A
et Z π
sin(nt)
dt = 0.
0 cha + cos t
(c) Calculer : Z π
1
I= dt.

R
0 (cha + cos t)2
29. Soient x ∈ [0; +∞[, f : R −→ R, 2π-périodique, telle que f (π) = 0 et :

∀t ∈] − π; π[, f (t) = sh xt.

(a) Vérifier f ∈ CM2π et calculer les coefficients de Fourier (trigonomé-

D
triques) de f .
(b) Étudier la convergence de la série de Fourier de f , et montrer :
+∞
X 2(−1)n+1 n sh πx
∀t ∈] − π; π[, sh xt = sin(nt).
n=1 π(n2 + x2 )
(c) En déduire :
Z +∞
cos(xt) π
dt = .
0 ch t 2 ch πx
2

30. Déterminer l’ensemble des applications f : R −→ C, 2π-périodiques, de

classe C ∞ , telles qu’il existe M ∈ R+ tel que :
∀n ∈ N, ∀x ∈ R, |f (n) (x)| ≤ M.

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 28

31. Soit α ∈]1; +∞[.

(a) Montrer :
1 1
Z +∞
dt Z 1 α−1
u + uα
α = du.
0 tα + 1 0 1+u
(b) En déduire :
+∞
X (−1)n+1 2
Z +∞
dt
α = α + 1 .
0 tα + 1 n=1 n − α2
2

T
(c) Établir :
+∞
X 2(−1)n+1 x 1 1
∀x ∈]0; 1[, = − ,
n=1 π(n − x )
2 2 sin(πx) πx

F
en étudiant, pour x ∈]0; 1[ fixé, la fonction f : R −→ R, 2π-périodique,
telle que f (t) = cos(xt) si t ∈] − π; π[.
(d) Démontrer :
Z +∞
dt π
α
=  .
0 t +1 α sin π α

A
(e) En Rdéduire les valeurs des intégrales suivantes :
x −1
(1) 0+∞ t1+t dt, x ∈]0; 1[.
R +∞ x−2
(2) R0 t ln(1 + t) dt, x ∈]0; 1[.
+∞ eat 3
(3) −∞ bt +ect dt, (a, b, c) ∈ R , b < a < c.
R +∞ e eat
(4) −∞ ch(ct) dt, (a, c) ∈ R2 , |a| < c.

R
R +∞ ch(at)
(5) 0 ch(ct)
dt, (a, c) ∈ R2 , |a| < c.
32. Soit (αn )n≥0 une suite à termes dans R+ , convergeant vers 0.
P
(a) Montrer qu’il existe une extractrice σ telle que la série n≥0 ασ(n)
converge.

D
(b) En déduire qu’il existe f : R −→ R, 2π-périodique, continue, telle
que, en notant an (f ), bn (f ) (pour n ∈ N) les coefficients de Fourier
trigonométriques de f , il existe une infinité de n ∈ N tels que : |an (f )| +
|bn (f )| > αn . (Utiliser l’exercice 4.26.)

Exercice 5 : TRANSFORMÉE EN Z
1. On pose xn = cos(ωnT ), ω et T fixés. Calculons la transformée en Z de :
(a) Z(xn )
(b) Z(nxn ), pour k = {1; 3}
d k
Z(nxn ), Rappel : Z(nk xn ) = (−z dz ) Z(xn )

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 29

4z
2. Déterminer la suite x dont la transformée en Z est f (z) = 3z 2 −2z−1

z2
3. Déterminer la suite x dont la transformée en Z est f (z) = (z+3)2

4. On considère une suite x = (xn )n∈N . On définit la suite y = (yn )n∈N par
yn = nk=0 xk
P

T
(a) Déterminer une équation récurrente entre les suites x et y.
(b) En déduire la transformée en Z de y (en fonction de celle de x).

Exercice 6 : TRANSFORMÉE DE FOURIER

F
1. Résoudre l’équation suivante :

Z ∞ sin(α), si 0 ≤ α ≤ π
f (x)cos(αx)dx =
0 0, si α > π

A
2. Soit le signal x(t), défini par :

1 + t,

 si − 1 ≤ t ≤ 0
x(t) = 1 − t, si 0 ≤ t ≤ 1

R



0, sinon

(a) Représenter le signal x(t).

(b) Déterminer :
i. le nom du signal x(t)

D
ii. la période du signal x(t)
(c) Représenter les signaux suivants :
i. x(t + 1) ;
ii. x(t − 1) ;
iii. U(t)x(t)
(d) Déterminons la transformée de Fourier des signaux suivants :

i. x(t − 1/2))
ii. x(2t)
d2 x(t)
iii. dt2

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 30

(a) Représenter le signal x(t).

(b) Déterminons :
i. le nom du signal x(t) : c’est un signal triangulaire.
ii. la période du signal x(t) : sa période vaut T = 2.
(c) Représenter les signaux suivants :

T
i. x(t + 1),
(d) Déterminons la transformée de Fourier des signaux suivants :

i. x(t − 1/2) :

F
ii. x(2t) :
Pour a = 2,F ({x( 12 t)})(ν) = 1
|2|
sin2c ( πν
2
)
d2 x(t)
iii. dt2

A
3. On considère le signal suivant : x[n] = α|n| , n ∈ Z et α ∈ ]−1; 1[ , α ̸= 0
(a) Calculer xb(ν) la transformée de Fourier de x(n)
4. On considère le signal défini pour a > 0 par fa (t) = U(t)e−at et les signaux

R
ga et ha tels que :ga (t) = fa (t) + fa (−t) et ha (t) = fa (t) − fa (−t)
(a) Traçons les graphes de fa ,ga ,ha et calculons leur transformée de Fourier.
(b) En déduire la valeur des intégrales 0∞ cos(ωx) dx et 0∞ xsin(ωx)
R R
1+x2 1+x2
dx.
cos(ωx) xsin(ωx)
(c) En déduire la valeur des intégrales 0+∞ 2 dx et 0+∞ 2
R R
2
dx.
x +a x + a2

D
5. En considérant la transformée sinus et cosinus de la fonction f (t) = e−kt ,
k > 0,
(a) évaluer les deux intégrales
R ∞ cos(kα)
i. A = 0 k2 +α2
dα
R ∞ αsin(kα)
ii. B = 0 k2 +α2
dα
(b) Utiliser l’égalité de Parseval pour évaluer les deux intégrales suivantes :
R∞ 1
i. C = 0 (k2 +α2 )2
dα
R∞ α 2
ii. D = 0 (k2 +α2 )2 dα

(c) En déduire la transformée sinus et cosinus de la fonction g(t) = sin(pt)e−kt

et la fonction h(t) = cos(pt)e−kt , où p est un réel donnée.

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 31

t−τ 2
−π( )
6. Soit gτ σ (t) = e σ , t ∈ R, paramétrée par τ ∈ R et σ ∈ R∗+ . On
2
admettra le résultat suivant : R e−πx dx = 1. Dans un premier temps on
R

s’intéresse à la gaussienne centrée réduite g = g01 .

2
(a) Montrer simplement que sa transformée de Fourier s’écrit :G(f ) = 2 0+∞ e−πt cos(2πf t)dt
R

(b) En dérivant sous le signe intégrale par rapport à la fréquence puis en

intégrant par partie, déterminez une équation différentielle simple en
G(f ).

T
(c) Résolvez l’équation différentielle et déterminez les constantes d’intégra-
2
tion en utilisant R e−πx dx = 1.
R

(d) Déterminez finalement la transformée de Fourier Gτ σ de gτ σ

F
Exercice 7 :
Résoudre
′
1. (U(t)sin(t))

A
′′
2. (U(t)sin(t))

Exercice 8 : Automatique
1. Soit un système placé en boucle fermée (BF) à retour unitaire. Ce système

R
est modélisé par une fonction de transfert.
K
G(p) = .
p(p2 + p + 3)

D
(a) Dessiner le système en BF.
(b) Donner la fonction de transfert du système en BF.
2. Soit la fonction de transfert suivante :
(p + 1)(p + 100)
G(p) =
p + 10
Tracer les diagrammes de Bode.
3. Soit la fonction de transfert suivante :
1000(p + 1)
G(p) =
(p + 10)2
Tracer le diagramme de Nyquist.

Dr. Ghislain PANDRY

 Travaux Pratiques Matlab/Simulink page 32

4. On considère un système de fonction de transfert G(p) avec :

K
G(p) = avec T = 0,1 s.
1 + Tp

(a) Calculer l’expression précise de la pulsation de coupure à 0 dB définie

par :
G(ωc0 ) = 1.

T
(b) Montrer que si K ≫ 1, on a :

K
ωc0 ≈ .
T

F
(c) Calculer la valeur du gain K qui permet d’obtenir une pulsation ωc0 =
10 rad/s.
5. On considère un système de fonction de transfert G(p) avec :

K
G(p) = avec K > 1.

A
2
p 2ξp
2
+ +1
ωn ωn

(a) Calculer en fonction de K, ωn et ξ, dans les trois cas ξ < 1, ξ = 1 et

ξ > 1, l’expression de la pulsation de coupure à 0 dB définie par :

R
G(ωc0 ) = 1.

(b) Montrer que si K ≫ 1, on a, dans tous les cas :

√
ωc0 ≈ ωn K.

D Dr. Ghislain PANDRY