PSD-EngTel (página)

De MediaWiki do Campus São José
Revisão de 10h45min de 24 de fevereiro de 2015 por Moecke (discussão | contribs) (→‎Atividades extra)
Ir para navegação Ir para pesquisar

MURAL DE AVISOS E OPORTUNIDADES DA ÁREA DE TELECOMUNICAÇÕES


Registro on-line das aulas

Aula 1 (6 fev)
  • Revisão de Sinais e Sistemas no tempo discreto em Matlab:
Variação do Experimento 1.2
%% Experimento 1.2
fs = 10; Ts = 1/fs; fase = 0;
time = 0:Ts:(1-Ts);
f1 = 3; f2 = 7;
s_1 = cos(2*pi*f1*time+fase);
s_2 = cos(2*pi*f2*time+fase);
fsa = 1000; Tsa = 1/fsa;
time_aux = 0:Tsa:(1-Tsa);
figure(1);
stem(time,s_1,'ob');
hold on;
plot(time_aux, cos(2*pi*f1*time_aux+fase),'--k');
stem(time,s_2,'+r');
plot(time_aux, cos(2*pi*f2*time_aux+fase),'--m');
hold off;
legend('s_1 discreto','s_1 contínuo','s_2 discreto','s_2 contínuo')

Varie o valor da frequência de amostragem de 6 até 20 Hz e observe:

  1. Em qual frequência deixa de ocorrer recobrimento do sinal 2.
  2. O que ocorre quando a frequência é 6, 7, 14 Hz? Explique
  3. Qual deveria ser a frequência do sinal f_2 para que as amostras tomadas sejam coincidentes como o sinal f_1 para uma frequência de amostragem f_s? Reescreva a equação e verifique no Matlab.
Aula 2 (10 fev)
  • Revisão de Sinais e Sistemas no tempo discreto em Matlab:
Variação do Experimento 2.2
%% Experimento 2.2
% Resposta em frequencia usando a função freqz
N = 1;
num = [1 0 0 0];
den = poly([0.8 0.2])
%den = [1 0.6 -0.16];
% modo 1
%[H,w]=freqz(num,den,[0:pi/100:N*pi-pi/100]);
%plot(w/pi, abs(H));
% modo 2
%[H,w]=freqz(num,den);
%plot(w/pi, abs(H));
% modo 3
%[H,w]=freqz(num, den, 'whole');
%plot(w/pi, abs(H));
% modo 4
freqz(num, den, 'whole');
figure(2);
zplane(num,den);

%% Resposta em frequencia substituindo z -> e^(jw)
syms z
Hf(z) = symfun(z^2/(z-0.2)/(z+0.8),z);
pretty(Hf)
latex(Hf)
N = 1;
w = [0:pi/100:N*pi-pi/100];
plot(w/pi,abs(Hf(exp(1i*w))))
%title(['$' latex(Hf) '$'],'interpreter','latex')
text(0.2,2,['H(z) = ' '$$' latex(Hf) '$$'],'interpreter','latex')
xlabel(['w/' '$$' '\pi' '$$'],'interpreter','latex')
  1. Verifique a diferença entre os tipos de plots comentados no código.
  2. substitua o denominador de H(z) por dois polos em [-0.8 -0.8].
  3. verifique o que ocorre se forem utilizados polos complexos conjugados [0.3-0.4i 0.3+0.4i 0.1]
  4. verifique o que ocorre se forem utilizados polos complexos não conjugados [0.3-0.4i 0.3+0.8i]
  5. verifique o que ocorre se os polos estiverem fora do circulo unitário [1.2 -0.2]. Interprete este resultado
Aula 3 (13 fev)
  • Revisão de Sinais e Sistemas no tempo discreto em Matlab:
  • Filtros Digitais (Experimento 2.3)
  • Filtragem de Sinais (Experimento 3.1)
  • Análise de Sinais (Experimento 3.2)
  • Uso das funções zp2tf, tf2zp, fft, ifft, fftfilt, fftshift
  • Ver pag. 141 a 145 e 230 a 235 de [1]
Aula 4 (20 fev)
  • Filtros Digitais: Filtros IIR:
  • Aproximação de magnitude de filtros analógicos: Butterworth.
  • Ver pag. 186 a 193 de [2]

Atividades extra

