Réponse fréquentielle d'un système discret sous matlab

bonjour,

je veux visualiser la réponse fréquentielle d'un système discret défini par :
x[n+1] = Ax[n] + Bu[n]
y[n] = Cx[n] + Du[n]

sachant que le vecteur u est une force sinusoïdale avec laquelle on excite notre système.
dt=0.05;
Tmp=[0: dt:20];
n_pas=length(Tmp);
u=0.5*sin(2*pi*F*(Tmp+dt))

j'arrive à visualiser la réponse temporelle par la fonction "dlsim" pour une fréquence F constante
[y,x]=dlsim(A,B,C,D,u(:,1:n_pas));

je veux visualiser sa réponse fréquentielle pour une plage de la fréquence F
 

tout exemple e programme même simple est le bienvenue
merci