Buscar

Comunicação Digital - atividade computacional 4

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Você viu 3, do total de 5 páginas

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Faça como milhares de estudantes: teste grátis o Passei Direto

Esse e outros conteúdos desbloqueados

16 milhões de materiais de várias disciplinas

Impressão de materiais

Agora você pode testar o

Passei Direto grátis

Prévia do material em texto

PONTIFÍCIA UNIVERSIDADE CATÓLICA DE MINAS GERAIS 
Graduação em Engenharia Elétrica 
 
 
 
Mariele Cristina de Lima Souza 
 
 
 
 
 
 
 
ATIVIDADE COMPUTACIONAL 4 
 
 
 
 
 
 
 
 
 
Poços de Caldas 
2017 
1. CÓDIGO FONTE 
 
%Modulação PSK e QAM 
close all;clear all;clc; 
l = 2640; %input ('Quantidade de bits transmitido:'); 
m_psk = input('Ordem de modulação PSK:'); 
m_qam = input('Ordem de modulação QAM:'); 
n_psk = log2(m_psk); 
n_qam = log2(m_qam); 
snr = 20; %input('Relação Sinal ruido:'); 
 
for i=0:1:snr %Para criarmos um grafico de desempenho 
 
 %Fonte 
 BB_psk = (randi(m_psk,1,l))-1; 
 BB_qam = (randi(m_qam,1,l))-1; 
 
 %Modulação 
 y_psk = modulate(modem.pskmod(m_psk),BB_psk); 
 y_qam = modulate(modem.qammod(m_qam),BB_qam); 
 %%scatterplot(y) 
 
 
 %Pulse Sahping 
 yp_psk = rectpulse(y_psk,n_psk); 
 yp_qam = rectpulse(y_qam,n_qam); 
 
 
 %Canal 
 ruido_psk = awgn(yp_psk,i,'measured'); %Normaliza o ruído 
 ruido_qam = awgn(yp_qam,i,'measured'); %Normaliza o ruído 
 %ynoisy = yp+awgn 
 %%scatterplot(ruido) 
 
 %Filtro Int/Dump 
 yf_psk = intdump(ruido_psk,n_psk); 
 yf_qam = intdump(ruido_qam,n_qam); 
 
 %Demodulador 
 z_psk = demodulate(modem.pskdemod(m_psk),yf_psk); 
 z_qam = demodulate(modem.qamdemod(m_qam),yf_qam); 
 
 %Contador de Erros 
 [Erro_psk(i+1), BER_psk(i+1)] = symerr(BB_psk,z_psk); 
 [Erro_qam(i+1), BER_qam(i+1)] = symerr(BB_qam,z_qam); 
 vetorx(i+1)=i; 
 
 a = 10^(i/10); %Anti-log 
 ser_psk (i+1)= 2*qfunc(sqrt(2*a*log2(m_psk))*sin(pi/m_psk)); 
 ser_qam(i+1) = 4*(1-1/sqrt(m_qam))*qfunc(sqrt(3*a*log2(m_qam)/(m_qam-1))); 
end 
 semilogy(vetorx,ser_psk,'r') % Grafico da Equação 
 hold on 
 semilogy(vetorx,ser_qam) % Grafico da Equação 
 xlabel SNR 
 ylabel SER 
 hold on 
 grid on 
 
2. GRÁFICOS 
 
 
Gráfico 4PSK – 4QAM 
 
Gráfico 8-PSK e 16-QAM 
 
 
Gráfico 16-PSK 64-QAM 
 
3. CONCLUSÃO 
 
A relação sinal-ruído compara o nível de um sinal desejado com o nível do ruído de 
fundo. Quanto mais alta for a relação sinal-ruído, menor é o efeito do ruído de fundo 
sobre a SER. 
Em relação ao número de bits por símbolo (M) podemos dizer que é diretamente 
proporcional a SER, tendo em vista que quanto maior for o valor de M maior será SER. 
Como podemos analisar no gráfico para 4-PSK temos um valor final na casa de 10-40 já 
para o 16-PSK esse valor final passa para a casa de 10-6.

Outros materiais