for equation set ax
=p若矩陣a滿秩,則:
1. if a is a square matrix and has a−
1 , (正定): x=
a−1p
. 2. if a is not a square matrix,使用廣義逆矩陣.
a超定,x=
(ata
)−1a
tp;
a欠定(大多數情況下),x=
at(a
at)−
1p.
a欠定:在影象重建中,意味著 投影線個數 < 影象畫素數。在高解析度下是很容易出現的情況,這也是迭代重建需要提供其他約束的原因。若矩陣a不滿秩,不是方陣,方程組不相容,則:
不滿秩:各條投影線可能相關,比如:0度和180度的投影線為同一條,當然相關了;使用奇異值分解svd尋找廣**。 am不是方陣:投影線個數 != 影象畫素數,這也是很容易理解的;
不相容:係數矩陣的秩 < 擴充套件矩陣的秩,這在投影資料含雜訊的情況下幾乎是必然的。
×n=u
m×mς
vtn×
n vt
v=in
×n u
tu=i
m×m
σm×n
=[di
ag00
0]
廣義逆矩陣為: a+
=vς+
ut其中, σ+
m×n=
[dr0
00]
dr=d
iag
x=a+
p=vς
+utp
線性方程組的數值解法
1.向量與矩陣的範數 norms of vectors and matrices 為了研究線性方程組數值解法的誤差估計和迭代法的收斂 性,有必要引進向量範數和矩陣範數的概念。歐式範數 設?1,2,1,2,稱 累加1 n 的?為向量?與?的內積,稱非負實數 2 根號 累加1 n的 2的1 2次方為向量...
方程 方程組的種類
線性方程也稱一次方程式。指未知數都是一次的方程。其一般的形式是ax by cz d 0。線性方程的本質是等式兩邊乘以任何相同的非零數,方程的本質都不受影響。線性方程 非線性方程,就是因變數與自變數之間的關係不是線性的關係,這類方程很多,例如平方關係 對數關係 指數關係 三角函式關係等等。求解此類方程...
c 解方程組 秩與方程組
今天要講的是兩個結論,通過對這兩個結論的理解和認識可以將很多東西串起來,既算是乙個深化認識,也算是乙個總結。對於方程組ax b 1 如果a是行滿秩的矩陣,那麼方程組要麼有唯一解,要麼有無窮多解。如果a是行滿秩的矩陣,因為矩陣的列秩等於矩陣的行秩,所以矩陣的列秩等於矩陣的行數,所以矩陣的列向量的線性組...