Gráficos códigos PH CALDEIRA

Prévia do material em texto

1 
 
UNIVERSIDADE FEDERAL RURAL DO SEMI-ÁRIDO 
CENTRO DE ENGENHARIAS 
CURSO DE ENGENHARIA ELÉTRICA 
DISCIPLINA: ANÁLISE DE SINAIS E SISTEMAS 
 PROF. Msc. ISAAC BARROS TAVARES DA SILVA 
 ALUNO: PAULO HENRIQUE DOS REIS CALDEIRA 
 
 
 
 
 
 
 
 
PLOTANDO GRÁFICOS DE FUNÇÕES NO DOMINIO DA FREQUENCIA COM O 
SOFTWARE MATLAB® 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
MOSSORÓ 
2017 
2 
 
 
 
%Grafico do modulo e fase da função X(w) 
% X(w)= (1/root(a^2+w^2))*exp(-j*arctan(w/a)) 
 
clear; 
clc 
 
a = 2; 
w=(-pi:0.01:pi); 
for i=(1:length(w)); 
 x(i)= 1/(a+(j*w(i))); 
end 
subplot(2,1,01); 
plot(w,abs(x),'g','linewidth',2) 
title('Gráfico do Módulo |X(j\omega)| ') 
xlabel('Frequência (\omega)') 
grid on; 
 
subplot(2,1,02) 
plot(w,angle(x),'r','linewidth',2.5); 
title('Gráfico da Fase \angleX(j\omega)') 
xlabel('Frequência (\omega)') 
grid on 
 
3 
 
 
 
% GRAFICO DA FUNÇÃO 
% X(w)= T*sinc(w*t/2) 
% com w variando de -pi a pi 
 
clear; 
clc 
 
 
w=(-pi:0.01:pi) 
T=2 
for j=(1:length(w)); 
x(j)= T*sinc((w(j)*T)/2); 
end 
hold on; 
plot(w,abs(x),'b','linewidth',2) 
plot(w,angle(x),'g','linewidth',2.5); 
title('Gráfico função envelope |X(j\omega)| ') 
xlabel('Frequência (\omega)') 
grid on; 
 
 
4 
 
 
 
% GRAICO DA FUNÇÃO 
% sinc(T)= (sin(T))/T 
% com w variando de -pi a pi 
 
clear; 
clc 
 
w=(-pi:0.01:pi) 
T=2 
for j=(1:length(w)); 
x(j)= 2*pi*(rectangularPulse(w(j)/T)); 
end 
hold on 
plot(w,abs(x),'g','linewidth',2.5) 
title('Função sinc(T) ') 
xlabel('Frequência (\omega)') 
grid on 
 
5 
 
 
 
% GREAFICO DA FUMÇÃO 
% x(t)= 1/(a+t) 
% X(W)= 2*pi*exp (-a*w(j)); 
 
clear; 
clc 
 
t= (0:0.01:1) 
a=1 
w= (0:0.01:1) 
for j=(1:length(w)); 
x(j)= 2*pi*exp (-a*w(j)); 
end 
for j=(1:length(t)) 
 Y(j)= 1/(a+t(j)); 
end 
subplot(1,2,1) 
plot(w,abs(x),'r','linewidth',2.5) 
title('|X(j\omega)| ') 
xlabel('Frequência (\omega)') 
grid on 
subplot(1,2,2) 
plot(t,abs(Y),'g','linewidth',2.5) 
title('|x(t)| ') 
xlabel('tempo') 
grid on 
 
6 
 
 
 
 
 
7 
 
% GRFICO DA FUNÇÃO 
% X(W)= 1/(1+y*sen(w)-j*sen(w)) 
% com w variando de -pi a pi 
% com w variando de 0 a 4*pi 
 
clear; 
clc 
 
a=0.5 
w0=(-pi:0.01:pi) 
w1=(0:0.01:4*pi) 
T=2 
y=0.8 
for j=(1:length(w0)); 
 x0(j)= 1/(1-(y*cos(w0(j)))+1i*y*sin(w0(j))); 
end 
for j=(1:length(w1)); 
 x1(j)= 1/(1-(y*cos(w1(j)))+1i*y*sin(w1(j))); 
end 
figure (1) 
subplot(1,2,1) 
plot(w0,abs(x0),'r') 
title('Gráfico do Módulo |X(j\omega)| ') 
xlabel('Frequência (\omega)') 
grid on 
subplot(1,2,2) 
plot(w0,angle(x0),'g') 
title('Gráfico da Fase \angleX(j\omega) ') 
xlabel('Frequência (\omega)') 
grid on 
 figure (2) 
 subplot(1,2,1) 
 plot(w1,abs(x1),'r') 
 title('Gráfico do Módulo |X(j\omega)| ') 
 xlabel('Frequência (\omega)') 
 grid on 
 subplot(1,2,2) 
 plot(w1,angle(x1),'g') 
 title('Gráfico da Fase \angleX(j\omega) ') 
 xlabel('Frequência (\omega)') 
 grid on 
 
 
8 
 
 
 
 
 
% GRFICO DA FUNÇÃO 
% X(W)= (sen(4.5*w))/sen(0.5*w) 
% com w variando de -pi a pi 
 
clear; 
clc 
 
w=(-pi:0.01:pi) 
T=2 
for j=(1:length(w)); 
x(j)=(sin(4.5*w(j)))/(sin(0.5*w(j))); 
end 
hold on 
plot(w,abs(x),'r','linewidth',2.5) 
title('Gráfico do Módulo |X(j\omega)| ') 
xlabel('Frequência (\omega)') 
grid on