Question

我有一个程序可以通过w得到一个矩阵3(n+1)，3(n+1)。我有一个fbar的向量3(n+1)，1。我希望得到一个矩阵，当w乘以它时，会得到fbar。

在数学符号中，w * A = fbar。我有w和fbar，我想要A。

我试图用这个命令解决它：

fsolve({seq(multiply(w, A)[i, 1] = fbar[i, 1], i = 1 .. 3*(n+1))})

但我不理解Maple给出的回应：

fsolve({2.025881905 A1[2,1]+7.814009150 A1[3,1]+...  

-7.071067816 10^(-13) A1[3,1]-0.0004999999990 

A1[4,1]-0.0007071067294 A1[5,1]-0.0004999999990 A1[6,1]) 

A3[6,1]=0},{A1[1,1],A1[2,1],A1[3,1],A1[4,1],A1[5,1],A1[6,1],A\

2[1,1],A2[2,1],A2[3,1],A2[4,1],A2[5,1],A2[6,1],A3[1,1],A3[2,1]\

,A3[3,1],A3[4,1],A3[5,1],A3[6,1]})

这是什么意思，我怎样才能得到更有意义的答案？

Answer 1

如果LinearSolve和LinearAlgebra分别定义为矩阵和向量，则可以使用w包中的fbar功能直接执行此操作。下面的代码是一个可重现的例子。请注意，LinearSolve的解决方案应该等于x。

w := Matrix(<<1,2,3>|<4,5,6>|<7,8,10>>);
LinearAlgebra[ReducedRowEchelonForm](%); ## Full rank => 1 solution)
x := <1,2,3>;
fbar := w.x;


## Solve the equation w.x = fbar
LinearAlgebra[LinearSolve](w,fbar);

如何解决枫树中的矩阵方程

1 个答案: