分数阶效应下饱和非线性介质中艾里高斯光束传输仿真代码功能说明

Matlab光场调控的仿真代码（全套复现论文） 之前本科搞大创发了篇文章，纯搞光场调控的仿真，后来读研不做这个方向了，寻思卖了 Tips：本科生毕设，研究生搞理论的，领域为非线性光学的、光场调控做仿真的，这部分代码可能真的会帮到你一些呢，哈哈 这部分matlab代码写出了分布傅立叶算法（用来解非线性薛定谔方程）的数值求解过程，光场函数你买了可以随便你调，不难理解

一、代码整体概述

本套代码基于Matlab开发，是一套完整复现“分数阶效应下饱和非线性介质中艾里高斯光束传输与交互特性”研究的仿真工具集。代码以分步傅里叶法为核心数值计算方法，系统探究了单/双艾里高斯光束在饱和非线性介质中的传输演化规律，重点分析了分数阶Lévy指数、饱和非线性强度、高斯调制参数等关键物理参数对光束传输特性的影响。

Matlab光场调控的仿真代码（全套复现论文） 之前本科搞大创发了篇文章，纯搞光场调控的仿真，后来读研不做这个方向了，寻思卖了 Tips：本科生毕设，研究生搞理论的，领域为非线性光学的、光场调控做仿真的，这部分代码可能真的会帮到你一些呢，哈哈 这部分matlab代码写出了分布傅立叶算法（用来解非线性薛定谔方程）的数值求解过程，光场函数你买了可以随便你调，不难理解

整套代码包含15个功能文件，按用途可分为核心算法文件、单光束传输仿真文件、双光束交互仿真文件、参数敏感性分析文件四类，具备模块化程度高、参数可配置性强、可视化效果专业等特点，可直接用于非线性光学领域的学术研究与论文图表生成。

二、核心算法与物理模型

（一）核心数值方法：分步傅里叶法

光束在非线性介质中的传输遵循非线性薛定谔方程，代码采用分步傅里叶法（Split-Step Fourier Method, SSFM）求解该方程，将传输过程分解为线性传输和非线性作用两个独立步骤，交替在频域和空域中计算，具体流程如下：

1. 频域线性传输：对当前光场进行傅里叶变换（fft），引入分数阶色散相位因子exp(-i/2dz/2.(w.*w).^(alpha/2))，其中w为角频率，alpha为分数阶Lévy指数，实现分数阶效应下的线性色散演化；
2. 空域非线性作用：将频域演化后的光场逆傅里叶变换（ifft）回空域，通过饱和非线性相位因子exp(betaidz*abs(Fai).^2./(1+abs(Fai).^2))模拟介质的饱和非线性响应，其中beta为饱和非线性强度参数；
3. 二次频域线性传输：再次对光场进行傅里叶变换，施加剩余的线性传输相位因子，完成一个传输步长（dz）的演化。

该方法兼顾了计算精度与效率，能够准确捕捉分数阶效应与饱和非线性共同作用下光场的复杂演化行为。

（二）光场初始模型

1. 单艾里高斯光束

初始光场由艾里函数、指数衰减项、高斯调制项和相位调制项组成，表达式如下：

u = u0*airy(x).*exp(a*x).*exp(-g*x.^2).*exp(i*v*x);

各参数物理意义：

2. 双艾里高斯光束

用于模拟两光束交互特性，初始光场为两束单艾里高斯光束的叠加，支持同相位/反相位配置：

u = u0*(airy(x-d).*exp(a*(x-d)).*exp(-g*(x-d).^2).*exp(i*v*(x-d))+... exp(i*Q*pi).*airy(-x-d).*exp(-a*(x+d)).*exp(-g*(x+d).^2).*exp(-i*v*(x+d)));

新增参数说明：

（三）关键物理参数定义

三、文件结构与功能分工

整套代码按功能分为四类文件，各文件职责明确、相互独立，具体分工如下：

（一）核心算法文件

1. H.m：傅里叶级数辅助函数，生成特定余弦级数组合（Y=sum(cos(i*x))，i取1到N的奇数），可用于光场初始分布的辅助构建与验证，输出级数曲线。
2. SSFFTsuanfa.m：分步傅里叶法核心演示文件，简化版单光束传输仿真代码，直接硬编码参数（alpha=1、beta=2），可快速验证算法正确性，输出光场强度随传输距离的演化数据（U为光场复振幅，I为强度）。

（二）单光束传输仿真文件

用于生成单光束传输的二维强度分布图，探究不同参数组合下的光场演化结构：

1. fig1.m：分析分数阶Lévy指数（alpha=[1,2]）和饱和非线性强度（beta=[0,2]）对单光束传输的影响，输出4个子图（2×2布局）的强度分布伪彩图，坐标轴标注规范，支持论文直接使用。
2. fig2.m：固定alpha=1，重点探究不同饱和非线性强度（beta=[-2,1,3,5]）的影响，输出4个子图（2×2布局），优化了坐标轴颜色与边框样式，可视化效果更清晰。

（三）参数敏感性分析文件

聚焦关键参数对光束传输特征量（峰值强度、孤子宽度）的定量影响，生成量化分析曲线：

1. 分数阶Lévy指数（alpha）影响分析

• fig3a.m：固定beta=2，计算alpha=[1,1.5,2]时峰值强度随传输距离的变化曲线，输出单图多曲线对比图，标注不同alpha值对应的曲线。
• fig3b.m：计算alpha在1~2区间（步长0.05）内的平均峰值强度，输出alpha-平均峰值强度关系曲线，用于定量分析分数阶效应对光束强度的整体影响。
• fig3c.m：分析不同alpha下孤子宽度的振荡特性，通过插值计算半高全宽（FWHM），输出孤子宽度随传输距离的振荡曲线。
• fig3d.m：定量分析alpha对平均孤子宽度的影响，输出alpha-平均孤子宽度关系曲线，为孤子宽度调控提供数据支撑。

