| 查看: 1922 | 回复: 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   |
回复此楼» 猜你喜欢
不自信的我
已经有6人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
要不要辞职读博？
已经有3人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
自荐读博
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
上海交通大学医学院王戈林课题组招聘博士后和助理研究员
+1/267
-
双一流南京医科大学招计算机、AI、统计、生物信息等方向26年9月入学博士
+1/177
-
江汉大学招聘AI for Materials/电解液/锂金属/全固态电池等方面的博士或者博士后
+1/170
-
留学导师避雷——望传播
+1/85
-
上海师范大学生物医用材料方向招收2026级博士研究生
+1/84
-
上海大学昝鹏教授、军事医学研究院伯晓晨研究员/倪铭副研究员 课题组招聘博士生
+1/80
-
湖南师范大学医工交叉科研团队招收计算机博士生
+1/80
-
北京-89175-事业单位-诚征女友
+1/79
-
澳门大学智慧城市物联网国重“结构智能感知、健康监测与无损检测”研究方向博士后招聘
+1/77
-
中科院长春光机所 招收计算材料学博士/硕士研究生（含机器学习辅助材料设计方向）
+1/74
-
山东征女友，坐标济南
+1/60
-
北京化工大学化学工程学院杨琪教授 邱介山教授，招收储能电池方向博士研究生
+1/44
-
中科院理化技术研究所张飞龙研究员/王树涛研究员团队招生（博士/硕士）
+1/35
-
征女友 @长安
+1/21
-
重庆医科大学-药学院-新靶标教育部医药基础研究创新中心-药物化学2026年博士招生
+1/20
-
中国科学院大学-杨晗课题组-诚聘-博士后、副研究员
+1/12
-
吉林大学材料物理本科生求问调剂信息
+1/11
-
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/8
-
液晶拓扑光子学博士招生（电子科技大学）
+1/6
-
代算！材料学理论计算
+1/3
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