Solucao da P1 de Modelagem de sistemas dinamicos

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

clear all
close all 
clc
% --------------------------------------------------------------
fmax =600; % Definido a partir do enunciado
tamos = 1/(10*fmax);
t = [0:tamos:10]; t=t';
f = logspace(0,3,100); %Vetor de Frequências
w = 2*pi*f;% Frequência em Rad/s
G = [];
 for k = w 
 u = sin(k*t);
 A = [sin(k*t) cos(k*t)];
 y = circuito3(u,tamos);%Subistitui o Circuito para ver o gráfico
 theta = (A'*A)\(A'*y);
 ymod = A*theta;
 ganhov = sqrt(theta(1)^2+theta(2)^2);%Ganho em V/V
 G = [G 20*log10(ganhov)]; %Ganho em db
 end
semilogx(f,G)
grid on
disp('O Circuito ideal é o 3, pois trata-se de um filtro passa Faixa que é o mais adequado a situação.')
% vv= exp(-14/20)%Ganho em V/V
% fmax = 1000; %Definindo a frequencia máxima
% tamos = 1/(fmax*10);
% 
% f = logspace(-1,3,100); f=f'; % vetor de frequências para visualizar o sinal;
% w = 2*pi*f;
% tfinal = 10*(1/fmax); % 5 CICLOS DO SINAL DE ENTRADA
% tempo = 0:tamos:tfinal; tempo=tempo'; % vetor de tempo
% % 
% % ts = 0.0001;
% % t = [0:ts:10]';
% % 
% % freq = logspace(-1,3,100); %espaço para o diagrama de bode
% % 
% % 
% % for n= 2*pi*freq
% % u = sin(n*t);
% % A = [sin(n*t) cos(n*t)];
% % y = filtro(u,ts);
% % theta = (A'*A)\(A'*y);
% % ymod = A*theta;
% % gain = [gain 20*log10(sqrt(theta(1)^2+theta(2)^2))];
% % end
% % 
% % 
% % semilogx(freq,gain)
% % grid on
% 
% 
% A=[];
% for k = w
% 
% u = sin(w*tempo);
% A = [ sin(w*tempo) cos(w*tempo) ];
% [y] = circuito3(u,tamos);
% theta = (A'*A)\(A'*y);
% 
% G = sqrt(theta(1)^2 + theta(2)^2);
% F = atan2d(theta(2),theta(1));
% 
% plot(tempo, y, 'b')
% end
% 
% G=G';
% F=F';
% 
% ganhodb = 20*log(G);% Ganho em db
% fase = F*(180/pi); % Fase em rad/s
% figure(1)
% plot(f, ganhodb, 'k' )