analyse fréquentielle

[invite9b8a9ef4]

bonjour,

je veux réprésenter la réponse fréquentielle d'un système discret décrit par
x[n+1] = Ax[n] + Bu[n]
y[n] = Cx[n] + Du[n]
avec le vecteur u est de forme sinusoïdale
dt=0.05;
Tmp=[0:dt:20];
n_pas=length(Tmp);
u= 0.5*sin(2*pi*F*(Tmp+dt))

j'arrive à représenter la réponse temporelle pour une fréquence F donnée par
[y,x]=dlsim(A,B,C,D,u(:,1:n_pas));
je veux représenter x et y en fonction de la fréquence F variable
merci

-----

[JPL]

C'est des maths ou de la programmation ? Parce que si c'est de la programmation il faudrait dire avec quel outil tu veux programmer

[invite9b8a9ef4]

C'est des maths ou de la programmation ? Parce que si c'est de la programmation il faudrait dire avec quel outil tu veux programmer

c est un programme en matlab
en fait, avec plot(Tmp, x) j arrive à visualiser la réponse temporelle
mais je veux un prpgramme pour visualiser la réponse fréqquentielle

[inviteb9f49292]

Si tu veux visualiser le spectre, utilise la fonction fft...
Mais si j'ai bien compris, c'est plus compliqué que ça, tu veux prévoir le spectre pour n'importe quelle fréquences...
En utiliseant la linéarité de la transformée de Fourier, et en réfléchissant un peu à ce qu'est x[n+1] (et y[n+1]) en fonction uniquement de u[n], et sauf erreur de ma part, c'est assez facile, il n'y a pas besoin de MATLAB pour ça

[A voir en vidéo sur Futura]

[invite9b8a9ef4]

bonjour,

le système s écrit aussi :
x' = Ax + Bu
y = Cx + Du

je dois continuer de travailler avec matlab et pour la fréquence, je veux varier F de 1 à 100 par exemple.
et arriver enfin à visualiser x en fonction de F