做矩阵分块解析运算（证明也解决了，软件是不存在的，更新在3页）

orionsnow

唉今天又出新问题了，我发现我昨天忽略掉一个参数，这个参数每次simulation 都变的，这就是说每次 那个逆矩阵都要重新求。 我印象里头我们的simulation 原来是1-5 秒，这样如果算上求逆的200 秒就是205 秒了。

继续求可以符号计算矩阵运算的软件

熊猫羊

怪蜀黍的英语真强悍啊

熊猫羊

怪蜀黍，不是鄙视你，你懂的太多了，但一样都不精，standford网页上那个东东是国内大一线数的内容
实在搞不定的话，短我吧，matlab程序，自己写的，200欧，保讲懂，要c++的话，400欧

庄十三爷

怪蜀黍，不是鄙视你，你懂的太多了，但一样都不精，standford网页上那个东东是国内大一线数的内容实在搞不定的话，短我吧，matlab程序，自己写的，200欧，保讲懂，要c++的话，400欧
熊猫羊 发表于 2010-11-1 14:39

400 欧太多了，叔叔还是短我吧， 8euro 88 cent 就可以了 就可以了

熊猫羊

庄13爷，您对一个连 矩阵分解 和 矩阵分块 都没搞清的冤大头也太仁慈了吧？

庄十三爷

不少了，你想阿，如果是在柏林踩三轮车给人家帮忙搬家，一个小时十快钱最多了，写个程序怎么也要一天8个小时吧。所以换算成欧元差不多就是这个价。

运气好还可以卖人家个沙发，赚个50 大欧， 不过这种事情不是每次都能碰上的。 再说还有被人家人肉的风险，我可不认识版主，弹压不下去。

农大

2000*2000 用 “解析" 是不合适的
而且所有所谓的数学软件matlab之类，对待这样的矩阵也不会用所谓的”解析“解
lz 要是会c或者c++ 随便找个solver 解一下就行。。。

orionsnow

上来汇报下。 花了一个多月， 两个超大矩阵在特殊情况下的解析解都被我给推导出来了。

然后我老板在网上搜了几天，居然找到了第一个矩阵解法的推导，  那人最后把特例朝他的一般推导里头一代入，就得到了我推出来的解析解。 根据这个事实和数值算法的验证，就是前边无数同学提到的， solver（） inverse（） 这类的函数。 确认第一个矩阵的解法正确。 然后老板自己验证了下第二个公式也是对的。

我现在在模仿那人的一般推导推第二个公式。

论文在这里，感兴趣的可以看了以后讨论下， 公式41-46 是主要证明， 结论在第5章，公式50多。

Searle SR (1996) The matrix handling of BLUE and BLUP in the mixed linear
model. Linear Algebra Appl 264:291-311 232

这个人挺牛的， 普通解推导用的矩阵计算， 特解带入用的矩阵分块运算，主要证明就两页纸。

我为了和老板说明我的思路前后写了40页纸。 现在正在好好学习他的矩阵分块算法的推导。

矩阵分块计算的推导软件应该是不存在的， 我查了这个人paper 的被引用情况，没有人提到更深入的解法了。

估计现在计算机技术还没有那么发达，只能求数值解。 如果计算机能做推导，那估计数学家可以改行了。

orionsnow

再顶起来一下， 下一篇文章还是关于这类公式推导的。现在老板都快把这个当产业了，把以前的老算法都拿出来， 设某个参数为特殊值，然后加上点解析化简，就变成新文章了。