Discussion :

[pelotudo]

Bonjour,

je voudrais faire le réponse impulsionelle de la fonction h() suivante (avec W=13):

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
function h = h(i,W)
 
if (i>= -W & i<=0)
    h=-f(-i);
elseif (i>= 1 & i<=W)
    h=f(i);
end
 
function f = f(x)
 
A=0.2208;
K=[1.583 1.468 -0.078 -0.036 -0.872 -0.56];
s=0.5385;
 
f=exp(A*x)*(K(1)*sin(A*x)+K(2)*cos(A*x))+exp(-A*x)*(K(3)*sin(A*x)+K(4)*cos(A*x))+K(5)+K(6)*exp(s*x);

Je ne sais pas comment faire. Normalement, on a une fonction de transfert bien définie pour tracer la réponse impulsionelle. Là on a une fonction h() qui dépend d'une condition...

[rostomus]

Salut,

Tu veux faire comme ça?

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
function h = rep_imp(i,W)
 
h=zeros(size(i));
h(i>= -W & i<=0)=-f(-i(i>= -W & i<=0));
h(i>=0 & i<=W)=f(i(i>= 0 & i<=W));
 
 
 
function f = f(x)
 
A=0.2208;
K=[1.583 1.468 -0.078 -0.036 -0.872 -0.56];
s=0.5385;
 
f=exp(A*x).*(K(1)*sin(A*x)+K(2)*cos(A*x))+exp(-A*x).*(K(3)*sin(A*x)+K(4)*cos(A*x))+K(5)+K(6)*exp(s*x);

et par exemple:

Code :

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

1
2
3
4
 
t=-15:.1:15;
h=rep_imp(t,13);
plot(t,h)

[pelotudo]

Oui c'est ce que je trouve aussi avec mon code.
Mais, c'est la réponse impulsionnelle de h() ça? On ne devrait pas faire un système et l'utiliser avec la foncion step() pour avoir la réponse impulsionnelle?
En fait, je suis sensé trouver une réponse du type de la figure adjointe.