第30章 ，计算周期解的数值方法（1/1）

“没关系，慢慢来改，一天两天改不完很正常，不过你要记住我说的要求，可不能瞎改，要不然我肯定会让你重改的...”

隔天一早，王多鱼回到主编办公室内，就主动跟周老太太汇报工作，结果老太太就对他说了这么一句。

如此就更好，他就可以更加心安理得地磨洋工，反正每一天都有一块钱的工资，过段时间就直接将那一版结局交给老太太审核就是了。

中午的时候，王多鱼去吃午饭时碰巧看到了傅吉祥，便询问对方有没有多余的副食品粮票，对方还真有多余的票，而且很友好地换给了他。

“你回东北的时候，如果想要买这些糕点的话，到时候缺票的时候可以找我，我提前帮你换...”

“好嘞，谢谢傅队长，太感谢了！”

有了这些副食品粮票，王多鱼就可以去购买桂香村食品店的糕点了，一个是他也嘴馋，二是想要还掉钱大妈前天晚上那碗米饭的人情。

时间一晃而过，转眼到了傍晚，跟昨天一样，王多鱼还是溜达去了西单，今天晚上准备吃这边非常出名的烤鸭。

听说烤鸭最好吃的店是便宜坊，也有说是全聚德，但是对王多鱼来说都差不多，因为这时代的饭菜都是真材实料，不用心的门店还比较少见。

反正时间充裕，就在西单这一带随便吃就是了，等过段时间交稿之后有时间了，再去便宜坊或者全聚德试一试。

“嚯，真满足！”

一只烤鸭全被他吃进肚子里了，皮是全都片出来了，还剩下的鸭架子一半油炸，另一半做成汤，就这么烤鸭三吃，能不满足么？

二十二岁的年轻身体，胃口是真的好，肚子也是真的大，上辈子他才三十岁的时候，每天小半碗米饭，菜也不多吃，因为怕体重超标，而且晚饭还经常忘了吃。

就这，他BMI指数还是超了一点，体检报告上面说他的体重就是肥胖，表面上看起来他就是壮，而不是胖。

吃饱喝足之后，王多鱼这才漫步溜达着出来，在西单这边逛了一会儿，又买了几本书，然后来到桂香村食品店门口排队购买糕点。

喜欢吃糕点的人很多，得排队购买才行。

“同志，麻烦帮我拿一盒椒盐三角酥，两盒绿豆糕，两盒太师饼，再来一盒杏仁酥就好了。”

等到王多鱼的时候，他一开口，顿时让前后其他顾客和售货员都张大了嘴巴，别人买糕点都是几块几块地买，他可倒好，直接来进货啊。

总共六盒呢，那可真是不少了。

“同志你有这么多票么？哦哦，我数一下，行，我再跟你确认一下是三角酥和杏仁酥各一盒，绿豆糕和太师饼各两盒，对么？”

售货员从王多鱼手中接过钱和票之后，还再次确认了一遍，排在他身后的其他顾客都伸长了脖子，眼巴巴地看着，其中几个七八岁的小朋友看到这一幕，扁着嘴巴跟他们父母委屈地表达了他们的担忧：那个叔叔把糕点都买完了！

还别说，这个时间点已经很晚了，食品店也都差不多快要下班了，货架上摆放的糕点确实没有很多，可是队伍还很长。

王多鱼没有听到小朋友们的怨念和哀嚎，他提着六盒糕点回了大方胡同三号院，进了院子之后就直奔钱大妈家，直接给了两盒糕点给对方。

一礼还一礼，相信对方应该能明白他的意思。

钱大妈不想要这两盒糕点，但王多鱼举着他手里的另外四盒糕点，表示他买多了，自己一个人也吃不完，现在可是夏天，糕点放久了容易变味。

这个理由无懈可击，钱大妈就算不乐意也没辙。

张大妈、鲁大妈等邻居们看到这一幕，都不由暗乐不已：钱大妈也算是吃瘪了吧？以为人家王作家年轻好拿捏，嘿，没想到人家挺会来事儿的，滑不溜秋，看来是甭想使唤他咯。

关上门，王多鱼瞥了一眼桌子的稿件，昨天是怎么样的，现在还是怎么样，他根本没有动过，韩崇谦的文稿都不需要看，绝对是浪费时间。

蜡烛、煤油灯和电灯全都用上，王多鱼开始了今晚的工作，继续昨天晚上未完成的数学论文：

由于解初值问题时，相对于初值的选取具有内在的不稳定性，因而不妨采用多级打靶(Multiplicate shooting)，即将原区域(0，1)分隔成许多子区间，在这些子区间上求解初值问题其效果会更好些。

最笨拙的方法是直接利用求解常微分方程初值问题的任何一种数值积分方法，如最简单的是Euler折线法，从任何一个不是平衡点的初值出发，一个时间步长接着一个时间步长地往前计算，若能在屏幕上动态显示每一步的计算结果为最好。

一旦发现周期轨线出现时，就输出其计算结果，仔细观察其数值大小。近似相等的数值重复地出现时，就可以确定出周期解和周期的大小。这种做法显然带有一定的盲目性，会造成计算时间的浪费，还需要一定的人工干预。虽然笨抽一点，但能有效的计算出结果，当然它不属于我们要研究的方法之列。

王多鱼写的这篇论文是关于系统叙述计算常微分方程自治系统周期解的数值方法，因为在物理、化学、生物、电子学和工程技术等问题中经常可以遇到周期震荡现象，并且这种现象在自然界同样普遍存在着。

而在数学中求解微分方程是周期解是一个古老且有困难的问题，早在一九零五年时，戴维希尔伯特第16问题就是有关于微分方程周期解存在性的判定问题。

然而至今依然没有什么进展，并不是在过去半个多世纪来没有数学家去拆解，而是难度之大远超想象。

特别是在六十年代之后，计算机的普及和计算机技术的发展，大家的研究思路和方法已经发生改变，在实际问题的研究中，比如生物化学中的布鲁塞尔振子，大家的兴趣并不在于周期解存在性的理论证明，更多是关心周期解的位置、形状及周期的大小。

所以王多鱼编写的关于计算方法的研究就显得尤为重要了。