无线电自学教程 03:电路基础——从欧姆定律到 LC 谐振

0 阅读15分钟

系列导航:目录 | 上一章:数学工具箱 | 下一章:时域与频域

本章目标:读完本章,你将能够

  • 熟练运用欧姆定律、基尔霍夫定律计算电压、电流和功率,并推导最大功率传输定理;
  • 说明电容和电感的工作原理,算出 RC、RL 电路的时间常数,用 Wheeler 公式估算空心线圈的电感;
  • 用复阻抗分析交流电路,区分有功功率和视在功率,设计一个 RC 低通滤波器;
  • 推导串联与并联谐振的谐振频率、品质因数 Q 和带宽,理解谐振回路的“电压放大”;
  • 从能量角度理解 Q,用数值方法求解 RLC 电路的自由振荡;
  • 为一台中波矿石收音机设计调谐回路:选定可变电容,算出线圈要绕多少匝。

前置知识:第 02 章的复数、相量和分贝。第二、三、七节用到一阶和二阶常微分方程,只需会验证一个函数是否满足方程。

收音机上那个“调台”的旋钮,转动的是什么?

在老式收音机里,它转动的是一个可变电容。电容和一个线圈组成谐振回路:在成千上万个同时到达天线的电台信号中,只有频率等于回路谐振频率的那一个能“通过”,其他的都被挡在外面。转动旋钮改变电容,就改变了谐振频率,也就换了电台。

这就是本章的主线。我们从最基本的电阻电路讲起,逐步引入电容和电感,用第 02 章的复数工具分析它们在交流电路中的行为,最终来到无线电里最重要的电路之一:LC 谐振回路。最后,我们用学到的知识为第 05 章要做的矿石收音机设计调谐回路。

一、电压、电流与电阻

1. 基本物理量

物理量符号单位含义
电荷qq库仑(C)一个电子带 −1.602×10−19-1.602\times10^{-19} C
电流i=dq/dti = dq/dt安培(A)单位时间流过某截面的电荷
电压vv伏特(V)两点之间的电势差,单位电荷从一点移到另一点所做的功
功率p=vip = vi瓦特(W)单位时间消耗或输出的能量

2. 欧姆定律

对于电阻,电压和电流成正比:

v=iRv = iR

比例系数 RR 就是电阻,单位欧姆(Ω)。电阻消耗的功率是

p=vi=i2R=v2Rp = vi = i^2 R = \frac{v^2}{R}

这些功率全部变成热。一台 100 W 的短波电台,用 13.8 V 直流电源供电时,发射状态下电流约 20 A。电源线只要有 0.04 Ω 的电阻,就会产生 0.8 V 的压降,白白发热 16 W(练习 1)。所以电台的电源线都很粗。

3. 基尔霍夫定律

电路分析的两条基本定律:

  • 基尔霍夫电流定律(KCL):流入任一节点的电流之和等于流出的电流之和,因为电荷不会在节点上堆积;
  • 基尔霍夫电压定律(KVL):沿任一闭合回路,各段电压的代数和为零,因为电势是单值的。

由它们可以推出:串联电阻 R=R1+R2R = R_1 + R_2;并联电阻 1R=1R1+1R2\dfrac{1}{R} = \dfrac{1}{R_1} + \dfrac{1}{R_2};串联电阻的分压公式 v2=v R2R1+R2v_2 = v\,\dfrac{R_2}{R_1 + R_2}。这些在第四节换成复阻抗之后,仍然完全成立。

4. 最大功率传输定理

任何实际的信号源都有内阻。把信号源看成一个电压为 VsV_s 的理想电压源串联一个内阻 RsR_s,接上负载 RLR_L:

fig1_max_power.png

图 1:左:有内阻的信号源接负载。右:负载功率和效率随负载电阻的变化。负载等于内阻时功率最大。

回路电流 I=VsRs+RLI = \dfrac{V_s}{R_s + R_L},负载功率

PL=I2RL=Vs2 RL(Rs+RL)2P_L = I^2 R_L = \frac{V_s^2\,R_L}{(R_s + R_L)^2}

RLR_L 太小,电流大但负载上的电压小;RLR_L 太大,电压大但电流小。中间某处功率最大。对 RLR_L 求导:

dPLdRL=Vs2 (Rs+RL)2−2RL(Rs+RL)(Rs+RL)4=Vs2 Rs−RL(Rs+RL)3\frac{dP_L}{dR_L} = V_s^2\,\frac{(R_s + R_L)^2 - 2R_L(R_s + R_L)}{(R_s + R_L)^4} = V_s^2\,\frac{R_s - R_L}{(R_s + R_L)^3}

令它等于零,得到

 RL=Rs⟹PL,max⁡=Vs24Rs \boxed{\ R_L = R_s\quad\Longrightarrow\quad P_{L,\max} = \frac{V_s^2}{4R_s}\ }

负载电阻等于信号源内阻时,负载得到的功率最大,这叫阻抗匹配。附录代码以 Vs=100V_s = 100 V、Rs=50R_s = 50 Ω 为例:

RLR_L12.5 Ω25 Ω50 Ω100 Ω200 Ω
负载功率32.00 W44.44 W50.00 W44.44 W32.00 W
效率20.0 %33.3 %50.0 %66.7 %80.0 %

注意两点:

  1. 匹配时效率只有 50%,一半功率消耗在内阻上。所以电力系统从来不追求匹配(那里要的是效率)。射频系统追求匹配,是因为信号往往很微弱,每一点功率都很宝贵:天线收到的信号要尽可能多地送进接收机。
  2. 实际的发射机输出级并不真是“一个内阻 50 Ω 的电压源”。它被设计成接 50 Ω 负载时输出额定功率、失真最小,所以负载偏离 50 Ω 时,功率会下降,严重时还会损坏功放管。这就是业余电台要求天线系统“驻波比低”的原因,第 08 章讲传输线时详细讨论。

如果信号源内阻是复数 Zs=Rs+jXsZ_s = R_s + jX_s,最大功率传输的条件变成 ZL=Zs∗=Rs−jXsZ_L = Z_s^* = R_s - jX_s,称为共轭匹配:负载的电抗要正好抵消信号源的电抗(练习 10 是一个例子)。

二、电容

1. 定义与平行板电容

两块彼此绝缘的导体,加上电压 vv 后分别带上 +q+q 和 −q-q 的电荷。电荷与电压成正比,比例系数就是电容:

C=qvC = \frac{q}{v}

单位法拉(F)。1 F 非常大,射频电路里常用的是 pF(10−1210^{-12} F)和 nF(10−910^{-9} F)。对于面积为 AA、间距为 dd、中间填充介电常数 ε=εrε0\varepsilon = \varepsilon_r\varepsilon_0 的平行板电容器:

C=εAdC = \frac{\varepsilon A}{d}

附录代码算出:面积 10 cm²、间距 1 mm 的空气电容只有 8.854 pF。要得到更大的电容,就要增大面积、减小间距或者使用介电常数高的材料。

电容并联时电容相加,串联时倒数相加(与电阻正好相反):100 pF 和 220 pF 并联是 320 pF,串联是 68.75 pF。

2. 电压与电流的关系

对 q=Cvq = Cv 求导:

i=C dvdti = C\,\frac{dv}{dt}

