Universidad avanzadaComputationPlanificado
Numerical analysis
Numerical linear algebra (LU and QR factorisations, conditioning, least squares, eigenvalue algorithms), numerical solution of ODEs (Euler and Runge–Kutta methods, stiffness), and convergence and stability analysis.
math.comp.numerical-analysis
La lección de este tema todavía no se ha escrito. El esquema siguiente forma parte del mapa del plan de estudios y no ha sido revisado.
Se basa en
Requisitos previos obligatorios, hasta el fondo: 41
Aprenderás a
- Solve linear systems by LU factorisation with partial pivoting and analyse their conditioning.
- Solve least-squares problems via QR factorisation and compute eigenvalues with the QR algorithm.
- Implement Euler and Runge–Kutta methods and analyse their order and stability, including for stiff problems.
Definición
- The condition number κ(A) = ‖A‖ ‖A⁻¹‖ measures how sensitive the solution of Ax = b is to perturbations of the data.
- An algorithm is backward stable if it computes the exact answer to a slightly perturbed problem.
- A one-step method for ODEs has order p if its global error is O(hᵖ) for step size h.
Teorema
- The relative error in the solution of Ax = b is bounded, to first order, by κ(A) times the relative error in the data.
- Gaussian elimination with partial pivoting is backward stable in practice, although its worst-case growth factor is 2ⁿ⁻¹.
- The classical Runge–Kutta method has order 4.
- Explicit methods applied to stiff problems need very small steps, whereas some implicit methods are A-stable.
Errores frecuentes
- Computing A⁻¹ in order to solve Ax = b.
- Assuming a smaller step size always reduces the total error; rounding error eventually grows.
Descubrimientos opcionales
Itinerarios de aprendizaje
Lecturas complementarias
- Lloyd N. Trefethen and David Bau III, Numerical Linear Algebra (1997)
- Endre Süli and David F. Mayers, An Introduction to Numerical Analysis (2003)