登入首頁
收藏 0
返回課程
Week5 Unit5 QR Decomposition
  • [自學課程]數值分析Numerical Analysis
  • 第一週:數值分析課程簡介、數值誤差的來源
  • Week1 Unit0 電腦的暴力
  • Week1 Unit1 授課方式介紹
  • Week1 Unit2 數值誤差的來源
  • Week1 Unit3 Python程式簡介
  • Week1 Unit4 數值誤差的來源
  • 第一週測驗
  • 第二週:基礎Python 程式入門
  • Week2 Unit0 Python程式安裝教學
  • Week2 Unit1 Python初步入門
  • Week2 Unit2 Python 常用的型別
  • Week2 Unit3 Python 流程控制
  • Week2 Unit4 Python 模組引用與常用函數
  • Week2 Unit5 Python 函數的定義方式
  • Week2 Unit6 Python 檔案的操作
  • Week2 Unit7 Python 繪圖
  • 基礎Python 程式入門
  • 第三週:非線性方程式的根、高維度非線性方程式的根
  • Week3 Unit00 Newton Method
  • Week3 Unit01 Genetic method
  • Week3 Unit02 Bisection method
  • Week3 Unit03 Simple iteration and python
  • Week3 Unit04 brouwer's fixed point theorem
  • Week3 Unit05 constraction mapping
  • Week3 Unit06 Newton's method
  • Week3 Unit07 Recall Newton's method
  • Week3 Unit08 newton's method 2nd convergence
  • Week3 Unit09 big O
  • Week3 Unit10 Mooc's 收斂速度定義
  • Week3 Unit11 Secant Method
  • 非線性方程式的根、高維度非線性方程式的根
  • 第四週:求解線性方程組
  • Week4 Unit0 當週課程介紹
  • Week4 Unit1 矩陣與線性方程組
  • Week4 Unit2 row operation
  • Week4 Unit3 Gauss elimination
  • Week4 Unit4 LU decomposition
  • Week4 Unit5 solving triangular system
  • 求解線性方程組
  • 第五週:正交化與特徵值
  • Week5 Unit0 當週課程介紹
  • Week5 Unit1 inverse power method
  • Week5 Unit2 Power method
  • Week5 Unit3 Orthogonalization of vectors
  • Week5 Unit4 Gram Schmidt
  • Week5 Unit5 QR Decomposition
  • 正交化與特徵值
  • 第六週:插值法
  • Week6 Unit0 當週課程介紹
  • Week6 Unit1 introduction to interpolation
  • Week6 Unit2 Lagrange interpolation
  • Week6 Unit3 Hermite interpolation
  • Week6 Unit4 Neville's algorithm
  • Week6 Unit5 basis of polynomial
  • Week6 Unit6 uniqueness of polynomial interpolation
  • Week6 Unit7 Newton interpolation
  • Week6 Unit8 B spline
  • 插值法
  • 第七週:數值微分
  • Week7 Unit0 當週課程介紹
  • Week7 Unit 1 數值微分:Forward, backward and central method
  • Week7 Unit 2 Richardson extrapolation
  • Week7 Unit 3 effect of rounding error
  • Week7 Unit 4 high order differential
  • 第八週:數值積分
  • Week8 Unit0 當週課程介紹
  • Week8 Unit1 中點法,梯形法和辛普森法
  • Week8 Unit2 數值積分誤差分析
  • Week8 Unit3 adaptive quadrature
  • Week8 Unit4 Gauss quadrature
  • 第九週:數值微分方程
  • Week9 Unit0 當週課程介紹
  • Week9 Unit1 Euler Method
  • Week9 Unit2 Heun's method
  • Week9 Unit3 Runge Kutta method
  • Week9 Unit4 shooting method
  • 索引
  • 重點
  • 討論 (4)
  • 筆記
長度: 17:25, 發表時間 : 2015-09-08 15:23
觀看次數 : 769
  • 17:25
    1. Week5 Unit5 QR Decomposition
Uploaded by NCCU Moocs on 2014-12-16.
附件
長度: 17:25, 發表時間 : 2015-09-08 15:23
觀看次數 : 769
Uploaded by NCCU Moocs on 2014-12-16.
原始資料來自 : http://moocs.nccu.edu.tw/media/148
附件
討論功能僅開放給課程成員,請先加入課程
最新的回應 ...more
2019-12-04
瑋辰 (asd415305@yahoo.com.tw) : 如果你想要得到長方形的由來,動手算一遍大致上就知道了,可以分成m>n and m<n
瑋辰 (asd415305@yahoo.com.tw)
並且他實際上在SQR有隱藏一個小技巧會跟你動手計算有不一樣的方法
for i in range(n):
for j in range(i):
R[j,i]=Q[:,j].T*A[:,i] ##不同之處,這是動手計算的方法
Q[:,i]-=R[j,i]*Q[:,j] ***********
if la.norm(Q[:,i])>tol:
R[i,i]=la.norm(Q[:,i])
Q[:,i]=Q[:,i]/R[i,i]
return Q,R


for i in range(n):

for j in range(i):
R[j,i]=Q[:,j].T*Q[:,i]這是老師的方法
Q[:,i]-=R[j,i]*Q[:,j] *********
if la.norm(Q[:,i])>tol:
R[i,i]=la.norm(Q[:,i])
Q[:,i]=Q[:,i]/R[i,i]
return Q,R
為何兩者結論相同,Q[:,i]-=R[j,i]*Q[:,j] =A[:,i] -\sum (R[1,i]*Q[:,1]+R[2,i]*Q[:,2]+......R[j,i]*Q[:,j])
下一次的迴圈 j+1
R[j+1,i]=Q[:,j+1].T*Q[:,i]=Q[:,j+1].T*[A[:,i] -\sum (R[1,i]*Q[:,1]+R[2,i]*Q[:,2]+......R[j,i]*Q[:,j])]
藉由前面得到的都是orthogonal set 所以會等於
Q[:,j+1].T*[A[:,i] over
2019-12-04 1 樓
瑋辰 (asd415305@yahoo.com.tw)
很精彩,尤其是設計長方形QR,以及當遇到l.d時如何避免。
2019-12-03 1 樓
瑋辰 (asd415305@yahoo.com.tw)
如果你想要得到長方形的由來,動手算一遍大致上就知道了,可以分成m>n and m<n
2019-12-04 2 樓
瑋辰 (asd415305@yahoo.com.tw)
在計算Q的同時計算出R
2019-12-03 1 樓
承翰 (chiha8888@gmail.com)
咳嗽聲...
2019-08-12 1 樓
筆記功能僅開放給課程成員,請先加入課程
Prev
Week5 Unit4 Gram Schmidt
Next
正交化與特徵值