| 查看: 775 | 回复: 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 |
回复此楼» 猜你喜欢
材料与化工（0856）304求 B区 调剂
已经有3人回复
-
268求调剂
已经有8人回复
-
265求调剂
已经有14人回复
-
350 本科985求调剂，求老登收留
已经有3人回复
-
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有9人回复
-
初试 317
已经有3人回复
-
279分求调剂 一志愿211
已经有15人回复
-
286求调剂
已经有4人回复
-
279求调剂
已经有4人回复
-
求调剂
已经有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
