平方根最佳逼近的计算机实现.pdf

第 18智 第2期 高 师理 科 学 刊 Vo1．18 No．2 journalotScMlceofT~=hers CoIl~ andUniversity ]998年 6月 June 1998 平方根最佳逼近的计算机实现 李希文 徐肇玉 (二)2钾． (一 ， — — 一 一 (齐齐暗尔大学啦学幕) 摘要 利用垒 PeLI序列 ．给 出了对非完全平方数N(N为整数)开平方的快速最佳有 理逼近算法． ． 美羹词Pell序列竖 蛆拱现 劲枨 最位 定义11令·(者 竺 三方根逼近定理本·文将给出计算机实现方法· 理 z y)是 PeIl方程 ～ ，一1 ( 为非完全平方数) (1) -a ． 一 的基本解．记 一玉 (xl ,／-NYD N Pell √ 一 + ，则 {只j称为 的全 序列． _ 定理 I。。关于非平方正整数N的全Pell序列 {尸_}及其子列是√ 的最佳过剩有理逼 近序列． 为方便，我们给出 定理2跚 若( ) 是(1)的任一组正整数解，则(∞ )_)必是f1)的一组正整数解，且 宅，都是全Pell序列的某项．其中 ‰ ．毫+Nd． (2) 现在应用定理 1，定理2，给出如下计算、／ (Ⅳ为非完全平方数)的最佳逼近算法 1．求出(1)的任一组正整数解．方法有二． 查表法：例如查啪． 计算法：取 =1，2，3，…，计算z(=√ 研 )，得出一组正整数解(z ． 2．令 -r， 一Y．应用公式(2)．(3)，计算 ，‰ 3．反复进行 ：令 一 ． 一yz 4一竺即为√Ⅳ的最佳有理逼近．其精度估值为： I 一√ =—— L 0 j +-／~yi 瑰给出此方法的计算机参考程序． 奉文1998年 4月 15日收到． 尊 2期 李希文辱：平 7 DmfnASDD1IBI．E xASDD1瑚 E．YASD0 7BI．E DIM PASDOUBLE．aASDOUBLE．bASDOUBLE．plASDOIYBLE．p2ASDOLrBLE a0王AR：CLS TNPI丁1-。PIea8einputanumber(O)：：‘n p~SQR缸) IFp=1wr )THEN [注1] PR rr n‘= ‘：SQR缸) 日 范 y