Modulaçao AM, codigo matlab

UNIVERSIDADE FEDERAL DE OURO PRETO
INSTITUTO DE CIÊNCIAS EXATAS E APLICADAS
DEPARTAMENTO DE CIÊNCIAS EXATAS E APLICADAS

Luiz Gustavo Silva

Trabalho de fundamentos de telecomunicação
Modulação AM

João Monlevade
06 de agosto de 2013

Introdução
Esse trabalho tem como objetivo simular o processo de modulação e demodulação de AM no ambiente computacional Matlab. As modulaçoes devem ser realizadas com os respectivos coeficientes u1 = 0.5, u2 = 1 e u3 = 1.5.

Desenvolvimento
Começamos pela programaçao em matlab para modulaçao e demodução:
%% Modulação AM_DSB_SC clc close all clear all f1 = 10; f2 = 20; f3 = 30; f4 = 40; f5 = 50; f6 = 60; f7 = 70; f8 = 80; f9 = 90; fo = 10*f9; u1=0.5; u2 = 1; u3 = 1.5; t = 0: 1/1000 : 1 - 1/1000; x = 10 + 5*cos(2*pi*f1.*t) + 7*cos(2*pi*f2.*t) + 3*cos(2*pi*f3.*t) +
6*cos(2*pi*f4.*t) + 2*cos(2*pi*f5.*t) + 1*cos(2*pi*f6.*t) +
4*cos(2*pi*f7.*t) + 2*cos(2*pi*f8.*t) + 9*cos(2*pi*f9.*t);
X = fft(x);
X_mag = abs(X);
Xmin = min(x); xo = cos(2*pi*fo.*t);
Xo = fft(xo);
Xo_mag = abs(Xo);
Xomin = min(xo); xm = x.*xo; %modulação
Xm_mag = abs(fft(xm)); r = 1000/length(xo); f = 0: r : 999;
%% Demodulação_AM_DSB_SC fc = 90; o=2; fs=1000; %intervalo wn=2*fc/fs; [b,a]=butter(o,wn,'low'); sinal= 2*xm.*xo; %Demodulação sinalF=filtfilt(b,a,sinal); %Recuperação

figure(1); subplot(5,2,1), plot(t,x), title('Sinal Modulante no Dominio do tempo'), xlabel('Tempo(s)'), ylabel('Amplitude'), subplot(5,2,2), plot(t,xo), title('Portadora Dominio do tempo'), xlabel('Tempo s'), ylabel('Amplitude'), subplot(5,2,3), plot(t,xm), title('Sinal Modulado no Dominio do tempo'), xlabel('Tempos'), ylabel('Amplitude'), subplot(5,2,4), plot(t,sinal), title('Sinal Demodulado no Domínio do tempo'), xlabel('Tempo s'), ylabel('Amplitude'), subplot(5,2,5), plot(t,sinalF), title('recuperaçao no Dominio do tempo'), xlabel('Tempos'), ylabel('Amplitude'), figure(2); subplot(5,2,1),