En análisis inverso, siempre calcular el número de condición. Si es alto, se necesita regularización. 4. Ejercicio 3: Regularización de Tikhonov Enunciado: Un experimento entrega el sistema mal condicionado: [ \beginpmatrix 1 & 1 \ 1 & 1.01 \endpmatrix \mathbfx = \beginpmatrix 2 \ 2.01 \endpmatrix ] Pero las mediciones tienen ruido. La solución por mínimos cuadrados da ( \mathbfx = (1, 1)^T ) exacto. Si añadimos ruido ( \mathbfb = (2, 2.02)^T ), la solución se vuelve ( \mathbfx = (0, 2)^T ). Aplique regularización de Tikhonov con ( \lambda = 0.1 ) para estabilizar.
El estimador de mínimos cuadrados minimiza ( | \mathbfF - \mathbfA k |^2 ). La solución normal: [ k = (\mathbfA^T \mathbfA)^-1 \mathbfA^T \mathbfF ] [ \mathbfA^T \mathbfA = 1^2+2^2+3^2+4^2 = 30 ] [ \mathbfA^T \mathbfF = 1\cdot2.1 + 2\cdot4.0 + 3\cdot6.2 + 4\cdot7.9 = 2.1+8.0+18.6+31.6 = 60.3 ] [ k = \frac60.330 = 2.01 \ \textN/m ] analisis inverso ejercicios resueltos
Hemos resuelto un problema inverso lineal bien condicionado. La matriz ( \mathbfA^T \mathbfA ) es invertible. Si los datos tuvieran ruido, el resultado cambiaría ligeramente. Este es el caso más simple de inversión. 3. Ejercicio 2: SVD y análisis del condicionamiento Enunciado: Considere un sistema ( \mathbfA \mathbfx = \mathbfb ) donde [ \mathbfA = \beginpmatrix 1 & 1 \ 1 & 1.0001 \endpmatrix, \quad \mathbfb = \beginpmatrix 2 \ 2.0001 \endpmatrix ] Resuelva el sistema y analice su sensibilidad al ruido. Luego, añada un pequeño error ( \delta \mathbfb = (0, 0.0001)^T ) y observe el cambio en la solución. En análisis inverso, siempre calcular el número de
