• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            CG@CPPBLOG

            /*=========================================*/
            隨筆 - 76, 文章 - 39, 評論 - 137, 引用 - 0
            數(shù)據(jù)加載中……

            我的SICP習(xí)題答案(1.7)

            對于正文中的 good-enough? 謂詞,設(shè)所求 x 的真實平方根為 xt,那么我們的依據(jù)是當(dāng) guess 從 1 趨向與 xt 的平方差 小于 c(0.001) 時,guess 近似于 xt。實際當(dāng) xt^2 也就是 x 足夠小時, guess 會 逐漸穩(wěn)定在 √c 附近。從下面實驗可以看出:

            > (sqrt (square 0.1))
            0.10032578510960607
            > (sqrt (square 
            0.05))
            0.054237622808967656
            > (sqrt (square 
            0.01))
            0.03230844833048122
            > (sqrt (square 
            0.005))
            0.031515954454847304
            > (sqrt (square 
            0.001))
            0.031260655525445276
            > (sqrt (square 
            0.0001))
            0.03125010656242753
            >
             
            因為 guess^2 < c + x , 當(dāng) x < c 時,guess 就更多的依賴于 c 了。

            對于特別大的數(shù),由于浮點數(shù)在特別大時,離散性非常明顯,相鄰的兩個數(shù)之間的差距會非常大,導(dǎo)致 |guess^2 - x| 始終 大于 c,計算便進入了 無限循環(huán)。

            比如計算 (sqrt 1e300)

            使用變化率改進后的代碼如下:

            (define (sqrt-new x)
              (sqrt-iter-new x 
            1.0 x))

            (define (sqrt-iter-new s1 s2 x)
              (if (enuf-new? s1 s2)
                  s2
                  (sqrt-iter-new s2 (improve s2 x) x)))

            (define (enuf-new? s1 s2)
              (< (/ (abs (- s1 s2)) s1) 
            0.001))

            可以測算幾個數(shù),驗證結(jié)果還是比較好的。

            > (sqrt-new (square 1e150))
            1.0000000000084744e+150
            > (sqrt-new (square 
            1e-150))
            1.0000000000084744e-150


            posted on 2008-03-11 00:07 cuigang 閱讀(1650) 評論(2)  編輯 收藏 引用 所屬分類: Lisp/Scheme我的SICP答案

            評論

            # re: 我的SICP習(xí)題答案(1.7)  回復(fù)  更多評論   

            我用的DrScheme (#lang planet neil/sicp)
            我定義的名字是SQRT,程序和你的殊途同歸
            比較奇怪的是
            如果我用(SQRT 1e8)沒問題
            如果我用(SQRT 100000000)就不行了,卡住了
            如果我用(SQRT 100000000.0)又可以了
            這種現(xiàn)象有點沒頭緒,不知道你的系統(tǒng)會不會這樣
            2010-02-01 12:08 | mxj

            # re: 我的SICP習(xí)題答案(1.7)  回復(fù)  更多評論   

            @mxj

            (define (sqrt-new x)
            (sqrt-iter-new x 1.0 x))

            這句估計你寫成了
            (define (sqrt-new x)
            (sqrt-iter-new x 1 x))

            用了整數(shù) 1 而不是浮點數(shù) 1.0
            2010-02-19 22:21 | cuigang
            久久国产免费观看精品3| 久久国产精品久久国产精品| 色综合久久88色综合天天 | 久久人人爽人人爽人人片AV不| 亚洲va国产va天堂va久久| 99久久精品国产综合一区| 亚洲午夜无码AV毛片久久| 久久精品a亚洲国产v高清不卡| 精品久久久无码中文字幕| 五月丁香综合激情六月久久| 国产AV影片久久久久久| 久久超碰97人人做人人爱| 久久亚洲中文字幕精品一区| 狠狠色丁香久久综合五月| 无码人妻久久一区二区三区蜜桃| 国产亚洲美女精品久久久久狼| 亚洲精品无码专区久久同性男| 久久精品男人影院| 久久精品www人人爽人人| 精品久久久久久久久免费影院| 国产精品一区二区久久精品无码 | 亚洲精品乱码久久久久久按摩| 国产亚洲美女精品久久久| 国产日产久久高清欧美一区| 日韩久久久久久中文人妻| 伊人久久国产免费观看视频| 青青草国产成人久久91网| 久久亚洲私人国产精品| 亚洲午夜久久久久久久久久| 久久久久久精品免费免费自慰| 日韩十八禁一区二区久久| 久久综合九色综合欧美就去吻| 久久99久久成人免费播放| 久久99精品久久久久久水蜜桃| 女人香蕉久久**毛片精品| 日韩欧美亚洲综合久久影院d3| 青青草原综合久久| 久久精品国产黑森林| 久久综合九色综合欧美就去吻| 性做久久久久久久久久久| 中文字幕精品久久|