| 查看: 999 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab 三重数值积分 已有1人参与
|
||
|
自己按这个公式写了一个三重积分，用的int函数，但一直warning。 Z1u=sqrt(2*Rs*r1*cos(theta1)-r1^2); r1l = sqrt(R(j)^2-(Rp^2*sin(theta2)^1)-Rp*cos(theta1); r1u = 2*Rs*cos(theta1); theta1u = acos(sqrt((R(j)^2-Rp^2)/(4*Rs*(2*Rp)))); V(j) = 4*int(int(int(1,z1,0,Z1u)*r1,r1,r1l,r1u),theta1,0,theta1u); 希望大神能给一个三重数值积分的代码。 多谢！  IMG_1282(20170312-161908).jpg |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大神指点怎么用MATLAB求这个三重积分
已经有0人回复
4楼2017-03-12 19:46:52
tuburfransth
金虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2380.1
- 散金: 68
- 红花: 5
- 帖子: 384
- 在线: 124.8小时
- 虫号: 3425987
- 注册: 2014-09-18
- 性别: GG
- 专业: 模式识别
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2017-03-12 22:06:45
jindan8865: 金币+5, ★有帮助 2017-03-22 09:03:01
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2017-03-12 22:06:45
jindan8865: 金币+5, ★有帮助 2017-03-22 09:03:01
|
syms Rs r1 theta1 Rn Rp z1u = sqrt(2*Rs*r1*cos(theta1) - r1^2); r11 = sqrt(Rn^2 - Rp^2*(sin(theta1)^2)) - Rp*cos(theta1); r1u = 2*Rs*cos(theta1); theta1u = acos(sqrt((Rn^2 - Rp^2)/(4*Rs*(Rs + Rp)))); Vn = int(int(int(r1, r1, 0, z1u), Rs, r11, r1u), Rn, 0, theta1u) #-------------------------------------------------------------------- Warning: Explicit integral could not be found. Vn = int(- (Rp*cos(theta1) - (Rn^2 - Rp^2*sin(theta1)^2)^(1/2))*(r1^2/2 + (cos(theta1)*(Rp*cos(theta1) - (Rn^2 - Rp^2*sin(theta1)^2)^(1/2))*r1)/2) - 2*Rs*cos(theta1)*(r1^2/2 - Rs*r1*cos(theta1)^2), Rn == 0..acos(((Rn^2 - Rp^2)/(4*Rs*(Rp + Rs)))^(1/2))) #---------------------------------------------------------------------- 看结果，貌似无法得到显式积分结果 |
2楼2017-03-12 18:02:37
3楼2017-03-12 19:44:42
200