%modulazione AM
t=[0:0.001:0.1];
ma=0.8;
xm=cos(2*pi*50*t);
xam=cos(2*pi*200*t)+ma*(xm.*cos(2*pi*200*t));
subplot(221);
plot(t,xam);
title('modulazione di ampiezza');
xlabel('tempo(sec)');pause;

Xs=fft(xam,512);
X=fftshift(Xs);
f=[0:255]*500/256;
subplot(222);
plot(f,abs(X(257:512)));
title('spettro della AM');
xlabel('frequenza(Hz)');
pause;

%modulazione FM
t=[0:0.001:0.1];
mf=0.2;
xm=sin(2*pi*50*t);
xfm=cos(2*pi*200*t+mf*xm);
subplot(223);
plot(t,xfm);
title('modulazione di frequenza');
xlabel('tempo(sec)');pause;

Xs=fft(xfm,512);
X=fftshift(Xs);
f=[0:255]*500/256;
subplot(224);
plot(f,abs(X(257:512)));
title('spettro della FM');
xlabel('frequenza(Hz)');
pause;