laba6

Fd = 1000;
t = 0:(1/Fd):(2*pi)/5;
Fn = 100;

%x = 4*cos(2*pi*t)+2*sin(4*pi*t).*sin(6*pi*t);
x = 2*cos(2*pi*t)+0.5*cos(-2*pi*t)-0.5*cos(10*pi*t);

figure(1);
plot(t,x,'k');
hold on;
grid on;
title('Модулируемый сигнал');

Msd = ammod(x,Fn,Fd);
furd = fft(Msd);
furd = abs(furd(1:length(furd)/2+1));
M = max(furd);
furd = abs(furd(1:length(furd)/2+1))/(0.5*M);
frd = [0:(length(furd)-1)]*Fd/length(furd)/2;

figure(2);
plot(frd,furd,'g');
hold on;
grid on;
title ('Спектр сигнала после двухполосной амплитудной модуляции');

Mss = ssbmod(x,Fn,Fd);
furs = fft(Mss);
furs = abs(furs(1:length(furs)/2+1));
M1 = max(furs);
furs = abs(furs(1:length(furs)/2+1))/(0.25*M1);
frs = [0:length(furs)-1]*Fd/length(furs)/2;

figure(3);
plot(frs,furs,'b');
hold on;
grid on;
title('Спектр сигнала после однополосной амплитудной модуляции');