三维高斯积分,用两个坐标积分点和三个坐标积分点的代码
1条回答 默认 最新
- joel_1993 2021-10-20 14:46关注
你好,代码供参考
%两个积分点 P = [-1/sqrt(3), 1/sqrt(3)];%积分点 W = [1,1];%权函数 F = @(x,y,z)x^2+y^2+z^2;%比如说将F的xyz都从-1积分到1 s = 0; for i = 1:numel(P) for j = 1:numel(P) for k = 1:numel(P) s = s + F(P(i), P(j), P(k))*W(i)*W(j)*W(k); end end end s
结果:s =
8.000000
%三个积分点 P = [-sqrt(0.6), 0, sqrt(0.6)];%积分点 W = [5/9, 8/9, 5/9];%权函数 F = @(x,y,z)x^2+y^2+z^2;%比如说将F的xyz都从-1积分到1 s = 0; for i = 1:numel(P) for j = 1:numel(P) for k = 1:numel(P) s = s + F(P(i), P(j), P(k))*W(i)*W(j)*W(k); end end end s
结果:s =
8.000000
答题不易,有帮助望采纳呢本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 无用评论 打赏 举报
悬赏问题
- ¥15 c语言怎么用printf(“\b \b”)与getch()实现黑框里写入与删除?
- ¥20 怎么用dlib库的算法识别小麦病虫害
- ¥15 华为ensp模拟器中S5700交换机在配置过程中老是反复重启
- ¥15 java写代码遇到问题,求帮助
- ¥15 uniapp uview http 如何实现统一的请求异常信息提示?
- ¥15 有了解d3和topogram.js库的吗?有偿请教
- ¥100 任意维数的K均值聚类
- ¥15 stamps做sbas-insar,时序沉降图怎么画
- ¥15 买了个传感器,根据商家发的代码和步骤使用但是代码报错了不会改,有没有人可以看看
- ¥15 关于#Java#的问题,如何解决?