| 查看: 3511 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助龙格库塔法解激光器调Q速率方程
|
||
|
在欧攀的《高等光学仿真（MATLAB）版》——光波导、激光中，有关于被动调Q的matlab解法。 附上程序 %被动调Q速率方程数值求解 clc clear close all T0 = 0.7; R = 0.8; Rp = 2e28; y0 = [1;0;0]; tspan=[0 0.05]; tic [t,y] = ode23('rate_eq',tspan,y0,[],Rp,T0,R); 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})'); % T(t)=T0 * (t - floor(t*freqQ)/freqQ<tQ) subplot(3,1,3); plot(t,y(:,3)); xlabel('时间(s)'); ylabel('基态粒子数密度(m^{-3})'); %%%%%%%%%%%%%%%%%%%%%%% 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)]; 请问如何将这段程序改为声光调Q适用的？书作者回复： 1. 声光调Q属于主动调Q，与本书中介绍的是被动调Q是有所差异的，不过也可以用 MATLAB 来求解。 关键在于声开关透过率函数要写对，比如说声光开光的重复频率是 freqQ，打开时长是 tQ，打开时的透过率为T0，关闭时的透过率为0，则在时刻 t 的声开关透过率函数可以写成： T(t)=T0 * (t - floor(t*freqQ)/freqQ<tQ) 将该函数带入激光器的调Q速率方程，编写好 MATLAB 中的常微分方程组，然后给定初值调用 ode45()求解 实在基础太差，连改哪里都不知道。请高手帮忙，多谢！ |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！！对于调用龙格库塔法如何对结果判断反馈
已经有6人回复
- 求助用matlab编一个程序，采用龙格库塔发求求解动力学方程的参数问题
已经有5人回复
- Matlab中ode45出现大误差，关于用ODE45解偏微分方程的疑问，有待解决，大家共同讨论！
已经有4人回复
- 急！求大神用Matlab四阶龙格库塔解个方程！
已经有5人回复
- 求高手帮我解一阶常微分方程组，非常感谢！！
已经有5人回复
- 请问什么是分段龙格库塔的方法？需要用matlab软件采用这种方法求二阶微分方程的解
已经有3人回复
- mathcad怎样用四阶龙格库塔法求解二阶微分方程？
已经有8人回复
- 龙格-库塔（runger-kuta）求解微分方程
已经有3人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 一阶微分方程的解析解(伯努利方程带常数项)
已经有12人回复
- 求matlab 用四阶龙格-库塔法求解微分方程
已经有4人回复
- 龙格库塔法 matlab
已经有6人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- matlab解微分方程组
已经有15人回复
- 一个微分方程组，求解
已经有18人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
- 【求助】请教非齐次常微分方程组的解析解法
已经有4人回复
- 【求助】Matlab中利用四阶龙格-库塔法求解微分方程！！！！
已经有9人回复
- 【讨论】【【四阶龙格-库塔法求解微分方程】】
已经有8人回复
- 【求助】求解一个二阶非线性微分方程x''+Ax-B/(x^3)=C
已经有6人回复
- 【求助】急请微分方程高手看下这个方程组！
已经有14人回复
- 【求助】请教二阶微分方程的特解（含初始条件）
已经有8人回复
- 【求助】四阶龙格-库塔方法！！！！
已经有17人回复
- 【求助】20金币求一个微积分方程
已经有15人回复
3楼2018-03-22 15:04:23
2楼2016-03-31 20:58:55
