Fs=40000;  % Frecuencia de muestreo
t=[0:1/Fs:40000/Fs-1/Fs]';
F0=200;  
y=cos(2*pi*F0*t);  % señal de mensaje 
Fc=10000;   % Frecuencia de la portadora
ka=0.5; % Factor de modulacion
Ac=5;
ymod=Ac*(1+ka*y).*cos(2*pi*Fc*t);  % señal modulada
figure(1)
plot(t,ymod)
xlabel('Tiempo [seg]')
title('Señal modulada AM-DSB')
N=2^14;
Ymod=fft(ymod,N);
F=[-Fs/2:Fs/N:Fs/2-Fs/N]';
figure(2)
plot(F,abs(fftshift(Ymod)))
xlabel('Frecuencia [Hz]')
title('Espectro de la Señal modulada AM-DSB')
%==================================================
% Run Simulink model demodulador_AM.mdl
%==================================================
tt=[0:1/Fs:length(ydemod)/Fs-1/Fs]';
figure(3)
plot(tt,ydemod)
xlabel('Tiempo [seg]')
title('Señal demodulada AM-DSB')
figure(4)
plot(tt,ydemodfilt)
xlabel('Tiempo [seg]')
title('Señal demodulada AM-DSB filtrada')
N=2^14;
Ydemod=fft(ydemod,N);
Ydemodfilt=fft(ydemodfilt,N);
F=[-Fs/2:Fs/N:Fs/2-Fs/N]';
figure(5)
plot(F,abs(fftshift(Ydemod)))
xlabel('Frecuencia [Hz]')
title('Espectro de la Señal demodulada AM-DSB')
figure(6)
plot(F,abs(fftshift(Ydemodfilt)))
xlabel('Frecuencia [Hz]')
title('Espectro de la Señal demodulada AM-DSB filtrada')