PONTIFÍCIA UNIVERSIDADE CATÓLICA DE MINAS GERAIS
Por: Mariele Souza • 20/3/2019 • Trabalho acadêmico • 395 Palavras (2 Páginas) • 125 Visualizações
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
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
GRÁFICOS
[pic 1]
Gráfico 4PSK – 4QAM
[pic 2]
Gráfico 8-PSK e 16-QAM
[pic 3]
Gráfico 16-PSK 64-QAM
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.
...