qq_46129398的博客

原创 Mamyshev振荡器原理详解

Mamyshev振荡器是一种基于Mamyshev再生器级联的全光纤锁模激光器，利用自相位调制(SPM)光谱展宽和偏置滤波实现高强度脉冲

原创 混合双环结构（阴阳腔）+ 碳纳米管可饱和吸收体（SWCNT-SA）

摘要：高功率全光纤脉冲激光器的混合腔结构设计，结合非线性光纤环镜(NOLM)、非线性放大环镜(NALM)和碳纳米管可饱和吸收体(SWCNT-SA)，

原创 超快NOLM光纤激光器仿真模型的代码全解与物理内涵(3)

本文提出了一种层次化、模块化的光纤激光器仿真模型，采用自适应步长控制算法确保计算精度。模型包含线性算符构造、非线性算符计算、有源增益模块等核心组件，通过交互图像法和四阶Runge-Kutta方法实现脉冲演化求解。系统支持完整的激光腔结构模拟，包括光纤传播、NOLM非线性镜等关键器件，并提供可视化模块用于分析脉冲时频特性。该设计具有明确的物理意义和良好的可扩展性，能有效模拟激光器从初始扰动到稳态锁模的全过程。

原创 超快NOLM光纤激光器仿真模型的代码全解与物理内涵(2)

本文详细分析了2×2光纤耦合器的数学模型与物理特性。基于耦合模理论，推导了同步耦合情况下的传输矩阵表达式，建立了功率耦合比ρ与耦合长度的关系。重点讨论了四端口耦合器的相位特性和能量守恒验证，分析了不同耦合比下的特殊情形（单端口输入、50/50耦合器等）。文章还探讨了耦合器在非线性光纤环形镜中的应用原理，以及实际器件中的偏振无关性、频率依赖性和损耗模型等工程实现问题。通过严格的数学推导，证明了该模型严格满足能量守恒定律，为光纤系统的设计与分析提供了理论基础。

原创 超快NOLM光纤激光器仿真模型的代码全解与物理内涵(1)

本文提出了一种基于相互作用表象龙格库塔算法和守恒量误差法(CQEM)的超快光纤激光器数值仿真模型。该模型通过将线性算符和非线性算符分离处理，在频域积分非线性项，实现了高效稳定的脉冲传播模拟。核心求解器IP_CQEM_FD.m采用自适应步长控制，完整考虑了光纤色散、非线性效应及增益动态。模型采用模块化设计，包含耦合器、滤波器、放大器等组件，适用于锁模激光器和超连续谱产生等非线性光学现象的研究。通过光子数守恒原理实现步长自适应调整，保证了计算精度和稳定性。

原创 光纤激光器放大器内部脉冲演化：简化的奥秘

本文探讨了光纤激光器放大器模型中脉冲演化的简化方法。研究通过忽略次要损耗因素，将掺镱光纤中的皮秒脉冲视为准连续光处理，提出基于平均功率的增益饱和公式，大幅降低了计算复杂度。该方法将复杂的速率方程组简化为静态增益参数，并通过分离变量处理增益谱线型，实现了计算效率与精度的平衡。文中还展示了配套的MATLAB框架结构，包含光纤参数定义、主调用函数和核心求解器等模块，为光纤激光器设计提供了高效模拟工具。

原创 # 1智能体的任务环境与运行机制

本文阐述了智能体的任务环境特征及其核心运行机制。

原创 用SHAP解释机器学习模型预测：原理与实践指南

随着AI在关键决策中扮演越来越重要的角色，可解释性不再是"可有可无"的功能，而是必须满足的需求。SHAP为我们提供了一套强大的工具，让黑箱模型变得透明、可信、可审计。最好的模型不仅要准确，更要能解释清楚为什么准确。SHAP正是连接模型复杂性与人类理解的关键桥梁。

原创 KAN神经网络进行工程问题建模与优化

