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

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

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

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日

悬赏问题

  • ¥15 nginx反向代理获取ip,java获取真实ip
  • ¥15 eda:门禁系统设计
  • ¥50 如何使用js去调用vscode-js-debugger的方法去调试网页
  • ¥15 376.1电表主站通信协议下发指令全被否认问题
  • ¥15 物体双站RCS和其组成阵列后的双站RCS关系验证
  • ¥15 复杂网络,变滞后传递熵,FDA
  • ¥20 csv格式数据集预处理及模型选择
  • ¥15 部分网页页面无法显示!
  • ¥15 怎样解决power bi 中设置管理聚合,详细信息表和详细信息列显示灰色,而不能选择相应的内容呢?
  • ¥15 QTOF MSE数据分析