| 查看: 2351 | 回复: 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   |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有14人回复
-
范进中举一文的中心思想
已经有5人回复
-
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
散金币
+1/929
-
面上第三年申请，改了一遍又一遍，希望能中啊！散金祈福！！！
+5/520
-
每年这时候我都再仔细读一遍《范进中举》
+1/472
-
上海交通大学张航课题组招聘博士后（电化学储能）
+1/89
-
中国药科大学博士，博后招聘（高分子，金属有机方向），欢迎大家交流推荐
+1/84
-
那年春暖花开，今朝如梦初醒 —— 北漂十五载，相亲十年，归乡孑然一身
+1/67
-
焦虑
+1/64
-
面上基金祈福
+1/42
-
强耀生物金秋开学季活动，2026年8月18日-10月31日，肽库筛选低至2W5
+3/37
-
瑞典-博后-固态电池-粘结剂-锂离子电池方向
+1/25
-
南京理工大学优青团队催化化学方向招收推免研究生
+1/22
-
上海大学王亮课题组诚招博士后（碳材料、光催化、电催化方向）
+2/20
-
流放岭南浮肿了求救
+1/20
-
北理工集成电路杰青团队 | 诚招科助理
+1/18
-
中科院大连化学物理研究所 招聘 催化剂研发方向 科研助理 2名
+1/5
-
南京大学 智能驱动与感知材料实验室 诚招申请考核博士生/科研助理/博士后
+1/4
-
祈祷老公中国基
+1/4
-
出全新蔡司 Plan?APOCHROMAT 20×/0.8 平场复消色差物镜
+1/2
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/1
-
【产品测评分享】体外转录产量不理想、dsRNA 副产物高？T7 RNA 聚合酶实测体验
+1/1
5楼2011-03-22 23:24:52
2楼2011-03-22 21:18:04
:tiger01
*^__^*) 嘻嘻…… 不会呢！！我帮你转转吧！！。。。 3楼2011-03-22 21:25:12
4楼2011-03-22 22:30:25
