%%% ECA004 - Sistemas de Controle Clssico% Laboratrio 10 - Resposta em Frequncia% Clara D.%% Exerccio 01% Obtenha o Diagrama de Bode de G(s) = 20s+2/s+4s+25num = [0 20 2]; % vetor do numerador da funo de transfernciaden = [1 4 25]; % vetor do denominador da funo de transfernciasys = tf(num,den); % funo de transfernciaw = logspace(-3,4); % especifica a faixa de frequnciafigure(1) % abre janela para plotar o grficobode(num, den, w) % funo do Diagrama de Bodetitle('Diagrama de Bode de G(s) = 20s+2/s+4s+25') % titula a figura(1)%% Exerccio 02% Para o sistema do Exemplo 6.1 acrescente um sensor com valor de 0,02% na realimentao. obtenha o diagrama de Bode para a funo de% transferncia de malha fechada. Obtenha tambm o pico de ressonncia, e a% frequncia de ressonncia.% H = 0.02% Sistema com H:num = [0 0 0 1]; % vetor do numerador da funo de transfernciaden = [0.5 1.5 1 0.02]; % vetor do denominador da funo de transfernciasys1 = tf(num,den); % funo de transfernciasys = feedback(sys1,1);w = logspace(-1,1); % especifica a faixa de frequnciafigure(2) % abre janela para plotar o grficobode(sys,w) % % funo do Diagrama de Bodetitle('Diagrama de Bode de G(s) = 1/0.5s+1.5s+s+0.02') % titula a figura(2)[mag,phase,w] = bode(sys,w);%incorpora os pontos de referencia especificado[Mp,k] = max(mag); % valor de mxima aplituderesonant_peak=20*log10(Mp) % pico de ressonnciaresonant_frequency=w(k) % frequncia de ressonncia