【今日份GPT废话大赏——最费滚轮的一集】
为什么GPT生成的东西总是—— 换行。 说一句话,没说完。 然后。 就再换一行啊!!!
文章目录
- 第一章 多普勒频移到底改变了什么:从一束光到一个运动光源
-
- 1.1 光源不动时,观察者看到的频率是什么
- 1.2 为什么光的多普勒效应不能直接照搬声波公式
- 1.3 远离和靠近的符号一定要先定义清楚
- 第二章 低速多普勒为什么看起来很简单:工程模型与误差来源
-
- 2.1 低速情况下最常用的线性公式
- 2.2 为什么低速运动特别适合连续波激光测量
- 2.3 低速时最容易犯的一个错误:把“频率变化”误认为“固定频移”
- 第三章 高速运动进入相对论区间:时间膨胀和传播时延到底谁在起作用
-
- 3.1 相对论 Doppler 不是简单的“加一个时间膨胀”
- 3.2 传播时延才是第二个不能忽略的因素
- 3.3 为什么高速下经典公式会明显失效
- 第四章 光频梳遇到多普勒:为什么问题突然从“一个频率”变成“整个频谱”
-
- 4.1 光频梳不是一根谱线,而是一排高度规则的谱线
- 4.2 理想情况下,多普勒会让整个光频梳一起“缩放”
- 4.3 为什么不同纵模看起来会产生不同的绝对频移
- 4.4 光频梳的重复频率也会动,这意味着什么
- 第五章 光频梳高速运动的真正难点:时间、频率和脉冲相位会一起变化
-
- 5.1 光频梳的 Doppler 不只有“梳齿移动”
- 5.2 加速运动为什么会产生光学 chirp
- 5.3 光频梳中不同纵模会出现不同的 Doppler 相位斜率
- 第六章 双光梳中的 Doppler:为什么它比连续激光复杂很多
-
- 6.1 双光梳系统为什么要把光学频率搬到 RF
- 6.2 运动光频梳和静止光频梳做拍频会发生什么
- 6.3 这也是为什么光频梳可以拿来测速
- 第七章 异步采样和运动 Doppler:为什么速度会干扰时间测量
-
- 7.1 异步光学采样本来就是一个时间扫描器
- 7.2 多普勒会同时进入“频率轴”和“时间轴”
- 7.3 这就是时延-多普勒耦合为什么会出现
- 第八章 TPFC 为什么在运动目标场景有意思:从“扫描”变成“跟踪”
-
- 8.1 普通双梳和 TPFC 的思路差异
- 8.2 TPFC 锁定目标后,不意味着 Doppler 消失
- 8.3 TPFC 的时序反馈量本身就带着运动信息
- 第九章 光源低速、高速、线性和非线性运动应该怎么理解
-
- 9.1 低速恒速:最简单的情况
- 9.2 低速线性加速:Doppler 变成一条直线
- 9.3 高速线性加速:速度线性不意味着 Doppler 线性
- 9.4 非线性速度:STFT 会出现更加复杂的曲线
- 第十章 相对论时间膨胀对光频梳到底意味着什么
-
- 10.1 频率是时间的倒数,所以时间尺度变化必然进入频梳
- 10.2 不要把时间膨胀和传播 Doppler 分成两个互不相关的补丁
- 第十一章 连续激光和光频梳的 Doppler,到底应该怎么对比
- 第十二章 一个特别容易踩的坑:光频梳“整体平移”和“整体伸缩”不是一回事
-
- 12.1 固定频移和比例缩放有什么区别
- 12.2 为什么这个区别对双光梳特别重要
- 第十三章 高速运动下的频谱图应该怎么看
-
- 13.1 连续激光的频谱
- 13.2 加速运动时,单次 FFT 为什么会变宽
- 13.3 光频梳的 STFT 会出现什么
- 第十四章 为什么“运动速度越高,光频梳越适合做 Doppler 测量”并不完全正确
-
- 14.1 高速度意味着信号变大,也意味着模型变复杂
- 14.2 高速运动会让“低速近似”越来越危险
- 第十五章 对双光梳时频传递而言,真正需要跟踪的不是一个频移
-
- 15.1 单纯测量 Doppler 和完成时间同步不是同一个目标
- 15.2 TPFC 可以把“跟踪运动目标”变成一个闭环问题
- 第十六章 如果用数学模型描述整个光频梳运动过程,可以建立什么框架
- 第十七章 工程仿真中应该怎样分层,否则很容易把模型写乱
-
- 17.1 第一层:运动学模型
- 17.2 第二层:传播模型
- 17.3 第三层:光学 Doppler 模型
- 17.4 第四层:相干相位模型
- 第十八章 哪些情况下可以用简化公式,哪些情况下必须上相对论模型
- 第十九章 对光频梳来说,真正值得关注的几个 Doppler 观测量
- 第二十章 一个非常重要的认识:光频梳让 Doppler 从“一个数字”变成“一个结构”
- 第二十一章 为什么“不同纵模伸缩”是一个很值得继续研究的问题
- 第二十二章 和你的双光梳时频传递研究结合起来应该怎么理解
- 第二十三章 为什么单纯看 Doppler 频率还不够
- 第二十四章 最后把整个问题压缩成几句话
- 第二十五章 从仿真角度,最值得做的不是“画一张 Doppler 曲线”
在光通信、激光测距、激光雷达、光学时频传递、双光梳测距以及自由空间光频梳同步系统里,“多普勒频移”几乎是绕不开的问题。
很多资料把它简单写成:
f
′
=
f
(
1
−
v
c
)
f'=f\\left(1-\\frac{v}{c}\\right)
f′=f(1−cv),然后告诉你“运动速度越大,频率变化越大”。这个公式在低速情况下没有问题,但真正进入光频梳、高速运动、动态时频传递以及相对论速度范围以后,仅仅记住这一条公式远远不够。
尤其到了光频梳系统,问题会从“一个激光频率偏了多少”进一步变成:
光源运动
→
光学载频变化
→
梳齿整体伸缩
→
f
r
e
p
变化
→
脉冲到达时间变化
→
双梳拍频变化
→
时延-多普勒耦合
\\text{光源运动} \\rightarrow \\text{光学载频变化} \\rightarrow \\text{梳齿整体伸缩} \\rightarrow f_{\\mathrm{rep}}\\text{变化} \\rightarrow \\text{脉冲到达时间变化} \\rightarrow \\text{双梳拍频变化} \\rightarrow \\text{时延-多普勒耦合}
光源运动→光学载频变化→梳齿整体伸缩→frep变化→脉冲到达时间变化→双梳拍频变化→时延–多普勒耦合
也就是说,对连续激光来说,多普勒首先表现为频率问题;对光频梳来说,多普勒同时表现为频率、脉冲时序、重复频率、载波相位和拍频结构的问题。
这也是为什么在双光梳时频传递系统中,不能简单地把多普勒频移当成一个“额外的 MHz 偏移量”扔进仿真程序。
本文不讨论代码,而是从物理图像出发,把这个问题从连续波激光一直梳理到光频梳、双光梳和高速相对论运动。重点放在**“到底发生了什么、为什么发生、仿真时应该怎么理解”**。
第一章 多普勒频移到底改变了什么:从一束光到一个运动光源
1.1 光源不动时,观察者看到的频率是什么
先从最简单的连续波激光开始。
假设有一个理想单频激光器,发出的电场可以写成:
E
(
t
)
=
E
0
cos
(
2
π
f
0
t
+
ϕ
0
)
E(t)=E_0\\cos(2\\pi f_0t+\\phi_0)
E(t)=E0cos(2πf0t+ϕ0)
其中:
f
0
f_0
f0是激光器本身的光学频率。
以常见的 1550 nm 激光器为例:
f
0
=
c
λ
0
f_0=\\frac{c}{\\lambda_0}
f0=λ0c,大约为:
f
0
≈
193.5
T
H
z
f_0\\approx193.5\\ \\mathrm{THz}
f0≈193.5 THz
所以当我们说“1550 nm 激光”,实际上对应的是一个接近 193 THz 的超高频电磁振荡。
如果激光源和接收机都静止,并且传播介质也不考虑变化,那么接收端仍然会测到接近: 193.5的光学频率。
真正的问题从光源开始运动以后出现。假设光源沿着观察方向运动。如果光源远离观察者,那么后续发出的每一个波峰都会比前一个波峰稍微多走一点距离,导致观察者接收到的波峰时间间隔变长。
波峰时间间隔变长意味着:
f
o
b
s
<
f
0
f_{\\mathrm{obs}}<f_0
fobs<f0
这就是红移。
反过来,如果光源朝观察者运动,那么连续发出的波峰之间传播距离逐渐缩短:
f
o
b
s
>
f
0
f_{\\mathrm{obs}}>f_0
fobs>f0
这就是蓝移。
所以从物理图像来看,多普勒效应其实可以先理解成一句非常直观的话:
运动改变了观察者收到连续波峰的时间间隔。
频率变化只是时间间隔变化在频域中的表现。
1.2 为什么光的多普勒效应不能直接照搬声波公式
很多初学者接触多普勒效应时,首先学到的是声波。
例如声源向外运动时,经常写成:
f
′
=
f
1
+
v
/
v
s
f'=\\frac{f}{1+v/v_{\\mathrm{s}}}
f′=1+v/vsf
这里的传播速度是声速。
光波和声波最大的不同,是光不需要介质传播,而且真空中的光速对所有惯性观察者都相同。
因此光学多普勒效应必须建立在狭义相对论框架下,而不能直接把声波公式中的声速替换成光速。
低速时,两种思路的结果看起来很接近,这是因为:
v
c
≪
1
\\frac{v}{c}\\ll1
cv≪1
此时高阶相对论修正非常小。
所以工程领域经常看到:
f
o
b
s
≈
f
0
(
1
−
v
c
)
f_{\\mathrm{obs}}\\approx f_0\\left(1-\\frac{v}{c}\\right)
fobs≈f0(1−cv)
这个公式不是一个完全独立的“光学多普勒定律”,而是相对论 Doppler 公式在低速条件下的一阶近似。
1.3 远离和靠近的符号一定要先定义清楚
这是做仿真时最容易出现的问题之一。
本文统一采用v>0 表示光源远离观察者。
定义:
β
=
v
c
\\beta=\\frac{v}{c}
β=cv
那么远离情况下,相对论 Doppler 因子为:
D
=
1
−
β
1
+
β
D=\\sqrt{\\frac{1-\\beta}{1+\\beta}}
D=1+β1−β
因此:
f
o
b
s
=
D
f
0
f_{\\mathrm{obs}}=Df_0
fobs=Df0
由于:
D
<
1
D<1
D<1
所以得到红移。
如果光源朝观察者运动,则:
v
<
0
v<0
v<0
代入同一个公式自然得到 D>1,因此:
f
o
b
s
>
f
0
f_{\\mathrm{obs}}>f_0
fobs>f0
得到蓝移。
在工程仿真中,最推荐直接使用有符号速度
β
\\beta
β,然后统一使用一个公式。这样比程序中到处写“如果速度大于零就减、否则就加”更不容易出错。
第二章 低速多普勒为什么看起来很简单:工程模型与误差来源
2.1 低速情况下最常用的线性公式
当:
β
=
v
c
≪
1
\\beta=\\frac{v}{c}\\ll1
β=cv≪1
时,可以展开:
1
−
β
1
+
β
\\sqrt{\\frac{1-\\beta}{1+\\beta}}
1+β1−β
得到:
D
≈
1
−
β
+
β
2
2
−
β
3
2
+
⋯
D\\approx1-\\beta+\\frac{\\beta^2}{2}-\\frac{\\beta^3}{2}+\\cdots
D≈1−β+2β2−2β3+⋯
因此第一阶近似就是:
D
≈
1
−
β
D\\approx1-\\beta
D≈1−β
于是:
f
o
b
s
≈
f
0
(
1
−
β
)
f_{\\mathrm{obs}}\\approx f_0(1-\\beta)
fobs≈f0(1−β)
也可以写成:
Δ
f
=
f
o
b
s
−
f
0
\\Delta f=f_{\\mathrm{obs}}-f_0
Δf=fobs−f0
所以:
Δ
f
≈
−
f
0
v
c
\\Delta f\\approx-f_0\\frac{v}{c}
Δf≈−f0cv
这就是激光测量中最常出现的 Doppler 公式。
这里有一个非常重要的工程现象:
虽然速度
v
/
c
v/c
v/c 非常小,但光学载频
f
0
f_0
f0 极高,所以最后得到的频移仍然可以非常明显。
例如对于 1550 nm:
f
0
≈
193.5
T
H
z
f_0\\approx193.5\\ \\mathrm{THz}
f0≈193.5 THz
如果速度是:
v
=
1
m
/
s
v=1\\ \\mathrm{m/s}
v=1 m/s
那么:
v
c
≈
3.34
×
10
−
9
\\frac{v}{c}\\approx3.34\\times10^{-9}
cv≈3.34×10−9
看起来非常小。
但因为:
193.5
T
H
z
×
3.34
×
10
−
9
193.5\\ \\mathrm{THz}\\times3.34\\times10^{-9}
193.5 THz×3.34×10−9
最后仍然能得到约数百 kHz 量级的频移。
这就是光学 Doppler 测量非常灵敏的原因之一。
2.2 为什么低速运动特别适合连续波激光测量
对于连续波激光,原始光场频率极高,但是工程系统往往并不直接测量 193 THz。
系统通常把被测激光和本地参考激光进行光学外差:
E
1
(
t
)
=
A
1
cos
(
2
π
f
1
t
)
E_1(t)=A_1\\cos(2\\pi f_1t)
E1(t)=A1cos(2πf1t)
E
2
(
t
)
=
A
2
cos
(
2
π
f
2
t
)
E_2(t)=A_2\\cos(2\\pi f_2t)
E2(t)=A2cos(2πf2t)
两者在探测器中混合后,会出现:
f
b
e
a
t
=
∣
f
1
−
f
2
∣
f_{\\mathrm{beat}}=|f_1-f_2|
fbeat=∣f1−f2∣
如果两个光源频率非常接近,那么原本几百 THz 的变化会被搬移到 MHz、GHz 甚至更低的电子频率范围。
这就是为什么光学 Doppler 测量能够被普通高速电子学系统处理。
在实际系统中,多普勒频移本身通常不是直接“看到的 193 THz 搬动”,而是通过光学频率差 转化成可电子处理的拍频这一步对后面的光频梳系统也极其重要。
2.3 低速时最容易犯的一个错误:把“频率变化”误认为“固定频移”
假设光源速度恒定:
v
(
t
)
=
v
0
v(t)=v_0
v(t)=v0
那么:
Δ
f
D
=
c
o
n
s
t
a
n
t
\\Delta f_D=\\mathrm{constant}
ΔfD=constant
这时候 Doppler 只是一个固定频偏。
但很多实际运动目标是:
v
=
v
(
t
)
v=v(t)
v=v(t)
例如:
v
(
t
)
=
a
t
v(t)=at
v(t)=at
那么:
Δ
f
D
(
t
)
=
−
f
0
a
t
c
\\Delta f_D(t)=-f_0\\frac{a t}{c}
ΔfD(t)=−f0cat
Doppler 频率本身就随时间变化。
这时候系统看到的就不再是一个固定频率,而是:
f
D
=
f
D
(
t
)
f_D=f_D(t)
fD=fD(t)
也就是一个chirp。
因此:
恒速运动 → 固定 Doppler 偏移;
加速运动 → Doppler 随时间变化;
非线性运动 → Doppler 曲线本身也呈非线性。
如果把这种信号直接送进 FFT,那么得到的频谱通常不是一根漂亮的窄线,而会出现扩展、展宽甚至明显的频率轨迹。
因此对于运动光源,STFT 往往比单次 FFT 更有信息量。
第三章 高速运动进入相对论区间:时间膨胀和传播时延到底谁在起作用
3.1 相对论 Doppler 不是简单的“加一个时间膨胀”
这是理解高速光源运动最重要的一部分。
很多资料会说:
高速运动时,因为时间膨胀,所以频率变低。
这句话只对了一半。
如果仅仅考虑运动光源自身的时间膨胀,那么光源固有时间
d
τ
d\\tau
dτ和实验室坐标时间
d
t
e
dt_e
dte满足:
d
τ
=
d
t
e
γ
d\\tau=\\frac{dt_e}{\\gamma}
dτ=γdte
其中:
γ
=
1
1
−
β
2
\\gamma=\\frac{1}{\\sqrt{1-\\beta^2}}
γ=1−β2
1
因此在实验室中看来,光源内部的振荡已经变慢。
如果光源自身频率是
f
0
f_0
f0,那么仅从时间膨胀角度看,已经有:
f
l
a
b
=
f
0
γ
f_{\\mathrm{lab}}=\\frac{f_0}{\\gamma}
flab=γf0
但是这还不是最终观测到的 Doppler 频率。
因为光是传播过来的,而且光源正在运动。
3.2 传播时延才是第二个不能忽略的因素
设光源在发射时刻
t
e
t_e
te 距离观察者:
R
(
t
e
)
R(t_e)
R(te)
一束光从光源发出后,到达观察者的时间为:
t
o
b
s
=
t
e
+
R
(
t
e
)
c
t_{\\mathrm{obs}}=t_e+\\frac{R(t_e)}{c}
tobs=te+cR(te)
如果光源正在远离:
d
R
d
t
e
=
v
\\frac{dR}{dt_e}=v
dtedR=v
于是:
d
t
o
b
s
d
t
e
=
1
+
v
c
1
+
β
\\frac{dt_{\\mathrm{obs}}}{dt_e}= 1+\\frac{v}{c} 1+\\beta
dtedtobs=1+cv1+β
这个因素意味着:
光源发射事件之间的时间间隔,在传播到观察者之后还会被额外拉长。
因此最终观察者看到的频率不能只写成:
f
0
/
γ
f_0/\\gamma
f0/γ
而应该进一步除以传播造成的时间拉伸:
f
o
b
s
=
f
0
γ
(
1
+
β
)
f_{\\mathrm{obs}}=\\frac{f_0}{\\gamma(1+\\beta)}
fobs=γ(1+β)f0
于是:
f
o
b
s
=
f
0
1
−
β
1
+
β
f_{\\mathrm{obs}}=f_0\\sqrt{\\frac{1-\\beta}{1+\\beta}}
fobs=f01+β1−β
这就是严格的一维径向相对论 Doppler 公式。
所以可以非常清楚地把它拆成:
相对论 Doppler
=
时间膨胀
+
传播时延变化
\\boxed{\\text{相对论 Doppler}= \\text{时间膨胀} + \\text{传播时延变化} }
相对论 Doppler=时间膨胀+传播时延变化
这是很多简单教材容易略过,但在高速运动仿真中必须明确的地方。
3.3 为什么高速下经典公式会明显失效
经典近似是:
f
c
l
a
s
s
i
c
a
l
=
f
0
(
1
−
β
)
f_{\\mathrm{classical}} =f_0(1-\\beta)
fclassical=f0(1−β)
严格公式是:
f
r
e
l
=
f
0
1
−
β
1
+
β
f_{\\mathrm{rel}} =f_0\\sqrt{\\frac{1-\\beta}{1+\\beta}}
frel=f01+β1−β
当:
β
→
0
\\beta\\rightarrow0
β→0
二者几乎完全一致。
但是随着速度越来越接近光速,两者开始快速分离。
例如:
β
=
0.1
\\beta=0.1
β=0.1
相对论 Doppler 因子约为:
D
≈
0.9045
D\\approx0.9045
D≈0.9045
而经典公式给出:
D
c
l
a
s
s
i
c
a
l
=
0.9
D_{\\mathrm{classical}}=0.9
Dclassical=0.9
此时已经出现可见差异。
到了:
β
=
0.5
\\beta=0.5
β=0.5
相对论结果:
D
=
0.5
1.5
≈
0.577
D=\\sqrt{\\frac{0.5}{1.5}} \\approx0.577
D=1.50.5
≈0.577
经典公式却只有:
D
=
0.5
D=0.5
D=0.5
差别已经不能忽略。
再到:
β
=
0.8
\\beta=0.8
β=0.8
严格结果:
D
=
0.2
1.8
=
1
3
D=\\sqrt{\\frac{0.2}{1.8}} =\\frac13
D=1.80.2
=31
而经典近似仍然会给出:
D
=
0.2
D=0.2
D=0.2
两者已经完全不是一个量级的模型了。
所以工程上可以这样理解:
m/s、km/s、普通车辆速度、机械运动速度,通常完全可以使用线性 Doppler 模型。
当研究问题明确进入高能粒子、相对论运动或者超高速航天器场景时,必须切换到相对论表达式。
第四章 光频梳遇到多普勒:为什么问题突然从“一个频率”变成“整个频谱”
4.1 光频梳不是一根谱线,而是一排高度规则的谱线
连续波激光可以近似写成:
f
0
f_0
f0
一个频率。
而光频梳具有:
f
n
=
n
f
r
e
p
+
f
c
e
o
f_n=nf_{\\mathrm{rep}}+f_{\\mathrm{ceo}}
fn=nfrep+fceo
其中: n是梳齿编号,
f
r
e
p
f_{\\mathrm{rep}}
frep是重复频率,
f
c
e
o
f_{\\mathrm{ceo}}
fceo是载波包络偏移频率。
因此光频梳的频谱并不是“一条线”,而是一列规律分布的纵模。
例如:
f
−
2
,
f
−
1
,
f
0
,
f
1
,
f
2
,
…
f_{-2},f_{-1},f_0,f_1,f_2,\\ldots
f−2,f−1,f0,f1,f2,…
相邻梳齿之间间隔基本为:
f
r
e
p
f_{\\mathrm{rep}}
frep
如果
f
r
e
p
=
200
M
H
z
f_{\\mathrm{rep}}=200\\ \\mathrm{MHz}
frep=200 MHz,那么光频梳可以理解成在几百 THz 光学载频附近,每隔 200 MHz 放置一根非常稳定的光谱线。
因此,一旦整个光源发生 Doppler 运动,问题就变成:
这一整排梳齿应该怎么移动?
4.2 理想情况下,多普勒会让整个光频梳一起“缩放”
如果光频梳作为一个整体运动,并且不考虑额外的色散、介质效应以及光梳自身结构变化,那么每一个纵模都会乘上同一个 Doppler 因子:
f
n
′
=
D
f
n
f_n'=D f_n
fn′=Dfn,即:
f
n
′
=
D
(
n
f
r
e
p
+
f
c
e
o
)
f_n'=D(nf_{\\mathrm{rep}}+f_{\\mathrm{ceo}})
fn′=D(nfrep+fceo),进一步整理:
f
n
′
=
n
(
D
f
r
e
p
)
+
D
f
c
e
o
f_n'=n(Df_{\\mathrm{rep}}) + Df_{\\mathrm{ceo}}
fn′=n(Dfrep)+Dfceo
因此:
f
r
e
p
′
=
D
f
r
e
p
f
c
e
o
′
=
D
f
c
e
o
\\boxed{f_{\\mathrm{rep}}'=Df_{\\mathrm{rep}}} \\boxed{f_{\\mathrm{ceo}}'=Df_{\\mathrm{ceo}}}
frep′=Dfrepfceo′=Dfceo
这意味着:
Doppler 不是简单地给整套光频梳加一个固定频率偏移,而是把整个频谱按比例缩放。
这就是光频梳 Doppler 与连续波激光 Doppler 最大的差别之一。
4.3 为什么不同纵模看起来会产生不同的绝对频移
假设某两个梳齿
f
1
f_1
f1和
f
2
f_2
f2都受到同一个 Doppler 因子
D
D
D,那么:
Δ
f
1
=
f
1
(
D
−
1
)
Δ
f
2
=
f
2
(
D
−
1
)
\\Delta f_1=f_1(D-1) \\\\ \\Delta f_2=f_2(D-1)
Δf1=f1(D−1)Δf2=f2(D−1)
因为:
f
1
≠
f
2
f_1\\neq f_2
f1=f2
所以:
Δ
f
1
≠
Δ
f
2
\\Delta f_1\\neq\\Delta f_2
Δf1=Δf2
也就是说:不同纵模的绝对 Doppler 频移不同
但是如果看相对变化:
Δ
f
n
f
n
=
D
−
1
\\frac{\\Delta f_n}{f_n}=D-1
fnΔfn=D−1
则所有纵模相同。
所以更准确的说法不是:
“不同纵模受到不同的相对论 Doppler 因子。”
而应该说:
“在理想共动情况下,不同纵模具有相同的相对频率缩放,但因为各自的绝对光学频率不同,所以绝对频移量不同。”
这个区别在光频梳 Doppler 仿真里非常重要。
4.4 光频梳的重复频率也会动,这意味着什么
这是很多只研究连续波激光的人容易忽略的一点。
因为:
f
r
e
p
′
=
D
f
r
e
p
f_{\\mathrm{rep}}'=Df_{\\mathrm{rep}}
frep′=Dfrep
所以重复频率也发生变化。
如果原来的:
f
r
e
p
=
200
M
H
z
f_{\\mathrm{rep}}=200\\ \\mathrm{MHz}
frep=200 MHz
那么运动后:
f
r
e
p
′
<
200
M
H
z
f_{\\mathrm{rep}}'<200\\ \\mathrm{MHz}
frep′<200 MHz
对于远离观察者的运动。
对应到时域,就是:
T
r
e
p
=
1
f
r
e
p
T_{\\mathrm{rep}}=\\frac1{f_{\\mathrm{rep}}}
Trep=frep1
变成:
T
r
e
p
′
=
1
D
f
r
e
p
T_{\\mathrm{rep}}'=\\frac1{Df_{\\mathrm{rep}}}
Trep′=Dfrep1
也就是:
T
r
e
p
′
=
T
r
e
p
D
T_{\\mathrm{rep}}'=\\frac{T_{\\mathrm{rep}}}{D}
Trep′=DTrep
由于远离时:
D
<
1
D<1
D<1
于是:
T
r
e
p
′
>
T
r
e
p
T_{\\mathrm{rep}}'>T_{\\mathrm{rep}}
Trep′>Trep
这意味着:
在频域里看到的是频梳整体向低频侧缩放,在时域里看到的则是光脉冲之间的间隔变长。
所以频域和时域看到的是同一个 Doppler 现象的两种表现。
第五章 光频梳高速运动的真正难点:时间、频率和脉冲相位会一起变化
5.1 光频梳的 Doppler 不只有“梳齿移动”
光频梳的脉冲序列可以写成:
E
(
t
)
=
∑
k
A
(
t
−
k
T
r
e
p
)
E(t)=\\sum_k A(t-kT_{\\mathrm{rep}})
E(t)=k∑A(t−kTrep)
这里每一个k代表一个脉冲。
当光源运动后,观察者看到的:
T
r
e
p
T_{\\mathrm{rep}}
Trep
也在变化。
所以:
k
T
r
e
p
kT_{\\mathrm{rep}}
kTrep
不再是一个固定线性时间轴。
如果速度是:
v
=
v
(
t
)
v=v(t)
v=v(t)
那么 Doppler 因子也是:
D
=
D
(
t
)
D=D(t)
D=D(t)
于是:
f
r
e
p
=
f
r
e
p
(
t
)
f_{\\mathrm{rep}}=f_{\\mathrm{rep}}(t)
frep=frep(t)
最终:
T
r
e
p
=
T
r
e
p
(
t
)
T_{\\mathrm{rep}}=T_{\\mathrm{rep}}(t)
Trep=Trep(t)
这意味着整个脉冲列发生动态伸缩。
所以对于光频梳:
Doppler
⇒
频梳伸缩
⇒
脉冲间隔变化
\\boxed{ \\text{Doppler} \\Rightarrow \\text{频梳伸缩} \\Rightarrow \\text{脉冲间隔变化} }
Doppler⇒频梳伸缩⇒脉冲间隔变化
5.2 加速运动为什么会产生光学 chirp
假设一个光源的速度是v(t),那么:
β
(
t
)
=
v
(
t
)
c
\\beta(t)=\\frac{v(t)}{c}
β(t)=cv(t)
相对论 Doppler 因子为:
D
(
t
)
=
1
−
β
(
t
)
1
+
β
(
t
)
D(t)=\\sqrt{\\frac{1-\\beta(t)}{1+\\beta(t)}}
D(t)=1+β(t)1−β(t)
连续波光源观察到:
f
o
b
s
(
t
)
=
f
0
D
(
t
)
f_{\\mathrm{obs}}(t) =f_0D(t)
fobs(t)=f0D(t)
所以如果:
d
v
d
t
≠
0
\\frac{dv}{dt}\\neq0
dtdv=0
就有:
d
f
o
b
s
d
t
≠
0
\\frac{df_{\\mathrm{obs}}}{dt}\\neq0
dtdfobs=0
于是收到的光信号具有一个随时间变化的瞬时频率。
对应相位:
ϕ
(
t
)
=
2
π
∫
0
t
f
o
b
s
(
t
′
)
d
t
′
\\phi(t)=2\\pi\\int_0^t f_{\\mathrm{obs}}(t')dt'
ϕ(t)=2π∫0tfobs(t′)dt′
而相对于本地参考:
ϕ
D
(
t
)
=
2
π
∫
0
t
Δ
f
D
(
t
′
)
d
t
′
\\phi_D(t)=2\\pi\\int_0^t\\Delta f_D(t')dt'
ϕD(t)=2π∫0tΔfD(t′)dt′
所以运动产生的 Doppler 并不只是一个频率变化问题,它同时对应一个相位积分问题。
这就是为什么在相干光学系统里,即使频率误差非常小,也可能通过长时间积分积累出巨大的相位误差。
5.3 光频梳中不同纵模会出现不同的 Doppler 相位斜率
第
n
n
n 根梳齿:
f
n
′
=
D
f
n
f_n'=Df_n
fn′=Dfn
其相位:
ϕ
n
(
t
)
=
2
π
∫
0
t
D
f
n
(
t
′
)
d
t
′
\\phi_n(t)=2\\pi\\int_0^tDf_n(t')dt'
ϕn(t)=2π∫0tDfn(t′)dt′
所以:
ϕ
n
(
t
)
=
2
π
f
n
∫
0
t
D
(
t
′
)
d
t
′
\\phi_n(t) =2\\pi f_n\\int_0^tD(t')dt'
ϕn(t)=2πfn∫0tD(t′)dt′
因为:
f
n
f_n
fn
不同,所以不同梳齿的相位累积速率也不同。
这就是为什么高速动态 Doppler 对光频梳的影响不能仅仅用“统一加一个频偏”模拟。
特别是在双光梳系统中,如果不同梳齿分别参与外差,最终 RF 域的拍频会出现与模式编号相关的变化。
这类现象进一步就会进入:
Doppler–delay coupling
\\boxed{\\text{Doppler–delay coupling}}
Doppler–delay coupling
和:
Doppler–phase coupling
\\boxed{\\text{Doppler–phase coupling}}
Doppler–phase coupling
第六章 双光梳中的 Doppler:为什么它比连续激光复杂很多
6.1 双光梳系统为什么要把光学频率搬到 RF
传统双光梳系统通常有两套光频梳:
f
r
e
p
,
1
f_{\\mathrm{rep},1}
frep,1
和:
f
r
e
p
,
2
f_{\\mathrm{rep},2}
frep,2
两者略有不同:
Δ
f
r
e
p
=
f
r
e
p
,
2
f
r
e
p
,
1
\\Delta f_{\\mathrm{rep}} = f_{\\mathrm{rep},2} f_{\\mathrm{rep},1}
Δfrep=frep,2frep,1
于是光学域中相距几十 MHz、几百 MHz 的梳齿,经过双梳外差之后,会映射到 RF 域。
这相当于:
optical spectrum
→
RF spectrum
\\text{optical spectrum} \\rightarrow \\text{RF spectrum}
optical spectrum→RF spectrum
因此双梳系统最大的优势之一,就是把难以直接电子采集的光学结构转化成能够处理的射频结构。
但是 Doppler 一旦进入其中,就会改变这张“映射关系”。
6.2 运动光频梳和静止光频梳做拍频会发生什么
假设参考梳:
f
n
(
A
)
=
n
f
r
e
p
,
A
+
f
c
e
o
,
A
f_n^{(A)} =nf_{\\mathrm{rep},A}+f_{\\mathrm{ceo},A}
fn(A)=nfrep,A+fceo,A
运动梳:
f
n
(
B
)
=
n
f
r
e
p
,
B
+
f
c
e
o
,
B
f_n^{(B)}=nf_{\\mathrm{rep},B}+f_{\\mathrm{ceo},B}
fn(B)=nfrep,B+fceo,B
运动以后:
f
n
(
B
)
′
=
D
f
n
(
B
)
f_n^{(B)\\prime} =D f_n^{(B)}
fn(B)′=Dfn(B)
那么第
n
n
n 根梳齿对应的拍频为:
f
b
e
a
t
,
n
=
∣
D
(
n
f
r
e
p
,
B
+
f
c
e
o
,
B
)
−
(
n
f
r
e
p
,
A
+
f
c
e
o
,
A
)
∣
f_{\\mathrm{beat},n}=\\left|D\\left(nf_{\\mathrm{rep},B}+f_{\\mathrm{ceo},B}\\right)- \\left(nf_{\\mathrm{rep},A}+f_{\\mathrm{ceo},A}\\right)\\right|
fbeat,n=∣D(nfrep,B+fceo,B)−(nfrep,A+fceo,A)∣
整理之后:
f
b
e
a
t
,
n
=
∣
n
(
D
f
r
e
p
,
B
−
f
r
e
p
,
A
)
+
(
D
f
c
e
o
,
B
−
f
c
e
o
,
A
)
∣
f_{\\mathrm{beat},n} =\\left| n \\left( Df_{\\mathrm{rep},B} -f_{\\mathrm{rep},A} \\right) + \\left( Df_{\\mathrm{ceo},B} -f_{\\mathrm{ceo},A} \\right) \\right|
fbeat,n=∣n(Dfrep,B−frep,A)+(Dfceo,B−fceo,A)∣
这个公式非常值得记住。
因为它直接告诉我们:
Doppler 不仅会改变某一根 RF 拍频,它还会改变整个 RF comb 的斜率。
也就是说:
f
r
e
p
f_{\\mathrm{rep}}
frep
变化会直接反映到 RF 模式间隔。
而:
f
c
e
o
f_{\\mathrm{ceo}}
fceo
的变化会反映到整体 RF offset。
6.3 这也是为什么光频梳可以拿来测速
连续波激光利用:
Δ
f
D
\\Delta f_D
ΔfD
测速。
光频梳则可以利用:
Δ
f
r
e
p
\\Delta f_{\\mathrm{rep}}
Δfrep
和:
Δ
f
c
e
o
\\Delta f_{\\mathrm{ceo}}
Δfceo
以及不同梳齿对应的拍频变化共同估计运动状态。
对于双梳而言,多根梳齿相当于同时提供了很多“观察通道”。
因此从信息量来看:
单频激光
<
单组光频梳
<
双光频梳
\\text{单频激光} < \\text{单组光频梳} < \\text{双光频梳}
单频激光<单组光频梳<双光频梳
但是与此同时,模型复杂度也是不断增加的。
第七章 异步采样和运动 Doppler:为什么速度会干扰时间测量
7.1 异步光学采样本来就是一个时间扫描器
传统双光梳测量常利用:
Δ
f
r
e
p
\\Delta f_{\\mathrm{rep}}
Δfrep
让两个脉冲列逐渐错开。
它本质上是在做一个“慢速时间扫描”。
可以把两组脉冲想象成:
Comb A:
∣
−
−
−
∣
−
−
−
∣
−
−
−
∣
−
−
−
∣
−
−
−
∣
\\text{Comb A: } |—|—|—|—|—|
Comb A: ∣−−−∣−−−∣−−−∣−−−∣−−−∣
Comb B:
∣
−
−
−
∣
−
−
−
∣
−
−
−
∣
−
−
−
∣
−
−
−
∣
\\text{Comb B: } \\ \\ |—|—|—|—|—|
Comb B: ∣−−−∣−−−∣−−−∣−−−∣−−−∣
随着时间推移,两者相对位置不断变化。
这就是异步采样的时间游标。
但是如果其中一组光频梳本身正在高速运动,那么其:
f
r
e
p
f_{\\mathrm{rep}}
frep
也发生变化。
这意味着本来固定的:
Δ
f
r
e
p
\\Delta f_{\\mathrm{rep}}
Δfrep
实际上变成:
Δ
f
r
e
p
(
t
)
\\Delta f_{\\mathrm{rep}}(t)
Δfrep(t)
于是时间扫描速度也开始变化。
7.2 多普勒会同时进入“频率轴”和“时间轴”
这是双梳系统里最值得注意的一件事情。
Doppler 改变:
f
r
e
p
f_{\\mathrm{rep}}
frep
而:
T
r
e
p
=
1
f
r
e
p
T_{\\mathrm{rep}}=\\frac1{f_{\\mathrm{rep}}}
Trep=frep1
因此它必然也改变:
pulse arrival time
\\text{pulse arrival time}
pulse arrival time
也就是说:
Doppler
→
f
r
e
p
→
T
r
e
p
→
time position
\\boxed{ \\text{Doppler} \\rightarrow f_{\\mathrm{rep}} \\rightarrow T_{\\mathrm{rep}} \\rightarrow \\text{time position} }
Doppler→frep→Trep→time position
所以如果你只看频谱,很容易觉得 Doppler 就是“频率偏移”。
但对于双光梳时间测量而言,它同时会变成一个:
timing error
\\boxed{\\text{timing error}}
timing error
7.3 这就是时延-多普勒耦合为什么会出现
假设目标距离:
R
(
t
)
R(t)
R(t)
那么传播时延:
τ
(
t
)
=
R
(
t
)
c
\\tau(t)=\\frac{R(t)}{c}
τ(t)=cR(t)
如果目标运动:
R
(
t
)
=
R
0
+
v
t
R(t)=R_0+vt
R(t)=R0+vt
那么:
τ
(
t
)
=
R
0
c
+
v
t
c
\\tau(t)=\\frac{R_0}{c}+\\frac{vt}{c}
τ(t)=cR0+cvt
所以时延本身随时间变化。
与此同时,光学频率也变化:
f
D
(
t
)
∼
f
0
v
c
f_D(t) \\sim f_0\\frac{v}{c}
fD(t)∼f0cv
于是系统同时观察到:
τ
(
t
)
\\tau(t)
τ(t)
和:
f
D
(
t
)
f_D(t)
fD(t)
这两个量彼此相关。
对于高速动态目标,如果仍然把它们当作两个完全独立的变量,就很容易产生参数估计偏差。
因此实际高精度双梳测量经常需要:
delay-Doppler joint estimation
\\boxed{\\text{delay-Doppler joint estimation}}
delay-Doppler joint estimation
或者:
delay-Doppler decoupling
\\boxed{\\text{delay-Doppler decoupling}}
delay-Doppler decoupling
第八章 TPFC 为什么在运动目标场景有意思:从“扫描”变成“跟踪”
8.1 普通双梳和 TPFC 的思路差异
你提供的 TPFC 材料对传统双光梳系统有一个非常清楚的总结:传统方案依赖两套光频梳重复频率的刻意失配形成游标式扫描,这会带来等效探测死区,并使测量速度、灵敏度和分辨率之间出现权衡。
TPFC 则增加了一个非常重要的自由度:
X
(
t
)
X(t)
X(t)
也就是直接控制脉冲时序的位置。
因此它不必始终依靠:
Δ
f
r
e
p
\\Delta f_{\\mathrm{rep}}
Δfrep
去慢慢扫描整个时间空间。
这就像普通双梳是拿着尺子不断来回扫:
scan
→
scan
→
scan
\\text{scan}\\rightarrow\\text{scan}\\rightarrow\\text{scan}
scan→scan→scan
而 TPFC 更接近:
find
→
lock
→
follow
\\text{find}\\rightarrow\\text{lock}\\rightarrow\\text{follow}
find→lock→follow
这种方式尤其适合动态目标。
8.2 TPFC 锁定目标后,不意味着 Doppler 消失
这一点一定不能误解。
如果目标在运动:
v
≠
0
v\\neq0
v=0
那么远端光脉冲仍然存在:
f
D
≠
0
f_D\\neq0
fD=0
TPFC 锁定做的是:
X
T
P
F
C
(
t
)
→
X
s
i
g
n
a
l
(
t
)
X_{\\mathrm{TPFC}}(t) \\rightarrow X_{\\mathrm{signal}}(t)
XTPFC(t)→Xsignal(t)
也就是让本地可控脉冲持续跟踪远端信号的到达位置。
因此:
TPFC 是在跟踪 Doppler 导致的时序变化,而不是把 Doppler 从物理世界里消掉。
这一点非常重要。
换句话说:
锁定
≠
没有 Doppler
\\boxed{ \\text{锁定} \\neq \\text{没有 Doppler} }
锁定=没有 Doppler
而是:
锁定
=
让本地系统持续跟着 Doppler 变化走
\\boxed{ \\text{锁定} =\\text{让本地系统持续跟着 Doppler 变化走} }
锁定=让本地系统持续跟着 Doppler 变化走
8.3 TPFC 的时序反馈量本身就带着运动信息
你提供的 TPFC 材料给出了一个非常直接的系统实现:跟踪模式下,时序鉴别器输出时序误差,反馈到 PII 控制器形成 TPFC 时序轨迹;同时载波环路修正解调数控振荡器以补偿回波 Doppler。也就是说,时序误差和载波频率误差是同时被测量和控制的。
这其实非常适合运动目标。
因为:
X
(
t
)
X(t)
X(t)
告诉你目标在时间上的位置变化。
而:
θ
(
t
)
\\theta(t)
θ(t)
的变化又携带 Doppler 和运动信息。
论文材料还明确给出了:
Δ
X
=
Δ
θ
0
C
2
π
N
f
r
e
p
\\Delta X =\\frac{\\Delta\\theta_0^C} {2\\pi Nf_{\\mathrm{rep}}}
ΔX=2πNfrepΔθ0C
并指出在跟踪模式中保持指定梳齿相位并调节时序,可以由有效载波频率提取相对运动信息。
这意味着 TPFC 不只是“一个会移动的本振”,而是可以同时成为:
timing tracker
+
phase tracker
\\boxed{ \\text{timing tracker} + \\text{phase tracker} }
timing tracker+phase tracker
第九章 光源低速、高速、线性和非线性运动应该怎么理解
9.1 低速恒速:最简单的情况
首先:
v
(
t
)
=
v
0
v(t)=v_0
v(t)=v0
那么:
β
=
v
0
c
\\beta=\\frac{v_0}{c}
β=cv0
是常数。
相对论 Doppler:
D
=
1
−
β
1
+
β
D=\\sqrt{\\frac{1-\\beta}{1+\\beta}}
D=1+β1−β
也是常数。
所以:
f
o
b
s
=
D
f
0
f_{\\mathrm{obs}}=Df_0
fobs=Df0
表现为一根固定频移的谱线。
对于光频梳:
f
n
′
=
D
f
n
f_n'=Df_n
fn′=Dfn
整个梳状结构保持规则,只是整体缩放。
这是最好处理的情况。
9.2 低速线性加速:Doppler 变成一条直线
如果:
v
(
t
)
=
v
0
+
a
t
v(t)=v_0+at
v(t)=v0+at
且:
v
≪
c
v\\ll c
v≪c
那么:
Δ
f
D
(
t
)
≈
−
f
0
c
(
v
0
+
a
t
)
\\Delta f_D(t) \\approx -\\frac{f_0}{c}(v_0+at)
ΔfD(t)≈−cf0(v0+at)
因此:
d
Δ
f
D
d
t
≈
−
f
0
a
c
\\frac{d\\Delta f_D}{dt} \\approx -\\frac{f_0a}{c}
dtdΔfD≈−cf0a
Doppler 频率随时间线性变化。
因此 STFT 图上会出现接近直线的频率轨迹。
这种情况就是最典型的 Doppler chirp。
9.3 高速线性加速:速度线性不意味着 Doppler 线性
如果仍然:
v
(
t
)
=
v
0
+
a
t
v(t)=v_0+at
v(t)=v0+at
但是:
v
∼
c
v\\sim c
v∼c
那么:
f
D
(
t
)
=
f
0
1
−
v
(
t
)
/
c
1
+
v
(
t
)
/
c
f_D(t) =f_0 \\sqrt{ \\frac{1-v(t)/c} {1+v(t)/c} }
fD(t)=f01+v(t)/c1−v(t)/c
由于根号中的关系是非线性的,因此:
f
D
(
t
)
f_D(t)
fD(t)
不再是直线。
即使速度是匀加速:
v
(
t
)
=
a
t
v(t)=at
v(t)=at
频率轨迹依然会发生弯曲。
这就是高速情况下非常直观的一个“相对论非线性 Doppler”现象。
9.4 非线性速度:STFT 会出现更加复杂的曲线
如果速度是:
v
(
t
)
=
v
max
sin
(
Ω
t
)
v(t)=v_{\\max}\\sin(\\Omega t)
v(t)=vmaxsin(Ωt)
那么:
f
D
(
t
)
=
f
0
1
−
β
max
sin
(
Ω
t
)
1
+
β
max
sin
(
Ω
t
)
f_D(t) =f_0 \\sqrt{ \\frac{1-\\beta_{\\max}\\sin(\\Omega t)} {1+\\beta_{\\max}\\sin(\\Omega t)} }
fD(t)=f01+βmaxsin(Ωt)1−βmaxsin(Ωt)
这时候频率变化也会呈现明显的非正弦性质。
特别是当:
β
max
\\beta_{\\max}
βmax
比较大时,波形的上下半周期不再具有简单的镜像关系。
因此对高速周期运动,不能简单地说:
“速度是正弦的,所以 Doppler 也是正弦的。”
低速时可以近似这样理解。
高速时则不能。
第十章 相对论时间膨胀对光频梳到底意味着什么
10.1 频率是时间的倒数,所以时间尺度变化必然进入频梳
频率和周期满足:
f
=
1
T
f=\\frac1T
f=T1
如果运动造成:
d
t
→
d
τ
dt\\rightarrow d\\tau
dt→dτ
发生变化,那么所有基于时间尺度定义的量都必须重新解释。
光频梳的重复频率:
f
r
e
p
f_{\\mathrm{rep}}
frep
本质上就是:
1
T
r
e
p
\\frac1{T_{\\mathrm{rep}}}
Trep1
因此时间尺度发生相对论变化以后:
f
r
e
p
f_{\\mathrm{rep}}
frep
自然也发生变化。
所以对于高速运动的光频梳,所谓“时间膨胀”并不是一个可以独立放在系统角落里的修正项。
它最终会进入:
f
r
e
p
f_{\\mathrm{rep}}
frep
再进入:
T
r
e
p
T_{\\mathrm{rep}}
Trep
最后表现为:
pulse timing
\\text{pulse timing}
pulse timing
变化。
10.2 不要把时间膨胀和传播 Doppler 分成两个互不相关的补丁
一个常见但不够严谨的仿真方法是:
先算:
f
1
=
f
0
γ
f_1=\\frac{f_0}{\\gamma}
f1=γf0
然后再人为乘一个:
1
−
v
c
1-\\frac vc
1−cv
这种处理很容易重复计算或者漏项。
更好的方式,是直接从发射时刻与到达时刻的映射出发:
t
o
b
s
=
t
e
+
R
(
t
e
)
c
t_{\\mathrm{obs}}=t_e+\\frac{R(t_e)}{c}
tobs=te+cR(te)
再从源的固有相位:
ϕ
e
=
2
π
f
0
τ
\\phi_e=2\\pi f_0\\tau
ϕe=2πf0τ
出发。
因为:
d
τ
=
d
t
e
γ
d\\tau=\\frac{dt_e}{\\gamma}
dτ=γdte
最终得到:
f
o
b
s
=
1
2
π
d
ϕ
e
d
t
o
b
s
f_{\\mathrm{obs}} =\\frac1{2\\pi} \\frac{d\\phi_e}{dt_{\\mathrm{obs}}}
fobs=2π1dtobsdϕe
自然就是:
f
o
b
s
=
f
0
γ
(
1
+
β
)
f_{\\mathrm{obs}} =\\frac{f_0} {\\gamma(1+\\beta)}
fobs=γ(1+β)f0
对于远离观察者的情况。
这种方法的好处是:
时间膨胀和传播时延并不是两个“拼接项”,而是同一个时空映射中的两个部分。
第十一章 连续激光和光频梳的 Doppler,到底应该怎么对比
| 连续波激光 | 光学频率整体偏移 | 相位持续积累 |
| 脉冲激光 | 脉冲到达时间改变 | 脉冲间隔变化 |
| 光频梳 | 整个梳谱按 Doppler 因子缩放 |
f r e p f_{\\mathrm{rep}} frep、 f c e o f_{\\mathrm{ceo}} fceo、脉冲时序同时变化 |
| 双光梳 | RF 拍频结构变化 | 时延-多普勒耦合 |
| TPFC 跟踪系统 | 时序控制量持续变化 | 可从时序和载波相位同时提取运动信息 |
这个表背后其实是一条很清楚的物理链:
连续波
→
频率
\\boxed{ \\text{连续波} \\rightarrow \\text{频率} }
连续波→频率
脉冲列
→
时间间隔
\\boxed{ \\text{脉冲列} \\rightarrow \\text{时间间隔} }
脉冲列→时间间隔
光频梳
→
频率+时间
\\boxed{ \\text{光频梳} \\rightarrow \\text{频率+时间} }
光频梳→频率+时间
双光梳
→
频率+时间+RF映射
\\boxed{ \\text{双光梳} \\rightarrow \\text{频率+时间+RF映射} }
双光梳→频率+时间+RF映射
TPFC
→
频率+时间+闭环跟踪
\\boxed{ \\text{TPFC} \\rightarrow \\text{频率+时间+闭环跟踪} }
TPFC→频率+时间+闭环跟踪
因此系统越往后,Doppler 的表现形式就越丰富。
第十二章 一个特别容易踩的坑:光频梳“整体平移”和“整体伸缩”不是一回事
12.1 固定频移和比例缩放有什么区别
假设原来的两个频率是:
f
1
=
100
f_1=100
f1=100
f
2
=
110
f_2=110
f2=110
如果给它们统一加:
+
5
+5
+5
那么得到:
105
,
115
105,\\quad115
105,115
两者间隔仍然是:
10
10
10
这是整体平移。
但如果进行比例缩放:
f
′
=
1.01
f
f'=1.01f
f′=1.01f
那么:
101
,
111.1
101,\\quad111.1
101,111.1
间隔变成:
10.1
10.1
10.1
这才叫整体伸缩。
而理想径向 Doppler 对光频梳更接近后者:
f
n
′
=
D
f
n
f_n'=Df_n
fn′=Dfn
因此:
f
r
e
p
′
=
D
f
r
e
p
f_{\\mathrm{rep}}'=Df_{\\mathrm{rep}}
frep′=Dfrep
12.2 为什么这个区别对双光梳特别重要
如果你错误地把 Doppler 写成:
f
n
′
=
f
n
+
Δ
f
D
f_n'=f_n+\\Delta f_D
fn′=fn+ΔfD
那么:
f
n
+
1
′
−
f
n
′
=
f
n
+
1
−
f
n
f_{n+1}'-f_n' =f_{n+1}-f_n
fn+1′−fn′=fn+1−fn
也就是说:
f
r
e
p
f_{\\mathrm{rep}}
frep
完全不变。
这实际上隐含了一个物理假设:
所有纵模整体刚性平移。
然而真正的径向 Doppler 关系是:
f
n
′
=
D
f
n
f_n'=Df_n
fn′=Dfn
于是:
f
n
+
1
′
−
f
n
′
=
D
(
f
n
+
1
−
f
n
)
f_{n+1}'-f_n' =D(f_{n+1}-f_n)
fn+1′−fn′=D(fn+1−fn)
得到:
f
r
e
p
′
=
D
f
r
e
p
f_{\\mathrm{rep}}'=Df_{\\mathrm{rep}}
frep′=Dfrep
因此如果你的目标是研究:
高速运动下光频梳不同纵模的变化
那么必须使用比例缩放模型,而不是简单给每根梳齿添加同一个频移。
第十三章 高速运动下的频谱图应该怎么看
13.1 连续激光的频谱
恒速运动时:
f
o
b
s
=
c
o
n
s
t
a
n
t
f_{\\mathrm{obs}}=\\mathrm{constant}
fobs=constant
所以频谱仍然是一根窄线:
Power
│
│ │
│ │
│ │
──┼─────────────┼──────── Frequency
fD
只不过谱线位置从:
f
0
f_0
f0
移动到:
D
f
0
Df_0
Df0
13.2 加速运动时,单次 FFT 为什么会变宽
如果:
f
D
=
f
D
(
t
)
f_D=f_D(t)
fD=fD(t)
那么整个观测窗口里频率一直在变化。
FFT相当于把整个时间窗口压缩成一个频谱统计结果。
因此一个不断变化的频率轨迹会被摊开成一个频率范围。
所以:
FFT 看到的是“扩展后的频谱”;STFT 看到的是“频率如何随时间运动”。
这就是为什么运动目标仿真中,一般建议:
FFT + STFT
\\boxed{ \\text{FFT + STFT} }
FFT + STFT
一起画。
13.3 光频梳的 STFT 会出现什么
对于频率固定的光频梳,STFT 会出现:
f
1
,
f
2
,
f
3
,
…
f_1,f_2,f_3,\\ldots
f1,f2,f3,…
一组平行水平线。
如果光源在恒速运动,那么所有梳齿都会发生比例变化。
如果速度随时间变化,那么这些线会同步向上或者向下运动。
而且不同梳齿:
f
n
′
=
D
f
n
f_n'=Df_n
fn′=Dfn
所以它们的绝对频率移动量不同。
于是高阶模式的轨迹变化更加明显。
这就是非常典型的:
Doppler-induced comb stretching
\\boxed{\\text{Doppler-induced comb stretching}}
Doppler-induced comb stretching
第十四章 为什么“运动速度越高,光频梳越适合做 Doppler 测量”并不完全正确
14.1 高速度意味着信号变大,也意味着模型变复杂
从最简单的公式:
Δ
f
≈
−
f
0
v
c
\\Delta f\\approx-f_0\\frac vc
Δf≈−f0cv
可以看到:
v
v
v
越大:
∣
Δ
f
∣
|\\Delta f|
∣Δf∣
越大。
这似乎意味着测量越来越容易。
但高速同时会产生:
f
r
e
p
f_{\\mathrm{rep}}
frep
变化、
f
c
e
o
f_{\\mathrm{ceo}}
fceo
变化、
T
r
e
p
T_{\\mathrm{rep}}
Trep
变化、
ϕ
(
t
)
\\phi(t)
ϕ(t)
快速累积、
Δ
f
D
(
t
)
\\Delta f_D(t)
ΔfD(t)
非线性变化。
所以高速并不是简单地“信号更强”。
更准确地说是:
信号更明显,但模型也更复杂
\\boxed{ \\text{信号更明显,但模型也更复杂} }
信号更明显,但模型也更复杂
14.2 高速运动会让“低速近似”越来越危险
如果你的仿真只覆盖:
v
=
1
∼
100
m
/
s
v=1\\sim100\\ \\mathrm{m/s}
v=1∼100 m/s
那么:
β
∼
10
−
7
∼
10
−
6
\\beta\\sim10^{-7}\\sim10^{-6}
β∼10−7∼10−6
相对论修正通常完全不是主导因素。
这时候最重要的反而是:
- 光路长度变化;
- 振动;
- 激光线宽;
- 光梳相位噪声;
- 探测器噪声;
- 环路带宽。
但如果你为了展示相对论现象,把速度提高到:
0.1
c
0.1c
0.1c
甚至:
0.5
c
0.5c
0.5c
那么相对论 Doppler 就会成为不可忽略的物理效应。
因此不要为了“看起来高级”而在一个本来只有几 m/s 的实验里强行加入
0.5
c
0.5c
0.5c。
更合理的研究方法是:
低速模型用于工程真实性,高速模型用于验证相对论算法。
第十五章 对双光梳时频传递而言,真正需要跟踪的不是一个频移
15.1 单纯测量 Doppler 和完成时间同步不是同一个目标
在时频传递中,我们真正关心的是:
Δ
t
A
B
=
t
B
−
t
A
\\Delta t_{\\mathrm{AB}} =t_B-t_A
ΔtAB=tB−tA
而运动目标会引入传播时间:
T
A
B
(
t
)
=
R
(
t
)
c
T_{AB}(t) =\\frac{R(t)}{c}
TAB(t)=cR(t)
如果:
R
(
t
)
=
R
0
+
v
t
R(t)=R_0+vt
R(t)=R0+vt
那么:
T
A
B
(
t
)
=
R
0
c
+
v
t
c
T_{AB}(t) =\\frac{R_0}{c} + \\frac{vt}{c}
TAB(t)=cR0+cvt
与此同时,又会出现:
f
D
(
t
)
f_D(t)
fD(t)
因此实际测量里经常出现:
clock offset
+
propagation delay
+
Doppler
\\boxed{ \\text{clock offset} + \\text{propagation delay} + \\text{Doppler} }
clock offset+propagation delay+Doppler
三者同时变化的情况。
这也是双向光学时频传递的重要原因。
15.2 TPFC 可以把“跟踪运动目标”变成一个闭环问题
你提供的资料中,TPFC被设计成可以动态调节:
X
(
t
)
X(t)
X(t)
和:
θ
(
t
)
\\theta(t)
θ(t)
并通过两套 PLL 对时序和相位实施控制。其系统能够在捕获后进入跟踪模式,让 TPFC 脉冲持续保持与输入脉冲序列的时间重合。
进一步,在跟踪模式中,系统同时使用时序误差信号和载波相位误差信号:时序控制维持脉冲重合,载波环路补偿 Doppler,二者联合用于提取距离和速度。
这类架构的价值,不在于“把 Doppler 消掉”,而在于:
把 Doppler 从一个扰动变成一个可跟踪、可估计的动态量
\\boxed{ \\text{把 Doppler 从一个扰动变成一个可跟踪、可估计的动态量} }
把 Doppler 从一个扰动变成一个可跟踪、可估计的动态量
这对未来的运动目标光学时频传递很有意义。
第十六章 如果用数学模型描述整个光频梳运动过程,可以建立什么框架
对于一个理想光频梳:
f
n
=
n
f
r
e
p
+
f
c
e
o
f_n=nf_{\\mathrm{rep}}+f_{\\mathrm{ceo}}
fn=nfrep+fceo
定义运动:
v
(
t
)
v(t)
v(t)
那么:
β
(
t
)
=
v
(
t
)
c
\\beta(t)=\\frac{v(t)}{c}
β(t)=cv(t)
相对论 Doppler 因子:
D
(
t
)
=
1
−
β
(
t
)
1
+
β
(
t
)
D(t) =\\sqrt{ \\frac{1-\\beta(t)} {1+\\beta(t)} }
D(t)=1+β(t)1−β(t)
于是:
f
n
o
b
s
(
t
)
=
D
(
t
)
[
n
f
r
e
p
+
f
c
e
o
]
f_n^{\\mathrm{obs}}(t) =D(t) \\left[ nf_{\\mathrm{rep}}+f_{\\mathrm{ceo}} \\right]
fnobs(t)=D(t)[nfrep+fceo]
进一步:
f
r
e
p
o
b
s
(
t
)
=
D
(
t
)
f
r
e
p
f_{\\mathrm{rep}}^{\\mathrm{obs}}(t) =D(t)f_{\\mathrm{rep}}
frepobs(t)=D(t)frep
以及:
f
c
e
o
o
b
s
(
t
)
=
D
(
t
)
f
c
e
o
f_{\\mathrm{ceo}}^{\\mathrm{obs}}(t) =D(t)f_{\\mathrm{ceo}}
fceoobs(t)=D(t)fceo
这三个式子实际上已经构成了一个非常完整的理想运动光频梳模型。
如果进一步加入传播距离:
R
(
t
)
R(t)
R(t)
那么到达时间满足:
t
o
b
s
=
t
e
+
R
(
t
e
)
c
t_{\\mathrm{obs}} =t_e+\\frac{R(t_e)}c
tobs=te+cR(te)
如果目标还发生加速度:
a
(
t
)
=
d
v
d
t
a(t)=\\frac{dv}{dt}
a(t)=dtdv
那么:
D
=
D
(
t
)
D=D(t)
D=D(t)
从而:
f
n
=
f
n
(
t
)
f_n=f_n(t)
fn=fn(t)
再继续积分:
ϕ
n
(
t
)
=
2
π
∫
f
n
(
t
′
)
d
t
′
\\phi_n(t) =2\\pi \\int f_n(t')dt'
ϕn(t)=2π∫fn(t′)dt′
于是你就可以从:
v
(
t
)
v(t)
v(t)
一路建立到:
f
n
(
t
)
→
ϕ
n
(
t
)
→
E
n
(
t
)
→
heterodyne
→
f
R
F
(
t
)
f_n(t) \\rightarrow \\phi_n(t) \\rightarrow E_n(t) \\rightarrow \\text{heterodyne} \\rightarrow f_{\\mathrm{RF}}(t)
fn(t)→ϕn(t)→En(t)→heterodyne→fRF(t)
这就是从运动学一直走到光电探测器的完整链路。
第十七章 工程仿真中应该怎样分层,否则很容易把模型写乱
17.1 第一层:运动学模型
先定义:
R
(
t
)
,
v
(
t
)
,
a
(
t
)
R(t),\\quad v(t),\\quad a(t)
R(t),v(t),a(t)
这个阶段只研究目标怎么运动。
例如:
R
(
t
)
=
R
0
+
v
0
t
+
1
2
a
t
2
R(t)=R_0+v_0t+\\frac12at^2
R(t)=R0+v0t+21at2
不要一开始就把 PLL、噪声、FFT、光梳全塞进去。
17.2 第二层:传播模型
由:
R
(
t
)
R(t)
R(t)
得到:
T
(
t
)
=
R
(
t
)
c
T(t)=\\frac{R(t)}c
T(t)=cR(t)
如果考虑相对论,就使用发射时间到接收时间的严格映射。
这样可以得到真实的:
t
o
b
s
t_{\\mathrm{obs}}
tobs
17.3 第三层:光学 Doppler 模型
然后计算:
D
(
t
)
D(t)
D(t)
得到:
f
l
a
s
e
r
(
t
)
f_{\\mathrm{laser}}(t)
flaser(t)
或者:
f
n
c
o
m
b
(
t
)
f_n^{\\mathrm{comb}}(t)
fncomb(t)
这一步是连续激光与光频梳真正产生分叉的地方。
17.4 第四层:相干相位模型
频率是相位的导数:
f
(
t
)
=
1
2
π
d
ϕ
d
t
f(t)=\\frac1{2\\pi}\\frac{d\\phi}{dt}
f(t)=2π1dtdϕ
所以:
ϕ
(
t
)
=
2
π
∫
f
(
t
)
d
t
\\phi(t)=2\\pi\\int f(t)dt
ϕ(t)=2π∫f(t)dt
不要只改变频率而不改变相位。
特别是模拟:
- 外差拍频;
- PLL;
- 光学相干积分;
- 双梳相位;
- Doppler chirp;
时,相位连续性比单纯频率曲线更加重要。
第十八章 哪些情况下可以用简化公式,哪些情况下必须上相对论模型
| 实验室平台机械运动 | 一阶 Doppler 通常足够 |
| 激光雷达常规目标 | 一阶 Doppler 通常足够 |
| m/s 到 km/s 级目标 | 一阶模型通常已经很好 |
| 亚米每秒级高精度测速 | 应关注传播时延、环境和系统噪声 |
| 高速航天器 | 需要评估相对论修正 |
| 粒子束 | 必须使用相对论 Doppler |
|
v ∼ 0.1 c v\\sim0.1c v∼0.1c 以上 |
不建议使用线性近似 |
| 相对论高速光梳仿真 | 直接使用 Doppler factor |
| 研究时间膨胀 | 必须显式建立固有时与接收时间关系 |
这里最容易犯的错误,就是把“物理上正确”与“工程上必要”混为一谈。
相对论模型当然更加完整,但如果你的实际系统速度只有几十 m/s,却把所有噪声和执行器都简化掉,只强调时间膨胀,最后得到的仿真反而没有实际意义。
第十九章 对光频梳来说,真正值得关注的几个 Doppler 观测量
如果以后你研究光频梳运动,不建议只画:
Δ
f
\\Delta f
Δf
一条曲线。
更推荐同时观察:
f
n
(
t
)
f_n(t)
fn(t)
也就是各纵模频率。
再观察:
f
r
e
p
(
t
)
f_{\\mathrm{rep}}(t)
frep(t)
看频梳间隔如何变化。
再观察:
f
c
e
o
(
t
)
f_{\\mathrm{ceo}}(t)
fceo(t)
看载波包络偏移如何变化。
再观察:
T
r
e
p
(
t
)
T_{\\mathrm{rep}}(t)
Trep(t)
看脉冲间隔怎么变。
再观察:
ϕ
n
(
t
)
\\phi_n(t)
ϕn(t)
看不同纵模的相位积累。
最后再观察:
f
b
e
a
t
,
n
(
t
)
f_{\\mathrm{beat},n}(t)
fbeat,n(t)
看这些变化怎样映射到 RF 拍频。
这样画出来的结果才真正能够说明:
Doppler 是如何从运动学一路传到光频梳和电子频域的。
第二十章 一个非常重要的认识:光频梳让 Doppler 从“一个数字”变成“一个结构”
连续激光的 Doppler 测量可能最终告诉你:
Δ
f
=
2.3
M
H
z
\\Delta f=2.3\\ \\mathrm{MHz}
Δf=2.3 MHz
这是一个数字。
但光频梳能够告诉你:
Δ
f
n
\\Delta f_n
Δfn
随着:
n
n
n
如何变化。
例如:
Δ
f
n
=
(
D
−
1
)
f
n
\\Delta f_n=(D-1)f_n
Δfn=(D−1)fn
因此:
Δ
f
n
∝
f
n
\\Delta f_n \\propto f_n
Δfn∝fn
如果你把:
Δ
f
n
\\Delta f_n
Δfn
画成模式编号
n
n
n 的函数,那么在理想共动模型下,会得到近似线性关系。
这实际上提供了一种非常漂亮的物理判据:
如果不同梳齿的 Doppler 绝对频移严格按光学频率比例变化,那么说明系统主要表现为整体 Doppler 缩放。
如果出现额外弯曲,则说明可能存在:
色散
\\text{色散}
色散
或者:
模式相关噪声
\\text{模式相关噪声}
模式相关噪声
或者:
传播介质效应
\\text{传播介质效应}
传播介质效应
或者:
光梳内部动力学
\\text{光梳内部动力学}
光梳内部动力学
等额外因素。
第二十一章 为什么“不同纵模伸缩”是一个很值得继续研究的问题
对理想真空径向 Doppler:
f
n
′
=
D
f
n
f_n'=Df_n
fn′=Dfn
所有模式的相对伸缩相同。
但真实光频梳不是一个数学上的无穷窄谱线集合。
它存在:
f
r
e
p
f_{\\mathrm{rep}}
frep
和:
f
c
e
o
f_{\\mathrm{ceo}}
fceo
同时还存在:
- 腔内色散;
- 模式相位;
- 光谱包络;
- 非线性效应;
- 载波包络相位;
- 环境扰动;
- PLL残差。
因此在高速、宽带、动态运动下,真正的实验系统并不一定表现为一个完全理想的“统一缩放器”。
这时候就会出现一个非常有价值的研究方向:
相对论 Doppler
+
频梳内部动力学
\\boxed{ \\text{相对论 Doppler} + \\text{频梳内部动力学} }
相对论 Doppler+频梳内部动力学
也就是把:
D
(
t
)
D(t)
D(t)
作为外部物理输入,然后研究:
f
r
e
p
f_{\\mathrm{rep}}
frep
和:
f
c
e
o
f_{\\mathrm{ceo}}
fceo
以及不同模式相位怎样共同响应。
对于未来双梳高速运动测量、空间光钟网络、自由空间激光时频传递,这种模型会比简单的单频 Doppler 模型更有意义。
第二十二章 和你的双光梳时频传递研究结合起来应该怎么理解
你前面一直在讨论:
双光梳、同步采样、TPFC、时间频率传递以及运动 Doppler。
把现在的知识串起来,其实可以得到一条非常漂亮的技术路线:
目标运动
\\boxed{ \\text{目标运动} }
目标运动
↓
R
(
t
)
,
v
(
t
)
\\boxed{ R(t),v(t) }
R(t),v(t)
↓
T
A
B
(
t
)
\\boxed{ T_{AB}(t) }
TAB(t)
↓
D
(
t
)
\\boxed{ D(t) }
D(t)
↓
f
n
(
t
)
,
f
r
e
p
(
t
)
,
f
c
e
o
(
t
)
\\boxed{ f_n(t),f_{\\mathrm{rep}}(t),f_{\\mathrm{ceo}}(t) }
fn(t),frep(t),fceo(t)
↓
pulse timing
\\boxed{ \\text{pulse timing} }
pulse timing
↓
dual-comb heterodyne
\\boxed{ \\text{dual-comb heterodyne} }
dual-comb heterodyne
↓
f
R
F
(
t
)
\\boxed{ f_{\\mathrm{RF}}(t) }
fRF(t)
↓
timing error + phase error
\\boxed{ \\text{timing error + phase error} }
timing error + phase error
↓
clock offset / velocity / delay
\\boxed{ \\text{clock offset / velocity / delay} }
clock offset / velocity / delay
这其实就是你目前研究方向里“运动 Doppler”和“时频传递”真正连接起来的地方。
第二十三章 为什么单纯看 Doppler 频率还不够
假设你最终得到:
f
D
(
t
)
f_D(t)
fD(t)
你可以从它计算:
v
(
t
)
v(t)
v(t)
但是时间频率传递还需要知道:
Δ
t
(
t
)
\\Delta t(t)
Δt(t)
而:
Δ
t
(
t
)
\\Delta t(t)
Δt(t)
又受到传播距离:
R
(
t
)
R(t)
R(t)
的影响。
所以仅仅知道:
f
D
(
t
)
f_D(t)
fD(t)
不能自动知道:
Δ
t
(
t
)
\\Delta t(t)
Δt(t)
这就是为什么在高精度时频传递中:
frequency transfer
≠
time transfer
\\boxed{ \\text{frequency transfer} \\neq \\text{time transfer} }
frequency transfer=time transfer
两者虽然使用同一个光学载体,但估计目标不同。
光频梳之所以特别有吸引力,是因为它天然把时间尺度和光学频率尺度绑定起来。
这也是光频梳从单纯“频率标尺”向“时间标尺”发展的原因。
你提供的 TPFC 材料也正是沿着这个思路展开:TPFC能够输出可调时序和相位,用于时间间隔基准,并进一步拓展到自由空间光频梳时频传递。
第二十四章 最后把整个问题压缩成几句话
如果只记住几个物理结论,可以把全文压缩成下面这套逻辑。
连续波激光受到 Doppler 后:
f
o
b
s
=
D
f
0
f_{\\mathrm{obs}}=Df_0
fobs=Df0
其中低速时:
D
≈
1
−
v
c
D\\approx1-\\frac vc
D≈1−cv
高速时必须使用:
D
=
1
−
β
1
+
β
D= \\sqrt{ \\frac{1-\\beta}{1+\\beta} }
D=1+β1−β
而且这个相对论结果来自:
时间膨胀
+
运动传播时延
\\boxed{ \\text{时间膨胀} + \\text{运动传播时延} }
时间膨胀+运动传播时延
对于光频梳:
f
n
=
n
f
r
e
p
+
f
c
e
o
f_n=nf_{\\mathrm{rep}}+f_{\\mathrm{ceo}}
fn=nfrep+fceo
理想径向运动下:
f
n
′
=
D
f
n
f_n'=Df_n
fn′=Dfn
所以:
f
r
e
p
′
=
D
f
r
e
p
f_{\\mathrm{rep}}'=Df_{\\mathrm{rep}}
frep′=Dfrep
以及:
f
c
e
o
′
=
D
f
c
e
o
f_{\\mathrm{ceo}}'=Df_{\\mathrm{ceo}}
fceo′=Dfceo
因此:
不同纵模的绝对 Doppler 频移不同,但理想情况下它们的相对频率缩放相同。
如果速度恒定:
f
D
=
c
o
n
s
t
a
n
t
f_D=\\mathrm{constant}
fD=constant
如果速度变化:
f
D
=
f
D
(
t
)
f_D=f_D(t)
fD=fD(t)
于是产生:
Doppler chirp
\\text{Doppler chirp}
Doppler chirp
对于双光梳:
f
r
e
p
f_{\\mathrm{rep}}
frep
和:
f
c
e
o
f_{\\mathrm{ceo}}
fceo
都会受到影响,因此 Doppler 最终会映射到:
RF beat note
\\text{RF beat note}
RF beat note
以及:
pulse timing
\\text{pulse timing}
pulse timing
最后就会出现:
delay-Doppler coupling
\\boxed{ \\text{delay-Doppler coupling} }
delay-Doppler coupling
而对于 TPFC:
它不是把 Doppler 消除,而是持续跟踪 Doppler 导致的时序和相位变化
\\boxed{ \\text{它不是把 Doppler 消除,而是持续跟踪 Doppler 导致的时序和相位变化} }
它不是把 Doppler 消除,而是持续跟踪 Doppler 导致的时序和相位变化
这也是为什么 TPFC 可以从单纯的“光频梳标尺”变成一种具有闭环跟踪能力的时域、频域光学振荡器。你提供的材料已经在实验上展示了这一思路,并明确把跟踪时序和载波相位结合起来用于动态测距和速度提取。
第二十五章 从仿真角度,最值得做的不是“画一张 Doppler 曲线”
如果后面准备把这个方向继续做成论文级仿真,我更建议把模型按照以下逻辑建立:
运动模型
→
传播模型
→
相对论 Doppler
→
光频梳模型
→
双梳拍频
→
时序估计
\\boxed{ \\text{运动模型} \\rightarrow \\text{传播模型} \\rightarrow \\text{相对论 Doppler} \\rightarrow \\text{光频梳模型} \\rightarrow \\text{双梳拍频} \\rightarrow \\text{时序估计} }
运动模型→传播模型→相对论 Doppler→光频梳模型→双梳拍频→时序估计
然后至少观察:
v
(
t
)
v(t)
v(t)
R
(
t
)
R(t)
R(t)
f
D
(
t
)
f_D(t)
fD(t)
f
r
e
p
(
t
)
f_{\\mathrm{rep}}(t)
frep(t)
f
c
e
o
(
t
)
f_{\\mathrm{ceo}}(t)
fceo(t)
ϕ
(
t
)
\\phi(t)
ϕ(t)
以及:
f
b
e
a
t
,
n
(
t
)
f_{\\mathrm{beat},n}(t)
fbeat,n(t)
最后再通过 FFT 和 STFT 去看:
Spectrum
\\text{Spectrum}
Spectrum
和:
Time-Frequency Distribution
\\text{Time-Frequency Distribution}
Time-Frequency Distribution
这样得到的结果才不是“画了几张很好看的 Doppler 图”,而是一条完整的物理链。
特别是你现在已经在做双光梳 O-TWTFT、同步采样和 TPFC,那么下一步最有价值的模型其实是:
运动目标
+
双光梳
+
Doppler
+
TPFC跟踪
+
时延估计
\\boxed{ \\text{运动目标} + \\text{双光梳} + \\text{Doppler} + \\text{TPFC跟踪} + \\text{时延估计} }
运动目标+双光梳+Doppler+TPFC跟踪+时延估计
这个体系可以进一步研究一个非常有意思的问题:
传统异步双光梳在高速运动条件下,Doppler 会如何改变
Δ
f
r
e
p
\\Delta f_{\\mathrm{rep}}
Δfrep 和等效时间轴;而 TPFC 主动跟踪以后,这些 Doppler-induced timing error 能否被控制环路直接吸收,并通过反馈量反推出目标速度。
这就从“研究 Doppler”真正走到了“研究 Doppler 对光学时频传递系统性能的影响”。
如果你在实际仿真中遇到过“明明理论 Doppler 频移正确,但双梳拍频位置还是不对”“
f
r
e
p
f_{\\mathrm{rep}}
frep 该不该跟着 Doppler 一起变”“不同梳齿到底该加同一个频移还是乘同一个比例”“高速运动时相对论修正到底应该放在哪里”等问题,欢迎继续交流。把这些地方真正理顺之后,后面的双光梳运动目标仿真会清楚很多。

