核心想法：
简谐振动 \(\cos\) 不是“算出来的”，而是一个匀速圆周运动在直径上的投影。
1. 几何构造
画一个圆，半径 \(A\)，圆心在原点。
一个质点以角速度 \(\omega\) 在圆上匀速转动，它的位置角为：
\[ \theta(t)=\omega t+\varphi \]
质点在圆上的坐标是：
\[ X=A\cos(\omega t+\varphi) \]
\[ Y=A\sin(\omega t+\varphi) \]
现在只看 \(X\) 坐标，也就是质点在 \(x\) 轴上的投影。
这个投影就是：
\[ \boxed{x(t)=A\cos(\omega t+\varphi)} \]
所以，简谐振动就是匀速圆周运动在一条直径上的投影。
2. 为什么这个投影满足 \(F=-kx\)
圆周运动的质点有向心加速度：
\[ \vec a_c=-\omega^2\vec r \quad \text{是矢量，方向指向圆心} \]
方向指向圆心。
把这个加速度投影到 \(x\) 轴上：
\[ a_x=-a_c\cos\theta=-\omega^2 A\cos\theta \]
但：
\[ x=A\cos\theta \]
所以：
\[ a_x=-\omega^2 x \]
两边乘 \(m\)：
\[ F_x=ma_x=-m\omega^2 x \]
令：
\[ k=m\omega^2 \]
就得到：
\[ \boxed{F_x=-kx} \]
所以，匀速圆周运动的投影，自动满足弹簧振子的方程。
3. 几何图像
你可以这样想：
- 圆上质点匀速转；
- 它在 \(x\) 轴上的影子来回移动；
- 影子在最右端时，速度为零，加速度指向圆心，即指向左；
- 影子在最左端时，速度为零，加速度指向右；
- 影子在中间时，速度最大，加速度为零。
这正好就是弹簧振子的运动：
- 端点：速度为零，加速度最大；
- 平衡位置：速度最大，加速度为零。
4. 用几何直接读出 \(x,v,a\)
设：
\[ \theta=\omega t+\varphi \]
圆上质点位置：
\[ \vec r=(A\cos\theta,\ A\sin\theta) \]
速度：
\[ \vec v=(-A\omega\sin\theta,\ A\omega\cos\theta) \]
加速度：
\[ \vec a=(-A\omega^2\cos\theta,\ -A\omega^2\sin\theta) \]
只看 \(x\) 分量：
\[ x=A\cos\theta \]
\[ v_x=-A\omega\sin\theta \]
\[ a_x=-A\omega^2\cos\theta=-\omega^2 x \]
所以：
\[ \boxed{a_x=-\omega^2 x} \]
这就是 \(F=-kx\) 的几何来源。
5. 为什么是 \(\cos\) 而不是别的
因为：
- 圆周运动是匀速的，角度随时间线性增加：\(\theta=\omega t+\varphi\)；
- 投影到 \(x\) 轴，自然出现 \(\cos\theta\)；
- 所以 \(x(t)=A\cos(\omega t+\varphi)\)。
不是先假设 \(\cos\)，而是圆周运动的投影天然就是 \(\cos\)。
6. 一句话总结
\[ \boxed{\text{简谐振动 }x=A\cos(\omega t+\varphi)\text{ 就是匀速圆周运动在直径上的投影。}} \]
几何上：
- 圆半径 \(A\) 是振幅；
- 角速度 \(\omega\) 是角频率；
- 初始角 \(\varphi\) 是初相位；
- 投影位置就是位移 \(x\)；
- 投影加速度满足 \(a=-\omega^2 x\)，即 \(F=-kx\)。
所以，不用解微分方程，画一个圆就能看出为什么解是 \(\cos\) 形式。
严格来说，欠阻尼受迫振动的稳态解并不衰减；衰减的是瞬态解。
总解 = 瞬态衰减振动 + 稳态等幅振动。
不过，如果你的 \(x(t)\) 是一个衰减振动形式，例如：
\[ x(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
那么求外力 \(F(t)\) 的方法仍然是：
\[ \boxed{F(t)=m\ddot{x}(t)+c\dot{x}(t)+kx(t)} \]
下面直接对一般衰减振动形式求一遍。
1. 设已知位移为
\[ x(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
其中：
\[ \beta=\frac{c}{2m},\qquad \omega_d=\sqrt{\omega_0^2-\beta^2},\qquad \omega_0=\sqrt{\frac{k}{m}} \]
2. 求一阶导数
\[ \dot{x}(t)
A\left[-\beta e^{-\beta t}\cos(\omega_d t+\varphi) -e^{-\beta t}\omega_d\sin(\omega_d t+\varphi)\right] \]
即：
\[ \dot{x}(t)
-Ae^{-\beta t}\left[\beta\cos(\omega_d t+\varphi) +\omega_d\sin(\omega_d t+\varphi)\right] \]
3. 求二阶导数
对 \(\dot{x}\) 再求导：
\[ \ddot{x}(t)
Ae^{-\beta t} \left[ (\beta^2-\omega_d^2)\cos(\omega_d t+\varphi) +2\beta\omega_d\sin(\omega_d t+\varphi) \right] \]
4. 代入受迫振动方程
受迫振动方程为：
\[ m\ddot{x}+c\dot{x}+kx=F(t) \]
所以：
\[ F(t)=m\ddot{x}+c\dot{x}+kx \]
代入 \(x,\dot{x},\ddot{x}\)。
注意：
\[ c=2m\beta \]
\[ k=m\omega_0^2 \]
并且：
\[ \omega_0^2=\omega_d^2+\beta^2 \]
先算 \(m\ddot{x}\)：
\[ m\ddot{x}
mAe^{-\beta t} \left[ (\beta^2-\omega_d^2)\cos(\omega_d t+\varphi) +2\beta\omega_d\sin(\omega_d t+\varphi) \right] \]
再算 \(c\dot{x}\)：
\[ c\dot{x}
2m\beta\cdot \left( -Ae^{-\beta t} \left[ \beta\cos(\omega_d t+\varphi) +\omega_d\sin(\omega_d t+\varphi) \right] \right) \]
即：
\[ c\dot{x}
-2m\beta Ae^{-\beta t} \left[ \beta\cos(\omega_d t+\varphi) +\omega_d\sin(\omega_d t+\varphi) \right] \]
再算 \(kx\)：
\[ kx
m\omega_0^2 Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
5. 合并余弦项和正弦项
把 \(F(t)\) 写成：
\[ F(t)=mAe^{-\beta t}\left[P\cos(\omega_d t+\varphi)+Q\sin(\omega_d t+\varphi)\right] \]
其中余弦项系数：
\[ P=(\beta^2-\omega_d^2)-2\beta^2+\omega_0^2 \]
因为：
\[ \omega_0^2=\omega_d^2+\beta^2 \]
所以：
\[ P=\beta^2-\omega_d^2-2\beta^2+\omega_d^2+\beta^2=0 \]
正弦项系数：
\[ Q=2\beta\omega_d-2\beta\omega_d=0 \]
所以：
\[ \boxed{F(t)=0} \]
6. 结论
如果
\[ x(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
并且满足：
\[ \beta=\frac{c}{2m},\qquad \omega_d=\sqrt{\frac{k}{m}-\beta^2} \]
那么：
\[ \boxed{F(t)=0} \]
也就是说：
这种衰减振动是欠阻尼自由振动，不是受迫振动。它不需要外力，靠初始条件维持。
7. 如果你想表达的是“受迫振动总解中的瞬态部分”
受迫振动总解一般为：
\[ x(t)=x_{\text{瞬态}}(t)+x_{\text{稳态}}(t) \]
其中瞬态部分：
\[ x_{\text{瞬态}}(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
它对应的是齐次方程的解，所以代入齐次方程：
\[ m\ddot{x}+c\dot{x}+kx=0 \]
因此对这部分而言，外力为零。
而稳态部分才对应外力：
\[ x_{\text{稳态}}(t)=B\cos(\omega t-\delta) \]
此时：
\[ F(t)=F_0\cos(\omega t) \]
8. 如果你已知的 \(x(t)\) 不是标准欠阻尼自由振动
那就要老老实实代入：
\[ \boxed{F(t)=m\ddot{x}(t)+c\dot{x}(t)+kx(t)} \]
你可以把具体的 \(x(t)\) 发给我，我直接帮你算出 \(F(t)\)。
1. 方程
你定义的力：
\[ F=F_{\text{驱}}-c\dot{x} \]
按前面的约定，运动方程是：
\[ m\ddot{x}+kx=F \]
如果：
\[ F=C \]
则：
\[ m\ddot{x}+kx=C \]
这是一个无阻尼、受常力的弹簧振子方程。
2. 求解
方程：
\[ m\ddot{x}+kx=C \]
先找平衡位置。令 \(\ddot{x}=0\)：
\[ kx=C \]
所以平衡位置移到：
\[ x_0=\frac{C}{k} \]
令：
\[ y=x-\frac{C}{k} \]
则：
\[ \ddot{y}=\ddot{x} \]
代入：
\[ m\ddot{y}+k\left(y+\frac{C}{k}\right)=C \]
\[ m\ddot{y}+ky+C=C \]
\[ m\ddot{y}+ky=0 \]
所以 \(y\) 是简谐振动：
\[ y(t)=A\cos(\omega_0 t+\varphi) \]
其中：
\[ \omega_0=\sqrt{\frac{k}{m}} \]
因此：
\[ \boxed{ x(t)=\frac{C}{k}+A\cos(\omega_0 t+\varphi) } \]
3. 物理解释
- 常数 \(C\) 的作用只是把平衡位置从 \(0\) 移到 \(C/k\)。
- 振动频率不变，仍是： \[ \omega_0=\sqrt{\frac{k}{m}} \]
- 振幅 \(A\) 和初相位 \(\varphi\) 由初始条件决定。
4. 如果原来的系统有阻尼
原来的受迫振动方程是：
\[ m\ddot{x}+c\dot{x}+kx=F_{\text{驱}} \]
你定义的力是：
\[ F=F_{\text{驱}}-c\dot{x} \]
如果 \(F=C\)，则：
\[ F_{\text{驱}}-c\dot{x}=C \]
也就是：
\[ F_{\text{驱}}=C+c\dot{x} \]
此时驱动力不再是常数，而是：
\[ F_{\text{驱}}(t)=C+c\dot{x}(t) \]
也就是说，要让“你定义的那个力”保持常数 \(C\)，驱动力必须随速度变化，去补偿阻尼项。
5. 如果只考虑无阻尼情形
如果你直接把 \(F\) 当成弹簧振子受到的一个常数力，那么：
\[ m\ddot{x}+kx=C \]
解就是：
\[ \boxed{ x(t)=\frac{C}{k}+A\cos\left(\sqrt{\frac{k}{m}},t+\varphi\right) } \]
即：
常数力只改变平衡位置，不改变振动频率和振动形式。
6. 对比：如果 \(F\) 是简谐力
若：
\[ F=F_0\cos(\omega t) \]
则方程：
\[ m\ddot{x}+kx=F_0\cos(\omega t) \]
解为：
\[ x(t)=A\cos(\omega_0 t+\varphi)+\frac{F_0}{k-m\omega^2}\cos(\omega t) \]
当：
\[ \omega=\omega_0=\sqrt{\frac{k}{m}} \]
发生共振，振幅趋于无穷（无阻尼理想情形）。
7. 总结
| \(F\) 的形式 | 方程 | 解的特征 |
|---|---|---|
| \(F=0\) | \(m\ddot{x}+kx=0\) | 简谐振动，平衡在 \(x=0\) |
| \(F=C\) | \(m\ddot{x}+kx=C\) | 简谐振动，平衡移到 \(x=C/k\) |
| \(F=F_0\cos(\omega t)\) | \(m\ddot{x}+kx=F_0\cos(\omega t)\) | 受迫振动，\(\omega=\omega_0\) 时共振 |
所以，常数力 \(C\) 不会改变振动频率，只会把平衡位置平移 \(C/k\)。
\[ F=F_{\text{驱}}-c\dot{x} \]
运动方程是：
\[ m\ddot{x}+kx=F \]
现在取：
\[ F=ax,\qquad a>k \]
代入：
\[ m\ddot{x}+kx=ax \]
移项：
\[ m\ddot{x}+kx-ax=0 \]
\[ m\ddot{x}+(k-a)x=0 \]
1. 方程类型
因为 \(a>k\)，所以：
\[ k-a<0 \]
令：
\[ k-a=-\lambda^2,\qquad \lambda>0 \]
即：
\[ \lambda^2=a-k \]
方程变成：
\[ m\ddot{x}-\lambda^2 x=0 \]
也就是：
\[ \boxed{\ddot{x}-\frac{\lambda^2}{m}x=0} \]
这是指数增长/衰减型方程，不是简谐振动。
2. 通解
方程：
\[ \ddot{x}=\frac{\lambda^2}{m}x \]
令：
\[ \gamma=\sqrt{\frac{\lambda^2}{m}}=\sqrt{\frac{a-k}{m}} \]
则：
\[ \ddot{x}=\gamma^2 x \]
通解为：
\[ \boxed{ x(t)=Ae^{\gamma t}+Be^{-\gamma t} } \]
或者写成双曲函数：
\[ \boxed{ x(t)=C\cosh(\gamma t)+D\sinh(\gamma t) } \]
其中：
\[ \gamma=\sqrt{\frac{a-k}{m}}>0 \]
3. 物理意义
- 当 \(a>k\) 时，力 \(F=ax\) 的方向与位移同向。
- 位移越大，力越大，继续把振子往同一方向推。
- 所以不会来回振动，而是偏离平衡位置越来越快。
- 其中 \(e^{\gamma t}\) 项表示指数发散，\(e^{-\gamma t}\) 项表示指数衰减。
- 实际运动由初始条件决定；一般情况下会包含发散项。
4. 若初始条件为
\[ x(0)=x_0,\qquad \dot{x}(0)=v_0 \]
则：
\[ A+B=x_0 \]
\[ \gamma(A-B)=v_0 \]
解得：
\[ A=\frac{1}{2}\left(x_0+\frac{v_0}{\gamma}\right) \]
\[ B=\frac{1}{2}\left(x_0-\frac{v_0}{\gamma}\right) \]
所以：
\[ \boxed{ x(t)
\frac{1}{2}\left(x_0+\frac{v_0}{\gamma}\right)e^{\gamma t} + \frac{1}{2}\left(x_0-\frac{v_0}{\gamma}\right)e^{-\gamma t} } \]
其中：
\[ \gamma=\sqrt{\frac{a-k}{m}} \]
5. 若 \(a<k\)
则：
\[ k-a>0 \]
方程变为：
\[ m\ddot{x}+(k-a)x=0 \]
即：
\[ \ddot{x}+\frac{k-a}{m}x=0 \]
这是简谐振动：
\[ \boxed{ x(t)=A\cos(\omega t+\varphi) } \]
其中：
\[ \omega=\sqrt{\frac{k-a}{m}} \]
6. 若 \(a=k\)
则：
\[ m\ddot{x}=0 \]
\[ \ddot{x}=0 \]
解为匀速直线运动：
\[ \boxed{ x(t)=x_0+v_0t } \]
7. 总结
| 条件 | 方程 | 运动形式 |
|---|---|---|
| \(a<k\) | \(m\ddot{x}+(k-a)x=0\) | 简谐振动 |
| \(a=k\) | \(m\ddot{x}=0\) | 匀速直线运动 |
| \(a>k\) | \(m\ddot{x}+(k-a)x=0\)，即 \(m\ddot{x}-\lambda^2x=0\) | 指数增长/衰减，不发散就是衰减 |
所以，当 \(a>k\) 时，系统不再振动，而是出现指数增长或指数衰减的解。
\[ F=-kx-cv \]
且
\[ F=ma \]
得到：
\[ ma=-kx-cv \]
1. 先写成微分方程
因为
\[ v=\dot{x},\qquad a=\ddot{x} \]
所以：
\[ m\ddot{x}=-kx-c\dot{x} \]
移项：
\[ m\ddot{x}+c\dot{x}+kx=0 \]
即：
\[ \boxed{\ddot{x}+\frac{c}{m}\dot{x}+\frac{k}{m}x=0} \]
2. 求 \(v(t)\)
令：
\[ \beta=\frac{c}{2m},\qquad \omega_0=\sqrt{\frac{k}{m}} \]
则方程写为：
\[ \ddot{x}+2\beta\dot{x}+\omega_0^2x=0 \]
欠阻尼情形：
\[ \beta<\omega_0 \]
令：
\[ \omega_d=\sqrt{\omega_0^2-\beta^2} \]
位移解为：
\[ x(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
对 \(x(t)\) 求导，得速度：
\[ v(t)=\dot{x}(t) \]
\[ v(t)
-Ae^{-\beta t} \left[ \beta\cos(\omega_d t+\varphi) +\omega_d\sin(\omega_d t+\varphi) \right] \]
所以：
\[ \boxed{ v(t)
-Ae^{-\beta t} \left[ \beta\cos(\omega_d t+\varphi) +\omega_d\sin(\omega_d t+\varphi) \right] } \]
3. 求 \(a(t)\)
再对 \(v(t)\) 求导，得加速度：
\[ a(t)=\dot{v}(t)=\ddot{x}(t) \]
\[ a(t)
Ae^{-\beta t} \left[ (\beta^2-\omega_d^2)\cos(\omega_d t+\varphi) +2\beta\omega_d\sin(\omega_d t+\varphi) \right] \]
所以：
\[ \boxed{ a(t)
Ae^{-\beta t} \left[ (\beta^2-\omega_d^2)\cos(\omega_d t+\varphi) +2\beta\omega_d\sin(\omega_d t+\varphi) \right] } \]
4. 用 \(\omega_0\) 改写加速度
因为：
\[ \omega_d^2=\omega_0^2-\beta^2 \]
所以：
\[ \beta^2-\omega_d^2
\beta^2-(\omega_0^2-\beta^2)
2\beta^2-\omega_0^2 \]
因此：
\[ \boxed{ a(t)
Ae^{-\beta t} \left[ (2\beta^2-\omega_0^2)\cos(\omega_d t+\varphi) +2\beta\omega_d\sin(\omega_d t+\varphi) \right] } \]
5. 验证
由方程：
\[ m\ddot{x}+c\dot{x}+kx=0 \]
即：
\[ ma+cv+kx=0 \]
所以：
\[ a=-\frac{c}{m}v-\frac{k}{m}x \]
代入 \(v(t)\) 和 \(x(t)\)，也可以得到同样的 \(a(t)\)。
6. 总结
欠阻尼情形下：
\[ \boxed{ x(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) } \]
\[ \boxed{ v(t)
-Ae^{-\beta t} \left[ \beta\cos(\omega_d t+\varphi) +\omega_d\sin(\omega_d t+\varphi) \right] } \]
\[ \boxed{ a(t)
Ae^{-\beta t} \left[ (\beta^2-\omega_d^2)\cos(\omega_d t+\varphi) +2\beta\omega_d\sin(\omega_d t+\varphi) \right] } \]
其中：
\[ \beta=\frac{c}{2m},\qquad \omega_0=\sqrt{\frac{k}{m}},\qquad \omega_d=\sqrt{\omega_0^2-\beta^2} \]
\[ v(t)
-Ae^{-\beta t} \left[ \beta\cos(\omega_d t+\varphi) +\omega_d\sin(\omega_d t+\varphi) \right] \]
下面分几种“平均值”来说。
1. 对一个周期取平均（短时平均）
如果只看稳态振荡的一个周期，或者忽略衰减，取：
\[ T_d=\frac{2\pi}{\omega_d} \]
则在一个周期内：
\[ \langle \cos(\omega_d t+\varphi)\rangle=0 \]
\[ \langle \sin(\omega_d t+\varphi)\rangle=0 \]
所以：
\[ \boxed{\langle v(t)\rangle_{T_d}=0} \]
这是因为正弦、余弦在一个完整周期内的平均值为零。
2. 从 \(t=0\) 到 \(t\to\infty\) 的时间平均
\[ \langle v\rangle
\lim_{T\to\infty}\frac{1}{T}\int_0^T v(t),dt \]
由于 \(e^{-\beta t}\to0\)，且 \(\cos,\sin\) 有界，所以：
\[ \int_0^\infty v(t),dt \]
是有限值，而分母 \(T\to\infty\)，因此：
\[ \boxed{\langle v\rangle=0} \]
3. 均方根速度（有效值）
如果你要的是“平均速度大小”或“有效速度”，通常用均方根：
\[ v_{\text{rms}}
\sqrt{\langle v^2\rangle} \]
忽略衰减，只取一个周期平均：
\[ v(t)\approx -A\left[ \beta\cos\theta+\omega_d\sin\theta \right] \]
其中：
\[ \theta=\omega_d t+\varphi \]
则：
\[ v^2
A^2 \left[ \beta^2\cos^2\theta +\omega_d^2\sin^2\theta +2\beta\omega_d\sin\theta\cos\theta \right] \]
一个周期平均：
\[ \langle \cos^2\theta\rangle=\frac12 \]
\[ \langle \sin^2\theta\rangle=\frac12 \]
\[ \langle \sin\theta\cos\theta\rangle=0 \]
所以：
\[ \langle v^2\rangle
A^2\left( \frac{\beta^2}{2} +\frac{\omega_d^2}{2} \right)
\frac{A^2}{2}(\beta^2+\omega_d^2) \]
但：
\[ \beta^2+\omega_d^2=\omega_0^2=\frac{k}{m} \]
因此：
\[ \langle v^2\rangle
\frac{A^2\omega_0^2}{2} \]
所以：
\[ \boxed{ v_{\text{rms}}
\frac{A\omega_0}{\sqrt2}
A\sqrt{\frac{k}{2m}} } \]
4. 如果考虑完整衰减过程
考虑完整时间积分：
\[ \langle v^2\rangle
\frac{\int_0^\infty v^2(t),dt}{\int_0^\infty dt} \]
分母发散，所以不能这样定义。
通常改为：
\[ \int_0^\infty v^2(t),dt \]
这个积分是有限的，可以用来表示总能量耗散。
计算可得：
\[ \int_0^\infty v^2(t),dt
\frac{A^2}{4\beta}(\beta^2+\omega_d^2)
\frac{A^2\omega_0^2}{4\beta} \]
其中：
\[ \beta=\frac{c}{2m},\qquad \omega_0^2=\frac{k}{m} \]
所以：
\[ \boxed{ \int_0^\infty v^2(t),dt
\frac{A^2 k}{2c} } \]
5. 总结
| 平均方式 | 结果 |
|---|---|
| 一个周期平均 \(\langle v\rangle\) | \(0\) |
| 长时间平均 \(\langle v\rangle\) | \(0\) |
| 均方根速度 \(v_{\text{rms}}\) | \(\displaystyle \frac{A\omega_0}{\sqrt2}=A\sqrt{\frac{k}{2m}}\) |
| 速度平方时间积分 \(\int_0^\infty v^2dt\) | \(\displaystyle \frac{A^2 k}{2c}\) |
所以，速度的平均值通常是零，因为正负半周相互抵消；如果要看“大小”，就用均方根速度。
\[ x(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
来算。
1. 能量定义
弹簧振子的总机械能：
\[ E(t)=\frac12 m v^2+\frac12 k x^2 \]
其中：
\[ \beta=\frac{c}{2m},\qquad \omega_0^2=\frac{k}{m},\qquad \omega_d^2=\omega_0^2-\beta^2 \]
2. 位移和速度
\[ x(t)=Ae^{-\beta t}\cos\theta \]
\[ v(t)=-Ae^{-\beta t}\left[\beta\cos\theta+\omega_d\sin\theta\right] \]
其中：
\[ \theta=\omega_d t+\varphi \]
3. 动能
\[ K(t)=\frac12 m v^2 \]
\[ K(t)=\frac12 mA^2e^{-2\beta t} \left[\beta\cos\theta+\omega_d\sin\theta\right]^2 \]
展开：
\[ K(t)=\frac12 mA^2e^{-2\beta t} \left[ \beta^2\cos^2\theta +\omega_d^2\sin^2\theta +2\beta\omega_d\sin\theta\cos\theta \right] \]
4. 势能
\[ U(t)=\frac12 k x^2 \]
\[ U(t)=\frac12 kA^2e^{-2\beta t}\cos^2\theta \]
因为：
\[ k=m\omega_0^2=m(\omega_d^2+\beta^2) \]
所以：
\[ U(t)=\frac12 mA^2e^{-2\beta t}(\omega_d^2+\beta^2)\cos^2\theta \]
5. 总能量
\[ E(t)=K(t)+U(t) \]
把动能和势能相加：
余弦项：
\[ \frac12 mA^2e^{-2\beta t} \left[ \beta^2\cos^2\theta +(\omega_d^2+\beta^2)\cos^2\theta \right]
\frac12 mA^2e^{-2\beta t} \left[ \omega_d^2\cos^2\theta +2\beta^2\cos^2\theta \right] \]
正弦项：
\[ \frac12 mA^2e^{-2\beta t}\omega_d^2\sin^2\theta \]
交叉项：
\[ mA^2e^{-2\beta t}\beta\omega_d\sin\theta\cos\theta \]
合并后可以整理为：
\[ E(t)=\frac12 mA^2e^{-2\beta t} \left[ \omega_d^2 +2\beta^2\cos^2\theta +\beta\omega_d\sin2\theta \right] \]
6. 对一个周期取平均
在一个周期内：
\[ \langle \cos^2\theta\rangle=\frac12 \]
\[ \langle \sin2\theta\rangle=0 \]
所以：
\[ \langle E(t)\rangle
\frac12 mA^2e^{-2\beta t} \left[ \omega_d^2 +2\beta^2\cdot\frac12 \right] \]
\[ \langle E(t)\rangle
\frac12 mA^2e^{-2\beta t} (\omega_d^2+\beta^2) \]
但：
\[ \omega_d^2+\beta^2=\omega_0^2=\frac{k}{m} \]
所以：
\[ \boxed{ \langle E(t)\rangle
\frac12 kA^2e^{-2\beta t} } \]
也可以写成：
\[ \boxed{ \langle E(t)\rangle
E_0 e^{-2\beta t} } \]
其中：
\[ E_0=\frac12 kA^2 \]
是初始时刻一个周期内的平均能量。
7. 物理解释
-
能量按指数衰减： \[ E(t)\propto e^{-2\beta t}=e^{-\frac{c}{m}t} \]
-
阻尼越大，能量衰减越快。
-
衰减时间常数： \[ \tau=\frac{1}{2\beta}=\frac{m}{c} \]
-
品质因数： \[ Q=\frac{\omega_0}{2\beta}=\frac{\sqrt{km}}{c} \]
8. 如果考虑完整时间平均
若从 \(t=0\) 到 \(t\to\infty\) 做时间平均，分母发散，不能这样定义。
但可以求总耗散能量：
\[ \int_0^\infty \langle E(t)\rangle,dt
\frac12 kA^2\int_0^\infty e^{-2\beta t}dt
\frac{kA^2}{4\beta} \]
即：
\[ \boxed{ \int_0^\infty \langle E(t)\rangle,dt
\frac{kA^2 m}{2c} } \]
9. 总结
| 量 | 结果 |
|---|---|
| 瞬时总能量 | \(E(t)=\frac12 mv^2+\frac12 kx^2\) |
| 一个周期平均能量 | \(\langle E(t)\rangle=\frac12 kA^2e^{-2\beta t}\) |
| 初始平均能量 | \(E_0=\frac12 kA^2\) |
| 衰减规律 | \(E(t)\propto e^{-\frac{c}{m}t}\) |
| 时间常数 | \(\tau=\frac{m}{c}\) |
| 品质因数 | \(Q=\frac{\sqrt{km}}{c}\) |
所以，能量的期望值就是一个周期平均后的能量，它按 \(e^{-2\beta t}\) 衰减。
\[ x(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
的包络线是：
\[ \pm Ae^{-\beta t} \]
所以：
\[ \boxed{-Ae^{-\beta t}\le x(t)\le Ae^{-\beta t}} \]
1. 为什么不是 \(-A\) 到 \(A\)
因为多了一个衰减因子：
\[ e^{-\beta t} \]
它随时间减小：
\[ e^{-\beta t}\to 0 \]
所以 \(x(t)\) 的振幅也在减小。
-
\(t=0\) 时： \[ -A\le x(0)\le A \]
-
\(t>0\) 时： \[ -Ae^{-\beta t}<x(t)<Ae^{-\beta t} \]
-
\(t\to\infty\) 时： \[ x(t)\to 0 \]
2. 几何图像
\(x(t)\) 被夹在两条衰减曲线之间：
\[ x_{\max}=Ae^{-\beta t} \]
\[ x_{\min}=-Ae^{-\beta t} \]
这两条线就是包络线。
振子在两条包络线之间来回振荡，但每次摆动的幅度都更小。
3. 如果只看一个周期内的局部范围
在某个时刻 \(t_0\) 附近，一个周期内 \(x\) 的局部最大值约为：
\[ Ae^{-\beta t_0} \]
局部最小值约为：
\[ -Ae^{-\beta t_0} \]
所以局部范围也随时间缩小。
4. 总结
| 情况 | 范围 |
|---|---|
| 无阻尼简谐振动 | \(-A\le x(t)\le A\) |
| 欠阻尼自由振动 | \(-Ae^{-\beta t}\le x(t)\le Ae^{-\beta t}\) |
| 长时间后 | \(x(t)\to 0\) |
所以，欠阻尼时 \(x(t)\) 的范围不再是固定的 \(-A\) 到 \(A\)，而是逐渐缩小到零。
1. 基本设定
设信号为 \(x(t)\)，它的傅里叶变换为：
\[ \tilde{x}(\omega)=\int_{-\infty}^{\infty}x(t)e^{-i\omega t}dt \]
逆变换：
\[ x(t)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\tilde{x}(\omega)e^{i\omega t}d\omega \]
定义归一化功率谱：
\[ P(\omega)=\frac{|\tilde{x}(\omega)|^2}{\int_{-\infty}^{\infty}|\tilde{x}(\omega)|^2d\omega} \]
时间概率密度：
\[ \rho(t)=\frac{|x(t)|^2}{\int_{-\infty}^{\infty}|x(t)|^2dt} \]
2. 时间矩
\[ \langle t\rangle=\int_{-\infty}^{\infty}t\rho(t)dt \]
\[ \langle t^2\rangle=\int_{-\infty}^{\infty}t^2\rho(t)dt \]
\[ \Delta t^2=\langle t^2\rangle-\langle t\rangle^2 \]
3. 频率矩
\[ \langle \omega\rangle=\int_{-\infty}^{\infty}\omega P(\omega)d\omega \]
\[ \langle \omega^2\rangle=\int_{-\infty}^{\infty}\omega^2P(\omega)d\omega \]
\[ \Delta\omega^2=\langle \omega^2\rangle-\langle \omega\rangle^2 \]
4. 傅里叶变换的微分性质
关键公式：
\[ \mathcal{F}\left\{\frac{dx}{dt}\right\}=i\omega\tilde{x}(\omega) \]
所以：
\[ \int_{-\infty}^{\infty}\omega^2|\tilde{x}(\omega)|^2d\omega
\int_{-\infty}^{\infty}\left|\frac{dx}{dt}\right|^2dt \]
这是因为帕塞瓦尔定理：
\[ \int_{-\infty}^{\infty}|\tilde{x}(\omega)|^2d\omega
2\pi\int_{-\infty}^{\infty}|x(t)|^2dt \]
以及：
\[ \int_{-\infty}^{\infty}\omega^2|\tilde{x}(\omega)|^2d\omega
2\pi\int_{-\infty}^{\infty}\left|\frac{dx}{dt}\right|^2dt \]
所以：
\[ \boxed{ \langle \omega^2\rangle
\frac{\int_{-\infty}^{\infty}\left|\frac{dx}{dt}\right|^2dt} {\int_{-\infty}^{\infty}|x(t)|^2dt} } \]
这就是通过傅里叶变换求 \(\langle \omega^2\rangle\) 的核心公式。
5. 用你的 \(x(t)\) 算
你假设：
\[ x(t)=Ae^{-\beta t}e^{-i\omega_d t},\qquad t\ge0 \]
则：
\[ |x(t)|^2=A^2e^{-2\beta t} \]
时间矩：
\[ \langle t\rangle=\frac{1}{2\beta} \]
\[ \langle t^2\rangle=\frac{1}{2\beta^2} \]
\[ \Delta t^2=\frac{1}{4\beta^2} \]
频率矩：
\[ \frac{dx}{dt}
A(-\beta-i\omega_d)e^{-(\beta+i\omega_d)t} \]
\[ \left|\frac{dx}{dt}\right|^2
A^2(\beta^2+\omega_d^2)e^{-2\beta t} \]
所以：
\[ \int_0^\infty\left|\frac{dx}{dt}\right|^2dt
A^2(\beta^2+\omega_d^2)\cdot\frac{1}{2\beta} \]
\[ \int_0^\infty|x(t)|^2dt
\frac{A^2}{2\beta} \]
因此：
\[ \langle \omega^2\rangle
\frac{A^2(\beta^2+\omega_d^2)\cdot\frac{1}{2\beta}} {\frac{A^2}{2\beta}}
\beta^2+\omega_d^2 \]
即：
\[ \boxed{ \langle \omega^2\rangle=\beta^2+\omega_d^2=\omega_0^2 } \]
6. 得到不确定性乘积
\[ \Delta t^2=\frac{1}{4\beta^2} \]
\[ \Delta\omega^2=\langle \omega^2\rangle-\langle \omega\rangle^2 \]
对纯复指数，\(\langle \omega\rangle=\omega_d\)，所以：
\[ \Delta\omega^2
(\beta^2+\omega_d^2)-\omega_d^2
\beta^2 \]
因此：
\[ \boxed{ \Delta\omega^2\Delta t^2
\beta^2\cdot\frac{1}{4\beta^2}
\frac14 } \]
即：
\[ \boxed{ \Delta\omega,\Delta t=\frac12 } \]
7. 结论
通过傅里叶变换求 \(\langle \omega^2\rangle\) 的公式是：
\[ \boxed{ \langle \omega^2\rangle
\frac{\int\left|\frac{dx}{dt}\right|^2dt} {\int|x(t)|^2dt} } \]
对你假设的：
\[ x(t)=Ae^{-\beta t}e^{-i\omega_d t} \]
得到：
\[ \boxed{ \langle \omega^2\rangle=\beta^2+\omega_d^2=\omega_0^2 } \]
\[ \boxed{ \Delta\omega^2=\beta^2 } \]
\[ \boxed{ \Delta t^2=\frac{1}{4\beta^2} } \]
\[ \boxed{ \Delta\omega^2\Delta t^2=\frac14 } \]
刚好达到不确定性原理的下限。
1. 如果你用的是复数信号
\[ x(t)=Ae^{-\beta t}e^{-i\omega_d t},\qquad t\ge0 \]
它的傅里叶变换：
\[ \tilde{x}(\omega)
A\int_0^\infty e^{-\beta t}e^{-i\omega_d t}e^{-i\omega t}dt \]
\[
A\int_0^\infty e^{-[\beta+i(\omega+\omega_d)]t}dt \]
\[
\frac{A}{\beta+i(\omega+\omega_d)} \]
所以功率谱：
\[ |\tilde{x}(\omega)|^2
\frac{A^2}{\beta^2+(\omega+\omega_d)^2} \]
这是一个洛伦兹线型，中心在：
\[ \omega=-\omega_d \]
如果你定义频率变量为 \(\omega\)，那么：
\[ \langle \omega\rangle
\frac{\int_{-\infty}^{\infty}\omega|\tilde{x}(\omega)|^2d\omega} {\int_{-\infty}^{\infty}|\tilde{x}(\omega)|^2d\omega} \]
由于谱关于 \(\omega=-\omega_d\) 对称，积分得到：
\[ \langle \omega\rangle=-\omega_d \]
而不是 \(+\omega_d\)。
2. 如果你用的是实信号
\[ x(t)=Ae^{-\beta t}\cos(\omega_d t),\qquad t\ge0 \]
它的频谱在：
\[ \omega=+\omega_d \]
和：
\[ \omega=-\omega_d \]
各有一个峰。
如果只取正频率部分，那么：
\[ \langle \omega\rangle=\omega_d \]
如果对全频率轴积分，两个峰对称，平均值为：
\[ \langle \omega\rangle=0 \]
3. 如果你用的是“解析信号”约定
有时把：
\[ x(t)=Ae^{-\beta t}e^{i\omega_d t} \]
作为解析信号，那么谱峰在：
\[ \omega=+\omega_d \]
此时：
\[ \langle \omega\rangle=\omega_d \]
4. 结论
\[ \langle \omega\rangle=\omega_d \]
是否成立，取决于你的约定：
| 约定 | \(\langle \omega\rangle\) |
|---|---|
| 复信号 \(e^{-i\omega_d t}\)，全频率积分 | \(-\omega_d\) |
| 复信号 \(e^{i\omega_d t}\)，全频率积分 | \(+\omega_d\) |
| 实信号 \(\cos(\omega_d t)\)，全频率积分 | \(0\) |
| 实信号，只取正频率 | \(\omega_d\) |
| 解析信号 \(e^{i\omega_d t}\) | \(\omega_d\) |
所以在你前面用的：
\[ x(t)=Ae^{-\beta t}e^{-i\omega_d t} \]
这个约定下，全频率积分得到的是：
\[ \boxed{\langle \omega\rangle=-\omega_d} \]
如果你只关心大小或只取正频率，才写成：
\[ \langle \omega\rangle=\omega_d \]
因此需要确认你用的是哪种频率约定。
1. 之前单边指数
\[ x(t)=Ae^{-\beta t},\qquad t\ge 0 \]
傅里叶变换：
\[ \tilde{x}(\omega)=\frac{A}{\beta+i\omega} \]
功率谱：
\[ |\tilde{x}(\omega)|^2=\frac{A^2}{\beta^2+\omega^2} \]
如果从 \(-\infty\) 到 \(+\infty\) 积分：
\[ \int_{-\infty}^{\infty}\omega^2\frac{A^2}{\beta^2+\omega^2}d\omega \]
被积函数当 \(|\omega|\to\infty\) 时趋于常数 \(A^2\)，所以积分发散：
\[ \Delta\omega^2=\infty \]
2. 双边指数
\[ x(t)=Ae^{-\beta|t|},\qquad -\infty<t<\infty \]
傅里叶变换：
\[ \tilde{x}(\omega)=\frac{2A\beta}{\beta^2+\omega^2} \]
功率谱：
\[ |\tilde{x}(\omega)|^2=\frac{4A^2\beta^2}{(\beta^2+\omega^2)^2} \]
从 \(-\infty\) 到 \(+\infty\) 积分：
\[ \int_{-\infty}^{\infty}\omega^2\frac{4A^2\beta^2}{(\beta^2+\omega^2)^2}d\omega \]
被积函数当 \(|\omega|\to\infty\) 时趋于：
\[ \frac{4A^2\beta^2}{\omega^2} \]
所以积分收敛，结果是有限的。
3. 关键区别
| 情况 | 积分范围 | 功率谱形状 | \(\Delta\omega^2\) |
|---|---|---|---|
| 单边指数 | \(0\) 到 \(\infty\) | 洛伦兹 \(1/(\beta^2+\omega^2)\) | \(\infty\) |
| 双边指数 | \(-\infty\) 到 \(\infty\) | 洛伦兹平方 \(1/(\beta^2+\omega^2)^2\) | \(\beta^2\) |
所以：
- 单边指数：频谱是纯洛伦兹，二阶矩发散；
- 双边指数：频谱是洛伦兹的平方，二阶矩收敛。
4. 结论
你说的对：
前面算单边指数时，积分范围取 \(0\) 到 \(\infty\)，但频谱是定义在 \(-\infty\) 到 \(+\infty\) 上的，所以真正对频率积分时应该从 \(-\infty\) 到 \(+\infty\)，结果发散。
而双边指数因为时间域对称，频谱衰减更快，二阶矩才收敛。
所以：
\[ \boxed{ \text{单边指数：}\Delta t^2\Delta\omega^2=\infty } \]
\[ \boxed{ \text{双边指数：}\Delta t^2\Delta\omega^2=\frac12 } \]
\[ \boxed{ \text{高斯：}\Delta t^2\Delta\omega^2=\frac14 } \]
只有高斯达到最小不确定度。
\[ m\ddot{x}+c\dot{x}+kx=0 \]
这是关于时间 \(t\) 的常微分方程，只有一个自变量 \(t\)，所以严格来说不能做空间分离变量，也没有 \(X(x)\)。
但如果你一定要写成“分离变量”的形式，只能设：
\[ x(t)=e^{rt} \]
代入：
\[ mr^2+cr+k=0 \]
解得：
\[ r=-\beta\pm i\omega_d \]
其中：
\[ \beta=\frac{c}{2m} \]
\[ \omega_d=\sqrt{\frac{k}{m}-\frac{c^2}{4m^2}} \]
于是：
\[ x(t)=Ae^{-\beta t}\cos(\omega_d t+\varphi) \]
如果你把阻尼方程中的空间变量 \(x\) 当成自变量，写成：
\[ m\frac{d^2X}{dx^2}+c\frac{dX}{dx}+kX=0 \]
那它和原方程数学同形，可以解出：
\[ X(x)=e^{-\frac{c}{2m}x} \left[ A\cos\left(\sqrt{\frac{k}{m}-\frac{c^2}{4m^2}},x\right) + B\sin\left(\sqrt{\frac{k}{m}-\frac{c^2}{4m^2}},x\right) \right] \]
即：
\[ \boxed{ X(x)=Ae^{-\beta x}\cos(\omega_d x+\varphi) } \]
其中：
\[ \beta=\frac{c}{2m} \]
\[ \omega_d=\sqrt{\frac{k}{m}-\frac{c^2}{4m^2}} \]
但这只是数学形式，物理上阻尼振动方程里 \(x\) 是位移，不是自变量。