频率响应复数两种形式互拆完整推导

目录

• 一、第1步：建立带阻尼单自由度时域动力学方程
• 二、第2步：推导频率响应函数 \(H(j\omega)\)
• 三、第3步：复数有理化，分离实部 \(U(\omega)\)、虚部 \(V(\omega)\)
  • 频响极坐标形式（幅值+相位）
  • 实部、虚部物理含义
• 四、第4步：从复数频响还原为时域稳态响应
  • 时域两项分量物理拆解
  • 三个典型频率时域特征
• 五、仿真参数与典型频率数值结果
  • 系统参数
• 六、Python完整可视化代码（频域+时域+Nyquist图）
• 七、全流程逻辑闭环总结

一、第1步：建立带阻尼单自由度时域动力学方程

受简谐激励的单自由度粘性阻尼系统力平衡方程：

\[m\ddot{x}(t)+c\dot{x}(t)+kx(t)=f(t) \]

• \(m\)：集中质量，\(c\)：粘性阻尼系数，\(k\)：结构刚度
• \(x(t)\)：位移输出响应，\(f(t)\)：外部激励力

采用复指数稳态激励法求解线性系统稳态解，设复数形式激励：

\[f(t)=F_0 e^{j\omega t} \]

线性系统稳态响应为同频率简谐振动，设位移复数解：

\[x(t)=X \cdot e^{j\omega t} \]

\(X\) 为复数振幅，同时包含振动幅值、相位偏移两类信息。

二、第2步：推导频率响应函数 \(H(j\omega)\)

对位移求导：

\[\dot{x}(t)=j\omega X e^{j\omega t},\quad \ddot{x}(t)=(j\omega)^2 X e^{j\omega t}=-\omega^2 X e^{j\omega t} \]

代入运动微分方程，约去公共项 \(e^{j\omega t}\)：

\[\left(k-m\omega^2 + jc\omega\right) X = F_0 \]

定义位移/力频率响应函数FRF：输出复数振幅与输入复数振幅的比值

\[\boldsymbol{H(j\omega)=\frac{X}{F_0}=\frac{1}{k-m\omega^2 + jc\omega}} \]

分母物理拆解：

1. 实项 \(k-m\omega^2\)：刚度力、惯性力，属于保守储能力，只在动能、势能间转换能量，不发热损耗；
2. 虚项 \(jc\omega\)：阻尼力，属于耗散力，和速度成正比，始终滞后位移90°，用虚数单位\(j\)表征相位正交关系。

三、第3步：复数有理化，分离实部 \(U(\omega)\)、虚部 \(V(\omega)\)

复数分离实虚部规则：分子分母同乘分母复数共轭
分母：\(D=k-m\omega^2+jc\omega\)，共轭：\(\bar{D}=k-m\omega^2-jc\omega\)

\[H(j\omega)=\frac{1}{k-m\omega^2+jc\omega} =\frac{(k-m\omega^2)-jc\omega}{\left(k-m\omega^2\right)^2+\left(c\omega\right)^2} \]

拆分后标准复数形式 \(H(j\omega)=U(\omega)+jV(\omega)\)：

\[\boxed{ \begin{aligned} \boldsymbol{U(\omega)} &= \frac{k-m\omega^2}{\left(k-m\omega^2\right)^2+\left(c\omega\right)^2} \quad \text{实部：同相储能分量}\\ \boldsymbol{V(\omega)} &= -\frac{c\omega}{\left(k-m\omega^2\right)^2+\left(c\omega\right)^2} \quad \text{虚部系数：90°滞后耗能分量} \end{aligned} } \]

频响极坐标形式（幅值+相位）

1. 幅频特性（位移放大倍数）：

\[A(\omega)=|H(j\omega)|=\sqrt{U^2(\omega)+V^2(\omega)} =\frac{1}{\sqrt{\left(k-m\omega^2\right)^2+\left(c\omega\right)^2}} \]

2. 相频特性（位移相对激励的滞后相位）：

\[\varphi(\omega)=\arctan2\left(V(\omega),U(\omega)\right) =-\arctan\left( \frac{c\omega}{k-m\omega^2} \right) \]

相位取值范围：\(0^\circ \sim -180^\circ\)，位移永远滞后激励，符合物理因果规律。

实部、虚部物理含义

1. 实部 \(U(\omega)\)（储能同相分量）
   由刚度力、惯性力合力决定，位移分量和激励力相位完全相同，能量仅在弹簧弹性势能、质量动能间来回转换，无能量损耗：

