| 查看: 2599 | 回复: 5 | ||
[求助]
MATLAB作不定积分并作图问题 已有1人参与
|
|
现有变量r1、r2，其中r1和r2关系如附图，f1是与r2和r1有关的表达式，f为积分式求z(r)并作图 程序如下： %%参数设定 RFL=3.153;w0=2.366;q=15;d=150;n=1.46071; syms r1 r f1 r2=RFL*(sqrt(1-(exp(-2*(r1/w0)^2)))/(1-(1-exp(-2*(r1/w0)^2))^(q/2))^(1/q)) f1=((n^2-1)+((n-1)*d/(r2-r1))^2)^(-0.5) r=0:0.5:5 I1=int(f1,r1,0,r) I1 = subs(I1) plot(r,I1) 出现警告及报错： Warning: Explicit integral could not be found. > In sym.int at 58 I1 = int(1/(2552806235672245/2251799813685248+23648162291281931128408331705169/4951760157141521099596496896/(3153/1000*(1-exp(-500000/1399489*r1^2))^(1/2)/(1-(1-exp(-500000/1399489*r1^2))^(15/2))^(1/15)-r1)^2)^(1/2),r1 = 0 .. 1/70704950625) ??? Error using ==> plot Conversion to double from sym is not possible. 此时I1仍是sym类型，所以plot不能作图？之前求助过有人说是不存在解析解的原因，那么如何求数值解？应该如何写这个程序？   公式1.jpg  公式2.jpg  公式3.jpg |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有87人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有86人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
2楼2015-10-09 10:06:17
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★
yujianmiao: 金币+1, ★★★★★最佳答案, 不好意思，新人金币少~ 2016-01-11 20:33:40
yujianmiao: 金币+1, ★★★★★最佳答案, 不好意思，新人金币少~ 2016-01-11 20:33:40
|
r不能等于0吧，会出现分母为0的错误： r I1 0.1 6.39433504405824E-5 0.2 0.000255044114661227 0.3 0.0005711252881483 0.4 0.0010085944724409 0.5 0.00156249693209275 0.6 0.00222658819034182 0.7 0.00299342456499068 0.8 0.00385446979521733 0.9 0.00480021603896343 1 0.00582031837637208 1.1 0.0069037437392736 1.2 0.00803893744028246 1.3 0.00921401173579675 1.4 0.0104169591403217 1.5 0.0116358874438529 1.6 0.0128592648683396 1.7 0.0140761561854483 1.8 0.015276427896934 1.9 0.0164509044546538 2 0.0175914663258592 2.1 0.018691090782268 2.2 0.0197438441653592 2.3 0.0207448384198801 2.4 0.0216901651267359 2.5 0.0225768183553321 2.6 0.0234026147591895 2.7 0.0241661164422961 2.8 0.0248665597374679 2.9 0.0255037913120838 3 0.0260782118988047 3.1 0.0265907273069974 3.2 0.0270427060627266 3.3 0.0274359429312615 3.4 0.0277726276059436 3.5 0.0280553179386943 3.6 0.0282869172001222 3.7 0.0284706549670885 3.8 0.0286100713299425 3.9 0.0287090041853041 4 0.0287715794329358 4.1 0.0288022039294651 4.2 0.0288084985085687 4.3 0.0288279728514223 4.4 0.0288633492488579 4.5 0.0289134770958946 4.6 0.0289677422953553 4.7 0.0290109719790368 4.8 0.0290698384095189 4.9 0.029085233650094 5 0.0291170612263809 |
3楼2015-10-09 21:34:58
4楼2015-10-19 15:34:36
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2015-10-19 16:12:49
6楼2015-10-19 17:10:58
30