求二维数组的第二个积分

编程问答 2022-05-30

问题描述

我有以下二维数组：

N = 6

M = np.zeros((N,N))
for i in range(N-1):
    for j in range(N-1):
        if i==0 or j==0:
            M[i,j] = 0
        else:
            M[i,j] = 1

基本上，它是边上的 0 和中间的 1 的数组。

我想计算以下积分：

enter image description here

所以，我使用下面的代码：

dx=dy=d=1
My,Mx = np.gradient(M,d,d)

Lx = 0
p = 1/(M.shape[0]*M.shape[1])
for i in range(N-1):
    for j in range(N-1):
        Lx += P*math.sqrt(1+(Mx[i,j])**2)*dx*dy

但是，我从计算中得到的值与应有的值不匹配： 我的答案：0.77 正确答案：0.98

我做错了什么吗？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

numerical-methodsnumpypython