Experimento 3.2 (Prazo de entrega 23/02/2015)
  1. Escreva um script com funções do Matlab para gerar as representações no domínio do tempo sinais da Figura 3.26 [1]. Use o subplot para inserir todas os sinais em uma única Figure do Matlab.
  2. Escreva um script com funções do Matlab para gerar as representações no domínio da frequência sinais da Figura 3.27 [1]. Use o subplot para inserir todas os sinais em uma única Figure do Matlab.
  3. Faça uma estimação melhor do espectro do sinal, segmentando o sinal x(n) e calculando a média dos espectros obtidos.
  4. É importante utilizar as escalas e legendas corretamente nos gráficos.
  5. Ver as funções do matlab legend, title, xlabel, ylabel.
  6. Escreva um pequeno relatório técnico mostrando os resultados obtidos e documentando o código escrito. Enviar para o meu email em formato pdf.
ESTUDE - Resposta em frequência a partir da posição dos pólos-zeros
  1. Estudar a Seção 5.6 pag. 485 a 488 do livro do Lathi [3].
  2. Use as funções zplane, pz2tf, freqz para realizar o estudo.
ESTUDE - Resposta em frequência em Sistemas de Tempo Discreto
  1. Estudar a Seção 5.5 pag. 474 a 475 do livro do Lathi [3].
ESTUDE - Desempenho do código de uma FFT
  1. Estudar o Exemplo M9.2 pag. 778 a 779 do livro do Lathi [3].
N = 10000;
%% tamanho do vetor de 1015 -> 5 x 7 x 29
N0 = 1015;
y1 = rand(N0,1);
tic; for k = 1:N; Y1 = fft(y1)/N0; end; toc

%% tamanho do vetor de 1019 -> primo
N_0 = 1019;
y2 = [y1; zeros(N_0-N0,1)];
tic; for k = 1:N; Y2 = fft(y2)/N_0; end; toc

%% tamanho do vetor de 1024 -> 2^10
N_0 = 1024; 
y3 = [y1; zeros(N_0-N0,1)];
tic; for k = 1:N; Y3 = fft(y3)/1024; end; toc
Projeto de Filtro Butterworth (Prazo de entrega 02/03/2015)
  • Projeto os seguintes filtros analógicos passa baixas do tipo Butterworth que atenda as seguintes especificações:
a) Frequência de corte (passband) f_s = XX Hz. Frequência de rejeição (stopband) f_p = XX Hz. Ganho de XX dB em w = 0, atenuação de XX dB em f_s e atenuação mínima de XX dB em f_p.
b) Frequência de corte (passband) f_s = XX Hz. Frequência de rejeição (stopband) f_p = XX Hz. Ganho de XX dB em w = 0, atenuação de XX dB em f_s e atenuação mínima de XX dB em f_p.
  • O projeto de cada filtro deve apresentar o cálculo da ordem do filtro, dos polos do filtro, a equação de H(p), H(s), |H(jw)|^2.
  • A resposta em frequência normalizada e real devem ser feitas usando o Matlab. Apresentar o gráfico do ganho em dB e da fase em cada caso com escalas corretas.
  • Deve ser apresentado o diagrama dos pólos do filtro H(p) e H(s).
  • Para fins de comparação dos filtros, o filtro 1 e o filtro 2 devem ser sobrepostos nesses gráficos.
  • Escreva um pequeno relatório técnico mostrando os resultados obtidos e documentando o código escrito. Enviar para o meu email em formato pdf com o Assunto: PSD29007 - Projeto de Filtro Butterworth.

Recursos necessários

  • O Software Matlab está disponível na maioria dos laboratórios do IFSC-campus São José em instalação local tanto em Windows como Linux. Adicionalmente de qualquer maquina do IFSC, pode-se fazer Acesso ao IFSC-CLOUD.
  • Adicionalmente de qualquer maquina do IFSC, pode-se fazer usar o IFSC-CLOUD para ter acesso a estes softwares.

Links auxiliares

Referências Bibliográficas

  1. 1,0 1,1 1,2 1,3 1,4 DINIZ, P. S. R., DA SILVA, E. A. B., e LIMA NETTO, S. Processamento Digital de Sinais: Projeto e Análise de Sistemas. 2. ed. Porto Alegre: Bookman, 2014. 976 p. ISBN 978-8582601235
  2. SHENOI, B. A. Introduction to Digital Signal Processing and Filter Design. 1.ed. New Jersey: John Wiley-Interscience, 2006. 440 p. ISBN 978-0471464822
  3. 3,0 3,1 3,2 LATHI, Bhagwandas P. Sinais e Sistemas Lineares. 2. ed. Porto Alegre: Artmed-Bookman, 2007. 856 p. ISBN 978-8560031139


Curso de Engenharia de Telecomunicações