| 查看: 2091 | 回复: 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   |
回复此楼» 猜你喜欢
博士招生
已经有3人回复
-
夜，静悄悄的
已经有8人回复
-
国自科送审了吗
已经有8人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有5人回复
-
研究生做的很差，你们会让毕业吗？
已经有6人回复
-
26年博士申请自荐-电催化
已经有5人回复
-
2026博士或科研助理转27年博士
已经有5人回复
-
考博
已经有6人回复
-
26年申博自荐-计算机视觉
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
东南大学废弃碳资源利用-杰青团队-27级博士生(提前入学)-及博后、科研助理招募
+1/176
-
国自然申报交流祈福
+4/134
-
东南大学废弃碳资源利用团队（杰青/海优团队）27级博士研究生和科研助理招生
+1/127
-
2026年江西师范大学人文地理与低碳经济方向拟招收1名“申请入学制”学术博士研究生
+2/110
-
中科大 冯伟 课题组招生-机器人界面控制、电学机械界面设计
+1/83
-
太行山徒步之凌水河（一）
+1/79
-
松山湖材料实验室-大连理工大学联合招收博士研究生
+1/38
-
松山湖材料实验室-大连理工大学联合招收2027年博士研究生调剂考生若干名
+1/37
-
西北工业大学物理科学与技术学院-诚聘【AI与光学交叉方向】博士后
+1/30
-
2026年江西师范大学药学院陈芬儿院士课题组招收智能药学博士生
+1/29
-
Google ai pro 会员分享
+1/28
-
坐标南京
+1/20
-
暨南大学中医学院邓力课题组博士后招聘【中西医结合/交叉学科方向】
+1/17
-
大连理工大学张硕课题组2026年秋季博士生额外招生（补招1人，有机合成/糖化学方向）
+1/14
-
又是一年基金申报，焦虑，某些学术大佬能不能少点干预
+1/10
-
李老师课题组急招收生物信息学方向博士研究生！
+1/8
-
成都中医药大学药学院 谢恬教授招收2026级药物化学，有机合成方向的博士生
+1/5
-
南京大学医学院附属鼓楼医院院士团队招聘博士后、科研助理
+1/2
-
长安大学隧道通风与爆破课题组招聘科研助理
+1/2
-
请问有推荐的博后流动站么 / 博士电子信息专业 /Serverless + AI Agent 领域
+1/2
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