24小时热门版块排行榜    

查看: 1066  |  回复: 7
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 MyDianaL 的 8 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

MyDianaL

金虫 (正式写手)

[交流] matlab奇异矩阵求解

Ax=y,x=A/y已试过,效果不好,有没有其他方法?
回复此楼

» 猜你喜欢

一只伟大而神奇的猫!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yourmav

铁虫 (初入文坛)


小木虫(金币+0.2):抢了个小板凳,给个红包
MyDianaL(金币+6):谢谢讨论!!! 2010-03-11 10:26
如果就是解一个Y=AX的方程呢
A是32*59的,该怎么判断呢?
引用回帖:
Originally posted by javeey at 2010-03-08 23:11:15:
是近似奇异矩阵还是完全奇异矩阵?
近似奇异矩阵让matlab很为难,因为逆矩阵数量级太大,也使解的扰动太大,很可能与真实解相差太远。如果A中元素行之间的数量级很大,可考虑数据标准化,否则按完全奇异对待。
...

my girl
3楼2010-03-08 23:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

javeey

荣誉版主 (职业作家)

力拔山兮气盖世

优秀版主优秀版主


小木虫(金币+0.5):恭喜抢沙发,给个红包
MyDianaL(金币+6):谢谢讨论! 2010-03-11 10:26
是近似奇异矩阵还是完全奇异矩阵?
近似奇异矩阵让matlab很为难,因为逆矩阵数量级太大,也使解的扰动太大,很可能与真实解相差太远。如果A中元素行之间的数量级很大,可考虑数据标准化,否则按完全奇异对待。
如果是完全奇异矩阵,其要么无解,要么有无穷多个解。就是看秩(A|y)是否等于秩(A),相等则有无穷多个解,或者说确定x的方程不够,要增加样本;不相等的话,就用最小二乘法确定,这种情况是最好的,条件多可以避免实验过程中的误差。做回归的话,一般方程数大于未知数个数的5倍才好。
按你所说,据我所知,若方程为实验所得,最好是增加方程个数,即再多弄些实验数据来,否则matlab也很难处理。
早起的鸟儿有虫吃,早起的虫儿被鸟吃
2楼2010-03-08 23:11:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

javeey

荣誉版主 (职业作家)

力拔山兮气盖世

优秀版主优秀版主

引用回帖:
Originally posted by yourmav at 2010-03-08 23:28:01:
如果就是解一个Y=AX的方程呢
A是32*59的,该怎么判断呢?




很明显,方程个数是32,未知数个数是59,不定方程,有无穷多个解,但能求出解得一般形式(稍难,对应用来说无多大意义)或者某一个解。
早起的鸟儿有虫吃,早起的虫儿被鸟吃
4楼2010-03-08 23:37:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yourmav

铁虫 (初入文坛)

引用回帖:
Originally posted by javeey at 2010-03-08 23:37:01:

很明显,方程个数是32,未知数个数是59,不定方程,有无穷多个解,但能求出解得一般形式(稍难,对应用来说无多大意义)或者某一个解。

那最小二乘法解,还有矩阵分解的方法有使用的限制吗?
my girl
5楼2010-03-08 23:41:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见