电容的电流取决于电压的变化率。直流电压不变,所以电容隔断直流;电压不能突变(否则电流无穷大)。电容储存的能量是

W=12Cv2W = \frac{1}{2}Cv^2

3. RC 电路的充放电

电压源 VV 通过电阻 RR 给初始不带电的电容充电。由 KVL,V=iR+vC=RC dvCdt+vCV = iR + v_C = RC\,\dfrac{dv_C}{dt} + v_C,这个一阶微分方程的解是

vC(t)=V(1−e−t/τ),τ=RCv_C(t) = V\left(1 - e^{-t/\tau}\right),\qquad \tau = RC

τ\tau 称为时间常数。代回方程即可验证:RC⋅Vτe−t/τ+V−Ve−t/τ=VRC\cdot\frac{V}{\tau}e^{-t/\tau} + V - Ve^{-t/\tau} = V。放电时则是 vC=Ve−t/τv_C = Ve^{-t/\tau}。

fig2_rc_transient.png

图 2:RC 电路的充电与放电。每过一个时间常数,离终值的差距就缩小到原来的 1/e。

附录代码算出(R=10R = 10 kΩ,C=100C = 100 nF,τ=1.0\tau = 1.0 ms):

时间τ2τ3τ5τ
充电程度63.2 %86.5 %95.0 %99.3 %

工程上通常认为经过 5τ 就充满了。无线电里时间常数无处不在:电源的滤波、自动增益控制(AGC)的响应速度、CW 电键波形的上升沿(第 28 章会看到,上升沿太陡会产生“键击”干扰)。

4. 容抗

对正弦信号,第 02 章已经由相量推出电容的阻抗

ZC=1jωC=−jXC,XC=1ωCZ_C = \frac{1}{j\omega C} = -jX_C,\qquad X_C = \frac{1}{\omega C}

XCX_C 称为容抗。频率越高,容抗越小:电容“隔直流、通交流”,而且对高频的阻碍比对低频小得多。因子 −j-j 意味着电容的电流比电压超前 90°。射频电路里的耦合电容(让信号通过、挡住直流偏置)和旁路电容(把电源线上的高频干扰短路到地)都是这个原理。

三、电感

1. 定义

电流流过线圈会产生磁场,磁通量与电流成正比。电流变化时,变化的磁通在线圈里感应出电压(法拉第电磁感应定律),方向总是阻碍电流的变化(楞次定律):

v=L didtv = L\,\frac{di}{dt}

比例系数 LL 就是电感,单位亨利(H),射频电路里常用 μH 和 nH。电感储存的能量是 W=12Li2W = \frac12 Li^2。和电容对偶:电感的电流不能突变。突然断开一个有电流的电感,di/dtdi/dt 极大,会产生很高的电压,这就是继电器线圈两端要并联续流二极管的原因。

2. 空心线圈的电感:Wheeler 公式

业余无线电爱好者经常要自己绕线圈。单层密绕空心线圈的电感可以用 Wheeler 在 1928 年给出的经验公式估算:

L [μH]=r2N29r+10 lL\,[\mu\mathrm{H}] = \frac{r^2 N^2}{9r + 10\,l}

其中 NN 是匝数,rr 是线圈半径、ll 是绕线长度,单位都是英寸(1 in = 25.4 mm)。当 l>0.8rl > 0.8r 时,误差约为 1%。

附录代码算出:直径 20 mm、长 20 mm、15 匝的线圈,L=3.055L = 3.055 μH。同一骨架上绕 30 匝、长度也加倍到 40 mm,电感是 7.231 μH,只增加到 2.4 倍,而不是 N2N^2 意味的 4 倍,因为线圈变长后分母也变大了。只有在长度不变、把线绕得更密时,电感才接近与 N2N^2 成正比。

3. RL 电路与感抗

电感与电阻串联接通电压源时,电流按 i=VR(1−e−t/τ)i = \frac{V}{R}\left(1 - e^{-t/\tau}\right) 上升,时间常数 τ=L/R\tau = L/R,规律和 RC 充电完全一样(图 2)。例如 100 μH 与 10 Ω,τ=10.0\tau = 10.0 μs。

对正弦信号,电感的阻抗

ZL=jωL=jXL,XL=ωLZ_L = j\omega L = jX_L,\qquad X_L = \omega L

XLX_L 称为感抗,频率越高越大:电感“通直流、阻交流”。因子 jj 意味着电感的电流比电压滞后 90°。射频扼流圈(RFC)就是一个电感,让直流供电通过,同时挡住高频信号。

fig3_reactance.png

图 3:左:1 μH 的感抗与 100 pF 的容抗随频率的变化(双对数坐标)。右:电感和电容上电压与电流的相位关系。

附录代码列出了几个频率下的电抗:

频率1 MHz7.1 MHz30 MHz145 MHz
XLX_L(1 μH)6.28 Ω44.61 Ω188.50 Ω911.06 Ω
XCX_C(100 pF)1591.55 Ω224.16 Ω53.05 Ω10.98 Ω

两条直线在 15.915 MHz 处相交,此时感抗和容抗都是 100 Ω。这个交点就是第五节的谐振频率。

四、复阻抗与交流电路

1. 阻抗

把电阻、电容、电感统一用复数阻抗描述:

Z=R+jXZ = R + jX

实部 RR 是电阻,虚部 XX 是电抗:X>0X > 0 为感性,X<0X < 0 为容性。相量形式的欧姆定律 V=ZI\mathbf{V} = Z\mathbf{I},以及串并联公式、分压公式,全都和直流电路一样,只是换成了复数运算。阻抗的倒数称为导纳 Y=1/Z=G+jBY = 1/Z = G + jB,分析并联电路时用导纳更方便。

2. 两个例子

50 Ω 电阻串联 1 μH 电感,在 7.1 MHz:

Z=50+jωL=50+j44.61 Ω=67.01∠41.74° ΩZ = 50 + j\omega L = 50 + j44.61\ \Omega = 67.01\angle 41.74°\ \Omega

50 Ω 电阻并联 100 pF 电容,在 7.1 MHz:

Z=50⋅(−j224.16)50−j224.16=47.63−j10.62 Ω=48.8∠−12.57° ΩZ = \frac{50\cdot(-j224.16)}{50 - j224.16} = 47.63 - j10.62\ \Omega = 48.8\angle -12.57°\ \Omega

fig4_impedance_plane.png

图 4:阻抗平面。横轴是电阻,纵轴是电抗;上半平面是感性,下半平面是容性。

阻抗平面是理解匹配网络和史密斯圆图(第 08 章)的基础。比如天线分析仪测出某副天线的阻抗是 50+j44.6150 + j44.61 Ω,你就知道它“偏感性”:串一个合适的电容抵消 j44.61j44.61,就能变成纯 50 Ω(练习 10)。

3. 交流功率

给上面的串联电路加 10 V(有效值)电压,电流 I=V/Z=0.1492∠−41.74°\mathbf{I} = \mathbf{V}/Z = 0.1492\angle -41.74° A,滞后电压 41.74°。附录代码算出:

  • 有功功率 P=∣I∣2R=1.114P = |I|^2 R = 1.114 W:只有电阻消耗能量;
  • 视在功率 S=∣V∣∣I∣=1.492S = |V||I| = 1.492 VA;
  • 功率因数 cos⁡φ=P/S=0.746\cos\varphi = P/S = 0.746。

