作者Jerome1129 (KK)
看板MATLAB
标题[讨论] 有关联立方程式运算
时间Tue Nov 18 23:00:12 2014
大大们好，小弟有27个方程式欲求解
想用矩阵的方式来解开，
矩阵如下:
a=[1 2 1 1 2 1 1 2 1 zeros(1,18);
1 -2 1 1 2 1 1 2 1 zeros(1,18);
1 2 1 1 2 1 1 2 1 4 -4 -4 zeros(1,15);
1 2 1 1 -2 1 1 2 1 zeros(1,18);
1 -2 1 1 -2 1 1 2 1 zeros(1,18);
1 2 1 1 -2 1 1 2 1 4 -4 -4 zeros(1,15);
1 2 1 1 2 1 1 2 1 0 0 0 4 -4 -4 zeros(1,12);
1 -2 1 1 2 1 1 2 1 0 0 0 4 -4 -4 zeros(1,12);
1 2 1 1 2 1 1 2 1 4 -4 -4 4 -4 -4 zeros(1,12);
1 2 1 1 2 1 1 -2 1 zeros(1,18);
1 -2 1 1 2 1 1 -2 1 zeros(1,18);
1 2 1 1 2 1 1 -2 1 4 -4 -4 zeros(1,15);
1 2 1 1 -2 1 1 -2 1 zeros(1,18);
1 -2 1 1 -2 1 1 -2 1 zeros(1,18);
1 2 1 1 -2 1 1 -2 1 4 -4 -4 zeros(1,15);
1 2 1 1 2 1 1 -2 1 0 0 0 4 -4 -4 zeros(1,12);
1 -2 1 1 2 1 1 -2 1 0 0 0 4 -4 -4 zeros(1,12);
1 2 1 1 2 1 1 -2 1 4 -4 -4 4 -4 -4 zeros(1,12);
1 2 1 1 2 1 1 2 1 0 0 0 0 0 0 4 -4 -4 zeros(1,9);
1 -2 1 1 2 1 1 2 1 0 0 0 0 0 0 4 -4 -4 zeros(1,9);
1 2 1 1 2 1 1 2 1 4 -4 -4 0 0 0 4 -4 -4 zeros(1,9);
1 2 1 1 -2 1 1 2 1 0 0 0 0 0 0 4 -4 -4 zeros(1,9);
1 -2 1 1 -2 1 1 2 1 0 0 0 0 0 0 4 -4 -4 zeros(1,9);
1 2 1 1 -2 1 1 2 1 4 -4 -4 0 0 0 4 -4 -4 zeros(1,9);
1 2 1 1 2 1 1 2 1 0 0 0 4 -4 -4 4 -4 -4 zeros(1,9);
1 - 1 1 2 1 1 2 1 0 0 0 4 -4 -4 4 -4 -4 zeros(1,9);
1 2 1 1 2 1 1 2 1 4 -4 -4 4 -4 -4 4 -4 -4 zeros(1,9)]
c=[171 51 75 155 35 59 227 107 131 147 27 51 131 11 35 203 83 107 155 35 59 139 19 43 211 91 115]'
想要用inv(a)*c 的方法来求到答案，
但是在inv(a)的时候得到所有的值都是Inf，
大概知道是因为此矩阵det(a)=0，为奇异矩阵的关系，
想请问一下真的没办法解了吗?
感谢大大们拨空看这个问题，谢谢你们。
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.138.176.192
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/MATLAB/M.1416322815.A.E2D.html
1F：推 YoursEver: 你的a矩阵总共有九个column的0; underdetermined. 11/19 00:00
2F：→ Jerome1129: Y大你好，我有试过把後面九个为0的删掉变成18x18 11/19 00:18
3F：→ Jerome1129: 的矩阵，但结果还是一样，所以这个方程式无解吗? 11/19 00:19
4F：→ Jerome1129: 非常感谢你的回覆喔。 11/19 00:19
5F：→ rex0707: 试试自己写叠代法求解 11/19 23:22
6F：推 jack10313: 你的矩阵元素a(26,2) 没有打出来 不过不影响计算 11/20 23:04
7F：推 jack10313: 矩阵没满秩 可以使用虚反矩阵做计算(苏都因佛斯) 11/20 23:11
8F：→ jack10313: 简单来说 pinv(a)*c 要知道原理请查阅linear algebra 11/20 23:13
9F：推 tn00364361: 关键字：Moore-Penrose pseudo inverse跟SVD 11/21 15:45