• 低频段 \(\omega\ll\omega_n\)：\(U\approx \dfrac{1}{k}\)，弹簧刚度主导，位移与激励几乎同步；
• 共振点 \(\omega=\omega_n\)：\(k-m\omega^2=0\)，\(U=0\)，惯性力与刚度力相互抵消，储能分量完全消失；
• 高频段 \(\omega\gg\omega_n\)：\(U<0\)，惯性力主导，位移与激励接近180°反相。

2. 虚部 \(V(\omega)\)（正交耗能分量）
   完全由阻尼决定，位移分量滞后激励90°，激励输入的能量通过阻尼摩擦转化为热能耗散：

• 全程\(V(\omega)<0\)，表征位移始终滞后于激励；
• 共振点处虚部绝对值达到最大值，输入能量几乎全部用于克服阻尼做功，系统共振振幅最大；
• 低频、高频两端\(V\)趋近于0，阻尼几乎不消耗能量。

四、第4步：从复数频响还原为时域稳态响应

实际物理激励为实数简谐信号，取复激励的实部：

\[f(t)=F_0\cos\omega t=\mathrm{Re}\left[F_0 e^{j\omega t}\right] \]

系统复数位移响应：

\[\tilde{x}(t)=H(j\omega)F_0 e^{j\omega t}=F_0\left[U(\omega)+jV(\omega)\right]\left(\cos\omega t+j\sin\omega t\right) \]

展开后取实部，得到时域位移表达式：

\[\boldsymbol{x(t)=F_0 \cdot \Big[ U(\omega)\cos\omega t + V(\omega)\sin\omega t \Big]} \]

时域两项分量物理拆解

1. \(U(\omega)\cos\omega t\)：同相储能项
   和激励波形完全同步，仅由弹簧、质量的保守力贡献，只吞吐能量不消耗能量；共振时该项直接归零。
2. \(V(\omega)\sin\omega t\)：正交滞后耗能项
   由三角变换 \(\sin\omega t=\cos(\omega t-90^\circ)\)，该分量滞后激励90°，全部来自阻尼损耗；共振时该项幅值最大，是共振大振幅的根本原因。

三个典型频率时域特征

1. 低频 \(\omega\ll\omega_n\)
   \(U\approx \dfrac{1}{k},\ V\approx0\)，\(x(t)\approx \dfrac{F_0}{k}\cos\omega t\)，位移与激励几乎重合，弹簧形变储能，几乎无能量损耗。
2. 共振 \(\omega=\omega_n=\sqrt{\dfrac{k}{m}}\)
   \(U=0\)，仅保留滞后项：\(x(t)=F_0 V(\omega)\sin\omega t\)，位移严格滞后激励90°，激励能量全部被阻尼耗散，振幅达到峰值。
3. 高频 \(\omega\gg\omega_n\)
   \(U<0\)，惯性主导，位移与激励近似反向，振动过快导致质量几乎静止，阻尼耗散减弱，共振振幅快速衰减。

五、仿真参数与典型频率数值结果

系统参数

\(m=1\ \text{kg},\ k=100\ \text{N/m},\ c=2\ \text{Ns/m}\)，固有角频率 \(\omega_n=\sqrt{\dfrac{k}{m}}=10\ \text{rad/s}\)

工况	角频率 \(\omega\)	实部\(U\)	虚部\(V\)	幅值\(A\)	相位\(\varphi\)
低频	\(0.1\ \text{rad/s}\)	\(0.0100\)	\(\approx0\)	\(0.0100\)	\(-0.11^\circ\)（近乎同相）
共振	\(10.0\ \text{rad/s}\)	\(0.0000\)	\(-0.0500\)（负向峰值）	\(0.0500\)（振幅峰值）	\(-90^\circ\)（标准90°滞后）
高频	\(15.0\ \text{rad/s}\)	\(-0.0076\)	\(-0.0018\)	\(0.0078\)	\(-166.50^\circ\)（近乎反相）

六、Python完整可视化代码（频域+时域+Nyquist图）

import numpy as np
import matplotlib.pyplot as plt

# 1. 系统参数
m = 1
k = 100
c = 2
wn = np.sqrt(k / m)

