Simulación De Una Modulación Digital 8-Psk Y 256-Qam En Matlab 2010A

clc;clear all;close all; bit = randi([0 1],240000,1); % Creando el tren de bits nsamp=4; mod2 = modem.qammod('M', 256, 'PhaseOffset', 0,'InputType','bit'); % Creando el modulador QAM mod3 =modem.pskmod('M', 8, 'PhaseOffset', 0,'InputType','bit'); % Creando el modulador QAM y = modulate(mod2,bit); % Señal modulada , Modulacion 256QAM p = modulate(mod3,bit); % Señal modulada , Modulacion 8PSK%% Filtro N = 40; rolloff = 0.25; % Orden del filtro % Roll-off

filtdef = fdesign.pulseshaping(nsamp,'Square Root Raised Cosine','N,Beta',N,rolloff);% Definición del filtro filtrcl =design(filtdef); % Diseño el filtro

filtrcl.Numerator = filtrcl.Numerator * sqrt(nsamp);% Normalizar el filtro

%% Aplicación del filtro

ymues = upsample(y, nsamp);

% Sobremuestrear la señal

ytrans =filter(filtrcl, ymues);% Aplicar el filtro de raiz de coseno levantado pmues = upsample(p, nsamp); % Sobremuestrear la señal

ptrans = filter(filtrcl, pmues);% Aplicar el filtro de raiz de cosenolevantado

EbN0dB=-6:2:12; t=1; for l=EbN0dB SNR1 = EbN0dB + 10*log10(8) - 10*log10(nsamp); SNR2 = EbN0dB + 10*log10(3) - 10*log10(nsamp); xns = awgn(ytrans,SNR1(t),'measured');% Señal recibida conruido en dB wns = awgn(ptrans,SNR2(t),'measured');% Señal recibida con ruido en dB

%% Recibida yrcv = filter(filtrcl, xns); % Aplicar el filtro de raiz de coseno levantado % Quitar sobremuestreoyrcv = downsample(yrcv, nsamp); ret = N/(nsamp*2);

% Retardo de grupo (cantidad de muestras de entrada)

yrcv = yrcv(2*ret+1:end-2*ret);

% Conteo de retardo

prcv = filter(filtrcl, wns);% Aplicar el filtro de raiz de coseno levantado % Quitar sobremuestreo

prcv = downsample(prcv, nsamp); ret = N/(nsamp*2);

% Retardo de grupo (cantidad de muestras de entrada) % Conteo de retardoprcv = prcv(2*ret+1:end-2*ret);

%% Demodulación y cálculo del BER

h = modem.qamdemod('M', 256 , 'PhaseOffset', 0,'OutputType','bit'); % Creando demodulador 256QAM h2 = modem.pskdemod('M',...