唉今天又出新问题了，我发现我昨天忽略掉一个参数，这个参数每次simulation 都变的，这就是说每次 那个逆矩阵都要重新求。 我印象里头我们的simulation 原来是1-5 秒，这样如果算上求逆的200 秒就是205 秒了。
继续求可以符号计算矩阵运算的软件
怪蜀黍的英语真强悍啊
怪蜀黍，不是鄙视你，你懂的太多了，但一样都不精，standford网页上那个东东是国内大一线数的内容
实在搞不定的话，短我吧，matlab程序，自己写的，200欧，保讲懂，要c++的话，400欧
本帖最后由 庄十三爷 于 2010-11-1 20:32 编辑
怪蜀黍，不是鄙视你，你懂的太多了，但一样都不精，standford网页上那个东东是国内大一线数的内容实在搞不定的话，短我吧，matlab程序，自己写的，200欧，保讲懂，要c++的话，400欧
熊猫羊 发表于 2010-11-1 14:39 http://www.dolc.de/forum/images/common/back.gif
400 欧太多了，叔叔还是短我吧， 8euro 88 cent 就可以了 就可以了
庄13爷，您对一个连 矩阵分解 和 矩阵分块 都没搞清的冤大头也太仁慈了吧？
不少了，你想阿，如果是在柏林踩三轮车给人家帮忙搬家，一个小时十快钱最多了，写个程序怎么也要一天8个小时吧。所以换算成欧元差不多就是这个价。
运气好还可以卖人家个沙发，赚个50 大欧， 不过这种事情不是每次都能碰上的。 再说还有被人家人肉的风险，我可不认识版主，弹压不下去。
本帖最后由 农大 于 2010-11-11 21:44 编辑
2000*2000 用 “解析" 是不合适的
而且所有所谓的数学软件matlab之类，对待这样的矩阵也不会用所谓的”解析“解
lz 要是会c或者c++ 随便找个solver 解一下就行。。。
本帖最后由 orionsnow 于 2011-1-12 19:59 编辑
上来汇报下。 花了一个多月， 两个超大矩阵在特殊情况下的解析解都被我给推导出来了。
然后我老板在网上搜了几天，居然找到了第一个矩阵解法的推导，那人最后把特例朝他的一般推导里头一代入，就得到了我推出来的解析解。 根据这个事实和数值算法的验证，就是前边无数同学提到的， 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 的被引用情况，没有人提到更深入的解法了。
估计现在计算机技术还没有那么发达，只能求数值解。 如果计算机能做推导，那估计数学家可以改行了。
再顶起来一下， 下一篇文章还是关于这类公式推导的。现在老板都快把这个当产业了，把以前的老算法都拿出来， 设某个参数为特殊值，然后加上点解析化简，就变成新文章了。