# 2. 频率轴
w = np.linspace(0.01, 20, 2000)
D = (k - m * w ** 2) ** 2 + (c * w) ** 2
U = (k - m * w ** 2) / D
V = -(c * w) / D
A = np.sqrt(U ** 2 + V ** 2)
phi = np.arctan2(V, U) * 180 / np.pi

# 3. 选定三个特征频率点
w_list = [0.1, 10, 15]
t = np.linspace(0, 6, 1000)
F0 = 1

# 4. 绘制频域四张子图：实部、虚部、幅频、相频
plt.figure(figsize=(16, 10))

# 实部曲线
plt.subplot(2, 2, 1)
plt.plot(w, U, label='实部 U(ω) 储能分量')
plt.axvline(wn, color='r', linestyle='--', label=f'固有频率 ωₙ={wn} rad/s')
plt.grid(True)
plt.xlabel('角频率 ω (rad/s)')
plt.ylabel('实部 U')
plt.legend()
plt.title('实部特性曲线')

# 虚部曲线
plt.subplot(2, 2, 2)
plt.plot(w, V, label='虚部 V(ω) 耗能分量', color='orange')
plt.axvline(wn, color='r', linestyle='--')
plt.grid(True)
plt.xlabel('角频率 ω (rad/s)')
plt.ylabel('虚部 V')
plt.legend()
plt.title('虚部特性曲线（共振处负向峰值）')

# 幅频特性
plt.subplot(2, 2, 3)
plt.plot(w, A, label='幅值 A(ω)')
plt.axvline(wn, color='r', linestyle='--')
plt.grid(True)
plt.xlabel('角频率 ω (rad/s)')
plt.ylabel('幅值')
plt.legend()
plt.title('幅频特性曲线（共振幅值峰值）')

# 相频特性
plt.subplot(2, 2, 4)
plt.plot(w, phi, label='相位 φ(°)', color='green')
plt.axvline(wn, color='r', linestyle='--')
plt.axhline(-90, color='k', linestyle=':', label='-90°')
plt.grid(True)
plt.xlabel('角频率 ω (rad/s)')
plt.ylabel('相位 (°)')
plt.legend()
plt.title('相频特性曲线（共振相位严格-90°）')

plt.tight_layout()
plt.show()

# 5. 三个频率点时域波形对比图
plt.figure(figsize=(15, 5))
labels = ['低频 ω=0.1 rad/s', '共振 ω=10 rad/s', '高频 ω=15 rad/s']
for idx, wi in enumerate(w_list):
    # 计算当前频率下的U、V
    di = (k - m * wi ** 2) ** 2 + (c * wi) ** 2
    ui = (k - m * wi ** 2) / di
    vi = -(c * wi) / di
    f_t = F0 * np.cos(wi * t)
    x_t = F0 * (ui * np.cos(wi * t) + vi * np.sin(wi * t))

    plt.subplot(1, 3, idx + 1)
    plt.plot(t, f_t, 'k', label='激励 f(t)')
    plt.plot(t, x_t, 'r', label='位移 x(t)')
    plt.title(labels[idx])
    plt.xlabel('时间 t')
    plt.grid(True)
    plt.legend()

plt.tight_layout()
plt.show()

# 6. Nyquist复平面图（实部-虚部轨迹）
plt.figure(figsize=(6, 6))
plt.plot(U, V, color='darkred', linewidth=2, label='Nyquist轨迹')
plt.scatter(1/k, 0, c='blue', label='低频起点')
res_idx = np.argmin(np.abs(w - wn))
plt.scatter(U[res_idx], V[res_idx], c='orange', s=60, label='共振点')
plt.grid(True)
plt.xlabel('实部 U')
plt.ylabel('虚部 V')
plt.axis('equal')
plt.legend()
plt.title('单自由度系统Nyquist半圆（下半平面）')
plt.show()

七、全流程逻辑闭环总结

1. 时域力平衡微分方程 → 复指数稳态假设，推导出有理形式频响函数 \(H(j\omega)=\dfrac{1}{k-m\omega^2+jc\omega}\)；
2. 通过复数共轭有理化，拆分得到实部（刚度+惯性·同相储能）、虚部（阻尼·90°正交耗能）；
3. 利用欧拉公式展开复数激励并取实部，将稳态位移拆解为「同相振动分量+滞后90°耗能分量」，实现频域参数与时域波形物理意义一一对应；
4. 低频由实部刚度主导、共振由虚部阻尼主导、高频由实部惯性主导，Nyquist下半圆轨迹可用于模态参数识别求解阻尼比。