2. 高斯调制参数（g）影响分析

• fig4a.m：固定alpha=1、beta=2，计算g=[0.001,0.01,0.1,1]时峰值强度随传输距离的变化，输出多曲线对比图，并通过箭头标注各曲线对应的g值。
• fig4b.m：计算g在0.001~1区间（20个采样点）内的平均峰值强度，输出g-平均峰值强度关系曲线，量化高斯调制对光束强度的抑制/增强效应。

3. 光场衰减系数（a）影响分析

• fig4c.m：固定g=0.01、alpha=1、beta=2，计算a=[0.05,0.25,0.5,0.8]时峰值强度随传输距离的变化，输出多曲线对比图，标注各曲线对应的a值。
• fig4d.m：计算a在0.01~0.99区间（12个采样点）内的平均峰值强度，输出a-平均峰值强度关系曲线，分析衰减系数对光束传输稳定性的影响。

（四）双光束交互仿真文件

探究双光束在不同相位差、间隔、参数组合下的交互特性，生成多组对比图：

1. fig5.m（函数名标注为fig7）：分析双光束间隔（d=[-6,-2,0,2,6]）、高斯调制参数（g=[0.01,1]）和相位差（Q=[0,1]）对交互的影响，输出20个子图（4行5列布局），涵盖不同参数组合的强度分布，支持同相吸引、反相排斥现象的可视化验证。
2. fig6.m（函数名标注为fig8）：固定g=0.01、d=-2，分析饱和非线性强度（beta=[-1,0,0.5,1,1.5]）和分数阶Lévy指数（alpha=[1,1.5]）对双光束交互的影响，输出20个子图（4行5列布局），重点展示非线性强度与分数阶效应的耦合作用。
3. fig8_ceshi.m：fig6.m的补充验证文件，聚焦负饱和非线性强度（beta=[-5,-4,-3,-2,-1]），分析强自散焦情况下双光束的交互特性，输出20个子图（4行5列布局）。

四、代码运行流程与可视化输出

（一）单文件运行流程

以任意仿真文件（如fig1.m）为例，运行流程如下：

1. 初始化环境：clear all; close all; clc; clf;清除工作区变量、关闭现有图形窗口、清空命令行；
2. 加载配置：加载颜色映射（load mycolor）、设置字体大小（Tfs/Afs/Lfs）、定义参数组合（alpha/beta/g等）；
3. 循环仿真：遍历所有参数组合，调用内部函数evolution1执行分步傅里叶迭代计算；
4. 结果可视化：通过pcolor绘制强度分布伪彩图（shading interp实现插值平滑），设置颜色范围（caxis）、坐标轴范围，添加子图标注（text）和图例；
5. 输出信息：显示传输距离进度（disp），计算运行时间（tic/toc）。

（二）可视化输出类型

1. 二维强度分布伪彩图：用于展示光束在“横向坐标（x）-传输距离（z）”平面的强度演化，颜色深浅表示强度大小，典型应用于fig1.m、fig5.m等，可直观观察光束聚焦、孤子形成、双光束吸引/排斥等现象；
2. 峰值强度演化曲线：以传输距离（z）为横轴，峰值强度为纵轴，展示强度随传输的振荡或衰减趋势，典型应用于fig3a.m、fig4a.m等；
3. 参数-特征量关系曲线：以物理参数（alpha/g/a）为横轴，平均峰值强度或平均孤子宽度为纵轴，用于定量分析参数对光束特性的影响，典型应用于fig3b.m、fig4b.m等。

五、代码特点与使用建议

（一）核心特点

1. 模块化设计：所有仿真文件共享evolution1内部函数，避免代码冗余，参数修改集中在文件头部，便于快速调整；
2. 专业可视化：采用Times/Times New Roman字体，坐标轴标注规范（斜体表示物理量），子图布局合理，颜色映射统一，输出图表符合学术论文要求；
3. 高可扩展性：支持灵活调整参数范围（如alpha、beta、g等），可通过增加参数组合、修改传输步长（dz）或采样点数（N）优化仿真精度；
4. 实时进度反馈：仿真过程中通过disp输出当前传输距离，便于监控计算进度，tic/toc记录运行时间，可评估计算效率。

（二）使用建议

1. 环境要求：需安装Matlab R2016b及以上版本，无需额外工具箱，直接运行.m文件即可；
2. 参数调整：修改参数时建议保持单一变量原则（如分析alpha影响时固定其他参数），避免多变量耦合导致结果难以解读；
3. 精度与效率平衡：传输步长dz越小、采样点数N越大，计算精度越高，但运行时间越长，默认参数（dz=0.1、N=2^12）已兼顾精度与效率；
4. 结果保存：可在代码中添加save命令保存仿真数据（如save('data.mat','x','zplot','I')），或通过print命令导出高质量图片（如print('-dpng','-r300','fig1.png')）。

六、典型应用场景

1. 学术研究：复现《基于分数阶效应的饱和非线性介质中艾里高斯光束的传输与交互》论文中的仿真结果，支撑非线性光学领域的学术研究；
2. 参数优化：通过调整alpha、beta、g等参数，探究呼吸孤子的形成条件与调控方法，为实验研究提供理论参考；
3. 教学演示：直观展示分数阶效应、饱和非线性、光束交互等物理现象，辅助非线性光学课程的教学演示。