Strumenti musicali

1
2
3
4
5
data=loadtxt('sunspots.txt')
time=data[:,0]
nspots=data[:,1]

plot(time,nspots)

1
2
3
4
5
6
7
tau=time[1]-time[0]
T=time[-1]
N=len(time)
freq=arange(N)/T
Transf = my_dft(nspots,tau)

plot(freq,abs(Transf)**2,'-d')

1
2
3
4
5
r=2*pi*time/T
han=0.5-0.5*cos(r)
Transf_w=my_dft(nspots*han)

plot(freq,abs(Hw)**2,'-d')

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
samples_piano=loadtxt('piano.txt')
samples_trumpet=loadtxt('trumpet.txt')

N=len(samples_piano)
tau=1./44100
T=N*tau
time=linspace(0,T,N)

plot(time,samples_piano)
plot(time,samples_trumpet)

1
2
3
4
5
6
7
Transf_piano=my_dft(samples_piano,tau)
Transf_trumpet=my_dft(samples_trumpet,tau)
freq=arange(N)/T

plot(freq,abs(Transf_piano)**2)
xlim(0,10000/T) # imposto il range da visualizzare lungo x
ylim(0,3E+16)   # e lungo y