你可能会足够近
雅可比矩阵
得到完整的输出:
信谊x y真实
%简单2 x 2的矩阵:
M = (sin (x + y) cos (x * y); cos (x / (x ^ 2 + y ^ 2)谭(x));
%应该返回:
% M =
%
% (sin (x + y), cos (x * y)]
% (cos (x / (x ^ 2 + y ^ 2)),谭(x))
%
%这将把2 x2矩阵M变成4 x1 column-array,雅可比矩阵
%可以处理:
J =雅可比矩阵(M (:), (x, y))
%应该返回
% J =
%
% (cos (x + y), cos (x + y)]
% (sin (x / (x ^ 2 + y ^ 2)) * (1 / (x ^ 2 + y ^ 2)——(2 * x ^ 2) / (x ^ 2 + y ^ 2) ^ 2)、- (2 * x * y * sin (x) / (x ^ 2 + y ^ 2))) / (x ^ 2 + y ^ 2) ^ 2]
% (- y * sin (x * y), - x * sin (x * y)]
% (tan (x) ^ 2 + 1, 0]
这里你得到所有4 x 2衍生品,你将不得不重新对行排序命令你期望。toy-model-example得到:
J3D =重塑(J [2 2 2]);
像你所说。在你的情况中这可能是一个相当大的数组的大小
n1 x x n_vars n2
——如果我理解你的问题。
HTH
