matlab四重积分

专栏收录该内容
7 篇文章 0 订阅

这是在用matlab做光学衍射模拟的时候衍生出的一个问题

clear all;
f = 'x^2+y^2+m^2+n^2'

M = 10;
N = 10;
A = zeros(M,N);
for m = 1:M
    for n = 1:N
        ftmp = f;
        ftmp = subs(ftmp,'m',num2str(m));
        ftmp = subs(ftmp,'n',num2str(n));
        f1 = @(x,y) eval(ftmp);
        A(m,n) = dblquad(f1,0,1,0,1)
    end
end

 

或者

clear all;
sum = 0;
f2 = 'x^2+y^2+m^2+n^2'
for n = 1:10
    f = f2;
    f = subs(f,'n',num2str(n));
    f1 = @(x,y,m) eval(f);
    sum = sum + triplequad(f1,0,1,0,1,0,1)
end

  • 1
    点赞
  • 0
    评论
  • 0
    收藏
  • 一键三连
    一键三连
  • 扫一扫,分享海报

相关推荐
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值