Portifólio - Processamento Digital de Sinais

Prévia do material em texto

............................................................................................................................... 
 
 
 
ENGENHARIA ELÉTRICA 
WAGNER VIEIRA ALVES – RA 268522014 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
PORTFÓLIO 
PROCESSAMENTO DIGITAL DE SINAIS 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
............................................................................................................................... 
Guarulhos 
2019 
 
 
 
 
WAGNER VIEIRA ALVES 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
PORTFÓLIO 
PROCESSAMENTO DIGITAL DE SINAIS 
 
 
 
 
 
 
 
 
 
Trabalho apresentado ao Curso Engenharia 
Elétrica da Faculdade ENIAC para a disciplina 
Processamento Digital de Sinais. 
 
Prof. Edilson Alexandre Camargo 
 
 
 
 
 
 
 
 
 
 
 
 
 
Guarulhos 
2019 
 
 
 
............................................................................................................. 
1: Gere uma senóide com frequência de 100 Hz e amplitude 6V, amostrada a 
uma frequência de amostragem de 2 kHz. Adicione a esta senóide ruído 
gaussiano com desvio padrão de 2V. Plote o sinal resultante. 
t=0:0.0005:0.1; 
x=6*sin(2*pi*100*t); 
plot(t,x,'-*'); 
p=0:0.0005:0.1; 
y=2*sin(2*pi*p); 
plot(p,y,'-*'); 
ruido=randn(1,length(p)); 
plot(p,ruido); 
y_ruido=x+ruido; 
plot(p,y_ruido,'-*'); 
plot(t,x,'r-',p,ruido); 
subplot(3,1,1); 
plot(t,x,'-*'); 
title('Seno de 100 Hz'); 
xlabel('Tempo(segundos)'); 
ylabel('Amplitude(volts)'); 
grid 
subplot(3,1,2); 
plot(p,ruido); 
title('Ruido de 2 Volt'); 
xlabel('Tempo(segundos)'); 
ylabel('Amplitude(volts)'); 
grid 
subplot(3,1,3); 
plot(t,x,'r-',p,y_ruido); 
title('Soma do ruido no Seno'); 
xlabel('Tempo(segundos)'); 
ylabel('Amplitude(volts)'); 
grid 
 
 
 
 
2: Usando o comando help square, gere uma onda quadrada de 40 Hz, 
amostrada a uma frequência de 1 KHz. 
 
[u1,t] = gensig('square',tau,Tf,Ts); 
[u2,t] = gensig('pulse',tau,Tf,Ts); 
figure 
subplot(2,1,1) 
plot(t,u1) 
tt = sprintf(['Onda quadrada de frequ^encia %dHz, dura¸c~ao de %ds'... 
'e amostrada a %dkHz'],sigfreq,Tf,sampfreq/1000); 
title(tt) 
subplot(2,1,2) 
plot(t,u2) 
tt = sprintf(['Pulso peri´odico de frequ^encia %dHz, dura¸c~ao de %ds' 
'e amostrada a %dkHz'],sigfreq,Tf,sampfreq/1000); 
grid on 
Tf = 5; 
T = Tf*(1/sigfreq); 
dt = 1/sampfreq; 
t = 0:dt:T-dt; 
x = sawtooth(2*pi*sigfreq*t); 
figure 
plot(t,x) 
tt = sprintf(['Onda dente de serra de frequ^encia %dHz, dura¸c~ao de %ds' 
 'e amostrada a %dkHz'],sigfreq,Tf,sampfreq/1000); 
title(tt) 
grid on 
 
 
 
3: Calcule e plote sua resposta ao degrau s[n] para n = -10,.., 100. E determine 
se o sistema é estável. Usar como base o programa 49 
 
 
n=-10:100; 
x1 = n.*(stepseq(0,-10,100)-stepseq(10,-10,100)); 
x2 = 10*exp(-0.3*(n-10)).*(stepseq(10,-10,100)-stepseq(100,-10,100)); 
x = x1+x2; 
stem(n,x);title('Sequência de Degraus');xlabel('n');ylabel('x[n]'); 
grid 
 
 
 
O sistema é estável pois a amplitude vai até 10V e depois cai mantendo-
se em zero até n=100. Ou seja, o sistema não oscila. 
 4: Considere as sequências: Usar como base o programa 47 
x = [3, 11, 7, 0 , -1, 4, 2], -3 < n < 3 
h = [2, 3 , 0, -5, 2, 1]; -1 < n < 4 
Plote as sequencias acima e a convolução das duas sequências. 
 
x=[3,11,7,0,-1,4,2]; 
h=[2,3,0,-5,2,1]; 
y=conv(x,h); 
subplot(3,1,1); 
stem(x); 
title('Sequência') 
xlabel('X[N]'); 
ylabel('Y[N]'); 
grid 
x=[3,11,7,0,-1,4,2]; 
h=[2,3,0,-5,2,1]; 
p=conv(x,h); 
subplot(3,1,2); 
stem(h); 
title('Sequência') 
xlabel('X[N]'); 
ylabel('Y[N]'); 
grid 
x=[3,11,7,0,-1,4,2]; 
h=[2,3,0,-5,2,1]; 
m=conv(x,h); 
subplot(3,1,3); 
stem(m); 
title('Sequência (Convolução)') 
xlabel('X[N]'); 
ylabel('Y[N]'); 
grid 
 
 
 5: Usando a imagem do palhaço do Matlab, Plote a imagem com diferentes 
características de cor e brilho. 
load clown 
whos 
image(X) 
 
 
Programa: 
load clown 
whos 
image(X) 
colormap(map) 
 
Programa: 
load clown 
whos 
image(X) 
colormap(map) 
brighten(.1) 
 
 
Programa: 
load clown 
whos 
image(X) 
colormap(map) 
brighten(-.6) 
 
 
Programa: 
load clown 
whos 
image(X) 
colormap(map) 
brighten(-.9) 
 
 
Programa: 
load clown 
whos 
image(X) 
colormap(map) 
brighten(.7) 
 
 
 
Conclusão 
Os exercicio foi de grande importância para desemvolver os principios do 
sofware Matlab, entender a plotagen de sinal e o seu comprotamento, sendo 
uma materia muito complexa, tenho a percepção que será necessário se 
aprofundar nos proximo desafios. Acredito que contribuiu imensamente com os 
meus conhecimentos acadêmicos.