Discussion :

[virtuadrack]

bonjour,

je possède l'excitation qui suit

3.5y'(t) + 5.8y(t) = -1.3x'(t) + 6.2x (t)

et j'ai trouvé la réponse impulsionnelle

HH=.43*exp(-0.603*t)*u -1.08*A;

(u représente la 'unit step function' et A représente la réponse impulsionnelle)

le problème c'est que je ne comprend pas pourquoi le graphique j'obtiens n'est pas a 100% identique a celui que je dois avoir...

l'image qui suit représente les 2 graphiques, en premier le graphe correcte (celui que je dois obtenir) et celui qui suit est celui que j'ai actuellement.



Voila mon code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
wmax=20;
dw=0.1;
w= [-wmax:dw:wmax]';
t=0:1/120:16;
 
if t == 0
   A = 1;
else 
   A = 0;
end
if t >= 0
    u = 1;
else
    u = 0;
end
 
HH=.43*exp(-0.603*t)*u -1.08*A;
 
%Impulse Response of H 
 
p= plot(t, HH,'k'); 
set (p, 'LineWidth',2);
grid on;
xlabel('Impulse Response, {\omega}', 'Fontsize', 18, 'FontName','Times');
ylabel('h(t)','Fontsize', 18, 'FontName','Times');

Merci

[Jerome Briot]

Ce ne serait pas plutôt ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
wmax=20;
dw=0.1;
w= [-wmax:dw:wmax]';
t=0:1/120:16;

A = double(t==0);
u = double(t>=0);

% if t == 0
%    A = 1;
% else 
%    A = 0;
% end
% if t >= 0
%     u = 1;
% else
%     u = 0;
% end

HH=.43*exp(-0.603*t).*u -1.08*A;

%Impulse Response of H 

p= plot(t, HH,'k'); 
set (p, 'LineWidth',2);
grid on;
xlabel('Impulse Response, {\omega}', 'Fontsize', 18, 'FontName','Times');
ylabel('h(t)','Fontsize', 18, 'FontName','Times');

Mais je ne suis pas un spécialiste, donc c'est peut être totalement faux

[FR119492]

Salut!
A vue de nez, tu dois avoir deux fois la mème fonction, mais avec des échelles différentes. Essaie de les mettre les deux sur le même dessin.
Jean-Marc Blanc

[Antonin08]

Salut

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
if t >= 0
    u = 1;
else
    u = 0;
end

Ce ne serait pas plutôt

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
if t > 0
    u = 1;
else
    u = 0;
end

alors :
A t= 0, u =0, A=1
HH(t=0) = 0.43*exp(0)*0 - 1.08*1 = -1.08 CQFD

[Jerome Briot]

Euhh y'a rien qui vous gêne dans le code ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
>> t=0:1:16
 
t =
 
     0     1     2     3     4     5     6     7     8     9    10    11    12    13    14    15    16
 
>> t>0
 
ans =
 
     0     1     1     1     1     1     1     1     1     1     1     1     1     1     1     1     1

[virtuadrack]

Salut



Ce ne serait pas plutôt

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
if t > 0
    u = 1;
else
    u = 0;
end

alors :
A t= 0, u =0, A=1
HH(t=0) = 0.43*exp(0)*0 - 1.08*1 = -1.08 CQFD

ca ne donne pas la bonne courbe, c'est plutot une ligne droite...

il n'y a pas d'alternative

je sais qu'il faut utiliser cette function y(t)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
function y = u(t)
	zro = t == 0 ; pos = t > 0 ; y = zro/2 + pos ;

mais je sais pas comment l'utiliser. des idees ?

[phryte]

Salut.

celui que je dois obtenir

Il faut vectoriser :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
t=0:1/120:16;
A=zeros(1,length(t));
A(1)=1;
 u=ones(1,length(t)); 
 u(1)=0;
HH=.43*exp(-0.603*t).*u -1.08*A;
 
%Impulse Response of H 
 
p= plot(t, HH,'k'); 
set (p, 'LineWidth',2);
grid on;
xlabel('Impulse Response, {\omega}', 'Fontsize', 18, 'FontName','Times');
ylabel('h(t)','Fontsize', 18, 'FontName','Times');
axis([-1 20 -2 1])