电感和电容不消耗能量。它们在半个周期里吸收能量,下半个周期又如数还回去,这部分来回流动的功率称为无功功率。所以理想的电感和电容不会发热。实际元件总有一点损耗,下面会用品质因数 Q 来衡量。

4. 第一个滤波器:RC 低通

fig5_rc_lowpass.png

图 5:RC 低通滤波器及其伯德图。截止频率处衰减 3 dB、相移 −45°;高于截止频率后,每十倍频程衰减 20 dB。

把电阻和电容串联,从电容两端取输出。由分压公式:

H(jω)=VoutVin=1/(jωC)R+1/(jωC)=11+jωRCH(j\omega) = \frac{\mathbf{V}_{out}}{\mathbf{V}_{in}} = \frac{1/(j\omega C)}{R + 1/(j\omega C)} = \frac{1}{1 + j\omega RC}

∣H∣=1/1+(ωRC)2|H| = 1/\sqrt{1 + (\omega RC)^2}。低频时 ∣H∣≈1|H| \approx 1,信号直接通过;高频时电容近似短路,信号被衰减。ωRC=1\omega RC = 1 时 ∣H∣=1/2|H| = 1/\sqrt2,即 −3.01 dB,相位 −45°,对应的频率称为截止频率:

fc=12πRCf_c = \frac{1}{2\pi RC}

附录代码取 R=1R = 1 kΩ、C=47C = 47 nF,fc=3386.3f_c = 3386.3 Hz:

频率338.6 Hz(fc/10f_c/10)3000 Hz3386.3 Hz(fcf_c)33863 Hz(10fc10f_c)
∣H∣\lvert H\rvert−0.04 dB−2.52 dB−3.01 dB−20.04 dB
相位−5.71°−41.54°−45.00°−84.29°

这样一个滤波器可以放在接收机的音频输出端,滤掉 3 kHz 以上的“咝咝”噪声。但你也看到了它的局限:衰减太平缓,3 kHz 时已经衰减了 2.5 dB,到 10 倍频率才衰减 20 dB。要得到陡峭的滤波器,需要电感和电容组合的高阶滤波器(第 16 章),或者数字滤波器(第 20 章)。

五、串联谐振

1. 谐振频率

电阻、电感、电容串联,总阻抗

Z=R+j(ωL−1ωC)Z = R + j\left(\omega L - \frac{1}{\omega C}\right)

感抗随频率上升,容抗随频率下降(图 3),总有一个频率两者正好相等、互相抵消,电抗为零:

ω0L=1ω0C⟹ f0=12πLC \omega_0 L = \frac{1}{\omega_0 C}\quad\Longrightarrow\quad \boxed{\ f_0 = \frac{1}{2\pi\sqrt{LC}}\ }

这就是谐振频率。谐振时阻抗最小,等于纯电阻 RR,电流最大,并且与电压同相。低于谐振频率时容抗占优,电路呈容性;高于谐振频率时感抗占优,电路呈感性。

以 40 米波段为例:取 L=2.2L = 2.2 μH,要谐振在 7.1 MHz,需要

C=1ω02L=228.40 pFC = \frac{1}{\omega_0^2 L} = 228.40\ \mathrm{pF}

2. 品质因数与带宽

谐振曲线有多“尖”,用品质因数 QQ 衡量:

Q=ω0LR=1ω0CR=1RLCQ = \frac{\omega_0 L}{R} = \frac{1}{\omega_0 C R} = \frac{1}{R}\sqrt{\frac{L}{C}}

L/C\sqrt{L/C} 是谐振时感抗(也是容抗)的大小,称为回路的特性阻抗。R=5R = 5 Ω 时,附录代码算出 L/C=98.14\sqrt{L/C} = 98.14 Ω,Q=19.63Q = 19.63。

fig6_series_resonance.png

图 6:串联 RLC 的谐振曲线(上)和电流相位(下)。R 越小,Q 越高,曲线越尖。

带宽定义为电流降到最大值的 1/21/\sqrt2(功率降到一半,即 −3 dB)的两个频率之差。这发生在 ∣Z∣=2R|Z| = \sqrt2 R,也就是电抗的大小等于电阻时:

ωL−1ωC=±R⟹ω2L∓Rω−1C=0\omega L - \frac{1}{\omega C} = \pm R\quad\Longrightarrow\quad \omega^2 L \mp R\omega - \frac{1}{C} = 0

取正根:

ω1,2=∓R+R2+4L/C2L\omega_{1,2} = \frac{\mp R + \sqrt{R^2 + 4L/C}}{2L}

两者相减,根号项抵消:

ω2−ω1=RL⟹ B=f2−f1=R2πL=f0Q \omega_2 - \omega_1 = \frac{R}{L}\quad\Longrightarrow\quad \boxed{\ B = f_2 - f_1 = \frac{R}{2\pi L} = \frac{f_0}{Q}\ }

对串联 RLC 电路,这个结果是精确的。两根相乘还可得 ω1ω2=1/(LC)=ω02\omega_1\omega_2 = 1/(LC) = \omega_0^2,即谐振频率是两个半功率频率的几何平均。附录代码在 6~8.2 MHz 之间取了 200 万个点数值搜索,得到半功率点 6.9214 MHz 和 7.2832 MHz,带宽 361.7 kHz,与 f0/Q=361.7f_0/Q = 361.7 kHz 完全一致。

QQ 越高,选择性越好:Q = 19.6 时带宽 362 kHz,几乎占满整个 40 米波段;Q = 98 时带宽只有约 72 kHz。

3. 电压放大

谐振时电流 I=Vs/RI = V_s/R,电感两端的电压

∣VL∣=I⋅ω0L=VsR ω0L=Q Vs|V_L| = I\cdot\omega_0 L = \frac{V_s}{R}\,\omega_0 L = Q\,V_s

电容两端也是 QVsQV_s,只是相位相反,两者互相抵消,所以总电压仍然是 VsV_s。附录代码算出:信号源只有 1 V 时,电感和电容上的电压各有 19.63 V。

这既是好事也是隐患。好处是谐振回路天然有“电压增益”,接收机前端的调谐回路能把天线上的微弱信号放大 Q 倍;隐患是高 Q 回路里的元件要承受很高的电压。业余爱好者常用的磁环天线(第 11 章)Q 值可达数百,接 100 W 发射时,调谐电容两端会有几千伏的高压(练习 6),必须用耐压足够的真空电容或者大间距的空气电容,调节时绝对不能用手碰。

六、并联谐振

1. 理想情况

电感和电容并联,总导纳

Y=jωC+1jωL=j(ωC−1ωL)Y = j\omega C + \frac{1}{j\omega L} = j\left(\omega C - \frac{1}{\omega L}\right)

在同一个 ω0=1/LC\omega_0 = 1/\sqrt{LC} 处 Y=0Y = 0,阻抗为无穷大。这与串联谐振正好相反:串联谐振时阻抗最小,并联谐振时阻抗最大。

2. 带损耗的并联回路

实际的电感线圈总有电阻 rr(导线电阻,以及高频下的趋肤效应)。电感支路 r+jωLr + j\omega L 与电容并联:

Z=(r+jωL)⋅1jωCr+jωL+1jωC=r+jωL1−ω2LC+jωCrZ = \frac{(r + j\omega L)\cdot\dfrac{1}{j\omega C}}{r + j\omega L + \dfrac{1}{j\omega C}} = \frac{r + j\omega L}{1 - \omega^2 LC + j\omega Cr}

