版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>基金申请 (7625)
>论坛更新日志 (4946)
>休闲灌水 (1389)
>导师招生 (334)
>文献求助 (326)
>虫友互识 (320)
>论文投稿 (305)
>有奖起名 (165)
>硕博家园 (130)
>考博 (115)
>健康生活 (89)
>考研 (42)
>博后之家 (41)
>留学DIY (39)
>绿色求助(高悬赏) (35)
>公派出国 (33)

返回列表

tang3117022

金虫 (小有名气)

应助: 6 (幼儿园)
金币: 1058.7
帖子: 295
在线: 376.4小时
虫号: 1809326
注册: 2012-05-10
专业: 交叉学科中的光学问题

[求助] Matlab多重积分数值计算已有1人参与

求助：

我想用matlab计算一个多重积分（目前是8重积分），用数值方式求解。

我想到一个思路，把一个8维的矩形空间区域均匀划分一下，分成许多个8维小区域，再计算求和。

但是不知道matlab有没有函数能快速划分多维区域。

另外，大家还有没有更好的思路？

谢谢！

回复此楼

» 猜你喜欢

有多少人是今天查系统知道结果的？已经有16人回复
基金不中，共勉已经有3人回复
怎么查啊已经有5人回复
范进中举一文的中心思想已经有10人回复
我不理解！已经有13人回复
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战已经有9人回复
2026国自然函评费到账已经有25人回复
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？已经有28人回复
能否退出参与的面上项目解除限项已经有26人回复
出来了已经有9人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-10-26 22:08:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wwhblue

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 197.4
帖子: 45
在线: 29.1小时
虫号: 685949
注册: 2009-01-01
性别: GG
专业: 结构工程

【答案】应助回帖

给你三重辛普森复化积分算法，8重积分依次自己补充吧

% 复化辛普森公式3D
function s = simpson3D(fname,a,b,h1,c,d,h2,e,f,h3)
fe = simpson2D(fname,a,b,h1,c,d,h2,e);
ff = simpson2D(fname,a,b,h1,c,d,h2,f);
s = fe - ff;
x = e;
n = ceil((f-e)/h3);
for k=1:n
x= x + h3/2; s = s + 4* simpson2D(fname,a,b,h1,c,d,h2,x);
x= x + h3/2; s = s + 2* simpson2D(fname,a,b,h1,c,d,h2,x);
end
s = s * h3 / 6;
return;

% 复化辛普森公式2D
function s = simpson2D(fname,a,b,h1,c,d,h2,zi)
fc = simpson1D(fname,a,b,h1,c,zi);
fd = simpson1D(fname,a,b,h1,d,zi);
s = fc - fd;
x = c;
n = ceil((d-c)/h2);
for k = 1:n
x= x + h2/2; s = s + 4* simpson1D(fname,a,b,h1,x,zi);
x= x + h2/2; s = s + 2* simpson1D(fname,a,b,h1,x,zi);
end
s = s * h2 / 6;
return;

% 复化辛普森公式1D
function s = simpson1D(fname,a,b,h1,yi,zi)
%fname为被积函数，a,b分别为下界和上界，h1步长
fa = feval(fname,a,yi,zi);
fb = feval(fname,b,yi,zi);
s = fa - fb;
x = a;
n = ceil((b-a)/h1);
for k = 1: n
x= x + h1/2; s = s + 4* feval(fname,x,yi,zi);
x= x + h1/2; s = s + 2* feval(fname,x,yi,zi);
end
s = s * h1 / 6;
return;

赞一下