Solução da prova de modelagem de sistemas Dinâmicos

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

clc
clear 
close
bebe = 1:5;
bebe=bebe';
w = logspace(3, 4, 100);%Vetor de Frêquencia em Rad/seg, vai de 10^0 à 10^4, com 100 pontos
w = w';
% f = w./(2*pi);
f = 50:50:6000; f=f';
for i=1:length(f)
 
 tamos = 1/(10*f(i));
 tfinal = (5/f(i))+0.1;
 temp = 0:tamos:tfinal; temp = temp';
 A = [ sin(2*pi*f(i)*temp) cos(2*pi*f(i)*temp) ];
 saida = gravador(bebe(2),tamos,tfinal);
 saida = saida';
 theta = (A'*A)\(A'*saida);
 G(i) = sqrt(theta(1)^2 + theta(2)^2);%Ganho
 
end
G=G';
ganhodb = 20*log(G); % Ganho em db
figure(2)
plot(f, G, 'b')
title('Diagrama de bode - Amplitude')
xlabel('Frequêcia (hz)')
ylabel('Ganho (db)')
grid on
 
PosMagnitude = find(G == max(G));%Posição onde está a maior Amplitude
disp('A frequência em (Hz) que possui a maior magnitude é: ' )
disp(f(PosMagnitude))