| 查看: 934 | 回复: 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 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大神指点怎么用MATLAB求这个三重积分
已经有0人回复
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
4楼2017-03-12 19:46:52
20
