Baixe o app para aproveitar ainda mais
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.
Compartilhar