| 查看: 2301 | 回复: 7 | |||
[交流]
【求助】matlab求解一个简单的积分（急！！！）
|
|
大家好，在matlab方面我是个很业余的选手。 对于以下的方程，有人用蒙特卡罗法解出来了，但我看看了蒙特卡罗的方法，感觉还是比较有难度的，所以最近一直在用matlab求解如下的方程，自己写出了求解的言语，虽然得出了结果，但是结果中有警告，不知结果对不对，请大家帮忙？式中β=0.0015，γ=0.001（程序中β用i代替，γ用r代替） syms rb H i k ts t t1 r a b c d e f g h DA DB A B G H=100 rb=0.15 i=rb/H k=1.2.*10.^(-6) ts=H.^2/(9.*k) t=24.*3600 t1=t./ts r=3./2.*sqrt(t1) a=sqrt((i.^2)+1) b=sqrt((i.^2)+4) c=exp(-r.^2.*((i.^2)+1)) d=exp(-r.^2.*((i.^2)+4)) e=exp(-r.^2.*i.^2) f=erfc(r.*i) g=erfc(r.*b) h=erfc(r.*a) DA=a.*h-i.*f-(c-e)./(i.*sqrt(3.1415)) DB=a.*h-0.5.*(r.*f+b.*g)-(c-0.5.*(e+d))/(i.*sqrt(3.1415)) A=quadl('erfc(0.0145.*z)./sqrt(z.^2-0.0015.^2)',0.0015,sqrt(0.0015.^2+1)) B=quadl('erfc(0.0145.*z)./sqrt(z.^2-0.0015.^2)',sqrt(0.0015.^2+1),sqrt(0.0015.^2+4)) G=A-B-DA-DB   |
回复此楼» 猜你喜欢
面上提前没消息，有中的吗
已经有16人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有25人回复
-
2026年国自然面上资助率
已经有22人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有12人回复
-
这种情况还有戏吗
已经有12人回复
-
微信指数没变化，科研之友没阅读
已经有19人回复
-
准备明年的基金了
已经有7人回复
-
2027年申博
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有7人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
坐标北京，92男，诚征结婚对象
+1/166
-
国家纳米科学中心段鹏飞课题组招收2027年入学博士研究生
+1/80
-
★★★中山大学（深圳）院士团队招收27届推免硕士生（低维材料与器件方向）
+1/78
-
华东理工大学 生物工程学院 生物反应器工程国家重点室 博士后招聘
+1/76
-
Postdoctoral fellow in Medicinal Chemistry, McGill University
+1/70
-
采用本质安全型氢气流量控制的化学合成 - Alicat ISMC本安防爆质量流量控制器
+2/56
-
西湖大学物理学、光学、电子信息方向博士生招生（长期有效）
+5/25
-
西湖大学拓扑光学、非厄米光学、太赫兹方向博士后招聘（长期有效）
+5/15
-
招生宣讲会倒计时！7月29日20：00，不见不散！
+1/9
-
北理工集成电路杰青团队 | 诚招科助理
+1/7
-
西安交通大学补亚忠课题组招收申请考核制博士生
+2/6
-
中科院动物所招收27届推免或报考硕士生（生物或化学背景，北京）
+1/6
-
北理工集成电路杰青团队 | 诚招科助理
+1/6
-
国家自然科学基金申报者发榜前心理状态调查，麻烦帮忙填一下哇！
+1/4
-
新西兰坎特伯雷大学招博士后
+2/3
-
准现货/原装全新超微整机8x NVIDIA HGX A100 80G
+1/3
-
2027年硕士推免生 招生--哈尔滨工业大学王磊
+1/3
-
（Scopus期刊）Biofunctional Materials《生物功能性材料》招募青年编委及征稿
+1/3
-
国自然发榜在即：此刻的你，是什么心情？
+1/1
-
致每一位等待国自然结果的科研人：我们想听听你的心声
+1/1
2楼2011-03-22 21:18:04
:tiger01
*^__^*) 嘻嘻…… 不会呢！！我帮你转转吧！！。。。 3楼2011-03-22 21:25:12
4楼2011-03-22 22:30:25
5楼2011-03-22 23:24:52
6楼2011-03-23 07:00:59
7楼2011-03-23 07:40:34
8楼2012-11-08 15:56:28