基于PGC相位生成载波调制及其解调算法实例分析

相较于其他类型的光纤传感器,目前人们普遍认为干涉型光纤传感器是最有发展前途的。干涉型光纤传感器属于相位调制型光纤传感器,它是以光纤中光的相位变化来表示被测信息,首先需要对光束进行干涉得到干涉信号,再对得到的干涉信号进行解调得到相位信息。无论此时传感器感知的被测信息是哪种物理现象:磁场、声场、还是温度、位移等,在理论上,干涉型光纤传感器都具有很高的灵敏度和分辨率,结构灵活多样,适用场合广泛。其基本传感原理为:当被测信息量作用于光纤,使得光纤内传播的光波的相位发生变化,再使用相位生成载波调制解调(phase generated carrier-PGC)这种相干解调技术将相位变化转换为光强变化,从而检测出待测的信息量。

一、PGC相位生成载波调制实现原理分析

在这里插入图片描述
如图所示外调制方式的PGC系统,干涉输出信号可以表示为
在这里插入图片描述
载波信号为 C cos ⁡ ω 0 t C\cos {{\omega }_{0}}t Ccosω0t φ ( t ) \varphi \left( t \right) φ(t)为待测信号和环境漂移共同引起的相位变化。 C cos ⁡ ω 0 t + φ ( t ) C\cos {{\omega }_{0}}t+\varphi \left( t \right) Ccosω0t+φ(t)处于相位上,故将其看作相位变化,事实上载波信号和待测信号处于同等地位,由于载波频率远远高于待测信号,故可认为它载着待测信号,这就体现了相位载波的含义。

对式(1)进行贝塞尔函数展开得
在这里插入图片描述

由式(2)可知,干涉输出信号I,当 φ ( t ) = 0 \varphi \left( t \right)=0 φ(t)=0时,信号I中只存在 ω 0 {{\omega }_{0}} ω0的偶次谐波项;当 φ ( t ) = π / 2    \varphi \left( t \right)={\pi }/{2}\; φ(t)=π/2时,信号I中只存在 ω 0 {{\omega }_{0}} ω0的奇次谐波项。 φ ( t ) \varphi \left( t \right) φ(t)可表示为

在这里插入图片描述
其中, D cos ⁡ ω s t D\cos {{\omega }_{s}}t Dcosωst为待测信号幅度, φ 0 ( t ) {{\varphi }_{0}}\left( t \right) φ0(t)为环境噪声引起的相位变化。同理, cos ⁡ ω ( t ) \cos \omega (t) cosω(t) sin ⁡ ω ( t ) \sin \omega \left( t \right) sinω(t)按贝塞尔函数展开

在这里插入图片描述
由式(4)可以看出,当 φ 0 ( t ) = k π ( k = 012 ⋯   ) {{\varphi }_{0}}\left( t \right)=k\pi \begin{matrix} {} & (k=012\cdots ) \\ \end{matrix} φ0(t)=(k=012)时,输出干涉信号频谱中,偶(奇)次角频率 ω s {{\omega }_{s}} ωs出现在偶(奇)次角频率 ω 0 {{\omega }_{0}} ω0的两侧;当 φ 0 ( t ) = k π + π / 2    ( k = 012 ⋯   ) {{\varphi }_{0}}\left( t \right)=k\pi +{\pi }/{2}\;\begin{matrix} {} & (k=012\cdots ) \\ \end{matrix} φ0(t)=+π/2(k=012)时,频谱上偶(奇)次角频率 ω s {{\omega }_{s}} ωs出现在奇(偶)次角频率 ω 0 {{\omega }_{0}} ω0的两侧,待测信号的信息包含在奇(偶)次角频率 ω 0 {{\omega }_{0}} ω0的两侧的边带频谱中,它们或以 ω 0 {{\omega }_{0}} ω0偶次角频率,或以 ω 0 {{\omega }_{0}} ω0奇次角频率为中心。从而实现了相位生成载波调制。

未加载波调制前,当 φ ( t ) = k π ( k = 012 ⋯   ) \varphi \left( t \right)=k\pi \begin{matrix} {} & (k=012\cdots ) \\ \end{matrix} φ(t)=(k=012)时, cos ⁡ ω t = ± 1 \cos \omega t=\pm 1 cosωt=±1 φ ( t ) = k π + π / 2    ( k = 012 ⋯   ) \varphi \left( t \right)=k\pi +{\pi }/{2}\;\begin{matrix} {} & (k=012\cdots ) \\ \end{matrix} φ(t)=+π/2(k=012)时, cos ⁡ ω t = 0 \cos \omega t=0 cosωt=0。此时干涉信号将发生消隐或畸变,无法从中解调将待测信号。根据以上分析,加入载波信号后,即使出现 φ ( t ) = k π \varphi \left( t \right)=k\pi φ(t)= φ ( t ) = k π + π / 2    \varphi \left( t \right)=k\pi +{\pi }/{2}\; φ(t)=+π/2也不会发生信号的消隐或畸变,从而实现抗相位衰落。这就是进行相位生成载波调制的原因和意义所在。

二、PGC解调算法原理分析及相关参数选取

传统的PGC解调算法分为微分交叉相乘(DCM)算法反正切(Arctan)算法两种。以下将分别对两种解调算法进行详细分析,并分别对两种算法中的一些参数的选取原则进行详细的理论分析,如最佳相位调制度的选取、系统最低采样频率及系统解调动态范围上限等。

在这里插入图片描述

2.1 微分交叉相乘DCM解调算法原理

