当前位置：首页 > news >正文

第【7】期--自由空间光通信(FSO)在Gamma-Gamma湍流信道下的BER性能仿真-maltab完整代码+报告

news 2026/6/8 1:33:38

文章目录

- 摘要
- 1 研究背景
- 2 理论基础
- - 2.1 FSO基础
  - 2.2 gamma-gamma 湍流模型
- 3 仿真流程
- 4 仿真结果
- 5 结论

摘要

自由空间光通信（FSO）凭借其高带宽、免许可部署和抗电磁干扰等优势，已成为无线接入与空间通信领域的研究热点。然而，大气湍流引起的光强闪烁会严重恶化链路误码率性能。本文基于 Gamma‑Gamma 湍流模型，利用 MATLAB 平台建立了蒙特卡洛仿真框架，系统评估了 BPSK、DPSK、16PSK、64PSK 和 16QAM 五种调制格式在弱、中、强三种湍流强度下的误码率特性。仿真中同时考虑了热噪声和背景光噪声，并采用无衰落 BPSK 作为性能基准。结果表明：湍流强度越大，所有调制的性能退化越显著；BPSK 抗衰落能力最强，DPSK 提供了无需载波同步的工程折衷；在同等频谱效率下，16QAM 优于 16PSK，而 64PSK 在 FSO 中实用性较低。结果可为 FSO 系统自适应调制与编码策略的设计提供定量参考。

1 研究背景

自由空间光通信（Free Space Optical Communication, FSO）利用激光在大气中传输数据，具有高带宽、无需频谱许可、抗电磁干扰、部署快速等优点，被视为解决“最后一公里”接入瓶颈和应急通信的关键技术之一。然而，大气湍流引起的光强闪烁会导致接收信号剧烈起伏，严重恶化链路误码率（BER）。

为了定量评估不同调制方案在湍流环境下的鲁棒性，本文基于Gamma-Gamma湍流模型，通过蒙特卡洛仿真比较了 BPSK、DPSK、16PSK、64PSK 和 16QAM 五种调制方式在弱、中、强三种湍流条件下的BER性能。所有代码均在MATLAB中实现，可供研究者复现和扩展。

2 理论基础

2.1 FSO基础

自由空间光通信（Free Space Optical Communication, FSO）是一种利用激光或LED在大气、真空或太空间直接传输数据的无线通信技术。它兼具光纤通信的大带宽（可达Gbps至Tbps级）和无线部署的灵活性，无需频谱许可，且具有高安全性和抗电磁干扰能力，特别适用于楼宇互联、应急通信、卫星链路及最后一公里接入等场景。

然而，FSO的性能严重受制于大气湍流（引起光强闪烁、光束漂移）和恶劣天气（浓雾、雨雪会大幅衰减信号）。为了克服这些挑战，常采用多孔径接收、自适应光学、信道编码以及混合FSO/RF备份链路等技术。尽管存在环境影响，FSO仍是5G/6G回传、星地激光通信和量子密钥分发等领域的关键技术方向。

FSO系统的工作原理：

发送端：将电信号调制到激光二极管或LED发出的光波上（常用调制方式：OOK、BPSK、QAM等）。
光学发射天线：对光束进行准直和扩束，使其具有较小的发散角。
大气信道：光束穿过大气（或自由空间）到达接收端。
接收端：光学天线汇聚光束，光电探测器（如APD、PIN）将光信号转换回电信号，经解调恢复原始数据。

典型工作波长：850 nm（短距离）、1310 nm / 1550 nm（长距离，人眼安全等级更高）。

2.2 gamma-gamma 湍流模型

Gamma-Gamma分布被广泛用于描述从弱到强湍流范围内的光强起伏。其概率密度函数为：

3 仿真流程

接收机噪声包括热噪声和背景光噪声（太阳+天空辐射）。总噪声功率谱密度为：

参数	数值
蒙特卡洛比特数	1×10⁶
符号速率	155 MHz
探测器响应度	1 A/W
负载电阻	50 Ω
环境温度	300 K
光源波长	850 nm
接收孔径面积	1 cm²
链路距离	1 km
折射率结构常数	0.75×10⁻¹⁴ m⁻²/³
视场角	0.6 rad
光滤波带宽	1 nm
平均辐照度范围	10⁻⁸ ~ 10⁻⁵ W
弱湍流参数 (α, β)	(11.6, 10.1)
中等湍流参数 (α, β)	(4.0, 1.9)
强湍流参数 (α, β)	(4.2, 1.4)

4 仿真结果

可看到：

弱湍流曲线最靠左，中等湍流右移，强湍流最靠右
-所有调制方式的 BER 均随接收光功率增加而下降，符合通信系统的一般规律。
无衰落 BPSK（曲线作为性能下界，位于最左侧，即同等 BER 下所需光功率最低；
高阶调制如 64PSK对相位噪声和信道估计误差更敏感，即使弱湍流下的微小幅度/相位扰动也会导致误码。
部分代码

N=1e6;% 总比特数 M_bpsk=2;% BPSK M_dpsk=2;% DPSK M_psk16=16;% 16PSK M_psk64=64;% 64PSK M_qam16=16;% 16QAMM=4;Rb=155e6;%符号速率R=1;%响应度M_ind=1;%调制指数A=1;%副载波信号的振幅RL=50;%负载电阻Temp=300;%环境温度wavl=850e-9;%光源波长 %背景噪声 %接收孔径为1cm^2sky_irra=1e-3;%在850nm波长处,单位W/cm~2-wm-srsun_irra=550e-4;%在850nm波长处,单位W/cm~2-wmFOV=0.6;%单位是弧度OBP=1e-3;%光滤波带宽(μm)Isky=sky_irra*OBP*(4/pi)* FOV^2;%天空辐照度Isun=sun_irra*OBP;%太阳辐照度 %Rytov方差Range=1e3;%链路长度(m)Cn=0.75e-14;%折射率结构参数Rhol=1.23*(Range^(11/6))*Cn*(2*pi/wavl)^(7/6);%对数辐照度方差(必须小于1)Varl=Rhol;%对数强度方差r=sqrt(Varl);%对数强度标准差 %物理常数E_c=1.602e-19;%电子电荷B_c=1.38e-23;%玻耳兹曼常数Pd=A^2/2;Ktemp=(4*B_c*Temp*Rb/RL)+(2*E_c*R*Rb*(Isun + Isky));K1=3*log2(M)*((R*M_ind)^2)*Pd/(2*(M-1)*Ktemp);%背景噪声与热噪声Io=logspace(-8,-5,20);%平均接收辐照度IodBm=10*log10(Io*1e3);%平均接收辐照度(dBm)SNR2=((R.*Io).^2)./(Ktemp);SNRdB=10*log10(SNR2);SNR=10.^(SNRdB/10);% Gamma-Gamma湍流参数 alpha_w=11.6;% 弱湍流 beta_w=10.1;alpha_m=4.0;% 中等湍流 beta_m=1.9;alpha_s=4.2;% 强湍流 beta_s=1.4;