在 ω0=1/LC\omega_0 = 1/\sqrt{LC} 处,分母的实部 1−ω02LC=01 - \omega_0^2 LC = 0,若 Q≫1Q \gg 1,分子中的 rr 也可忽略:

Z(ω0)≈jω0Ljω0Cr=LCr=Q2rZ(\omega_0) \approx \frac{j\omega_0 L}{j\omega_0 Cr} = \frac{L}{Cr} = Q^2 r

这个阻抗称为并联回路的谐振阻抗(动态电阻)。仍取 L=2.2L = 2.2 μH、C=228.40C = 228.40 pF、r=5r = 5 Ω,附录代码算出:公式给出 1926.4 Ω,数值计算的最大值是 1928.9 Ω,相差不到 0.2%。一个只有 5 Ω 损耗的线圈,和电容并联后在谐振频率上呈现出将近 2 kΩ 的阻抗。代码还显示,正好在 1/(2πLC)1/(2\pi\sqrt{LC}) 处,阻抗有 −2.92° 的相位,并不是纯电阻:有损耗时,并联回路“阻抗最大”“纯电阻”“1/LC1/\sqrt{LC}”三个条件不再严格重合。Q 越高,三者越接近。

fig7_parallel_tank.png 图 7:带损耗的并联 LC 回路的阻抗。谐振时阻抗最大,约为 L/(Cr);Q 越高,峰越尖。

3. 串联谐振与并联谐振对比

串联谐振并联谐振
谐振时阻抗最小,等于 RR最大,约 Q2rQ^2 r
被放大的量元件上的电压(QQ 倍)回路里循环的电流(QQ 倍)
偏离谐振时阻抗迅速增大阻抗迅速减小
典型用途串在信号通路里,只让谐振频率通过;陷波器接地跨接在信号与地之间,只在谐振频率呈高阻;振荡器的“槽路”、中频变压器、调谐回路

并联谐振回路也叫槽路(tank circuit),因为能量在电感的磁场和电容的电场之间来回倾倒,像水在槽里晃荡。业余天线里的陷波器(trap)就是一个并联 LC:把它串在偶极子的振子上,谐振频率上阻抗很高,相当于把天线“切断”,于是一副天线可以在两个波段上工作。

七、品质因数的能量意义:自由振荡

1. Q 的能量定义

Q 更本质的定义是

Q=2π×回路储存的能量每个周期损耗的能量Q = 2\pi\times\frac{\text{回路储存的能量}}{\text{每个周期损耗的能量}}

验证它与 ω0L/R\omega_0 L/R 一致:串联谐振时电流峰值为 IpI_p,储能在电感和电容之间转移,总量等于电流最大时电感的储能 12LIp2\frac12 L I_p^2;每个周期电阻消耗 Irms2R T=12Ip2R⋅2πω0I_{\mathrm{rms}}^2 R\,T = \frac12 I_p^2 R\cdot\frac{2\pi}{\omega_0}。两者之比乘以 2π2\pi,正好是 ω0L/R\omega_0 L/R。

Q 高,意味着每个周期只损失储能的一小部分,回路能“自己振荡”很久。

2. 自由振荡

把电容先充到 1 V,然后断开电源,让 RLC 回路自己振荡。由 KVL:

Ld2qdt2+Rdqdt+qC=0L\frac{d^2q}{dt^2} + R\frac{dq}{dt} + \frac{q}{C} = 0

当 Q>12Q > \frac12 时,它的解是衰减振荡:

vC(t)=e−αt(cos⁡ωdt+αωdsin⁡ωdt),α=R2L=ω02Q,ωd=ω01−14Q2v_C(t) = e^{-\alpha t}\left(\cos\omega_d t + \frac{\alpha}{\omega_d}\sin\omega_d t\right),\qquad \alpha = \frac{R}{2L} = \frac{\omega_0}{2Q},\quad \omega_d = \omega_0\sqrt{1 - \frac{1}{4Q^2}}

包络按 e−ω0t/2Qe^{-\omega_0 t/2Q} 衰减,每过一个周期,振幅乘以 e−π/Qe^{-\pi/Q}。所以振荡大约 Q 个周期之后,振幅衰减到 e−π≈4.3%e^{-\pi} \approx 4.3\%。这给了 Q 一个直观的含义:一个高 Q 的回路像一口好钟,敲一下能响很久。

fig8_ringdown.png

图 8:RLC 回路的自由振荡(7.1 MHz,Q = 19.63)。

3. 数值实验

附录代码不用上面的解析解,而是用四阶龙格-库塔(RK4)方法直接对微分方程数值积分,每个周期 400 步,然后从波形里“测量”Q:相邻两个正峰值之比取对数,称为对数衰减率 δ\delta,由 δ≈π/Q\delta \approx \pi/Q 反推 Q。结果:

量数值
对数衰减率 δ0.16010
由衰减测得的 Q = π/δ19.62
公式 ω0L/R\omega_0 L/R19.63
5 个周期后的峰值(按 e−nπ/Qe^{-n\pi/Q} 预测)0.4491 V(0.4492 V)
10 个周期后的峰值0.2017 V(0.2018 V)
20 个周期后的峰值0.0407 V(0.0407 V)

数值解和理论在四位有效数字上吻合。这种“写出微分方程、数值积分、再从结果中测量参数”的方法,在后面的章节(比如第 07 章用 FDTD 求解麦克斯韦方程)还会反复使用。

4. 常见元件的 Q

元件Q 的量级
HF 波段的空心线圈、磁环线圈100~300
优质陶瓷电容、空气可变电容1000 以上
石英晶体10410^4~10610^6
微波谐振腔10310^3~10410^4

LC 回路的 Q 通常由线圈决定。石英晶体的 Q 比 LC 回路高好几个数量级,所以用它做振荡器,频率极其稳定;用它做滤波器,带宽可以窄到几百赫兹,这正是接收 CW 和 SSB 信号所需要的。

八、实战:设计中波矿石机的调谐回路

第 05 章我们要动手做一台矿石收音机。它的核心就是一个并联 LC 调谐回路。现在来设计它。

要求:覆盖中波广播 530~1600 kHz。

第一步:电容范围。 选一个常见的电感值 L=240L = 240 μH。由 C=1/(ω02L)C = 1/(\omega_0^2 L),附录代码算出需要的电容从 41.2 pF(1600 kHz)到 375.7 pF(530 kHz)。两者之比 9.11,正好是 (1600/530)2(1600/530)^2:因为 f0∝1/Cf_0 \propto 1/\sqrt{C},频率覆盖比等于电容变化比的平方根。

第二步:选可变电容。 收音机常用的单联可变电容最大约 365 pF。它的最小电容加上线圈、接线的分布电容,总共在 15~25 pF 之间。附录代码算出,配 240 μH 的线圈,调谐范围是 538~2055 kHz(最小 25 pF 时)或 538~2653 kHz(最小 15 pF 时)。高端富余很多;低端的 538 kHz 稍高于 530 kHz,再加上约 10 pF 的分布电容就够了(练习 8)。

fig9_mw_tuning.png

图 9:240 μH 线圈配可变电容的谐振频率。频率覆盖比等于电容变化比的平方根。

