| 查看: 736 | 回复: 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 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求大神指点怎么用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