谐波饱和方程推导：以连续积分替代离散频率叠加的技术问询

嘿，你的思路完全没问题！把离散叠加扩展到连续频率范围的积分确实能得到更简洁的解析结果，我来一步步帮你推导清楚：

首先，我们先把离散的叠加推广到连续情况：假设我们考虑以F为中心，宽度为2ΔF的连续频率区间（也就是频率从F-ΔF到F+ΔF），每个频率分量的权重相同，那么归一化后的振幅可以用积分表示为：
$$A(t) = \frac{1}{2\Delta F} \int_{F-\Delta F}^{F+\Delta F} \sin(2\pi f t) df$$
这里除以2ΔF是为了和离散情况的“除以波的数量”对应，做归一化处理，保证中心频率F处的振幅不会随区间宽度变化。

接下来我们求解这个积分，利用基本的三角函数积分公式：$\int \sin(ax) dx = -\frac{\cos(ax)}{a} + C$，这里令a=2πt，代入上下限计算：
$$
\begin{align*}
\int_{F-\Delta F}^{F+\Delta F} \sin(2\pi f t) df &= \left. -\frac{\cos(2\pi f t)}{2\pi t} \right|_{F-\Delta F}^{F+\Delta F} \
&= \frac{ -\cos(2\pi(F+\Delta F)t) + \cos(2\pi(F-\Delta F)t) }{2\pi t}
\end{align*}
$$

然后用三角恒等式**$\cos A - \cos B = -2\sin\left(\frac{A+B}{2}\right)\sin\left(\frac{A-B}{2}\right)$**来简化分子：
令$A=2\pi(F-\Delta F)t$，$B=2\pi(F+\Delta F)t$，代入后可得：
$$\cos A - \cos B = 2\sin(2\pi F t)\sin(2\pi \Delta F t)$$

把这个结果代回积分式，再除以归一化因子2ΔF，最终得到：
$$A(t) = \sin(2\pi F t) \cdot \frac{\sin(2\pi \Delta F t)}{2\pi \Delta F t}$$

如果熟悉信号处理里的sinc函数（定义为$\text{sinc}(x) = \frac{\sin(\pi x)}{\pi x}$），我们还可以把式子写成更简洁的形式：
$$A(t) = \sin(2\pi F t) \cdot \text{sinc}(2\Delta F t)$$

结果解释 #

• 这个式子和你离散叠加的思路本质一致：离散叠加是用三个离散频率点近似连续积分，而连续积分得到的是被sinc函数调制的正弦波，sinc函数的形状决定了振幅的“饱和”特性——中心频率F处的振幅最大（t=0时sinc值为1），随着时间t增加，振幅会以sinc函数的规律衰减、振荡，最终趋近于0。
• 你可以调整ΔF的大小：ΔF越大，sinc函数的主瓣越窄，振幅衰减得越快；ΔF越小，主瓣越宽，振幅变化越平缓，和离散情况中邻近频率的间隔对应。

内容的提问来源于stack exchange，提问作者J Guerrero