被动调Q耦合速率方程组数值Matlab仿真
一、问题综述
参考各种文献以及实验条件,得到微晶片激光器各物理量系数的典型取值列于表5.1,其中n1和n2分别是Nd3+:YAG和Cr4+:YAG的折射率。在MATLAB中编程求解该微晶片激光器对应的速率方程组,并作出光子数密度、反转粒子数密度以及基态粒子数密度随时间变化线。
二、问题分析
对于描述连续抽运下的Q脉冲序列建立过程的耦合方程组,进行数值仿真,就可以得到不同参量下被调动Q脉冲输出的精确结果;在Matlab中调用常微分方程组初值求解函ode45(),将表中各参数值带入考虑饱和吸收体激发态吸收的连续抽运、被动调Q速率方程组,进行数值计算。
三、Mtalab数值计算
I.被调Qswitch.m主函数
clc
clear
close allT0=0.7; %可饱和吸体初始透过率
R=0.8; %输出镜反射率
Rp=2e28; %抽运速率y0=[1;0;0]; %设定初值
tspan=[0 0.05]; %设定计算时间范围
tic
[t,y] = ode45('rate_eq',tspan,y0,[],Rp,T0,R); %解耦合速率方程组rate_eq
toc
y(:,1)=max(y(:,1),1);figure
subplot(3,1,1);
plot(t,y(:,1));
xlabel('时间(s)');
ylabel('光子数密度 (m^{-3})');subplot(3,1,2);
plot(t,y(:,2));
xlabel('时间(s)');
ylabel('反转离子数密度 (m^{-3})');subplot(3,1,3);
plot(t,y(:,3));
xlabel('时间 (s)');
ylabel('基态粒子数密度 (m^{-3})');figure %将光子数密度和反转离子数密度随时间变化画于同一图中
[AX,H1,H2]=plotyy(t,y(:,1),t,y(:,2));
set(H2,'LineStyle','--')
xlabel('时间 (s)')
set(get(AX(1),'Ylabel'),'string','光子数密度 (m^{-3})')
set(get(AX(2),'Ylabel'),'string','反转粒子数密度 (m^{-3})')%各种曲线共同仿真并确定波形及频率、脉冲宽度
y1=y(:,1);%去光子密度的数值,并赋给向量y1
sdy1=sign(round(diff(y1)));
sdy2=sign(round(diff(sdy1)));
pos=find(sdy2<0) + 1;%找到光子密度最大数值点n=21;%用于计算的脉冲个数
N=50;%每个脉冲前后的数值点个数
%所选脉冲宽度、脉冲间隔,并将脉冲绘制在同一图形
figure
hold on
for i=2:ntp=[pos(i)-N:pos(i)+N]';dt=mean(diff(t(tp)));plot(-dt*N:dt:dt*N,y(tp,1))ymax(i)=y1(pos(i));pp=find(y1(tp)>=ymax(i)/2);pwidth(i)=t(pos(n)-N+pp(end))-t(pos(n)-N+pp(1));pspace(i)=t(pos(i+1))-t(pos(i));
end
frequency=1/mean(pspace(2:n));
width=mean(pwidth);
II.rate_eq.m函数
function Dy =rate_eq(t,y,flag,Rp,T0,R)
sigma=5.4e-23;
sigma_gs=8.7e-23;
sigma_es=2.2e-23;
N_T=1.68e26;
tao_a=750e-6;
tao_gs=3e-6;
n1=1.82;
n2=1.80;
delta=0.02;
l=0.001;
ls=0.001;
gamma=1;
c=2.997963e8;
lc=n1*l+n2*ls;tr=lc/c;%光子往返一周的时间
n0s= -log(T0)/(sigma_gs*ls);%饱和吸收粒子数密度
y(1)=max(y(1),1);%光子数密度的最小值
%被调用Q耦合方程组
Dy=[ y(1)*(2*sigma*y(2)*l-2*sigma_gs*y(3)*ls-2*sigma_es*...(n0s-y(3))*ls-(log(1/R)+delta))/tr;Rp*(1-y(2)/N_T)-gamma*sigma*c*y(1)*y(2)-y(2)/tao_a;(n0s-y(3))/tao_gs-sigma_gs*c*y(1)*y(3)];
四、图形绘制与分析
I.光子数密度、反粒子数密度、基态粒子数仿真
a.各粒子分别在时域分布曲线
b.各粒子具体图例分析
1)光子数密度每隔一段固定时间产生一个光脉冲:
2)增益反粒子数每隔一段固定时间缓慢上升然后急剧下降,再缓慢上升:
3)单峰放大显示,基态粒子数每隔一段时间急剧下降然后再急剧上升:
II.光子数密度和增益反转粒子数密度时域分布共同仿真
a.整体反转粒子数密度与光子数密度随时间变换关系
b.第四个单光子数脉冲与反转粒子随时间变换关系
c.[分析]在被调Q过程中,每产生一个光脉冲,增益反粒子数密度急剧减小。
III.多个激光脉冲波形共同仿真
a.多激光脉冲复合仿真图
b.复合波形的峰值及宽度显示
[注]:多个激光波形大致轮廓贴近。
IV.运行得到脉冲频率及脉冲宽度的结果
五、反思与总结
a.被动调Q主要是将表中各参数值带入考虑饱和吸收体激发态吸收的连续抽运、被动调Q速率方程组,再调用Matlab中常微分方程组初值求解函数ode45(),对被动调Q速率方程组进行数值计算,得到y:39401x3 double的矩阵,根据y(:,1),y(:,2),y(:,3)分别求得各粒子时间的曲线及对应关系。
b.最后对计算得到的第一列数组进行先微分、四舍五入,再取符号函数,进行两次;再循环21次处理求解频率、波宽及各激光复合曲线。
被动调Q耦合速率方程组数值Matlab仿真相关推荐
- matlab 调q速率方程,ode45求解被动调Q速率方程,一直显示输入参量太多,怎么回事?...
查看ode45帮助文档,你在使用ode45函数时错误 [t,y] = ode45(odefun,tspan,y0) [t,y] = ode45(odefun,tspan,y0,options) 个人改 ...
- matlab 调q速率方程,被动调Q激光器速率方程.pdf
被动调Q激光器速率方程 214 IEEE JOURNAL OF QUANTUM ELECTRONICS, VOL. 25, NO. 2, FEBRUARY 1989 Theory of the Opt ...
- [准连续泵浦被动调Q激光器速率方程求解
@[准连续泵浦被动调Q激光器速率方程求解) 有哪位大神会脉冲泵浦被动调Q速率方程的求解啊?主要是泵浦速率那一项的程序啊,在此先行谢过!!!!主程序如下: clc clear close all T0 ...
- matlab仿真转速波形为负,转速、电流双闭环直流调速系统的课程设计MATLAB仿真.docx...
转速.电流双闭环直流调速系统的课程设计MATLAB仿真 PAGE \* MERGEFORMAT25 任务书1.设计题目转速.电流双闭环直流调速系统的设计2.设计任务某晶闸管供电的双闭环直流调速系统,整 ...
- matlab交流调压调速课本,双闭环交流调速系统的设计及MATLAB仿真
内容简介: 毕业设计 双闭环交流调速系统的设计及MATLAB仿真,共44页,17796字. 摘 要 本文首先介绍了交流调速系统的基本概况,并对MATLAB做了简单介绍.同时提出了本文所要研究解决的问题 ...
- 码元速率 matlab,[转载]基于小波变换的移相键控信号符号速率估计(matlab仿真)...
西安电子科技大学 西电大宝 在电子技术迅猛发展的当代社会,空间中充满了各种各样不同频率.不同调制类型的通信信号.正 常通信条件下,发送方和接受方进行的是合作通信,即接受方预先知道发送信号的频率.调制类 ...
- m基于随机接入代价的异构网络速率分配算法matlab仿真
目录 1.算法描述 2.仿真效果预览 3.MATLAB核心程序 4.完整MATLAB 1.算法描述 无线接入技术发展迅速,异构网络并存的现象普遍存在:同时,随着终端用户数量的剧增.业务类型的多样化和高 ...
- fvdm 跟驰模型 matlab仿真_MATLAB数值计算在光学仿真和教学中的应用
摘要:在光学教学过程中从基本的物理概念出发,建立相应的理论模型,并将光学问题归纳为特征方程求根.积分求解.常微分方程求解等几类数值求解问题,结合MATLAB强大的数值计算和图形显示功能,完成光学问题的 ...
- [激光原理与应用-24]:《激光原理与技术》-10- 激光产生技术-调Q技术、Q开关、Q驱动器
目录 第1章 调Q技术概述 1.1 什么是Q 1.2 什么是调Q技术 1.3 调Q的目的 1.4 调Q的原理 第2章 主要的调Q技术分类 2.2 电光调Q 2.3 被动调Q 第3章 调Q激光器参数对比 ...
最新文章
- iptables白名单配置
- Oracle编程入门经典 第4章 新9i示例模式
- Python读取文本文档转化成列表
- CentOS系统安装Java
- UVA10679 I Love Strings!!【字符串匹配】
- jquery 插件 分析
- hx711压力传感器工作原理_压电式压力传感器原理,你了解吗?
- 海康威视监控使用html播放
- idea代码格式化详细总结(快捷键、格式化规则、格式化模板)
- C++开发坦克大战--补充(加入传送门)--附完整代码
- Java获取图像Exif信息
- php搜索银行所在支行,如何查询怎么查建设银行所属支行?
- UniAPP HBuilderX 运行到各个小程序开发工具
- ImageLoader 修改个人头像
- 普拉纳夫:第六感技术的惊异潜力
- 51单片机 | 并行I/O口扩展实例(74LS244/74LS373/4071)
- Linux上安装pstree命令(-bash: pstree: command not found)
- 技术总监的校招之旅——分享校招全攻略指南(万字长文)
- 离线环境安装python第三方库
- 步进电机和伺服电机优缺点
热门文章
- 2021年为什么这么多人劝退Android开发,学习正确姿势是什么?
- 联邦学习框架和数据隐私综述
- 计算机联用测定无机盐溶解热测试题,实验一 计算机联用测定无机盐溶解热.doc...
- 深入理解Java虚拟机:JVM高级特性与最佳实践(第3版)读书笔记
- (2)stata的基本使用--分类回归 logit
- 剥开区块链的坚果——FIBOS远航分享
- 北斗授时装置(定时服务)京准为您来解读
- 传智汇2018首秀:携手京东打造精品技术盛宴
- 解决xxx packages are looking for funding
- 网络IP检测框架的基本设计思路