版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

jw19890719

主管區(qū)長

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

[求助] 求高手指點關于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)];

回復此樓

» 猜你喜歡

天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有270人回復
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8 已經(jīng)有0人回復
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院已經(jīng)有28人回復
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院已經(jīng)有23人回復
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學微結構與器件課題組。已經(jīng)有0人回復

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

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

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關的精華帖子，希望有用哦~

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

點擊這里搜索更多相關資源

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

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