2019年2月19日 星期二

Julia 擷取矩陣的一部分

Julia 擷取矩陣的一部分

擷取矩陣向量積的一部分,當您將矩陣作為二維數組時,您可能希望提取其中的一部分。
例如,當您想要提取上面對角化的矩陣的特徵向量時。
如果檢查上面獲得的特徵值特徵向量是否是正確的特徵值特徵向量。

using LinearAlgebra
H = [1 3
3 2]
e,u = eigen(H)
println(e)
println(u)

u1 = u[:,1]
u2 = u[:,2]
println(u1)
println(u2)

r1 = H*u1 - e[1].*u1
println(r1)
r2 = H*u2 - e[2].*u2
println(r2)



輸出結果
[-1.54138, 4.54138]
[-0.76302 0.646375; 0.646375 0.76302]
[-0.76302, 0.646375]
[0.646375, 0.76302]
[2.22045e-16, 2.22045e-16]
[4.44089e-16, 8.88178e-16]

沒有留言:

張貼留言

Galois LFSR

  在密碼學與數位訊號處理中, LFSR(線性回饋移位暫存器,Linear Feedback Shift Register) 是用來產生偽隨機序列(也就是串流加密中所需的金鑰流)最核心的硬體架構。 LFSR 主要分為兩種實現架構: Fibonacci(斐波那契) 與 Galo...