 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 3391 | 回復(fù): 2 | ||
[求助]
求助龍格庫塔法解激光器調(diào)Q速率方程
|
|
在歐攀的《高等光學(xué)仿真（MATLAB）版》——光波導(dǎo)、激光中，有關(guān)于被動調(diào)Q的matlab解法。 附上程序 %被動調(diào)Q速率方程數(shù)值求解 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('光子數(shù)密度(m^{-3})'); subplot(3,1,2); plot(t,y(:,2)); xlabel('時間(s)'); ylabel('反轉(zhuǎn)粒子數(shù)密度(m^{-3})'); % T(t)=T0 * (t - floor(t*freqQ)/freqQ<tQ) subplot(3,1,3); plot(t,y(:,3)); xlabel('時間(s)'); ylabel('基態(tài)粒子數(shù)密度(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); %被動調(diào)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)]; 請問如何將這段程序改為聲光調(diào)Q適用的？書作者回復(fù)： 1. 聲光調(diào)Q屬于主動調(diào)Q，與本書中介紹的是被動調(diào)Q是有所差異的，不過也可以用 MATLAB 來求解。 關(guān)鍵在于聲開關(guān)透過率函數(shù)要寫對，比如說聲光開光的重復(fù)頻率是 freqQ，打開時長是 tQ，打開時的透過率為T0，關(guān)閉時的透過率為0，則在時刻 t 的聲開關(guān)透過率函數(shù)可以寫成： T(t)=T0 * (t - floor(t*freqQ)/freqQ<tQ) 將該函數(shù)帶入激光器的調(diào)Q速率方程，編寫好 MATLAB 中的常微分方程組，然后給定初值調(diào)用 ode45()求解 實在基礎(chǔ)太差，連改哪里都不知道。請高手幫忙，多謝！ |
回復(fù)此樓» 猜你喜歡
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有6人回復(fù)
-
318求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過
已經(jīng)有6人回復(fù)
-
一志愿武理材料工程348求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
275求調(diào)劑
已經(jīng)有6人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有7人回復(fù)
-
303求調(diào)劑
已經(jīng)有6人回復(fù)
-
306求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助��！對于調(diào)用龍格庫塔法如何對結(jié)果判斷反饋
已經(jīng)有6人回復(fù)
- 求助用matlab編一個程序，采用龍格庫塔發(fā)求求解動力學(xué)方程的參數(shù)問題
已經(jīng)有5人回復(fù)
- Matlab中ode45出現(xiàn)大誤差，關(guān)于用ODE45解偏微分方程的疑問，有待解決，大家共同討論！
已經(jīng)有4人回復(fù)
- 急！求大神用Matlab四階龍格庫塔解個方程！
已經(jīng)有5人回復(fù)
- 求高手幫我解一階常微分方程組，非常感謝��！
已經(jīng)有5人回復(fù)
- 請問什么是分段龍格庫塔的方法？需要用matlab軟件采用這種方法求二階微分方程的解
已經(jīng)有3人回復(fù)
- mathcad怎樣用四階龍格庫塔法求解二階微分方程？
已經(jīng)有8人回復(fù)
- 龍格-庫塔（runger-kuta）求解微分方程
已經(jīng)有3人回復(fù)
- 【求助】非線性方程組的求解問題
已經(jīng)有6人回復(fù)
- 一階微分方程的解析解(伯努利方程帶常數(shù)項)
已經(jīng)有12人回復(fù)
- 求matlab 用四階龍格-庫塔法求解微分方程
已經(jīng)有4人回復(fù)
- 龍格庫塔法 matlab
已經(jīng)有6人回復(fù)
- 如何求解描述振動的二階微分方程
已經(jīng)有21人回復(fù)
- matlab解微分方程組
已經(jīng)有15人回復(fù)
- 一個微分方程組，求解
已經(jīng)有18人回復(fù)
- LABVIEW四階龍格庫塔法求解一階微分方程組
已經(jīng)有12人回復(fù)
- 【求助】請教非齊次常微分方程組的解析解法
已經(jīng)有4人回復(fù)
- 【求助】Matlab中利用四階龍格-庫塔法求解微分方程！�。�！
已經(jīng)有9人回復(fù)
- 【討論】【【四階龍格-庫塔法求解微分方程】】
已經(jīng)有8人回復(fù)
- 【求助】求解一個二階非線性微分方程x''+Ax-B/(x^3)=C
已經(jīng)有6人回復(fù)
- 【求助】急請微分方程高手看下這個方程組！
已經(jīng)有14人回復(fù)
- 【求助】請教二階微分方程的特解（含初始條件）
已經(jīng)有8人回復(fù)
- 【求助】四階龍格-庫塔方法�。。�！
已經(jīng)有17人回復(fù)
- 【求助】20金幣求一個微積分方程
已經(jīng)有15人回復(fù)
wangyanb
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 96
- 散金: 20
- 帖子: 4
- 在線: 3.6小時
- 蟲號: 3921757
- 注冊: 2015-06-12
- 性別: GG
- 專業(yè): 激光
2樓2016-03-31 20:58:55
魔愛閻
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1315.5
- 帖子: 6
- 在線: 15.4小時
- 蟲號: 3927155
- 注冊: 2015-06-16
- 專業(yè): 光學(xué)
3樓2018-03-22 15:04:23
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 mia1129 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305分求調(diào)劑（食品工程） +3 | Sxy112 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+3
|
酥酥魚.. 2026-03-21 | 3/150 |
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +4 | 我愛生物生物愛?/a> 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 | 4/200 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|