任务目标

分析光源附加噪声的理论值与实验值的计算方法,分析附加噪声与系统哪些方面相关。

激光器光源

SLD光源

SLD:超辐射发光二极管
光谱宽度经验公式
Δλ≈λ2hcqkbT\Delta \lambda \approx \frac{\lambda ^2}{h c}qk_bT Δλ≈hcλ2​qkb​T
其中
qqq——与掺杂浓度、跃迁机构、吸收有关的因子,典型值为1.8
λ\lambdaλ——峰值工作波长
PSDRIN−SLD≈1ΔfPSD_{RIN-SLD}\approx\frac{1}{\Delta f} PSDRIN−SLD​≈Δf1​
单位 Hz−1\mathrm{Hz^{-1}}Hz−1
Δf=cΔλλ2\Delta f=c\frac{\Delta \lambda}{\lambda ^2} Δf=cλ2Δλ​

掺饵光纤光源

频谱不对称,定义等效的RIN,Δλ\Delta\lambdaΔλ通过一样的式子计算
一般比SLD大,由于掺铒光纤光源的功率更大,因此实际的强度噪声大的多

激光器噪声的来源

理论计算方法1

光纤陀螺仪,中文书
散粒噪声:不相关的粒子流
σN˙2=2N˙Δfbw\sigma_{\dot{N}}^2=2\dot{N}\Delta f_{bw} σN˙2​=2N˙Δfbw​
附加噪声:源于宽带光源谱宽范围内各个傅里叶分量之间的拍效应,在光电探测器的输出光电流中形成一种附加噪声。

σN˙exc2=N˙2ΔfbwΔf\sigma_{\dot{N}_{exc}}^2=\frac{\dot{N}^2\Delta f_{bw}}{\Delta_f} σN˙exc​2​=Δf​N˙2Δfbw​​
Δfbw\Delta f_{bw}Δfbw​——计数带宽(检测频带),计数时间的倒数
Δf\Delta fΔf——频宽

粒子数的附加噪声转换为功率的附加噪声则有
(σPhυ)2=P2Δfbw(hυ)2Δf(\frac{\sigma_{P}}{h\upsilon})^2=\frac{P^2\Delta f_{bw}}{(h\upsilon)^2\Delta_f} (hυσP​​)2=(hυ)2Δf​P2Δfbw​​
即:
σP2=P2ΔfΔfbw\sigma_{P}^2=\frac{P^2}{\Delta_f}\Delta f_{bw} σP2​=Δf​P2​Δfbw​
在上述理论计算式中,Δf\Delta_fΔf​可以通过激光光谱得到,P可以由光电探测测器探测平均功率得到,Δfbw\Delta f_{bw}Δfbw​与探测的时间相关。通过以上数据计算出激光的附加噪声。

理论计算方法2

光纤陀螺仪,英文书
散粒噪声(shot noise)

