首页 文章

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

提问于
浏览 次
0

我有一个程序做一些工作来获得矩阵 w ， 3(n+1) 由 3(n+1) . 我有 fbar 的矢量 fbar 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]})

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

maple

1 回答

  • 0

    如果 w 和 fbar 分别定义为矩阵和向量，则可以使用 LinearAlgebra 包中的 LinearSolve 函数直接执行此操作 . 下面的代码是一个可重现的例子 . 请注意， 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);
    回复于

相关问题