| 查看: 1979 | 回复: 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   |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有7人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
坐标广州，征女友
+2/140
-
天津科技大学海洋与环境学院殷焕顺团队招博士生1名---分析化学领域
+1/83
-
ChineseResearchLaTeX: 开源、免费的vibe coding辅助国自然写作
+1/81
-
国家级青年人才课题组招收2026级硕士研究生
+1/31
-
海南大学海洋技术与装备学院-科研助理招聘（可读博）--膜分离水处理方向
+1/30
-
海南大学海洋技术与装备学院-科研助理招聘（可读博）膜分离水处理方向
+1/28
-
清华大学深圳国际研究生院招聘-博士后（长期有效）
+1/28
-
2026年天津科技大学“新能源催化与膜材料团队”研究生招生
+1/18
-
太原理工大学集成电路学院招收2026年博士研究生
+1/14
-
英国南安普顿大学禅铎课题组诚招气候动力方向博士后
+1/11
-
上海交通大学-宁波东方理工大学联合培养博士生
+1/10
-
湖南大学-分析检测技术和生物柔性传感器-招收1名博士研究生 (2026年，第二批)
+1/7
-
墨尔本大学(QS13)招全奖博士、CSC资助博士/访问学者(生物医学材料/器官芯片等方向)
+1/6
-
中北大学冯瑞教授*开山大弟子*招募
+1/5
-
26申博自荐求博导-生物传感分析方向
+1/5
-
怎么发布了求助贴了， 一发就转到删除栏了
+1/4
-
澳科大招收2026年秋季药物递送/生物材料方向硕士研究生（3月5日18:00报名截止）
+1/3
-
中国地质大学（武汉）杨华明课题组刘磊研究员招收冶金固废资源高值化利用方向博士
+1/1
-
西电科大招博士和科研助理： 结构生物学/纳米孔/分子动力学
+1/1
-
上海大学生物有机电子材料及器件团队博士研究生招聘
+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