Abstract
Unstable linear boundary value problems can be solved by the method of Invariant Imbedding in a stable manner. Instead of integration of the system equations this method requires the integration of a matrix Riccati equation, which depends on the boundary values of the problem. The dimension of the Riccati equation is determined by a suitable decoupling of the system equations. Invariant Imbedding now fails, if this decoupling does not correspond with the boundary condition. In addition, the Riccati equation has to be solved once more for each new boundary condition. An extension algorithm is defined, which maps the boundary value problem into a problem of double dimension. This “extended” boundary value is solved by a modified Invariant Imbedding. The resulting “Extended (Dual) Invariant Imbedding” is always applicable and does not depend on the boundary conditions. The corresponding “extended” Riccati equation has to be integrated only once “offline”. If the boundary condition is changed, only systems of linear equations have to be solved “online”.
Zusammenfassung
Das Verfahren des Invarianten Einbettens kann instabile lineare Randwertaufgaben auf stabile Art lösen. Die dabei zu integrierende Matrixriccatigleichung hängt von den Randwerten des Problems ab. Ihre Dimension wird durch ein geeignetes Entkoppeln der Systemgleichungen bestimmt. Das Verfahren muß versagen, wenn dieselbe Entkopplung nicht auch mit den Randwerten durchführbar ist. Außerdem muß für neue Randwerte die Riccatigleichung nochmals gelöst werden. Es wird ein Erweiterungsalgorithmus definiert, der die Randwertaufgabe eindeutig auf einer doppelten Dimension abbildet, die dann durch modifiziertes Invariantes Einbetten gelöst wird. Dieses „Erweiterte (Duale) Invariante Einbetten” ist immer durchführbar. Die zugehörige Riccatigleichung hängt nicht mehr von den Randwerten ab, sie muß nur mehr einmal „offline” gelöst werden. Für neue Randwerte sind nur mehr lineare Gleichungssysteme „online” zu lösen.
Similar content being viewed by others
References
Allen, R. C., Wing, G. M.: An invariant imbedding algorithm for the solution of inhomogenous linear two-point boundary value problems. I. Comp. Phys.14, 40–58 (1974).
Bellman, R., Kagiwada, H. Kalaba, R.: Numerical results for the auxiliary equation of radiative transfer. J. Quant. Spectrosc. and Radiative Transfer6, 291–310 (1966).
Bellman, R., Wing, G. M.: An introduction to invariant imbedding. New York: J. Wiley 1975.
Breitenecker, F.: Extended invariant imbedding for the solution of unstable linear boundary value problems. Thesis. Technical University Vienna, 1979 (in German).
Busbridge, I.: The mathematics of radiative transfer. Cambridge: Cambridge University Press 1960.
Hahn, W.: Stability of motion. Berlin-Heidelberg-New York: Springer 1967.
Kalman, R. E.: Contributions to the theory of optimal control. Bol. Soc. Mat. Mexicana5, 102–119 (1960).
Knobloch, H. W., Kappel, F.: Gewöhnliche Differentialgleichungen. Stuttgart: B. G. Teubner 1974 (in German).
Meyer, G. H.: Initial value methods for boundary value problems. New York: Academic Press 1973.
Meyer, G. H.: On a general theory of characteristics and the method of invariant imbedding. SIAM J. Appl. Math.16, 488–509 (1968).
Reid, W. T.: Riccati differential equations. New York: Academic Press 1972.
Sagirow, P.: Satellitendynamik. (B. I. 719.) Mannheim-Wien-Zürich: Bibliographisches Institut (in German).
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Breitenecker, F. On the solution of linear boundary value problems by extended invariant imbedding. Computing 28, 333–343 (1982). https://doi.org/10.1007/BF02279816
Received:
Issue Date:
DOI: https://doi.org/10.1007/BF02279816