版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

jw19890719

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 161
帖子: 7
在線: 36.6小時
蟲號: 1795060
注冊: 2012-05-04
專業(yè): 光學(xué)

[求助] 求高手指點關(guān)于matlab數(shù)值積分問題

數(shù)值積分時在quadgk那兒出錯，但是一直不知道是為什么，求大家?guī)兔匆豢?br />
主程序：
y0 = [0.001;0.001];%y(1),y(2)初值
tspan=[0 0.005];%時間的范圍，單位為s
[t,y] = ode45('gauss_rq',tspan,y0);

matlab程序如下：
function dy = gauss_rq(t,y)
sigma = 7.6e-22;%單位m2
sigma_g = 4.3e-22;
sigma_e = 8.2e-23;
ns0 = 2e23;
tao = 90e-6;%單位s
tao_s = 3.2e-6;
wl = 200e-6;
n1 = 2.19;
n2 = 1.81;
l = 0.005;
ls = 0.001097;
lc=0.085;
wp = 330e-6;
L=0.08;
alpha = 532;%吸收系數(shù)，m-1
wg = 308e-6;
ws = 247e-6;

R=0.85;
T0 = 0.91;
Pin = 2.67;
hvp = 2.4585e-019;%808nm泵浦能量
c = 299792458;
tr = (l*n1+ls*n2+(lc-l-ls))/c;
syms r;
n00 = (log(1/(T0^2))+L)*(1+(wg/wp)^2)/(2*sigma*l);
Rin = Pin.*(1-exp(-alpha.*l)).*exp(-2.*r.^2./(wp).^2)./(hvp.*pi.*(wp).^2.*l)
%a=0.01;

xx = exp(sigma.*c.*((wl./wg).^2).*exp(-2.*r.^2./(wg.^2))).*quad(@(t)arrayfun(@(t)exp(y(2)),t),0,t)+t.^2./(2.*tao)
yy = exp(sigma_g.*c.*((wl./ws).^2).*exp(-2.*r.^2/(ws.^2))).*quad(@(t)arrayfun(@(t)exp(y(2)),t),0,t)+t.^2./(2.*tao_s)

n = exp(-sigma.*c.*((wl./wg).^2).*exp(-2.*r.^2./(wg.^2)).*y(2)-t./tao).*(Rin.*xx+n00.*exp(-2.*r^2/(wp^2)))
ns = exp(-sigma_g.*c.*((wl./ws).^2).*exp(-2.*r.^2./(ws.^2)).*y(2)-t/tao_s).*((ns0/tao_s).*yy+ns0)

zz = 2.*sigma.*n.*l.*(wl.^2./(wg.^2)).*exp(-2*r.^2/(wg.^2))-2*sigma_g.*ns.*ls.*(wl.^2/(ws.^2)).*exp(-2*r.^2./(ws.^2))-2*sigma_e.*(ns0-ns).*ls.*(wl.^2./(ws.^2)).*exp(-2.*r.^2./(ws.^2))-log(1./R).*((wl.^2./(ws.^2)).*exp(-2.*r.^2./(ws.^2))-L.*exp(-2.*r.^2./(wl.^2)).*r)
z=quadgk(@(r)zz,0,inf)

dy =[(4*y(1)/(wl^2*tr))*z;
y(1)];

回復(fù)此樓

» 猜你喜歡

請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有161人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有78人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)

1樓 2014-01-01 21:19:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

幫我畫出這個帶有積分的函數(shù)的圖形。已經(jīng)有7人回復(fù)
一個積分，總是算不出來已經(jīng)有4人回復(fù)
MATLAB高手幫我看看這個數(shù)據(jù)存的格式，謝謝已經(jīng)有3人回復(fù)
求助，用matlab解微分方程組，希望高手能給予指點已經(jīng)有17人回復(fù)
求高手解答Matlab求解微分方程的初始值和邊界值問題已經(jīng)有5人回復(fù)
請matlab高手看一下這個問題怎樣編程才好，附程序代碼，謝謝！已經(jīng)有9人回復(fù)
求Matlab高手解決線性方程組的迭代求解問題已經(jīng)有17人回復(fù)
關(guān)于Matlab畫三維曲面圖的問題，求助高手，感激感激已經(jīng)有3人回復(fù)
matlab數(shù)據(jù)擬合求助�。∠Ｍ呤謳兔懸幌耺atlab的程序謝謝了��！已經(jīng)有15人回復(fù)
用matlab求一個積分。我求不出來啊已經(jīng)有3人回復(fù)
高手進(jìn)來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
求助啊......matlab 已經(jīng)有10人回復(fù)
matlab 氣泡邊緣提取求助高手�。� 已經(jīng)有5人回復(fù)
請教matlab高手求定積分已經(jīng)有3人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab求數(shù)值積分？已經(jīng)有10人回復(fù)
【求助】拜求高手解釋下matlab的報錯-----“out of memory” 已經(jīng)有13人回復(fù)
【求助】請各位高手告訴我：這種想法能不能用matlab實現(xiàn)？200金幣求助已經(jīng)有16人回復(fù)
【求助】matlab求解一個簡單的積分（急�。。。� 已經(jīng)有7人回復(fù)
【求助】【求助】mathematica中的一種數(shù)值積分問題已經(jīng)有5人回復(fù)
【求助】請教Matlab高手求定積分近似值問題已經(jīng)有4人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 jw19890719 的主題更新