第三步:绕线圈。 在直径 50 mm 的骨架(比如一段 PVC 管)上,用 0.3 mm 的漆包线密绕,每匝约占 0.35 mm 长度。线圈长度 l=0.35Nl = 0.35N mm 也随匝数变化,把 Wheeler 公式整理成关于 NN 的一元二次方程:

r2N2−10 L p N−9 L r=0r^2 N^2 - 10\,L\,p\,N - 9\,L\,r = 0

其中 pp 是每匝长度,所有长度以英寸为单位,LL 以 μH 为单位。附录代码解出 N=66.9N = 66.9 匝,绕线长度 23.4 mm,代回 Wheeler 公式核算得 L=240.0L = 240.0 μH。实际制作时绕 67 匝即可。漆包线的实际粗细、绕制的松紧都会让电感有几个百分点的偏差,好在可变电容的调节范围有富余,足以覆盖这点误差。

这样,一个能收听整个中波广播段的调谐回路就设计好了。第 05 章我们会给它加上天线、地线、检波二极管和耳机,组成一台不需要电池的收音机。

小结

  1. 欧姆定律 v=iRv = iR,功率 p=i2R=v2/Rp = i^2R = v^2/R。基尔霍夫电流、电压定律是电路分析的基础。
  2. 最大功率传输:负载电阻等于信号源内阻时,负载功率最大,为 Vs2/4RsV_s^2/4R_s,此时效率 50%。复数情况下要求共轭匹配 ZL=Zs∗Z_L = Z_s^*。
  3. 电容 i=C dv/dti = C\,dv/dt,电压不能突变,隔直通交,容抗 XC=1/ωCX_C = 1/\omega C,电流超前电压 90°。平行板电容 C=εA/dC = \varepsilon A/d。
  4. 电感 v=L di/dtv = L\,di/dt,电流不能突变,通直阻交,感抗 XL=ωLX_L = \omega L,电流滞后电压 90°。单层空心线圈可用 Wheeler 公式 L=r2N2/(9r+10l)L = r^2N^2/(9r + 10l)(英寸,μH)估算。
  5. 时间常数:RC 电路 τ=RC\tau = RC,RL 电路 τ=L/R\tau = L/R;经过 τ 达到终值的 63.2%,经过 5τ 达到 99.3%。
  6. 复阻抗 Z=R+jXZ = R + jX 让交流电路可以像直流电路一样计算。只有电阻消耗有功功率;RC 低通的截止频率 fc=1/(2πRC)f_c = 1/(2\pi RC),之后每十倍频程衰减 20 dB。
  7. 谐振:f0=1/(2πLC)f_0 = 1/(2\pi\sqrt{LC})。串联谐振阻抗最小(= R),元件电压放大 Q 倍;并联谐振阻抗最大(≈ L/Cr = Q²r)。带宽 B=f0/QB = f_0/Q。
  8. 品质因数 Q=ω0L/R=2π×Q = \omega_0L/R = 2\pi\times储能/每周期损耗。自由振荡约 Q 个周期后振幅衰减到 e−πe^{-\pi}。
  9. 频率覆盖比等于调谐电容变化比的平方根;240 μH 的线圈配 365 pF 可变电容即可覆盖中波广播段。

下一章,我们换一个角度看信号:从“随时间变化的波形”转到“由哪些频率组成”,这就是傅里叶分析和频谱。

练习题

  1. (电源线) 一台短波电台输出 100 W 时,从 13.8 V 电源取用 20 A 电流。(a) 电台的整机效率是多少?(b) 电源线正负两根各 2 m,每米电阻 0.01 Ω,电源线上的压降和发热各是多少?
  2. (失配损失) 50 Ω 的信号源接 75 Ω 的负载(比如误用了电视用的 75 Ω 电缆),负载得到的功率是最大可得功率的百分之几?损失多少 dB?
  3. (感抗) 一个 10 μH 的电感,在什么频率下感抗等于 50 Ω?
  4. (串联谐振) 用 1.5 μH 的电感做一个谐振在 14.2 MHz 的串联谐振回路,需要多大的电容?回路的特性阻抗是多少?
  5. (选择性) 一个并联调谐回路谐振在 3.6 MHz,Q = 150。(a) 它的带宽是多少?(b) 相邻 10 kHz 处的另一个电台信号会被衰减多少 dB?仅靠这样一个回路能把两个电台分开吗?
  6. (电压放大) 一个 Q = 200 的串联谐振回路,用 10 V(有效值)的信号激励。谐振时电容两端的电压有效值和峰值各是多少?
  7. (编程) 把附录代码第 9 部分的电阻从 5 Ω 改成 1 Ω,并把振荡周期数增加到 120。测得的 Q 是多少?振荡 98 个周期后,峰值还剩多少?
  8. (调谐范围) 在第八节的设计中,如果可变电容的最大值 365 pF 再加上约 10 pF 的分布电容,调谐范围的低端是多少?
  9. (设计 RC 低通) 用一个 47 nF 的电容做截止频率 3000 Hz 的 RC 低通滤波器,电阻应取多大?
  10. (共轭匹配预览) 第四节的 50+j44.6150 + j44.61 Ω(7.1 MHz)要变成纯电阻 50 Ω,应串联一个多大的电容?

参考答案

  1. (a) 输入功率 13.8×20=27613.8 \times 20 = 276 W,效率 100/276=100/276 = 36.2%,其余约 176 W 变成了热,所以电台要有散热器和风扇。(b) 往返总长 4 m,电阻 0.04 Ω,压降 20×0.04=20 \times 0.04 = 0.8 V,发热 202×0.04=20^2 \times 0.04 = 16 W。电台实际只得到 13.0 V,可能已经低于它的额定工作电压。
  2. 负载功率与最大可得功率之比为 4RsRL(Rs+RL)2=4×50×751252=\dfrac{4R_sR_L}{(R_s + R_L)^2} = \dfrac{4\times50\times75}{125^2} = 96%,即损失 0.18 dB。轻微的失配损失并不大,严重失配(比如图 1 的两端)才真正有影响。
  3. f=XL2πL=502π×10−5=f = \dfrac{X_L}{2\pi L} = \dfrac{50}{2\pi\times10^{-5}} = 795.8 kHz。
  4. C=1(2π×14.2×106)2×1.5×10−6=C = \dfrac{1}{(2\pi\times14.2\times10^6)^2\times1.5\times10^{-6}} = 83.75 pF,特性阻抗 L/C=\sqrt{L/C} = 133.8 Ω。
  5. (a) B=3.6 MHz/150=B = 3.6\ \mathrm{MHz}/150 = 24 kHz。(b) 对高 Q 回路,偏离谐振 Δf\Delta f 处的相对响应约为 1/1+(2QΔf/f0)21/\sqrt{1 + (2Q\Delta f/f_0)^2}。2×150×10/3600=0.8332\times150\times10/3600 = 0.833,衰减 10lg⁡(1+0.8332)=10\lg(1 + 0.833^2) = 2.29 dB(按并联回路的完整公式计算是 2.29 dB)。只衰减了 2 dB 多,远远不够把两个电台分开。这就是为什么简单的调谐回路选择性很差,而超外差接收机要把信号变到固定的中频,再用晶体或陶瓷滤波器滤波(第 15 章)。
  6. VC=QVs=V_C = QV_s = 2000 V(有效值),峰值 20002=2000\sqrt2 = 2828 V。
  7. 测得 Q=Q = 98.14(公式值 ω0L/R=98.14\omega_0L/R = 98.14)。振荡 98 个周期后峰值剩 0.0434 V,约为 e−π=0.0432e^{-\pi} = 0.0432。
  8. Cmax⁡=375C_{\max} = 375 pF 时,f=12π240×10−6×375×10−12=f = \dfrac{1}{2\pi\sqrt{240\times10^{-6}\times375\times10^{-12}}} = 530.5 kHz,刚好覆盖中波段的低端。
  9. R=12πfcC=12π×3000×47×10−9=R = \dfrac{1}{2\pi f_c C} = \dfrac{1}{2\pi\times3000\times47\times10^{-9}} = 1129 Ω,取标准值 1.1 kΩ(截止频率约 3.08 kHz)或 1.2 kΩ(约 2.82 kHz)。
  10. 需要 XC=44.61X_C = 44.61 Ω,C=12π×7.1×106×44.61=C = \dfrac{1}{2\pi\times7.1\times10^6\times44.61} = 502.5 pF。串联后 Z=50+j44.61−j44.61=50Z = 50 + j44.61 - j44.61 = 50 Ω。注意这只在 7.1 MHz 成立,频率一变,感抗和容抗就不再相等。第 08 章会系统地讨论阻抗匹配。

