python---音频线性叠加

下面是音频的线性叠加代码：

import librosa
import numpy as np
from matplotlib import pyplot as plt

data, fs = librosa.load('C2_2_y.wav',sr=None)
s = np.arange(0,len(data),step=1)
time = s / fs
xmax = np.max(np.abs(data))
data = data / xmax    # 幅值归一化

y = np.random.randn(np.size(data))
ymax = np.max(np.abs(y))
y = y / ymax

z = data + y
zmax = np.max(np.abs(z))
z = z / zmax

plt.figure()

plt.subplot(3,1,1)
plt.plot(time,data)
plt.xlabel('time/s')
plt.ylabel('norm amp')
plt.title('original signal')


plt.subplot(3,1,2)
plt.plot(time,y)
plt.xlabel('time/s')
plt.ylabel('norm amp')
plt.title('random amp')

plt.subplot(3,1,3)
plt.plot(time,z)
plt.xlabel('time/s')
plt.ylabel('norm amp')
plt.title('linear add')

plt.show()

因为看的是赵力老师那本书，所以音频也用的是他本书对应的音频。画出来的图是：

np.size()函数是用来统计变量里面有多少元素，比如矩阵有多少个元素。

np.arange()--函数返回一个有终点和起点的固定步长的排列,其中,step代表跳跃的步长，但是属于左包右不包。比如取[3,9],也就是返回3，4，5，6，7，8.

np.linspace()返回一个有起点和终点的序列，但是它没有步长。它可以指定返回这个范围内多少个元素，有一个num的参数控制。