24小时热门版块排行榜    

CyRhmU.jpeg
查看: 243  |  回复: 1
当前主题已经存档。

乔丹1998

新虫 (初入文坛)

[交流] 【求助】有点特殊的线性超定方程组求解

小弟土木的,做谱方法在土木工程中应用的。
在算法实现中遇到下面的问题,请各位大侠帮忙了。工程类的数学基础不好,大侠们见笑了。
解一个线性超定方程组
假设AX=B
A是n+2行n列,B是n+2行1列
但由于实际问题所致
A矩阵前2行n列的元素都是个位数,而后n行n列的元素比较大达到几千的数量级。
而对应B向量的前2行也是个位数,但B向量的后n行都为0.举例形如如下表达式:
       A                            B
  2        3                       1
  4        1                       2
3000  1500                    0
1800  2500                    0
用matlab做最小二乘解,所得解为[1.73e-6,-1.24e-6], 在前两行上出现相当大的错误。
小弟想问,是不是这种方程组没有最合适的解?也就是说我的前面的理论推导过程错了。就不应该推导出这样的方程组。如果是这样,我得回去检查理论推导过程。如果能解,请大侠们指点。或是建议我看下什么书籍。
小弟在此谢过了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shanby3531

新虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
和matlab系统的浮点运算有很大关系,一般来说,我们的四则运算,尤其是加减法是有结合律的,可是在matlab里则不行,所以这也是matlab自己承认的一个小bug
2楼2009-05-26 11:16:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 乔丹1998 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见