参考资料

  1. W. H. Hayt, J. E. Kemmerly, J. D. Phillips, S. M. Durbin. Engineering Circuit Analysis. 9th ed. McGraw-Hill, 2018. 第 2~10 章、第 16 章(频率响应与谐振)。
  2. J. B. Hagen. Radio-Frequency Electronics: Circuits and Applications. 2nd ed. Cambridge University Press, 2009. 谐振电路与 Q 的相关章节。
  3. P. Horowitz, W. Hill. The Art of Electronics. 3rd ed. Cambridge University Press, 2015. 第 1 章。
  4. W. Hayward, R. Campbell, B. Larkin. Experimental Methods in RF Design. ARRL, 2003. 滤波器与谐振回路的相关章节。
  5. H. A. Wheeler. Simple Inductance Formulas for Radio Coils. Proceedings of the IRE, 16(10): 1398–1400, 1928.
  6. ARRL. The ARRL Handbook for Radio Communications. 最新版. “Electrical Fundamentals” 章节。

附录:完整代码

下面的代码复现了本章所有数值:最大功率传输、平行板电容与 RC 时间常数、Wheeler 公式、电抗、复阻抗与交流功率、RC 低通、串联谐振的 Q 与带宽(含数值核对)、并联谐振阻抗、用 RK4 求解自由振荡并测量 Q,以及中波矿石机调谐回路的设计,最后画出串联 RLC 的谐振曲线。只依赖 numpy 和 matplotlib,保存为 example.py 后运行 python example.py 即可。

"""无线电自学教程 03:电路基础 —— 附录代码

1. 欧姆定律与最大功率传输:50 Ω 信号源接不同负载;
2. 电容:平行板电容、RC 充电的时间常数;
3. 电感:单层空心线圈的 Wheeler 公式、RL 时间常数;
4. 电抗:1 μH 与 100 pF 在 7.1 MHz 的电抗,以及两者相等的频率;
5. 复阻抗:50 Ω 串 1 μH,50 Ω 并 100 pF,交流功率与功率因数;
6. RC 低通滤波器:截止频率、-3 dB 与 -45°;
7. 串联谐振:谐振频率、Q、带宽(数值核对)、电压放大;
8. 并联谐振:谐振阻抗 L/(Cr) 的数值核对;
9. 自由振荡:数值求解 RLC 电路的衰减振荡,用对数衰减率测 Q;
10. 中波矿石机的调谐回路:可变电容的范围与线圈匝数。
最后画出串联 RLC 的谐振曲线。
运行:python example.py
"""
import numpy as np
import matplotlib.pyplot as plt

plt.rcParams["font.sans-serif"] = ["Microsoft YaHei", "SimHei", "DejaVu Sans"]
plt.rcParams["axes.unicode_minus"] = False

EPS0 = 8.8541878128e-12        # 真空介电常数 [F/m]
INCH = 25.4e-3                 # 1 英寸 [m]


def polar_str(z, unit=""):
    return f"{abs(z):.4g}{unit}∠{np.degrees(np.angle(z)):.2f}°"


def rect_str(z, unit=""):
    sign = "+" if z.imag >= 0 else "-"
    return f"{z.real:.4g} {sign} j{abs(z.imag):.4g}{unit}"


def wheeler_uh(d, length, n):
    """单层空心线圈电感 [μH]。d:线圈直径 [m],length:绕线长度 [m],n:匝数。"""
    r, l = d / 2 / INCH, length / INCH
    return r**2 * n**2 / (9 * r + 10 * l)


def resonant_c(f, L):
    return 1 / ((2 * np.pi * f) ** 2 * L)


# ---------------- 1. 欧姆定律与最大功率传输 ----------------
print("=== 1. 最大功率传输:开路电压 100 V(有效值)、内阻 50 Ω ===")
VS, RS = 100.0, 50.0
for rl in (12.5, 25.0, 50.0, 100.0, 200.0):
    i = VS / (RS + rl)
    p_load = i**2 * rl
    eff = rl / (RS + rl)
    print(f"R_L = {rl:6.1f} Ω:电流 {i:.3f} A,负载功率 {p_load:6.2f} W,效率 {eff * 100:5.1f} %")
print(f"理论最大值 V^2/(4R_s) = {VS**2 / (4 * RS):.2f} W")

# ---------------- 2. 电容 ----------------
print("\n=== 2. 电容 ===")
A, d = 10e-4, 1e-3
print(f"平行板:面积 10 cm^2、间距 1 mm、空气:C = εA/d = {EPS0 * A / d * 1e12:.3f} pF")
R, Cc = 10e3, 100e-9
tau = R * Cc
print(f"RC 充电:R = 10 kΩ,C = 100 nF,τ = RC = {tau * 1e3:.1f} ms")
for k in (1, 2, 3, 5):
    print(f"  t = {k}τ:电容电压达到电源电压的 {(1 - np.exp(-k)) * 100:.1f} %")
c1, c2 = 100e-12, 220e-12
print(f"100 pF 与 220 pF:并联 {(c1 + c2) * 1e12:.1f} pF,串联 {c1 * c2 / (c1 + c2) * 1e12:.2f} pF")

# ---------------- 3. 电感 ----------------
print("\n=== 3. 电感 ===")
L_coil = wheeler_uh(20e-3, 20e-3, 15)
print(f"Wheeler 公式:直径 20 mm、长 20 mm、15 匝的空心线圈 L = {L_coil:.3f} μH")
print(f"同一线圈绕 30 匝(长度也加倍):L = {wheeler_uh(20e-3, 40e-3, 30):.3f} μH")
print(f"RL 电路:L = 100 μH,R = 10 Ω,τ = L/R = {100e-6 / 10 * 1e6:.1f} μs")

# ---------------- 4. 电抗 ----------------
print("\n=== 4. 电抗 X_L = ωL,X_C = 1/(ωC) ===")
L1, C1 = 1e-6, 100e-12
for f in (1e6, 7.1e6, 30e6, 145e6):
    w = 2 * np.pi * f
    print(f"f = {f / 1e6:6.1f} MHz:X_L(1 μH) = {w * L1:8.2f} Ω,X_C(100 pF) = {1 / (w * C1):8.2f} Ω")