将通过一个具体的工程案例，展示如何使用 KAN 神经网络进行系统建模与优化。面对的是一个具有4个输入参数和1个输出参数的系统。数据集来自两个 Excel 文件：inout.xlsx（输入数据）和 F_out.xlsx（输出数据）。目标是：1. 建立输入参数与输出参数之间的准确映射关系；2. 用建立的模型进行系统优化，找到最优的输入参数组合。

原创 DeepSeek AI 编程助手桌面应用程序 - 完整技术文档

DeepSeek AI 编程助手是一款基于 PySide6 构建的桌面应用程序，集成了 DeepSeek 大语言模型 API，为开发者提供智能代码分析、修正、补全、解释和优化功能。核心特性- 🚀 多编程语言支持（Python、JavaScript、Java、C++ 等）- 🎯 四种专业使用场景（代码修正、补全、解释、性能优化）- 🎨 现代化白色主题用户界面- ⚡ 异步 API 调用，避免界面冻结- 📊 可调节的 AI 创造性参数- 🔧 实时字体大小调节。

原创 深度解析主动锁模光纤激光器的腔体结构与核心公式

本文提出了一种基于双调制器结构的主动锁模光纤激光器仿真模型。该模型通过扩展非线性薛定谔方程，综合考虑了掺铒光纤增益、色散、非线性效应以及强度/相位调制器的作用，采用分布傅里叶法进行数值求解。仿真结果显示，系统能够产生稳定的锁模脉冲，脉冲宽度约5ps，光谱宽度约0.4nm，峰值功率达数十瓦。研究揭示了调制深度、色散管理等因素对脉冲特性的影响规律，为高性能主动锁模激光器设计提供了理论依据和优化方法。

原创 MATLAB .一元函数的图形

MATLAB .一元函数的图形

原创 一维插值的应用画图

利用MATLAB的插值函数进行插值，并画出你的手掌外形轮廓。

原创 MATLAB Quick Start

MATLAB 快速入门

原创 麦克斯韦速率分布

