Matematicas

%MODULACION Y DEMODULACION AM
%PROGRAMA PARA REALIZAR LA MODULACION Y DEMODULACION DE UNA SENAL
close all
clear all
fm=input('Introduzca frecuencia de muestreo: ' );
fs=input(' Introduzcafrecuencia de la señal: ' );
fc=input( ' Introduzca frecuencia de la señal portadora: ' );
%MODULACION
t=0:1/fm:5; %definiendo el tiempo real
L=length(t); %sacando la longitud del vector tiempof1=sin(2*pi*fs*t); %senal de informacion
subplot(2,2,1);
plot(t(1:500),f1(1:500), 'r ' );title(’senal de entrada’); grid on,xlabel(’tiempo’),ylabel(’m(t)’)
f2=cos(2*pi*fc*t); %creando cosenoidal paralograr la modulacion
subplot(2,2,2); plot(t(1:100),f2(1:100));title(’señal portadora’),grid on,xlabel(’tiempo’),ylabel(’p(t)’)
f3=f1.*f2; %multiplicacion de las dos funciones
subplot(2,2,3:4);plot(t(1:500),f3(1:500));title(’señal modulada’),grid on,xlabel(’tiempo’),ylabel(’m(t)*p(t)’)
hold on
plot(t(1:500),f1(1:500),’r')
pause
g1=abs(fftshift(fft(f1)))/L; %transformando en frecuencia f1k=[(-L/2)+1:L/2];
k1=(1/(0.795*2*pi))*k;
figure
subplot(3,1,1), plot(k1,g1),title(’espectro de magnitud de la senal de entrada’),grid on,xlabel(’frecuencia’),ylabel(’m(f)’)
g2=abs(fftshift(fft(f2)))/L;%transformando en frecuancia f2
subplot(3,1,2); plot(k1,g2), title(’espectro de magnitud de la senal portadora’),grid on,xlabel(’frecuencia’),ylabel(’p(f)’)
g3=abs(fftshift(fft(f3)))/L;%transformando en frecuencia f3
subplot(3,1,3); plot(k1,g3),title(’espectro de magnitud de la senal modulada’),grid on,xlabel(’frecuencia’),ylabel(’m(f) conv p(f)’)
pause
%DEMODULACION
%asumiendo ceroperdidas en el medio
f4=f3.*f2; %demodulando
figure %figura 3
subplot(2,2,1:2);plot(t(1:500),f4(1:500)),title(’senal demodulada’),grid on,xlabel(’tiempo’),ylabel(’R(t)’)
g4=abs(fftshift(fft(f4)))/L;%transformando en frecuencia la senal demodulada
subplot(2,2,3:4);plot(k1,g4),title(’espectro de magnitud de la senal demodulada’),grid on,xlabel(’frecuencia’),ylabel(’R(f)= m(f) conv p(f)’)
pause...