| 查看: 2448 | 回复: 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   |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
课题组会帮你通过考核么？
已经有4人回复
-
各位大佬，路线设计求助
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
2026年人文教育与社会发展国际会议（HESD 2026）
+1/888
-
求助DL 5009.1-2014
规程 第1部分：火力发电》
+1/487
-
爱丁顿实验
+1/37
-
高分辨质谱测试丨科研单位对外服务
+1/20
-
细胞检定样品，您用对了吗？
+1/3
-
委托检验的安全锁：质量保证协议，你签了吗？
+1/3
-
CAHSS（艺术、人文与社会科学国际会议）
+1/2
-
ICoEP（教育与心理学国际会议）
+1/2
-
ICoAAD（建筑、艺术与设计）
+1/2
-
GCHSD（全球文化遗产可持续发展大会）
+1/2
-
仪器验证全生命周期管理：检测数据精准与否的前提
+1/2
-
智能传感器网络实时数据处理：时序分析、DNN调度与MQTT边缘AI工作站配置
+1/1
-
清华大学-区域环境安全全国重点实验室招聘地下水科学/水文地质方向助理研究员/副研究�
+1/1
-
ICoFEE（林业、生态与环境）
+1/1
-
单细胞VDJ测序：从B细胞库中精准挖掘抗体药优势克隆
+1/1
-
CoIMCT（信息、媒体与通信技术）Informatics, Media and Communication Technology
+1/1
-
ECNLPIR（自然语言处理与信息检索）
+1/1
-
ECSFS（能源、气候与可持续发展峰会）
+1/1
-
CoEES（能源与地球科学）2027 International Conference on Energy and Earth Science
+1/1
-
新能源并网仿真：AI实现规程问答与试验报告生成，电磁暂态图形工作站硬件配置方案
+1/1
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