跳到主要內容

簡易檢索 / 詳目顯示

研究生: 陳揚敏
論文名稱: 縮基法初始值問題之數值研究
Numerical studies of reduced basis methos for initial value problems
指導教授: 林美佑
學位類別: 碩士
Master
系所名稱: 理學院 - 應用數學系
Department of Mathematical Sciences
論文出版年: 1990
畢業學年度: 78
論文頁數: 31
中文關鍵詞: 縮基法,投影法
外文關鍵詞: Reduced Basis Method, Projection
相關次數: 點閱:320下載:0
分享至:
查詢本校圖書館目錄 查詢臺灣博碩士論文知識加值系統 勘誤回報
  • 縮基法(RBM) 是對參數化的曲線求逼近解的一個方法,基本上乃使用投影法將解曲線投射到解空間的一子空間中,如此一來,可將原問題轉換成一較小的系統,並經由數值計算出小系統的解,來求得大系統的一逼近解。在本篇論文中主要的乃探討RBM在常微分方程組初始值問題上的應用,並發展一套含有誤差控制的演算法。

    本篇論文中所採用的ODE Solver 乃由Gordon 和Shampine 基於Adams PECE方法所發展的。在求解的過程中,對於計算解誤差的控制我們除了利用ODE Solver 的誤差估計,另外我們又發展對縮基解(reduced basis solution) 的後(Aposteriori)

    誤差估計,以確保數值計算解的準確性。我們所考慮使用的子空間有三種Taylor, Lagrange , Hermite 。同時為了要增加數值的特定性及簡化小系統的求解工作,我們先行將子空間的基底直交化。因此,除了誤差的控制外,我們也討論了roundoff error 對向量直交化及形成小系統時所造成的影響,並設立誤差標準以判別何時誤差過大到嚴重影響縮基解的準確度。

    本篇論文的目的是希望利用RBM發展出一套解常微分方程組初始值問題的求算法,以期計算解能在較短的時間內準確的被計算出來。


    The reduced basis method(RBM) is a scheme for approximating parametric solution curves. The basic technique of RBM is projection. By applying the method, we can find an approximate solution of the original system which satisfies a system of smaller size. In this paper, we mainly concern the applications of RBM for ODE initial value problems and develop an algorithm which contains a set of error controls.

    The ODE solver used in this paper is developed by Gordon and Shampine based on Adams PECE formulas. To assure the accuracy of the reduced basis approximation, we set up an appropriate automatic error control in calling GS solver and develop an a posteriori error estimate to keep the reduction error under control.

    The subspaces considered are Taylor, Lagrange and Hermite subspaces.In the meantime, in order to improve the numerical stability and simplify the computation of the reduced basis solution, we orthogonalize the generators of reduced subspaces. We also discuss the roundoff errors in the orthogonalization process and build up a criterion for identifying the case the accuracy of the reduced basis solution up a criterion for identifying the case the accuracy of the reduced basis solution is destroyed by the errors.

    The aim of this paper is to develop an algorithm to solve the ODE initial value problems efficiently.

    1. Introduction ......................... 1
    2. Reduced Basis Approximation . ......................... 2
    3. Ordinary Differential Equation Solver ......................... 3
    3.1 Operation Count of GS-solver ......................... 8
    4. Implementing RBM.......................... 9
    4.1 Selecting Subspaces ......................... 10
    4.2 Algorithm of Gram-Schmidt .................................... 14
    4.3 Roundoff Error......................... 16
    4.4 Order Control ......................... 22
    4.5 Operation Count of Overhead of RBM......................... 22
    5. Error Control ......................... 24
    5.1 A posteriori Error Estimate......................... 24
    5.2 A priori Error Estimate......................... 25
    6. Numerical Studies ......................... 26
    References......................... 30

    [1] N. N. Abdelmalek, Roundoff Error Analysis for Gram-Schmidt Method and Solution of Linear Least Squares Problems, BIT, 11(1971), pp.945-968.
    [2] B. O. Almroth, P. Stern and F. A. Brogan, Automatic Choice of Global Shape Functions in Structural Analysis, AIAA J., 16(1978), pp. 525-528.
    [3] E. Fehlberg, Klassische Runge-Kutta-Formeln vierter und niedrigerer Ordnung mit Schrittweiten-Kontrolle und ihre Anwendung auf Warmeleitungsprobleme, Computing, 6(1970), pp. 61-71.
    [4] J. P. Fink and W. C. Rheinboldt, On the Discretization Error of Parametrized Nonlinear Equations, SIAM J. Numer. Anal., 20(1989),pp. 792-746.
    [5] J. P. Fink and W. C. Rheinboldt, On the Error Behavior of the Reduced Basis Technique for Nonlinear Finite Element Approximations,Z. Angew. Math.Mech., 69(1989), pp. 21-28.
    [6] J . P. Fink and W. C. Rheinboldt, Local Error Estimates for Parametrized Nonlinear Equations, SIAM J. Numer. Anal., 22(1985),pp. 729-795.
    [7] C. W. Gear, Numerical Initial Value Problems in Ordinary Differential Equations, 1971, Prentice-Hall inc., Englewood Cliffs, New Jersey.
    [8] G. H. Golub and C. F. Van Loan, Matrix Computations, 1989.
    [9] M. K. Gordon and L. F. Shampine, Computer Solution of Ordinary Differential Equations, The Initial Value Problem, ~974, W. H. Freeman and Company, San Francisco.
    [10] J. D. Lambert, Computational Methods in Ordinary Differential Equations, 1989 J. W. Arrowsmith Ltd. Bristol.
    [11] M. Lin Lee, Estimation of the Error in the Reduced Basis Method Solution of Differential Algebraic Equation System, SIAM J. Numer. Anal., to appear.
    [12] M. Lin Lee, The Reduced Basis Method for Differential Algebraic Equation System, Technical Report ICMA-85-85, Inst. for Compo Math. And Appl., University of Pittsburgh, Pittsburgh, PA, July, 1985.
    [13] N. A. Nagy, Model Representation of Geometrically Nonlinear Behavior by the Finite Element Method, Computers and Structures,10(1977), pp. 689-688.
    [14] A. K. Noor and J. M. Peters, Reduced Basis Technique for Nonlinear Analysis of Structures, AIAA J., 18(1980), pp. 455-462.
    [15] A. K. Noor, C. M. Andersen and J. M. Peters, Reduced Basis Technique for Collapse of Shells, AIAA J., 19(1981), pp. 999-997.
    [16] T. A. Porsching, Estimation of the Error in the Reduced Basis Method Solution of Nonlinear Equations, Math. Comp., 45(1985), pp. 487-496.
    [17] T. A. Porsching and M. Lin Lee, The Reduced Basis Method For Initial Value Problems, SIAM J. Numer. Anal., 24(1987), pp. 1277-1287.
    [18] J. R. Rice, Experiments on Gram-Schmidt Orthogonalization, Math. Compo 20(1966), pp. 925-928.
    [19] J. Stoer and R. Bulirsch, Introduction to Numerical Analysis, 1980,Springer, New York, Heidelberg, Berlin.
    [20] G. W. Steward, Introduction to Matrix Computation, 1979, Academic Press, New York and London.
    [21] G. W. Steward, Perturbation Bounds for the QR Factorization of a Matrix, SIAM J, Numer. Anal., 14(1977), pp. 509-518.
    [22] J. S. Vandergraft, Introduction to Numerical Computation, 1989,Automated Sciences Group, Inc., Silver Spring, Maryland.
    [23] R. E. Williamson, Introduction to Differential Equations, ODE, PDE and Series.

    無法下載圖示 (限達賢圖書館四樓資訊教室A單機使用)
    QR CODE
    :::