| 查看: 2165 | 回复: 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   |
回复此楼» 猜你喜欢
宿州学院学报
已经有5人回复
-
青B发送上会通知了吗
已经有14人回复
-
博士申请
已经有5人回复
-
4,4二甲基联苯干啥用，有懂得吗
已经有4人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有7人回复
-
论文撤稿了
已经有9人回复
-
化学专业申博
已经有5人回复
-
招收2026级博士生
已经有5人回复
-
医学类期刊求推荐
已经有6人回复
-
26/27申博自荐
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 积分计算 算不出结果，求帮助！！！
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab求积分，出现奇点怎么办？
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
西安电子科技大学保研: 挑战高速IP设计NVMe，RDMA及NVMe-oF
+1/288
-
征结婚对象，男找女，坐标【北京】【北京】【北京】
+5/145
-
中科院二区期刊Gels 组稿
+1/81
-
中科院深圳先进技术研究院集成电路先进封装博士后招聘
+1/80
-
寻找做兼职的同学，长期的（长三角：血管介入医疗器械方面兼职）
+1/44
-
【急招】 北京工业大学 第二批博士申请-杰青课题组1-2名额，25日截止！！
+1/40
-
【通知】北京信息科技大学仪器科学与光电工程学院招收博士研究生（2026）
+2/36
-
【通知】北京信息科技大学仪器科学与光电工程学院招收博士研究生（2026）
+2/30
-
上海理工大学-赵斌教授课题组招收申请考核制博士【新能源材料】
+1/18
-
【26.9月入学】211大学补招学博1人
+1/15
-
郑州大学·安徽高等研究院 2026年博士招生（申请考核制）
+1/13
-
上海交通大学化学化工学院张智涛课题组诚聘博士后
+1/9
-
电子科技大学材料学院SFT创新中心招收准备考研和读博的科研助理 理工医交叉方向
+1/9
-
张兆威教授课题组长期招聘科研助理
+1/5
-
空间组学入门：GeoMx® DSP核心原理、实验流程与科研应用
+1/4
-
【皇家墨尔本理工大学】招大模型安全博士
+1/2
-
氢能专刊邀稿
+1/2
-
北航杭州国际校区招聘3D 打印，陶瓷，力学等方向博士后
+1/1
-
【有偿访谈招募】高才通来港后，你过得还好吗？
+1/1
-
【有偿访谈招募】高才通来港后，你过得还好吗？
+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