print(f"两者相等的频率 f = 1/(2π sqrt(LC)) = {1 / (2 * np.pi * np.sqrt(L1 * C1)) / 1e6:.3f} MHz")

# ---------------- 5. 复阻抗与交流功率 ----------------
print("\n=== 5. 复阻抗(f = 7.1 MHz) ===")
w = 2 * np.pi * 7.1e6
z_series = 50 + 1j * w * L1
print(f"50 Ω 串 1 μH:Z = {rect_str(z_series, ' Ω')} = {polar_str(z_series, ' Ω')}")
V = 10.0
I = V / z_series
print(f"  接 10 V(有效值):I = {polar_str(I, ' A')}")
print(f"  有功功率 P = |I|^2 R = {abs(I) ** 2 * 50:.3f} W,"
      f"视在功率 S = |V||I| = {V * abs(I):.3f} VA,功率因数 cos φ = {np.cos(np.angle(z_series)):.3f}")
zc = 1 / (1j * w * C1)
z_par = 50 * zc / (50 + zc)
print(f"50 Ω 并 100 pF:Z = {rect_str(z_par, ' Ω')} = {polar_str(z_par, ' Ω')}")

# ---------------- 6. RC 低通 ----------------
print("\n=== 6. RC 低通:R = 1 kΩ,C = 47 nF ===")
R6, C6 = 1e3, 47e-9
fc = 1 / (2 * np.pi * R6 * C6)
print(f"截止频率 f_c = 1/(2πRC) = {fc:.1f} Hz")
for f in (fc / 10, fc, 3e3, 10 * fc):
    H = 1 / (1 + 1j * 2 * np.pi * f * R6 * C6)
    print(f"  f = {f:8.1f} Hz:|H| = {20 * np.log10(abs(H)):7.2f} dB,相位 {np.degrees(np.angle(H)):7.2f}°")

# ---------------- 7. 串联谐振 ----------------
print("\n=== 7. 串联 RLC:L = 2.2 μH,谐振于 7.1 MHz,R = 5 Ω ===")
L7, R7, F0 = 2.2e-6, 5.0, 7.1e6
C7 = resonant_c(F0, L7)
w0 = 2 * np.pi * F0
Q7 = w0 * L7 / R7
print(f"C = 1/(ω0^2 L) = {C7 * 1e12:.2f} pF")
print(f"Q = ω0 L / R = {Q7:.2f},也等于 sqrt(L/C)/R = {np.sqrt(L7 / C7) / R7:.2f}")
print(f"特性阻抗 sqrt(L/C) = {np.sqrt(L7 / C7):.2f} Ω,带宽 f0/Q = {F0 / Q7 / 1e3:.1f} kHz")
f = np.linspace(6e6, 8.2e6, 2_000_001)
Z = R7 + 1j * (2 * np.pi * f * L7 - 1 / (2 * np.pi * f * C7))
Imag = 1 / np.abs(Z)
band = f[Imag >= Imag.max() / np.sqrt(2)]
print(f"数值核对:电流最大处 {f[np.argmax(Imag)] / 1e6:.4f} MHz,"
      f"半功率点 {band[0] / 1e6:.4f} ~ {band[-1] / 1e6:.4f} MHz,带宽 {(band[-1] - band[0]) / 1e3:.1f} kHz")
print(f"信号源 1 V(有效值):谐振时电流 {1 / R7:.2f} A,电感和电容上的电压各 {Q7:.2f} V")

# ---------------- 8. 并联谐振 ----------------
print("\n=== 8. 并联谐振:同样的 L、C,电感的损耗电阻 r = 5 Ω ===")
r8 = 5.0
zp = lambda f: (r8 + 1j * 2 * np.pi * f * L7) / (1 - (2 * np.pi * f) ** 2 * L7 * C7 + 1j * 2 * np.pi * f * C7 * r8)
Zp = zp(f)
k = np.argmax(np.abs(Zp))
print(f"公式 L/(C r) = {L7 / (C7 * r8):.1f} Ω,Q^2 r = {Q7**2 * r8:.1f} Ω")
print(f"数值:|Z| 最大 {abs(Zp[k]):.1f} Ω,出现在 {f[k] / 1e6:.4f} MHz;"
      f"在 1/(2π sqrt(LC)) 处 |Z| = {abs(zp(F0)):.1f} Ω,相位 {np.degrees(np.angle(zp(F0))):.2f}°")

# ---------------- 9. 自由振荡与 Q ----------------
print("\n=== 9. 自由振荡:电容先充到 1 V,然后让 RLC 回路自己振荡 ===")


def ring_down(L, C, R, cycles=40, steps_per_cycle=400):
    """RK4 求解 L di/dt + R i + q/C = 0,返回时间和电容电压。"""
    T0 = 2 * np.pi * np.sqrt(L * C)
    dt = T0 / steps_per_cycle
    n = int(cycles * steps_per_cycle)
    q, i = C * 1.0, 0.0
    vc = np.empty(n)

    def deriv(q, i):
        return i, -(R * i + q / C) / L

    for m in range(n):
        vc[m] = q / C
        k1 = deriv(q, i)
        k2 = deriv(q + dt / 2 * k1[0], i + dt / 2 * k1[1])
        k3 = deriv(q + dt / 2 * k2[0], i + dt / 2 * k2[1])
        k4 = deriv(q + dt * k3[0], i + dt * k3[1])
        q += dt / 6 * (k1[0] + 2 * k2[0] + 2 * k3[0] + k4[0])
        i += dt / 6 * (k1[1] + 2 * k2[1] + 2 * k3[1] + k4[1])
    return np.arange(n) * dt, vc


t9, vc9 = ring_down(L7, C7, R7)
peaks = [0] + [m for m in range(1, len(vc9) - 1)
               if vc9[m] > vc9[m - 1] and vc9[m] >= vc9[m + 1] and vc9[m] > 0]
pk = vc9[peaks]                                   # pk[n]:第 n 个周期结束时的峰值
delta = np.mean(np.log(pk[:-1] / pk[1:]))
print(f"相邻正峰值之比的对数(对数衰减率)δ = {delta:.5f}")
print(f"由衰减测得 Q = π/δ = {np.pi / delta:.2f},公式 ω0L/R = {Q7:.2f}")
for n in (5, 10, 20):
    print(f"振荡 {n:2d} 个周期后,峰值剩 {pk[n]:.4f} V(按 e^(-nπ/Q) 应为 {np.exp(-n * np.pi / Q7):.4f})")

# ---------------- 10. 中波矿石机的调谐回路 ----------------
print("\n=== 10. 中波调谐:L = 240 μH,覆盖 530 ~ 1600 kHz ===")
L10 = 240e-6
c_lo, c_hi = resonant_c(1600e3, L10), resonant_c(530e3, L10)
print(f"需要的电容:{c_lo * 1e12:.1f} pF(1600 kHz)~ {c_hi * 1e12:.1f} pF(530 kHz),"
      f"比值 {c_hi / c_lo:.2f} = (1600/530)^2 = {(1600 / 530) ** 2:.2f}")
