$\%$ system is limited by it's half of a sampling rate.

fmax $=$ ; $\%$ max frequency in the system simulation

nFFT $=2500$; $\%$ number of points in frequency domain

Xf $=$ fftshift$($fft$(,$nFFT$))$; $\%$ modulating signal spectrum

magXf $=$ abs$($Xf$)$; $\%$ absolute value of modulating spectrum

freq $=-$fmax:$($Fs$/$nFFT$)$:fmax$-($fmax$/$nFFT$)$; $\%$ double sided freq vector

Cf $=$ fftshift$($fft$(,$ nFFT$))$; $\%$ carrier signal spectrum

magCf $=$ abs$($Cf$)$; $\%$ absolute value of carrier signal spectrum

Yf $=$ fftshift$($fft$(,$ nFFT$))$; $\%$ modulated output signal spectrum

magYf $=$ abs$($Yf$)$; $\%$ absolute value of modulated output signal spectrum

figure;

plot$($freq$,20*$log$10($magXf$))$; hold on;grid on;

title$(\text{'}$Modulating signal magnitude spectrum'$)$;

xlabel$(\text{'}$frequency$\text{'})$; ylabel$(\text{'}$amplitude$\text{'})$;

plot$($freq$,20*$log$10($magCf$))$;

title$(\text{'}$Carrier signal magnitude spectrum'$)$;

xlabel$(\text{'}$frequency$\text{'})$; ylabel$(\text{'}$amplitude$\text{'})$;

plot$($freq$,20*$log$10($magYf$))$;

title$(\text{'}$Modulated output magnitude spectrum'$)$;

xlabel$(\text{'}$frequency$\text{'})$; ylabel$(\text{'}$amplitude$\text{'})$;