代码clc;clear;figure('numbertitle','off','name','麦克斯韦速率分布');global K T m;K=10;T=300;V=0:0.01:32;for m=50:50:250;subplot(1,2,1);FF1=4*pi*(m/(2*pi*K*T)).^(3/2);FF2=V.^2;FF3=exp(-m*V.^2/(2*K*T));% 坐标轴美化FF=FF1*FF2.*FF3;plot(V,FF,'linewidth',2.5

原创 2022数学建模A动态仿真

%% 2022数学建模A动态仿真clc;clear;figure('numbertitle','off','name','2022数学建模A动态仿真');%赋值global m M K C W C1 P Fm=4800;M=2400;K=80000;C=10000;W=1.8;C1=656;P=10000;%改变F=6000;%ode45求解绘图t1=[0:0.1:40];H1=[0;10;0;5];[t,H]=ode45(@dH,t1,H1);

原创 等势线matlab仿真

原创 有限差分法求解高阶导数

有限差分法的基本思想是把微分方程进行离散化，在运用边界条件，用一个一个不确定的差分方程进行逐级迭代去近似替换。有限差分法的实质是把特殊的偏微分方程改变成为有限差分形式迭代方程，把特殊的偏微分方程用线性代数微分方程组替换，以及把数学推导过程、改变成为数值计算，进而使连续特殊偏微分方程进行离散化。关于差分法主要讨论的问题有:差分格式的建立，把特殊的偏微分方程用差分方程进行代替，边界条件设置，以及差分矩阵相关概念理解与使用MATLAB编写相应数值分析求解的代码。在有限差分算法中，步长一般称为离散化之后的单元的

原创 拟合与插值的应用

原创 三个弹簧振子理论

三个弹簧振子理论建立如图 １ 所示的一维三弹簧振子的实验系 统:三个弹簧振子的质量均为ｍ1,m2,m3，弹簧的劲度系 数分别为k1,k2和k3 ，三个谐振子在微振动过程 中相对自身的平衡位置的位移分别为x1、x2 和 x3 ，取x1、x2 和 x3 广义坐标，则系统的动力学方程可表示为:代码弹簧摆摆动图clc;clear;global k1 k2 k3 M1 M2 M3k1=145000;k2=180000;k3=121000;M1=220;M2=240;M3...

原创 经纬度下温度变化

clc;clear;data1=xlsread('1998tmp.xlsx');aa1=xlsread('AAA1.xlsx');aa2=xlsread('AAA2.xlsx');aa3=xlsread('AAA3.xlsx');aa4=xlsread('AAA4.xlsx');aa5=xlsread('AAA5.xlsx');aa6=xlsread('AAA6.xlsx');aa9=xlsread('AAA9.xlsx');A1=data1(:,1);a1=A1(2:end);B1=d

原创 弹簧摆摆动轨迹

%% 弹簧摆摆动轨迹clc;clear;global L m k gtheta0=pi/10;m=1;k=80;g=9.8;L0=1;L=L0+m*g/k;t1=[0:0.05:15];H1=[L0 0 theta0 0];[t,H] = ode45(@dH,t1,H1);%求解动力学方程figure('numbertitle','off','name','弹簧摆摆动轨迹图');[y1,x1]=pol2cart(H(:,3),-H(:,1));plot(x1(1),y1(1)..

原创 小珠子在环内的凹槽中滚动仿真图

%% 小珠子在环内的凹槽中滚动仿真图clc;clear;figure('numbertitle','off','name','小珠子在环内的凹槽中滚动仿真图');global K RK=4;w=1.3;R=0.2;r=1;t1=[0,80];H1=[0;0.3];[t,H] = ode45(@dH,t1,H1);z=r-cos(H(:,1));x=sin(H(:,1)).*cos(w*t);y=sin(H(:,1)).*sin(w*t);plot3(x,y,z,'b',.

原创 Matlab相图

%% 水平动力学运动的相图clc;clear;figure('numbertitle','off','name','水平动力学运动的相图');global K RK=5;R=0.5;t1=[-10:0.1:10];for i=[-0.5 -1.5 -3 -4.5 0.5 1.5 3 4.5] H1=[0;i]; [t,H] = ode45(@dH,t1,H1); plot(H(1:end,1),H(1:end,2)); axis([-5 5 -4 4]...

原创 水平有摩擦力K的变化对运动曲线的影响

%% 水平有摩擦力K的变化对运动曲线的影响clc;clear;global K Rfigure('numbertitle','off','name','水平有摩擦力K的变化对运动曲线的影响');for K=0:0.5:3 R=0.25; t1=[0,80]; H1=[0;0.1]; [t,H] = ode45(@dH,t1,H1); plot(t,H(:,1),'LineWidth',1.5); hold on axis([0 60...

原创 串联弹簧动态物理仿真

%% 弹簧动图clc;clear;global k1 k2 k3 M1 M2 M3k1=145000;k2=180000;k3=121000;M1=220;M2=240;M3=350;t1=[0,5];H1=[4;0;8;0;12;0];[t,H] = ode45(@dH,t1,H1);%求解动力学方程figure('numbertitle','off','name','弹簧摆摆动轨迹图');R=0.3;xx1=20;xx2=40;xx3=60...

原创 撒克逊碗动态物理仿真

%% 撒克逊碗clc;clear;global M g p R r hd N N2 Lg=9.8;M=0.13581;R=0.05;p=1000;%千克%水的密度r=0.003;hd=0.002;N=0.0001;%碗的厚度%阻尼系数小孔的N2=1;L=0.1;%碗的阻力系数%碗的高度t1=[0,40];H1=[0;0;0;0][t,H] = ode45(@dH,t1,H1);%求解动力学方程%% 动画绘制figure('numbertitle','off','name','

原创 等势线动图

clc;clear;%%等势线动图for a=2:1:14;q1=-1.602e-19;q2=1.602e-19;radius = 2; %点电荷半径x = 25; y = 25; %x和y轴显示边界值k = 9e+9; %k = 1/4πεdelta = 1;%引入距离误差以消除点电荷处点位无限大的问题phi = linspace(-pi,pi,20); %电场线发射范围[X,Y]=meshgrid(-x:0.5:x,-y:0.5:y);%获取左半部分矩阵r1=sqr...