| 查看: 2246 | 回复: 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   |
回复此楼» 猜你喜欢
职称论文投稿
已经有9人回复
-
国自然申请五篇代表作大比拼，感觉这个是最重要的
已经有9人回复
-
无聊看看时间戳打发时间
已经有3人回复
-
基于自然哲学类比的风化壳型稀土矿
已经有14人回复
-
评委有多少概率知道其他专家手中有哪些人的本子？
已经有6人回复
-
求推荐期刊，重谢
已经有3人回复
-
中！中！中！
已经有4人回复
-
E0414, 我的本子有没有希望？
已经有17人回复
-
青A35岁以下通知答辩了吗
已经有4人回复
-
小城的小雨
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
代友征（上海）-江浙沪独女生，等一个年下的他
+3/411
-
荷兰格罗宁根大学F. Maresca教授新开放一个全奖PhD位置（沉淀强化+多尺度建模）
+1/98
-
温州医科大学李校堃院士团队宋林涛/黄志锋课题组诚聘博士后
+1/75
-
青年拔尖人才还有指标 25w*8年
+1/74
-
头部外企，招聘力学专业方向仿真岗
+1/41
-
深职大霍夫曼研究院吸附分离课题组诚聘博士后
+1/35
-
上海985专排第一的软工学院课题组接收AI生物计算方向的推免保研硕士
+1/32
-
四川大学周加境课题组招聘博士后/博士/研究助理（生物质与藻类资源利用/自组装材料）
+1/30
-
哈工大深圳-材料学院-招收申请审核制2027年春季/秋季入学博士生（2026年9月报名）
+1/28
-
中山大学智能工程学院【空间智能方向】招收2027年入学博士生
+2/24
-
澳大利亚西澳大学招收交通工程/智能交通方向博士生
+1/13
-
希望今年的面上能上哇
+2/8
-
东北师范大学荒漠与草地生态学方向诚聘博士后青年才俊
+1/8
-
广工-董华锋教授团队招收博士生（1学博-0-1专博）
+1/7
-
北理工集成电路杰青团队 | 诚招科助理
+1/5
-
从分子模拟到科学大模型：加入深圳理工大学新一代AI4Science团队
+1/5
-
密苏里大学生物材料合成生物学博士后招聘
+1/5
-
密苏里大学生物材料合成生物学博士后招聘
+1/4
-
密苏里大学生物材料合成生物学博士后招聘
+1/4
-
哈工大马樱教授招收2027级计算机类、集成电路类博士生
+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