臺灣博碩士論文加值系統

以有限元素法模擬的大型震動結構通常可以被離散化成二階動態系統的矩陣方程式。這些動態系統的解可以用公式表示成大型稀疏漢米爾頓矩陣特徵問題的解。傳統上採用的方法-舒而分解-所面臨的問題就是必須將整個矩陣存入電腦中，這種方法在電腦的記憶體不足時是無法運作的。而子空間逼近法通常會運用到矩陣的異動值和反矩陣來增加收斂速度，後者的運算量是相當貴的。並且實際上的困難是如何決定出這些異動值。
在這篇論文裡我們提議運用同倫分列式的方法，配合分割與合成的策略來計算異動值。並藉由主要奇子空間疊代法來逼近連續時間李卡提等式的解。論文中還列出了詳細的演算法和實驗的數據和經驗。最後我們針對實驗結果提出可能潛在的問題並討論之。

Large vibrating structures modeled by finite element methods are usually discretrized into second order dynamical systems of matrix equations. Solutions of such systems can be formulated as solving large sparse Hamiltonian eigenproblem. QR-type or Schur method suffer from fill-in's and are not practical for large problem. Subspace approaches usually require shift and invert to accelerate the convergence rate. It is, however, difficult in practice to determine the shift values. In this thesis we propose to use the divide-and-conquer homotopy-determinant algorithm for computing shift values, and dominant singular subspace updating method for approximating the solution of the underlined continuous-time Riccati equation. The detailed algorithms is presented, experiment results and observed difficulties is discussed.

1 Introduction 1
2 Feedback Stabilization and Optimal Control 3
2.1 State-Feedback Stabilization…………………………………………… 4
2.2 Continuous-time Linear Quadratic Regulator ………………………… 5
2.3 Existence and Uniqueness of the Stabilizing Solution of CARE……… 11
2.4 Invariant Subspace Methods for Riccati Equation …………………… 14
3 Numerical Methods for Large Space Problems 16
3.1 J-Lanczos Algorithm…………………………………………………… 16
3.2 Homotopy Determinant Method ……………………………………… 22
3.2.1 Hyman’s Method ……………………………………………………… 23
3.2.2 Following the Eigenvalue Paths ……………………………………… 23
3.3 Simple Subspace Iteration……………………………………………… 28
3.4 Recursive Calculation of Dominant Singular Subspaces……………… 29
4 Numerical Strategy for Large Sparse Problems 34
5 Numerical Experiments and Results 38
6 Conclusion 41
Reference 42

B.N. Datta, Numerical methods for linear control systems
design and analysis, preprint, 2000.
A.J. Laub, A Schur method for solving algebraic Riccati
equations, IEEE Trans AC, 1979.
W.R. Ferng, W-W. Lin, and C.S. Wang, The shift-inverted
$J$-Lanczos algorithm for the numerical sloutions of large algebraic Riccati equations, J. Comp. Math. Appl., 1997.
T.Y. Li and Z. Zeng, Homotopy-determinant algorithm for
solving nonsymmetric eigenvalue problems, Math. Comput., 1992.
B.T. Smith et al., Matrix eigensystem routines-EISPACK
guide, 2nd ed., Springer-Verlag, 1976.
G.H. Golub and C.F. Van Loan, Matrix Computation,
3rd ed., Johns Hopkins, London, 1996
Y. Chahlaoui, K. Gallivan, and P. Van Dooren, Recursive
calculation of dominant singular subspace}, preprint, 2000.