经典的散粒噪声是不相关的离散粒子流,特别是离散的电子组成的电流。而激光器基础的噪声是散粒噪声是光子流。若每秒钟粒子流为N˙\dot{N}N˙,测量时间是tmt_mtm​,则探测到的粒子数为
N=N˙tmN=\dot{N}t_m N=N˙tm​
不确定性为N\sqrt{N}N​
例如激光功率为13μ\muμW,波长1500nm,则光子流为101410^{14}1014 个每秒。若测量时间是1s(探测频率的带宽是Δfd=1Hz\Delta f_d=1HzΔfd​=1Hz),光子不确定性就是10710^{7}107 个,相对不确定性就是10−710^{-7}10−7 。
其分布满足正态分布,标准差为:
σN˙=2N˙Δfd\sigma_{\dot{N}}=\sqrt{2\dot{N}\Delta f_d} σN˙​=2N˙Δfd​​
对于电流,由离散的电荷q组成,标准差为;
σI=2qIΔfd\sigma_{I}=\sqrt{2qI\Delta f_d} σI​=2qIΔfd​​
对于能量P的光,由离散的光子组成,标准差为:
σPh=2hυPΔfd\sigma_{Ph}=\sqrt{2h\upsilon P\Delta f_d} σPh​=2hυPΔfd​​
相对噪声功率谱密度(ralative noise power spectral density, PSD):
PSDN˙=σN˙2/N˙2⋅Δfd=2/N˙PSDI=σI2/I2⋅Δfd=2q/IPSDPh=σPh2/P2⋅Δfd=2hυ/P\begin{aligned} &PSD_{\dot{N}}=\sigma_{\dot{N}}^2/\dot{N}^2\cdot\Delta f_d=2/\dot{N}\\ &PSD_{I}=\sigma_{I}^2/I^2\cdot\Delta f_d=2q/I\\ &PSD_{Ph}=\sigma_{Ph}^2/P^2\cdot\Delta f_d=2h\upsilon /P \end{aligned} ​PSDN˙​=σN˙2​/N˙2⋅Δfd​=2/N˙PSDI​=σI2​/I2⋅Δfd​=2q/IPSDPh​=σPh2​/P2⋅Δfd​=2hυ/P​
附加噪声(excess RIN)
激光器中原子或者离子自发向外辐射不相干的宽带光子,自发辐射光子有随机的相位和随机的频率成分,其产生的随机能量叫做附加相对强度噪声(RIN)。假设光谱半高全宽(FWHM)为ΔfFWHM=Δλ/λˉ2\Delta f_{FWHM}=\Delta \lambda/\bar{\lambda}^2ΔfFWHM​=Δλ/λˉ2
PSDRIN=1/ΔfPSD_{RIN}=1/\Delta f PSDRIN​=1/Δf
相对应的标准差为
σ=PSDRIN\sigma=\sqrt{PSD_{RIN}} σ=PSDRIN​​
自发辐射的光对于扩展光源是不相关的,很难较好耦合到单模光纤中去。

探测器光电子计算(实验测量)

Morkel, P.R, Laming, R.I, Payne, D.N. Noise characteristics of high-power doped-fibre superluminescent sources[J]. Electronics Letters, 26(2):96.

理论公式

偏振度为d的宽带光源,当测量时间远大于相干时间时总噪声方差为
σi2=2eiˉBe+(1+d2)iˉ2Be/Δυ\sigma^2_i=2e\bar iB_e+(1+d^2)\bar i^2B_e/\Delta \upsilon σi2​=2eiˉBe​+(1+d2)iˉ2Be​/Δυ
ddd为偏振度,d=1d=1d=1时表示偏振光,d=0d=0d=0时表示自然光;iii为探测器探测到的电流;BeB_eBe​为探测器测量带宽,为测量时间二倍的倒数,即B=12TB=\frac {1}{2T}B=2T1​;Δυ\Delta \upsilonΔυ为宽带光源的频谱宽度,定义如下
Δυ=[∫P(υ)dυ]2∫P2(υ)dυ\Delta \upsilon=\frac {[\int P(\upsilon )d \upsilon]^2}{\int P^2(\upsilon )d \upsilon}Δυ=∫P2(υ)dυ[∫P(υ)dυ]2​
P(υ)P(\upsilon )P(υ)代表光谱。其倒数τc=1/Δυ\tau_c=1/\Delta \upsilonτc​=1/Δυ表示宽带光源的相干长度

在上述公式中

  • 第一项为散粒噪声
  • 第二项为附加噪声

若测量附加噪声,且光源为非偏振光,则有
σi2=iˉ22TΔυ\sigma^2_i=\frac{\bar i^2}{2T\Delta \upsilon} σi2​=2TΔυiˉ2​

公式分析

系统信噪比(假设为自然光)
SNR=iˉ2σi2=iˉ2eBe+iˉBe/ΔυSNR=\frac{\bar i^2}{\sigma^2_i}=\frac{\bar i}{2eB_e+\bar iB_e/\Delta \upsilon} SNR=σi2​iˉ2​=2eBe​+iˉBe​/Δυiˉ​
探测电流的大小决定了系统两种误差相差的量级
当iˉ≪2eΔυ\bar i \ll 2e\Delta \upsiloniˉ≪2eΔυ,散粒噪声占主导
当iˉ≫2eΔυ\bar i \gg 2e\Delta \upsiloniˉ≫2eΔυ,附加噪声占主导

