[讨论] 有关联立方程式运算

楼主: Jerome1129 (KK)   2014-11-18 23:00:12
大大们好,小弟有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,为奇异矩阵的关系,
想请问一下真的没办法解了吗?
感谢大大们拨空看这个问题,谢谢你们。
作者: YoursEver (猪是妳)   2014-11-19 00:00:00
你的a矩阵总共有九个column的0; underdetermined.
楼主: Jerome1129 (KK)   2014-11-19 00:18:00
Y大你好,我有试过把后面九个为0的删掉变成18x18的矩阵,但结果还是一样,所以这个方程式无解吗?非常感谢你的回复喔。
作者: rex0707 (我就是我)   2014-11-19 23:22:00
试试自己写叠代法求解
作者: jack10313 (彩色笔)   2014-11-20 23:04:00
你的矩阵元素a(26,2) 没有打出来 不过不影响计算矩阵没满秩 可以使用虚反矩阵做计算(苏都因佛斯)简单来说 pinv(a)*c 要知道原理请查阅linear algebra
作者: tn00364361 (小氵斿)   2014-11-21 15:45:00
关键字:Moore-Penrose pseudo inverse跟SVD

Links booklink

Contact Us: admin [ a t ] ucptt.com