是否有一种简单的方法来填充矩阵f(K)*delo(K+1x,L)+f(K)*delo(K+1y,L)
其中K=(k1,k2), K+1x=(k1+1, k2), K+1y=(k1, k2+1), L=(l1,l2)是R ^ 2和
def delo(i,j):
if i==j:
return 1.0
else:
return 0.0？
在我有一个更简单的案例（例如f2(k1)*g(k2)*delo(K+1x,L)+f2(k1)*g(k2)*delo(K+1y,L)）之前，我可以将其分成np.kron(f2(k1)*delo(k1+1,l1),g(k2)*delo(k2,l2))+np.kron(f2(k1)*delo(k1+1,l1),g(k2)*delo(k2,l2))，但是现在f(K)是不可分割的。
我尝试制作一个for循环函数来填充矩阵，然后再次使用另一个for循环，用先前的矩阵填充新的矩阵项。从数学上讲，这就是解决问题的方式。
一旦delo矩阵被填充后，其结果将类似于定义函数
def G(x):
return np.kron(f(k1,k2,x)delo(k1+1,l1),delo(k2,l2))，其中x是任意变量？
抱歉，如果不清楚，我不知道如何在此处输入乳胶。