DCM 解调算法的思路是将干涉信号分别与单倍频和二倍频混频低通滤波后,得到一对相互正交的余弦项和正弦项,然后再经过微分交叉相乘相减、积分、高通滤波后实现待测信号的解调。具体过程如下:干涉信号分别与幅度分别为G、H,角频率为 ω 0 {{\omega }_{0}} ω0 2 ω 0 2{{\omega }_{0}} 2ω0的载波混频、低通滤波后得到

在这里插入图片描述
经过微分交叉相乘相减积分处理后的信号为
在这里插入图片描述
再通过高通滤波器得到了待测信号,即DCM算法解调输出为
在这里插入图片描述

待测信号被解调出来,只是幅值变化了一个系数 B 2 G H J 1 ( C ) J 2 ( C ) {{B}^{2}}GH{{J}_{1}}\left( C \right){{J}_{2}}\left( C \right) B2GHJ1(C)J2(C) 。为了减小输出结果对贝塞尔函数的依赖关系,通过选择适当的载波信号幅度即相位调制度C,使得 J 1 ( C ) J 2 ( C ) {{J}_{1}}\left( C \right){{J}_{2}}\left( C \right) J1(C)J2(C)取得极大值,且当C值稍有变化时系统解调输出幅值变化不大,再可以通过幅度补偿实现待测信号的完全解调。

2.2 Arctan 解调算法

Arctan 算法与 DCM 算法相同之处在于:二者都是分别与单倍频和二倍频混频低通滤波后得到两个相互正交的余弦项和正弦项。不同之处在于:反正切算法是将得到的两个正交项进行除法运算得到正切信号,然后对正切信号进行反正切算法,最后经过高通滤波实现信号的解调。

具体分析过程如下:干涉信号与载波混频低通后得到的正交项

在这里插入图片描述
对两式相除后得到
在这里插入图片描述求反正切再经高通滤波后得到
在这里插入图片描述
实现对待测信号的解调。

三、实例验证

3.1 Arctan反正切解调算法验证

首先分别模拟生成调制信号,与待测信号,并对待测信号添加噪声后进行调制,得到的调制信号如下图所示:
在这里插入图片描述

采用反正切解调算法进行解调,结果如下图所示:

在这里插入图片描述
该Matlab程序资源,可从下方获取

https://download.csdn.net/download/qq_36584460/85322776
下载后如有程序运行问题,请及时联系博主。

3.2 微分交叉相乘DCM解调算法实例演示

首先分别模拟生成调制信号,与待测信号,并对待测信号添加噪声后进行调制,得到的调制信号如下图所示:

在这里插入图片描述
采用微分交叉相乘DCM解调算法进行解调,结果如下图所示:

在这里插入图片描述
该Matlab程序资源,可从下方获取:

https://download.csdn.net/download/qq_36584460/85322781

下载后如有程序运行问题,请及时联系博主。

四、Labview用户界面

最后,开发了基于PGC相位生成载波调制及其解调的labview用户界面,如有需要,可私信博主。

在这里插入图片描述

五、参考文献

王燕. 干涉型光纤传感器及PGC解调技术研究, 天津理工大学, (2014).

  • 14
    点赞
  • 85
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
相位生成载波(Phase Generated Carrier,PGC调制是一种数字调制技术,它通过改变信号的相位来实现调制PGC调制的原理是将待调信号与一个参考载波进行相乘,使得待调信号的相位跟随参考载波相位变化。PGC调制的一个重要应用是数字通信中的频偏补偿。 PGC调制的步骤如下: 1. 将待调信号进行数字化处理,得到离散时间的信号(例如,使用matlab中的函数 `sampling`)。 2. 生成一个参考载波信号,频率与待调信号相同(例如,使用matlab中的函数 `cos`)。 3. 将待调信号与参考载波进行乘积运算,得到相位调制后的信号(例如,使用matlab中的点乘操作 `.*`)。具体公式为:调制信号 = 待调信号 * 参考载波。 4. 将相位调制后的信号发射出去。 相位生成载波解调的步骤如下: 1. 接收到相位调制后的信号。 2. 生成一个与发送端相位参考载波相同的本地载波信号(例如,使用matlab中的函数 `cos`)。 3. 将接收到的信号与本地载波进行乘积运算,得到解调后的信号(例如,使用matlab中的点乘操作 `.*`)。具体公式为:解调信号 = 接收信号 * 本地载波。 4. 对解调信号进行数字化处理,得到解调后的离散时间信号。 示例matlab代码如下: ```matlab % PGC调制 sampling_rate = 1000; % 采样率 time = 0:1/sampling_rate:1; % 时间序列 carrier_frequency = 20; % 载波频率 reference_carrier = cos(2*pi*carrier_frequency*time); % 参考载波信号 % 待调信号 message_signal = sin(2*pi*4*time); % 假设待调信号为正弦波 % PGC调制 modulated_signal = message_signal .* reference_carrier; % PGC解调 received_signal = modulated_signal; % 接收到的信号 local_carrier = reference_carrier; % 本地载波信号 demodulated_signal = received_signal .* local_carrier; % 显示结果 subplot(2,1,1); plot(time, modulated_signal); title('PGC Modulated Signal'); subplot(2,1,2); plot(time, demodulated_signal); title('PGC Demodulated Signal'); ``` 上述代码中,我们首先生成了一个参考载波信号和一个待调信号,然后进行PGC调制解调,最后绘制了调制解调信号的图形。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

简单光学

您的鼓励是我创作的最大动力!

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值