0

我有这个代码:

KM = [1 0; 0 1];    
syms x, syms y;
% definition of force field 
V = [x y]*KM*[x ; y]
Vdx = -diff(V,x);
Fx = @(x,y) eval(Vdx);

矩阵 KM 动态改变电位 V 的形状。我在 ode 求解器中使用的 Fx。我想要的是能够写出 Fx(A,B),其中 A,B 是矩阵。我需要将 V=(x^2+y^2) 的结果转换为 V=(x.^2+y.^2)。任何想法?

4

1 回答 1

0

我不太明白这个问题。不是Fx = @(x,y) -2*KM(1,1)*x-(KM(2,1)+KM(1,2))*y;当你手工计算的时候。

一般的 VV = KM(1,1)*x^2+(KM(1,2)+KM(2,1))*x*y+KM(2,2)*y^2只是一个二次形式。此外,如果是正定且对称的,您可以稍微简化一下。

于 2011-01-29T23:24:23.927 回答