x=2*pi/1025*[0:1024];
f=inline('0.5*sin(112*x)-cos(58*x)');
plot(x,f(x))
c=fft(f(x));
size(c)
z=abs(c);
plot(z)
find(z>1)
z(1)
s=load('signal1'); %% carica il segnale 
%%  http://www.dm.unibo.it/~achilles/calc/Programmi_CS4a/signal1
%%  che dovete prima salvare in una cartella dove Octave puo' trovarlo
plot(s)
plot(abs(fft(s)));
quit;
