第10章 疑惑
第二天，清晨，雨过天晴，阳光明媚，一缕柔和的阳光从窗外透了进来，照到了田中爱子的脸上，几只调皮的燕儿，似乎无比欢乐，在外面喳喳地叫个不停，将尚在梦中的她给惊醒了过来。睡眼惺忪的她睁眼一看，原来天早已经放晴了，连忙坐起身子，理了理有些蓬乱的头发，接着就急忙翻身下床，一边走出卧室，一边心里暗自懊恼：
完了完了，爱子啊爱子，你怎么这么不争气呢，第一天就睡懒床，你还要不要脸了？
看到田中爱子走了出来，孙月娥笑盈盈地说到：“闺女，怎么这就起来了，多睡会儿没事的。”
田中爱子听了，更加不自如起来，脸上微红，含羞地说到：“妈，不早了，天早就亮了，再睡我自己都不好意思见您了。”
“看，又跟妈见外了，自己家有什么不好意思的，闺女，既然起来了，那不如先去洗漱间洗洗，洗漱用品你爸一早起来就给你买好了。”
“好的，妈！”
“那妈领你去。”
看到洗漱台上摆好的洗漱用品，田中爱子心里不自觉的一甜，心想：真是上辈子修来的好福气，让我因祸得福，遇到了待我这么真心的爸妈，一定不能辜负了他们。
……
片刻功夫，一名翩若惊鸿，婉若游龙，荣耀秋菊，华茂春松的女子走了出来。孙月娥看到后震惊了，心里不自觉地暗自称赞个不停。
“妈，我爸呢，我怎么没看到他？”田中爱子轻问到。
孙月娥闻言，才回过神来，微笑地回到：“闺女，你爸他上班去了，他工作忙，任务重，不得不去，要妈今天在家陪你。”
“哦，也对，都9点了！都怪我起得太晚了，妈，我是不是给您添麻烦了，你也有工作的，哎……”田中爱子一脸懊恼的道。
“闺女，说到哪里去了，不要跟妈这么客气了，显得多生分……”说着说着，孙月娥竟伤心起来。
田中爱子看到孙月娥两眼微红，连忙往她怀里一扑，撒娇道：“妈，女儿知道了，女儿错啦，您别生气了。”
孙月娥被田中爱子这么一撒娇，破涕为笑的说：“嗨，我说乖女儿，你还真会看脸色，妈现在对你真是又疼又爱了，以后千万不要再和妈说生分的话了，好吗？”
“妈，女儿保证，再也不敢了。”说完田中爱子往孙月娥脸上亲了一口。
哈哈……一阵欢快的笑声游荡在客厅里……
此时的孙月娥，心里有说不出的欣慰与开心，她是好久没有这么发自内心的笑了，天真烂漫的田中爱子，让她从连失两子的悲痛中走了出来。
……
两人吃过早餐，孙月娥说到：
“闺女，你爸昨晚说要给你的财富，他后来找到了，给你留在书桌上的，我们母女是先出去散散心，还是……”
田中爱子闻言，连忙说到：“妈，我要看看爸给我留的财富，女儿好奇的很呢。”
听田中爱子这么一说，孙月娥哈哈地笑到：“哈哈……好吧，闺女，说真的，我也好奇，那我们先去书房看看去。”
说完两人手拉手往书房走去。
孙月娥一边开电脑，一边对田中爱子笑着说：
“闺女，你爸这人，哈哈，有时候妈真不知道拿什么好话说他了。”
“妈，爸这叫笃实好学，态度这么认真，办事这么稳重，天底下能有几个他这样的？妈，我们应该感到自豪，是不是？”田中爱子满脸崇拜地说到。
孙月娥见说得这么好，眉开眼笑的说到：
“闺女你这么说，妈也觉得有那么点意思了。来，坐下看看，你爸说你看到后一定会喜欢的。”
田中爱子问闻言，也不客气，坐下去好奇的看了起来。
这不看不知道，一看吓得她一跳。这真是一笔丰厚的知识财富，全是张起明在他所在航天动力学方面的知识积累和笔记。
然而更让田中爱子吃惊的是，原本以为，以她目前现有的知识量，不好好学个两三年，是看不懂多少的，可没想到啊，没想到啊，她大脑对这些知识似乎了若指掌，竟一看就懂了，看到有的地方甚至是茅塞顿开起来。
这太不可思议了，这不科学啊，自己明明记得没学过导航方面的知识啊，可我怎么还是看得懂，就连那些公式，脑子也可以飞速地推导验算。
只见：
以航天器的位置和速度作为状态变量，其状态方程如下所示:
x?(t)=f(x，t)+ω(t)(1)
其中，x(t)是状态变量;f(x，t)是状态动力学的非线性函数;ω(t)是噪声项．
考虑地球非球形引力摄动j2项，动力学方程如下式:
χ?=ν_x
y?=ν_y
z?=ν_z
(ν_x)?=-μr_x/r^3[1-j_2(r_e/r)^2(7.5(r_z^2)/r^2-1.5)]
(ν_y)?=-μr_y/r^3[1-j_2(r_e/r)^2(7.5(r_z^2)/r^2-1.5)]
(ν_z)?=-μr_z/r^3[1-j_2(r_e/r)^2(7.5(r_z^2)/r^2-4.5)]（2）
其中，(x，y，z)是航天器的位置矢量;(ν_x，ν_y，ν_z)是航天器相对于地球质心的速度;ｒe为地球半径．方程(2)的解可直接用四阶runge-kutta法数值计算得到．1．2测量方程在j2000.0地心赤道惯性坐标系，将航天器测得时间与地球质心的时间延迟量作为测量方程，在j2000.0地心赤道惯性坐标系，将航天器测得时间与地球质心的时间延迟量作为测量方程，有
?t=t_e-t_sc=1/cn?r+ν(3)
则得到观测方程为：
y=hx+v(4)
h=1/c[h^2，o_nx3]
h?=[n_1，n_2，?，n_n]^t
n_i=[n_1^i，n_2^i，n_3^i]^t
……
……
随着脑子里飞速的验算着这些陌生的公式，田中爱子额头上的汗也冒了出来，匪夷所思的让她再也不敢看下去了。
我这是怎么回事，我怎么突然凭空会这些我原本不会的知识了？
不可思议、不可捉摸的事，让田中爱子渐渐地陷入了深深的沉思……
(本章完)
完了完了，爱子啊爱子，你怎么这么不争气呢，第一天就睡懒床，你还要不要脸了？
看到田中爱子走了出来，孙月娥笑盈盈地说到：“闺女，怎么这就起来了，多睡会儿没事的。”
田中爱子听了，更加不自如起来，脸上微红，含羞地说到：“妈，不早了，天早就亮了，再睡我自己都不好意思见您了。”
“看，又跟妈见外了，自己家有什么不好意思的，闺女，既然起来了，那不如先去洗漱间洗洗，洗漱用品你爸一早起来就给你买好了。”
“好的，妈！”
“那妈领你去。”
看到洗漱台上摆好的洗漱用品，田中爱子心里不自觉的一甜，心想：真是上辈子修来的好福气，让我因祸得福，遇到了待我这么真心的爸妈，一定不能辜负了他们。
……
片刻功夫，一名翩若惊鸿，婉若游龙，荣耀秋菊，华茂春松的女子走了出来。孙月娥看到后震惊了，心里不自觉地暗自称赞个不停。
“妈，我爸呢，我怎么没看到他？”田中爱子轻问到。
孙月娥闻言，才回过神来，微笑地回到：“闺女，你爸他上班去了，他工作忙，任务重，不得不去，要妈今天在家陪你。”
“哦，也对，都9点了！都怪我起得太晚了，妈，我是不是给您添麻烦了，你也有工作的，哎……”田中爱子一脸懊恼的道。
“闺女，说到哪里去了，不要跟妈这么客气了，显得多生分……”说着说着，孙月娥竟伤心起来。
田中爱子看到孙月娥两眼微红，连忙往她怀里一扑，撒娇道：“妈，女儿知道了，女儿错啦，您别生气了。”
孙月娥被田中爱子这么一撒娇，破涕为笑的说：“嗨，我说乖女儿，你还真会看脸色，妈现在对你真是又疼又爱了，以后千万不要再和妈说生分的话了，好吗？”
“妈，女儿保证，再也不敢了。”说完田中爱子往孙月娥脸上亲了一口。
哈哈……一阵欢快的笑声游荡在客厅里……
此时的孙月娥，心里有说不出的欣慰与开心，她是好久没有这么发自内心的笑了，天真烂漫的田中爱子，让她从连失两子的悲痛中走了出来。
……
两人吃过早餐，孙月娥说到：
“闺女，你爸昨晚说要给你的财富，他后来找到了，给你留在书桌上的，我们母女是先出去散散心，还是……”
田中爱子闻言，连忙说到：“妈，我要看看爸给我留的财富，女儿好奇的很呢。”
听田中爱子这么一说，孙月娥哈哈地笑到：“哈哈……好吧，闺女，说真的，我也好奇，那我们先去书房看看去。”
说完两人手拉手往书房走去。
孙月娥一边开电脑，一边对田中爱子笑着说：
“闺女，你爸这人，哈哈，有时候妈真不知道拿什么好话说他了。”
“妈，爸这叫笃实好学，态度这么认真，办事这么稳重，天底下能有几个他这样的？妈，我们应该感到自豪，是不是？”田中爱子满脸崇拜地说到。
孙月娥见说得这么好，眉开眼笑的说到：
“闺女你这么说，妈也觉得有那么点意思了。来，坐下看看，你爸说你看到后一定会喜欢的。”
田中爱子问闻言，也不客气，坐下去好奇的看了起来。
这不看不知道，一看吓得她一跳。这真是一笔丰厚的知识财富，全是张起明在他所在航天动力学方面的知识积累和笔记。
然而更让田中爱子吃惊的是，原本以为，以她目前现有的知识量，不好好学个两三年，是看不懂多少的，可没想到啊，没想到啊，她大脑对这些知识似乎了若指掌，竟一看就懂了，看到有的地方甚至是茅塞顿开起来。
这太不可思议了，这不科学啊，自己明明记得没学过导航方面的知识啊，可我怎么还是看得懂，就连那些公式，脑子也可以飞速地推导验算。
只见：
以航天器的位置和速度作为状态变量，其状态方程如下所示:
x?(t)=f(x，t)+ω(t)(1)
其中，x(t)是状态变量;f(x，t)是状态动力学的非线性函数;ω(t)是噪声项．
考虑地球非球形引力摄动j2项，动力学方程如下式:
χ?=ν_x
y?=ν_y
z?=ν_z
(ν_x)?=-μr_x/r^3[1-j_2(r_e/r)^2(7.5(r_z^2)/r^2-1.5)]
(ν_y)?=-μr_y/r^3[1-j_2(r_e/r)^2(7.5(r_z^2)/r^2-1.5)]
(ν_z)?=-μr_z/r^3[1-j_2(r_e/r)^2(7.5(r_z^2)/r^2-4.5)]（2）
其中，(x，y，z)是航天器的位置矢量;(ν_x，ν_y，ν_z)是航天器相对于地球质心的速度;ｒe为地球半径．方程(2)的解可直接用四阶runge-kutta法数值计算得到．1．2测量方程在j2000.0地心赤道惯性坐标系，将航天器测得时间与地球质心的时间延迟量作为测量方程，在j2000.0地心赤道惯性坐标系，将航天器测得时间与地球质心的时间延迟量作为测量方程，有
?t=t_e-t_sc=1/cn?r+ν(3)
则得到观测方程为：
y=hx+v(4)
h=1/c[h^2，o_nx3]
h?=[n_1，n_2，?，n_n]^t
n_i=[n_1^i，n_2^i，n_3^i]^t
……
……
随着脑子里飞速的验算着这些陌生的公式，田中爱子额头上的汗也冒了出来，匪夷所思的让她再也不敢看下去了。
我这是怎么回事，我怎么突然凭空会这些我原本不会的知识了？
不可思议、不可捉摸的事，让田中爱子渐渐地陷入了深深的沉思……
(本章完)