测量装置

光纤
激光器
光谱分析仪
可调谐的光衰减器
光电探测器
低噪声的信号放大器
光功率探测器
平均光功率
理论值
射频谱线分析仪
实验值

原理:通过调节光的强度,逐步探测由弱光到强光时的信号,分析信号方差/标准差的波动。

光纤陀螺中的噪声

光纤陀螺中:
超辐射发光二极管,陀螺检测精度受到散粒噪声限制
掺饵光纤光源,陀螺检测精度受到相对强度噪声限制(光功率大)

RIN噪声

RIN定义

相对强度波动自相关函数的傅里叶变换,在给定的频率下,RIN表示1Hz带宽下的强度噪声功率和平均光功率对比值.

单位带宽内光强度噪声均方根与光功率平方的比值
田中成.相对强度噪声对微波光子链路性能影响规律研究
直流激光器在t时刻输出的激光功率:
P(t)=Pˉ+δP(t)P(t)=\bar{P}+\delta P(t) P(t)=Pˉ+δP(t)
Pˉ\bar{P}Pˉ为激光器输出的平均光功率,δP(t)\delta P(t)δP(t)为光功率的波动值
RIN(t)=[δP(t)]2Pˉ2=2[δP(t)]2PˉHz2BRIN(t)=\frac{[\delta P(t)]^2}{\bar{P}^2}=\frac{2[\delta P(t)]^2}{\bar{P}^2_{Hz}B} RIN(t)=Pˉ2[δP(t)]2​=PˉHz2​B2[δP(t)]2​
PˉHz2\bar{P}^2_{Hz}PˉHz2​单位带宽下的平均光功率,BBB噪声带宽因为RIN(t)的功率谱密度是单边带,因此带宽只考虑一半

实验值计算方法

方荣. 激光器相对强度噪声测量研究[D].北京邮电大学,2017.
RIN=1BΔP(ω)2P2RIN=\frac{1}{B}\frac{\Delta P(\omega)^2}{P^2}\\ RIN=B1​P2ΔP(ω)2​
其中BBB代表等效带宽,ΔP(ω)2\Delta P(\omega)^2ΔP(ω)2代表制定频率下的均方光强度波动,PPP为输出光强度平均值,RIN的单位一般用dBc/Hz\rm{dBc/Hz}dBc/Hz表示
RIN的主要来源
激光器本身:

  • 激光器的自发辐射(主要来源)
  • 载流子的涨落(主要来源)
  • 量子效率起伏
  • 外界环境波动
  • 泵浦源的电流波动
  • 多模激光器的模式分配噪声

光电探测器用光电流的方差表示
σRIN2=i2BΔυ\sigma^2_{RIN}=\frac{i^2B}{\Delta \upsilon} σRIN2​=Δυi2B​
其中Δυ=1RIN\Delta \upsilon=\frac{1}{RIN}Δυ=RIN1​,iii为探测器电流,BBB为带宽,因此RIN对整个系统附加噪声功率表示为
NRIN=1210RIN10i2BRLN_{RIN}=\frac{1}{2}10^{\frac{RIN}{10}}i^2BR_L NRIN​=21​1010RIN​i2BRL​
其中RLR_LRL​为负载阻抗。

