“电台”是文件激励源

from math import *
CARRIER1 = 1593E3
T = 0.25 / CARRIER1
PI2 = 2 * pi
CARRIER1 *= PI2
CARRIER2 = 639E3 * PI2
MIDC = 261.6 * PI2
MIDA = 440 * PI2
COEF = 0.5
for t in range(25000):
  t *= T
  s = COEF * sin(CARRIER1 * t) * sin(MIDC * t) + \
    COEF * sin(CARRIER2 * t) * sin(MIDA * t)
  print(t, s)

• 1593和639千周是中央人民广播电台的频率。千赫(兹)旧称千周。
• 261.6 Hz和440 Hz分别是钢琴的中央C和A。
• ω = 2πf
• 采样频率为1593 KHz的4倍。(25000*T)/(1/261.6)≈1
• 调幅是信号相乘。
• 系数0.5么，经实验0.01不行。
• 想把639 Khz检出来，LC≈6.2E-14，L取C的1000倍以上。

电路从One Transistor Radio抄的

国内网上有很多电路图，我为啥抄老外的呢？不是崇洋媚外而是：

1. NPN的电流方向我觉得顺眼
2. 他这个好像最简单
3. 简介了原理——虽然我只看懂了一点点。国内的么，刻薄地说有点“不想和外行多话；你看我头像牛B就行了”

结果

在Proteus里使用文件/音频激励源