| 查看: 631 | 回复: 0 | ||
[求助]
求高手指点关于matlab数值积分问题
|
|
数值积分时在quadgk那儿出错，但是一直不知道是为什么，求大家帮忙看一看 主程序： 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;%吸收系数，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)]; |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有180人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有112人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有69人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
找到一些相关的精华帖子，希望有用哦~
帮我画出这个带有积分的函数的图形。
已经有7人回复
- 一个积分，总是算不出来
已经有4人回复
- MATLAB高手帮我看看这个数据存的格式，谢谢
已经有3人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
已经有9人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 关于Matlab画三维曲面图的问题，求助高手，感激感激
已经有3人回复
- matlab数据拟合求助！！希望高手帮忙写一下matlab的程序 谢谢了！！
已经有15人回复
- 用matlab求一个积分。我求不出来啊
已经有3人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- 求助啊......matlab
已经有10人回复
- matlab 气泡边缘提取 求助高手！！
已经有5人回复
- 请教matlab高手 求定积分
已经有3人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab求数值积分？
已经有10人回复
- 【求助】拜求高手解释下matlab的报错-----“out of memory”
已经有13人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 【求助】【求助】mathematica中的一种数值积分问题
已经有5人回复
- 【求助】请教Matlab高手求定积分近似值问题
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
