13.9: Iterando
- Page ID
- 131052
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Ahora podemos usar las Ecuaciones 13.8.35a, b y obtener una mejor estimación de las relaciones triangulares. Los datos numéricos son
\(b_1 = 2/3 , \quad b_3 = 1/3 , \quad r_2 = 3.481 \ 33 , \)
\(τ_1 = t_3 − t_2 = 10\)días solares medios y días solares\(τ_3 = t_2 − t_1 = 5\) medios, pero recordemos que estamos expresando intervalos de tiempo en unidades de\(1/k\), que es días solares\(58.132 \ 440 \ 87\) promedio, y por lo tanto
\(τ_1 = 0.172 \ 021 \quad \text{and} \quad τ_3 = 0.086 \ 010.\)
Las ecuaciones 13.8.35 luego dan como resultado
\(a_1 = 0.666 \ 764, a_3 = 0.333 \ 411\).
Ahora podemos volver a la Ecuación 13.7.4 y comenzar de nuevo con nuestros nuevos valores para las relaciones triangulares — und so weiter − hasta que obtengamos nuevos valores para\(∆_1 , \ ∆_2 , \ ∆_3\) y\(r_2\). A continuación muestro en las dos primeras columnas las primeras estimaciones brutas (ya dadas anteriormente), en las 16 segundas dos columnas los resultados de la primera iteración, y, en las dos últimas columnas, los valores dados en las\(\text{IAU}\) efemérides publicadas.
\ begin {array} {l c c c}
&\ text {Primeras estimaciones crudas} &\ text {Primera iteración} &\ text {MPC}
\\\\
&\ quad r &\ quad r & □\ quad r\\\
1 & 2.72571\ 3.48532 & 2.65825\ 3.41952 & 2.644\ 3.406\
2 & 2.68160\ 3.48133 & 2.61558\ 3.41673 & 2.603\ 3.404\\
3 & 2.61073\ 3.47471 & 2.54579\ 3.41082 & 2.536\ 3.401\
\ final {array}
Vemos que hemos logrado una mejora sustancial, pero aún no estamos ahí. Ahora podemos calcular nuevos valores de\(a_1\) y\(a_3\) a partir de Ecuaciones 13.8.35a, b para obtener
\(a_1 = 0.666 \ 770 \quad a_3 = 0.333 \ 416\).
Podríamos (si así lo deseáramos) ahora volver a las Ecuaciones 13.7.4,5,6, e iterar de nuevo. Sin embargo, esto dará como resultado solo pequeños cambios a\(a_1\),\(a_3\),\(∆\) y\(r\), y tenemos que tener en cuenta que las Ecuaciones 13.8.35a, b son solo aproximaciones (a la orden\(τ^3\)). Por lo tanto, incluso si convergen iteraciones sucesivas, todavía no darán respuestas precisas y correctas para\(∆\) y\(r\).
Para anticiparnos, eventualmente llegaremos a algunas Ecuaciones exactas (Ecuaciones 13.12.25 y 13.12.26) que nos permitirán resolver el problema. Pero estas Ecuaciones no serán fáciles de resolver. Tienen que resolverse por iteración usando una primera suposición razonablemente buena. Nuestro objetivo actual es obtener una primera conjetura razonablemente buena para\(a_1\),\(∆\) y\(a_3\)\(r\), con el fin de prepararnos para la solución de las Ecuaciones exactas 13.12.25 y 13.12.26. Nuestros valores actuales de\(a_1\) y\(a_3\), si bien no son exactos, nos permitirán resolver exactamente las Ecuaciones 13.12.25 y 13.12.26, así que ahora, en lugar de volver de nuevo a las Ecuaciones 13.7.4,5,6, proceder directamente a las Secciones 13.11, 13.12 y 13.13.
Sin embargo, en la siguiente sección, proporcionamos (en las Ecuaciones 13.10.9 y 13.10.10), después de un esfuerzo considerable, expansiones de orden superior para\(a_1\) y\(a_3\). Estos pueden ser útiles, pero por razones explicadas en el párrafo anterior, puede ser más fácil saltarse por completo la Sección 13.10.