[問題] sparse矩陣如何使用pinv求解?!已回收
k=sparse(3,3);
k(1,1)=2;k(2,1)=1;k(2,2)=1;k(3,2)=1;k(3,3)=0;
a=[20 12 41]'
求解
pinv(k)*a
因為det(k)=0,也就是奇異矩陣,所以我用pinv來求解,
但好像是sparse不支援pinv的關係,會出現錯誤....
因為我的矩陣是k=144*144,所以上面用3*3舉例,
不用zeros寫是因為要減少迭代的時間,
故想請問有辦法讓這個值跑出來,或者修改哪邊讓它執行成功?!
感謝大大們的解惑~
感謝您!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 163.13.138.211
推
11/13 20:59, , 1F
11/13 20:59, 1F