d10, pitch = 50e-3, 0.35e-3                 # 线圈骨架直径、每匝占的长度(0.3 mm 漆包线密绕)
r_in, p_in, L_uh = d10 / 2 / INCH, pitch / INCH, L10 * 1e6
N = (10 * L_uh * p_in + np.sqrt((10 * L_uh * p_in) ** 2 + 36 * L_uh * r_in**3)) / (2 * r_in**2)
print(f"在直径 50 mm 的骨架上密绕:需要 {N:.1f} 匝,绕线长度 {N * pitch * 1e3:.1f} mm,"
      f"核算 L = {wheeler_uh(d10, N * pitch, N):.1f} μH")
for cmin in (15e-12, 25e-12):
    print(f"可变电容 {cmin * 1e12:.0f} ~ 365 pF 配 240 μH:调谐范围 "
          f"{1 / (2 * np.pi * np.sqrt(L10 * 365e-12)) / 1e3:.0f} ~ {1 / (2 * np.pi * np.sqrt(L10 * cmin)) / 1e3:.0f} kHz")

# ---------------- 画图:串联谐振曲线 ----------------
fig, ax = plt.subplots(figsize=(8, 4.5))
fx = np.linspace(6e6, 8.2e6, 2000)
for R_ in (20.0, 5.0, 1.0):
    Zx = R_ + 1j * (2 * np.pi * fx * L7 - 1 / (2 * np.pi * fx * C7))
    ax.plot(fx / 1e6, 20 * np.log10(R_ / np.abs(Zx)), label=f"R = {R_:g} Ω,Q = {w0 * L7 / R_:.0f}")
ax.set_xlabel("频率 / MHz")
ax.set_ylabel("电流(相对谐振时)/ dB")
ax.set_title("串联 RLC 谐振曲线:Q 越高,曲线越尖")
ax.set_ylim(-40, 2)
ax.grid(True, color="#e6e5e1")
ax.legend()
fig.tight_layout()
plt.show()

运行后的完整输出:

=== 1. 最大功率传输:开路电压 100 V(有效值)、内阻 50 Ω ===
R_L =   12.5 Ω:电流 1.600 A,负载功率  32.00 W,效率  20.0 %
R_L =   25.0 Ω:电流 1.333 A,负载功率  44.44 W,效率  33.3 %
R_L =   50.0 Ω:电流 1.000 A,负载功率  50.00 W,效率  50.0 %
R_L =  100.0 Ω:电流 0.667 A,负载功率  44.44 W,效率  66.7 %
R_L =  200.0 Ω:电流 0.400 A,负载功率  32.00 W,效率  80.0 %
理论最大值 V^2/(4R_s) = 50.00 W

=== 2. 电容 ===
平行板:面积 10 cm^2、间距 1 mm、空气:C = εA/d = 8.854 pF
RC 充电:R = 10 kΩ,C = 100 nF,τ = RC = 1.0 ms
  t = 1τ:电容电压达到电源电压的 63.2 %
  t = 2τ:电容电压达到电源电压的 86.5 %
  t = 3τ:电容电压达到电源电压的 95.0 %
  t = 5τ:电容电压达到电源电压的 99.3 %
100 pF 与 220 pF:并联 320.0 pF,串联 68.75 pF

=== 3. 电感 ===
Wheeler 公式:直径 20 mm、长 20 mm、15 匝的空心线圈 L = 3.055 μH
同一线圈绕 30 匝(长度也加倍):L = 7.231 μH
RL 电路:L = 100 μH,R = 10 Ω,τ = L/R = 10.0 μs

=== 4. 电抗 X_L = ωL,X_C = 1/(ωC) ===
f =    1.0 MHz:X_L(1 μH) =     6.28 Ω,X_C(100 pF) =  1591.55 Ω
f =    7.1 MHz:X_L(1 μH) =    44.61 Ω,X_C(100 pF) =   224.16 Ω
f =   30.0 MHz:X_L(1 μH) =   188.50 Ω,X_C(100 pF) =    53.05 Ω
f =  145.0 MHz:X_L(1 μH) =   911.06 Ω,X_C(100 pF) =    10.98 Ω
两者相等的频率 f = 1/(2π sqrt(LC)) = 15.915 MHz

=== 5. 复阻抗(f = 7.1 MHz) ===
50 Ω 串 1 μH:Z = 50 + j44.61 Ω = 67.01 Ω∠41.74°
  接 10 V(有效值):I = 0.1492 A∠-41.74°
  有功功率 P = |I|^2 R = 1.114 W,视在功率 S = |V||I| = 1.492 VA,功率因数 cos φ = 0.746
50 Ω 并 100 pF:Z = 47.63 - j10.62 Ω = 48.8 Ω∠-12.57°

=== 6. RC 低通:R = 1 kΩ,C = 47 nF ===
截止频率 f_c = 1/(2πRC) = 3386.3 Hz
  f =    338.6 Hz:|H| =   -0.04 dB,相位   -5.71°
  f =   3386.3 Hz:|H| =   -3.01 dB,相位  -45.00°
  f =   3000.0 Hz:|H| =   -2.52 dB,相位  -41.54°
  f =  33862.8 Hz:|H| =  -20.04 dB,相位  -84.29°

=== 7. 串联 RLC:L = 2.2 μH,谐振于 7.1 MHz,R = 5 Ω ===
C = 1/(ω0^2 L) = 228.40 pF
Q = ω0 L / R = 19.63,也等于 sqrt(L/C)/R = 19.63
特性阻抗 sqrt(L/C) = 98.14 Ω,带宽 f0/Q = 361.7 kHz
数值核对:电流最大处 7.1000 MHz,半功率点 6.9214 ~ 7.2832 MHz,带宽 361.7 kHz
信号源 1 V(有效值):谐振时电流 0.20 A,电感和电容上的电压各 19.63 V

=== 8. 并联谐振:同样的 L、C,电感的损耗电阻 r = 5 Ω ===
公式 L/(C r) = 1926.4 Ω,Q^2 r = 1926.4 Ω
数值:|Z| 最大 1928.9 Ω,出现在 7.1000 MHz;在 1/(2π sqrt(LC)) 处 |Z| = 1928.9 Ω,相位 -2.92°

=== 9. 自由振荡:电容先充到 1 V,然后让 RLC 回路自己振荡 ===
相邻正峰值之比的对数(对数衰减率)δ = 0.16010
由衰减测得 Q = π/δ = 19.62,公式 ω0L/R = 19.63
振荡  5 个周期后,峰值剩 0.4491 V(按 e^(-nπ/Q) 应为 0.4492)
振荡 10 个周期后,峰值剩 0.2017 V(按 e^(-nπ/Q) 应为 0.2018)
振荡 20 个周期后,峰值剩 0.0407 V(按 e^(-nπ/Q) 应为 0.0407)

=== 10. 中波调谐:L = 240 μH,覆盖 530 ~ 1600 kHz ===
需要的电容:41.2 pF(1600 kHz)~ 375.7 pF(530 kHz),比值 9.11 = (1600/530)^2 = 9.11
在直径 50 mm 的骨架上密绕:需要 66.9 匝,绕线长度 23.4 mm,核算 L = 240.0 μH
可变电容 15 ~ 365 pF 配 240 μH:调谐范围 538 ~ 2653 kHz
可变电容 25 ~ 365 pF 配 240 μH:调谐范围 538 ~ 2055 kHz