clear all, close all, clc

dt = 1e-4; % Tempo di Campionamento

f = 3;

x = -5:dt:5;

SinusoidalSignal = cos(2*pi*f*x);

[HF,F] = FourierTransform(SinusoidalSignal, dt); % Chiamo una funzione

figure,plot(F,abs(HF))

axis([-10 10 0 1.1*max(abs(HF))]),

RectangularSignal = (x >= -.5)&(x <= .5); % Costruisco un rettangolo con operazioni logiche

[HF,F] = FourierTransform(RectangularSignal, dt);

figure,plot(F,abs(HF))

axis([-10 10 0 1.1*max(abs(HF))]),

Moltiplicato = SinusoidalSignal.*RectangularSignal;

[HF,F] = FourierTransform(Moltiplicato, dt);

figure,plot(F,abs(HF))

axis([-10 10 0 1.1*max(abs(HF))]),