光源噪声及其计算方法相关推荐

  1. 接收机的噪声来源与噪声分析

    噪声分类 射频接收机中的噪声主要可以分为两类:内部噪声和外部噪声. 内部噪声 内部噪声主要来自于接收机内部的放大器.混频器.本振等元件所产生的噪声.根据不同的产生机制,内部噪声可以分为以下几类: a. ...

  2. 【学习笔记】《基于φ-OTDR的分布式扰动传感系统定位算法研究-北交-通信与信息系统-吴》重点笔记

    目录 一.绪论 1.1 引言 1.2 分布式光纤传感技术 1.2.1 基于瑞利散射的分布式光纤传感技术 1.2.2 基于拉曼散射的分布式光纤传感技术 1.2.3 基于布里渊散射的分布式光纤传感技术 1 ...

  3. 小波分析用于阈值去噪

    文章目录 目录 文章目录 前言 一.小波去噪概述 二.小波阈值去噪 介绍关于阈值去噪的方法 阈值的选取 优缺点 仅供参考 前言 在去噪领域中,小波理论由于其特殊的优点受到了许多学者的重视,他们应用小波 ...

  4. 基于分层卷积神经网络的皮肤镜图像分类方法

     [摘 要]针对皮肤镜图像数量不充足以及各类疾病之间影像数据不平衡的问题,提出一种融合类加权交叉熵损失函数和分层卷积神经网络的皮肤镜图像分类方法.首先对皮肤镜图像进行色彩恒常化处理,消除环境光源噪声: ...

  5. Towards Robust Deep Hiding Under Non-Differentiable Distortions for Practical Blind Watermarking论文阅读

    Towards Robust Deep Hiding Under Non-Differentiable Distortions for Practical Blind Watermarking 发表在 ...

  6. ASE光源 ASE噪声

    ASE光源 ASE的产生 自发辐射 高能级的电子不稳定,自发的掉到低能级,并维持能量守恒,释放一个光子.自发辐射光在频率.相位.偏振方向及传播方向都是随意的. 即非相干光. 受激辐射 高能级的电子不稳 ...

  7. 对于原始轨迹的噪声过滤方法和曲率计算方法

    对于地图采样获得的轨迹,通过cross entropy的方式,进行噪声过滤.曲率计算采用三点式计算圆半径获得. clc clear allload('D:\ccs\test.mat')% %%%%%% ...

  8. UnityShader28:噪声纹理

    一.噪声纹理 听名字挺蒙的,其实它的作用就是展现出一个事物"混沌"或"杂乱无章"的效果,就像雾气和雪地,它们绝对不会是均匀单一的,不然给人的感觉就像是屏幕变灰了 ...

  9. 视觉工程师必须知道选型知识与计算方法

    视觉工程师必须知道选型知识 1.工业相机的丢帧的问题是由什么原因引起的? 经常会有一些机器视觉工程师认为USB接口的工业相机会造成丢帧现象.一般而言,工业相机丢帧与工业相机所采用的传输接口是没有关系的 ...

最新文章

  1. phpcms_v9修改文章会提示 [hash]数据验证失败
  2. let const var 比较说明
  3. SQL Server 的历史与发展方向
  4. [BZOJ 1112] [POI2008] 砖块Klo 【区间K大】
  5. python查看dataframe数据类型_python pandas中DataFrame类型数据操作函数的方法
  6. 影响 oracle 性能的常见事件
  7. 夯实Java基础(十八)——泛型
  8. C#.NET常见问题(FAQ)-如何改变字符串编码
  9. 率辉c语言,C语言真题
  10. Javascript特效:五彩小球
  11. 实时汇率转换小程序(c++爬虫)
  12. Python 批量获取今日头条街拍美图
  13. 11个主流AI聊天机器人平台,你绝不能错过!
  14. 2021年危险化学品经营单位主要负责人免费试题及危险化学品经营单位主要负责人考试总结
  15. 数字逻辑综合工具实践-DC-08——静态时序分析(STA)
  16. java 接口重写_接口与继承系列教材 (三)- Java 重写方法
  17. 企业AAA信用等级认证好处
  18. NOIP-2011-J1-真题解析
  19. C语言计算级数fun,c语言编程 编写函数fun(),它的功能是:计算和输出下列级数的和....
  20. 往数据库中插入时间,与预期值相差8个小时(时区设置)

热门文章

  1. vf mysql_VF是什么?
  2. java8使用stream流将list转map
  3. MySQL JDBC编程
  4. matlab整理符号表达式,[2018年最新整理]MATLAB符号运算与符号方程求解.ppt
  5. 【开发问题系列1】无法定位程序输入点
  6. 什么是软件开发工具包(SDK)
  7. 基于hunt1.5的TCP会话劫持
  8. 求两个球体相交部分体积,公式及推导,C++实现
  9. iOS里很简单的推送,移动端。
  10. JS复制指定DIV文本内容