MATLAB:从矩阵转换为逐元素乘法
我有这样的代码:
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) 。 有什么想法吗?
I have this code:
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);
The matrix KM dynamically changes the shape of potencial V. The Fx I use in ode solver.
What I want is to be able write Fx(A,B), where A,B is matrix.
I need to convert result of V=(x^2+y^2) to V=(x.^2+y.^2).
Any idea?
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论
评论(1)
我不太明白这个问题。不是
Fx = @(x,y) -2*KM(1,1)*x-(KM(2,1)+KM(1,2))*y;
当您手工计算。一般V为
V = KM(1,1)*x^2+(KM(1,2)+KM(2,1))*x*y+KM(2,2)*y^2
这只是一个二次形式。另外,如果是正定且对称的,您可以稍微简化一下。I don't quite understand the question. Isn't
Fx = @(x,y) -2*KM(1,1)*x-(KM(2,1)+KM(1,2))*y;
when you do the math by hand.The general V is
V = KM(1,1)*x^2+(KM(1,2)+KM(2,1))*x*y+KM(2,2)*y^2
which is just a quadratic form. Also if is positive definite and symmetric you can simplify it a bit.