| 查看: 746 | 回复: 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 |
回复此楼» 猜你喜欢
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有7人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有15人回复
-
多组分精馏求助
已经有4人回复
-
免疫学博士有名额，速联系
已经有4人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有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