原文地址：https://dylanmeeus.github.io/posts/audio-from-scratch-pt11/
到目前为止，我们已经看到了如何生成正弦波，方波和三角波等纯信号。这些信号易于生成方便调试，但在现实世界中，仪器不会生成这种纯信号。例如，拔吉他弦时，吉他弦会沿着多个频率振动。这些不同的振动称为“谐波”，由“基本”频率+泛音频率组成。因此，我们听到的最终声音是沿这些不同频率振动的波的组合。
一种便于你想象理解的方式是用光进行类比。当你通过棱镜照射光时，它将分解为组成光的不同颜色（波长）。同样，你可以想象一个棱镜，通过该棱镜我们可以发送声波，将其分解为构成该特定声音的波长。
既然实际仪器组合了不同频率的波，那么问题就变成了，我们如何以编程方式生成这些波？
傅立叶加法
我们将通过应用傅立叶加法产生这些波。就像在上一篇帖子中一样，我们将生成的数据存储在一个表格中，我们可以用LookupOscillator来振荡。最后，我们将使用我们之前编写的函数normalize对查询表的幅度进行归一化。
func FourierTable(nharms int, amps []float64, length int, phase float64) []float64 {
table := make([]float64, length+2)
phase *= tau
for i := 0; i < nharms; i++ {
for n := 0; n < len(table); n++ {
amp := 1.0
if i < len(amps) {
amp = amps[i]
}
angle := float64(i+1) * (float64(n) * tau / float64(length))
table[n] += (amp * math.Cos(angle+phase))
}
}
return normalize(table)
}
在此函数中，我们通过nharms确定要生成的谐波次数，并通过 []float64设置不同的谐波幅度（如果未传，则默认为1.0）。此外，我们将为表格指定一个长度，以及一个开始阶段，该阶段实质上是波的偏移量。
基本思想是，对于每个谐波，我们都会遍历整个表格，并按照table[n]谐波的当前波幅增加波。随后的每个谐波都以N * fundamental频率振荡， 因此我们得到（harmonics = {1F, 2F ... NF}），其基波为“ 1F”。
我们可以根据生成的谐波来改变通过的角度math.Cos(..)来获得该值。
// Adjust angle for harmonic 'i' -> move N times forward in the wave
angle := float64(i+1) * (float64(n) * tau / float64(length))
// Add wave amplitude for harmonic 'i' to existing wave (table[n])
table[n] += (amp * math.Cos(angle+phase))
产生谐波
现在我们已经有了傅里叶加法，我们仍然必须使用此发生器来生成不同的波形。我们可以通过将一个float64切片传递给 FourierTable函数来更改谐波的外观，这对于生成我们想要的波谐波类型至关重要。
例如，要生成具有N个谐波的方波，三角波和锯齿波表：
func SquareTable(nharms, length int) []float64 {
amps := make([]float64, nharms)
for i := 0; i < len(amps); i += 2 {
amps[i] = 1.0 / float64(i+1)
}
return FourierTable(nharms, amps, length, -0.25)
}
func SawTable(nharms, length int) []float64 {
amps := make([]float64, nharms)
for i := 0; i < len(amps); i++ {
amps[i] = 1.0 / float64(i+1)
}
return FourierTable(nharms, amps, length, -0.25)
}
func TriangleTable(nharms, length int) []float64 {
amps := make([]float64, nharms)
for i := 0; i < nharms; i += 2 {
amps[i] = 1.0 / (float64(i+1) * float64(i+1))
}
return FourierTable(nharms, amps, length, 0)
}
要查看其外观和听起来的效果，我们可以用一个小的测试程序来生成此类波形。这样的程序包含在GoAudio示例中
要运行此程序以生成频率为440时长为4秒的方波，该方波具有6个谐波且最大幅度为0.8，我们将使用以下命令：
go run main.go -d 4 -s square -a 0.8 -h 6 -f 440 -o output.wav
输出效果如下
现在，我们具有开始使用GoAudio进行一些简单调音的基础知识，这也将是下一篇文章的主题。
