m0_61305569 2021-08-28 09:03 采纳率: 50%
浏览 81
已结题

三重累加求和数值积分怎么编程呢

最近遇到一个棘手的问题,面对这个三重的如何进行数值积分编程?
之前只进行过一重的数值积分代码如下

clear
clc
n=20;
d=0.1;
x=-0.2:0.0001:0.2;
N=length(x);
sum=0;
%trapz方法
for i=0:n
    for ii=1:N
        y=x(ii);
        f(ii)=0.9.*(%函数)
    end
    S=trapz(x,f);
    sum=sum+S;
end
sum

img

  • 写回答

1条回答 默认 最新

  • joel_1993 2021-08-28 10:39
    关注

    这个最好用integral3函数

    % 定义被积函数,注意用的是对应位置相乘的结果所以乘法、幂次和除法要加.号
    F = @(u,w,v,i,d,n) 0.9*(d+0.15*i/n-w)./((v-u).^2+(d+0.15*i/n-w).^2).^(3/2)...
    - 0.9*(d+0.15*i/n-w)./(0.025*0.025 + (v-u).^2 + (d+0.15*i/n-w).^2).^(3/2);
    u1 = -0.022; u2 = 0.022; % 积分区间
    w1 = -0.1; w2 = 0;% 积分区间
    v1 = -0.2; v2 = 0.2;% 积分区间
    n = 20;
    i_arr = 0:n;
    d = 0.1;
    S = 0;% 求和初值
    for i = 1:1:numel(i_arr)
    S = S + integral3(@(u,w,v)F(u,w,v,i_arr(i),d,n),u1,u2,w1,w2,v1,v2);
    end
    S
    

    结果

    
    S =
    
       0.012099286633777
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 9月5日
  • 已采纳回答 8月28日
  • 创建了问题 8月28日

悬赏问题

  • ¥20 数据排序,可选择排序方向
  • ¥15 修改一下代码,考虑进程到达时间不同的情况
  • ¥15 华为nova10pro ,关闭热动热点用量记录
  • ¥15 帮我看看这是个啥题,带解题过程和结果,条件如下FCF = 290471.33 g1 = 15% r = 8% g2 = 4% n = 5
  • ¥15 edem模拟颗粒不显示或者生成失败
  • ¥15 Python代码编写
  • ¥15 php 将rtmp协议转hls协议,无法播放
  • ¥20 python代码编写
  • ¥20 使用MPI广播数据遇到阻塞
  • ¥15 TinyMCE如何去掉自动弹出的“链接…”工具?