Scribd
Importer
7 vues10 pages

Analyse de la réponse d'un système RC

Ce document contient deux programmes qui analysent le comportement temporel de systèmes du premier ordre. Le premier programme détermine la réponse impulsionnelle, indicielle et à une rampe d'un système. Le deuxième programme étudie l'influence de la constante de temps sur ces réponses.

Transféré par

brbouchra091
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats RTF, PDF, TXT ou lisez en ligne sur Scribd
7 vues10 pages

Analyse de la réponse d'un système RC

Ce document contient deux programmes qui analysent le comportement temporel de systèmes du premier ordre. Le premier programme détermine la réponse impulsionnelle, indicielle et à une rampe d'un système. Le deuxième programme étudie l'influence de la constante de temps sur ces réponses.

Transféré par

brbouchra091
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats RTF, PDF, TXT ou lisez en ligne sur Scribd

programme 1 :

clear;

clc;

E=40;

R=50;

c=63*10^-6;

a=2;

tau=R*C ;

% Analyse temporelle d'un système de 1 ^ pi ordre %

t=0:0.001:0.1;

% Définir un vecteur de fin points linéairement %

% espacé d'un pas entre init et fin %

% Détermination de fonction de transfert %

num=[k];den=[tau 1];

% Création de la fonction de transfert %

sys = U(num, den) ;

% Réponse impulsionnelle d'un système de premier ordre %

y1=impulse(sys,t);

% Réponse à une impulsion unité %

figure(1)

% Ouverture d'une figure %

plot(t.y1)

% Traçage de la réponse impulsionnelle %

% Réponse indicielle d'un système de premier ordre %

y2=step(sys.t)
% Réponse indicielle du système %

figure(2)

plot(t,y2)

% Détermination de la constante de temps graphiquement (tau) à partir de la réponse indicielle %

hold on % Maintenir les graphes courants %

ym=max(y2); % Détermine la valeur maximale de y2([Link]) %

yt=0.63*ym; % Détermine la valeur de yt correspondante au temps de réponse %

max_idx=min(find(y2>=yt)); % find: permet de trouver l'indice de l'élément auquel y2>=yt %

tau=t(max_idx); % tau est le temps correspondant à cet indice dans le vecteur t %

plot([tau tau]‚[0 yt]‚’m--‘) % Tracer une ligne verticale à tr %

stf=sprintf(‘tau\n%3.2fs’,tau) % Créer l'étiquette volante %

gtext(str) %Permet de cliquer dans la fenêtre pour positionner le texte %

hold off

% Réponse à une rampe d'un système de premier ordre %

u=a*t;

% Détermination de la rampe u de pente a %

y3=lsim(sys.u,t)

% Réponse du système à une entrée u quelconque (t :temps) %

figure(3)

plot(t,y3);
programme 2 :

E = 40 ;

R = 50;

C=63*10^-6 ;

tau=R*C;

tau1=0.0016;

tau2=0.0032;

tau3=0.0064;

a =2;
k=1;

t=0:0.0001:0.1;num=[k];

den=[tau 1];

sys=tf(num, den);

num=[k]; den=[tau 2];

sys2=tf(num, den);

num =[k]; den=[tau 3];

sys3=tf(num, den);

figure (1)

y1=impulse(sys,t);

plot(t,y1)

y2= step(sys,t)

figure (2)

plot(t,y2);

hold on

ym=max(y2);

yt=0.63*ym;

max_idx=min(find(y2>=yt));

tau=t(max_idx);

plot([tau tau],[0 yt],'m--')

str-sprintf('tau\n%3.5f s', tau)

gtext(str)

hold off

tm=2.2*tau

trp=3.5*tau
u=a*t;

y3=lsim(sys,u,t);

figure (3)

plot(t,y3);

y4=step(sys,t);

y5=step(sys2,t);

y6=step(sys3,t)

figure (4)

hold on

plot(t,y4,'r');

plot(t,y5,'g');

plot(t,ye,'y');

hold

off

Vous aimerez peut-être aussi