scipy 小结

1from scipy.io import wavfile
2# read wav file
3filename = 'input.wav'
4rate,data = wavfile.read(filename)
5# write wav file
6filename = 'output.wav'
7wavfile.write(filename)

 1def sa(x):
 2    return np.sin(x)/x
 3def si(x):
 4    s = []
 5    for i in x:
 6        #防止被积函数分母出现0
 7        if i == 0:
 8            i = 1e-6
 9        s.append(quad(sa,1e-6,i)[0])
10        #return y,abserr
11    return s
12#plot sa函数与si函数的图像
13x = np.linspace(-10,10,100)
14f1 = sa(x)
15f2 = si(x)
16plt.plot(x,f1,label='si(x)')
17plt.plot(x,f2,label='sa(x)')
18plt.xlabel('t(s)')
19plt.ylabel('function')
20plt.legend()
21plt.show()

 1def pulse1(x):
 2    if x<=0:
 3        result = 0
 4    elif x>=1:
 5        result = 0
 6    else :
 7        result = x
 8    return result
 9t = np.linspace(0,2.5,50)
10#求卷积积分
11ulist = []
12for i in t:
13    f = lambda x:pulse1(x)*pulse1(i-x)
14    # quad产生0到i的等分矩阵X，设间距为t，求出此时f(X)的值，结果为t*f(X)的和
15    ulist.append(quad(f,0,i)[0])