Question

急急急！关于用matlab计算矩阵或编程计算矩阵。

一个矩阵是一个上三角矩阵与其转置矩阵的乘积,已知这个矩阵，然后如何用matlab编写函数或任意程序编写计算上三角矩阵。请大家注意是上三角矩阵与其转置矩阵而不是下三角。求高人帮忙啊

Answer 1

用 Maple 可以做如下计算：

1、设A为上三角矩阵，B为A 的转置矩阵；

2、继续，设C为已知矩阵并赋值，计算A.B-C（此例中，由该已知矩阵可以算得整数的结果）；

3、求解该矩阵方程，得abcdef的值。

4、最后该例解得共有8组解，如上。

5、换成其他的已知矩阵也可以解，但解的形式相应的要复杂很多，例如将C[2,2]换成59后，解如下：

需要注意，并不是任意的已知C都有结果。

Answer 2

受不了，MATLAB就是干这些的呀！这些都是最基本的了！
建议看下书吧，很简单的，先写好你的矩阵，然后用矩阵计算符很快就出来了！
平时用MATLAB编程就是因为它处理大量矩阵很方便！
建议楼主不要这样想走捷径！！！！

Answer 3

猜你应该懂英文的吧？

 

你的题目似乎和cholesky分解很像但其实和他没关系。

 

正确的是Kalman Filter有关。你可以参照这里：

"http://math.stackexchange.com/questions/110544/cholesky-factorisation-p-llt-rtr-vs-p-uut-ltl/110562#110562"，顺便说一句如果有问题去那问比较好，那里的人都比较厉害。

 

另外，你说的程序其实在上面的网址里提到了，我给你下载了下，这个包里面你只需要triag这个函数。

 

用的时候时和chol一起使用就可以解决问题。给你一个例子。假设你的矩阵是这样的

G =

     1    -1    -1    -1    -1
    -1     2     0     0     0
    -1     0     3     1     1
    -1     0     1     4     2
    -1     0     1     2     5，

你要的上三角阵用U=triag(chol(G)')得到是这样的 

U =

   -0.1078    0.7071    0.4194    0.3354    0.4472
         0       -1.4142         0             0             0
         0             0       -1.6394   -0.3354   -0.4472
         0             0             0       -1.7889   -0.8944
         0             0             0             0       -2.2361

你可以自己验证下。