SICP習題1.45 為什麼做average damp的次數需要大於等於log2n

阿新 • • 發佈：2018-11-13

簡單粗暴的標題，這個推導是我看了知乎的兩篇文章之後，結合自己的理解寫的，多畫了幾個影象，力求通俗易懂

知乎的參考文章如下

https://zhuanlan.zhihu.com/p/25601871

https://www.zhihu.com/question/28838814/answer/42283723

最初我們需要求解的fixed-point函式如下

$y=\frac{x}{y^{n-1}}$

經過m次average damp後，fixed-point函式變為

$y=\frac{1}{2^m}(\frac{x}{y{n-1}}+(2^m-1)y)$

每次fixed-point迭代計算，第k次迭代可以寫成如下形式，直到 $y_{k+1}\approx y_k$ 迭代結束

$y_{k+1}=\frac{1}{2^m}(\frac{x}{y_k{n-1}}+(2^m-1)y_k)$

以 $y_{k+1}$ 作為自變數，變換一下上面的公式，可以得到函式，這條直線的斜率為 $2^my_k^{n-1}$

$g(y_{k+1})=2^my_k^{n-1}y-x-(2^m-1)y_k^n$

每次迭代求 $y_{k+1}$ 就等同於求函式 $g(y_{k+1})$ 的零點

設函式 $f(y)=y^n-x$ ，該函式零點就是我們所要求的平方根，去掉 $g(y_{k+1})$ 的下標k+1，得到 $g(y)$

為了更好地理解，以一個例子代入，求解 $\sqrt[2]{16}$ ，初始值為1

1、不做average damp

兩個函式是這樣

$\left\{ \begin{aligned} f(y)&=y^2-16\\ g(y)&=y_ky-16 \end{aligned} \right.$

第一次迭代， $y_1=1$ ，得出 $y_2=16$

第二次迭代， $y_2=16$ ，得出 $y_3=1$

第三次迭代， $y_3=1$ ，得出 $y_3=16$

陷入了死迴圈，非常夭壽

2、做一次average damp

兩個函式是這樣

$\left\{ \begin{aligned} f(y)&=y^2-16\\ g(y)&=2y_ky-16-y_k^2 \end{aligned} \right.$

第一次迭代， $y_1=1$ ，得出 $y_2=8.5$

第二次迭代， $y_2=8.5$ ，得出 $y_3=5.19$

第三次迭代， $y_3=5.19$ ，得出 $y_3=4.14$

每次迭代的直線如下，可以看到直線的零點和曲線的零點越來越接近

對比一下上面的兩個圖的區別，不做average damp的時候，直線是在曲線零點的兩邊晃來晃去導致震盪不收斂，做了average damp後，直線一直在曲線的一側，慢慢貼近

如何保證直線不在曲線零點兩邊晃來晃去呢

看一下知乎大神給出的這個圖，藍色曲線表示 $f(y)$ ，其他直線表示 $g(y)$ ，紅點表示 $(y_k,y_k^n-x)$ 這個點，曲線和直線是必然通過這個點的，代入函式驗證一下就可以知道

現在看綠線，綠線的零點就是下一次的 $y_{k+1}$ ，和這次的 $y_k$ 分佈在了曲線零點的兩端，這兩個y值我用綠點標出來了，由上述分析可以知道，每次迭代y都分佈在曲線零點兩端就會引起震盪，所以綠線這種情況是不能收斂的，紅線和黃線的零點一直分佈在曲線零點的一側，可以收斂

很直觀可以看出來，只要每次在曲線直線交點處，直線的斜率大於等於曲線斜率即可，上面已經說了交點就是 $(y_k,y_k^n-x)$

曲線斜率為 $ny^{n-1}$ ，在交點處，斜率為 $ny_k^{n-1}$

直線斜率為 $2^my_k^{n-1}$

因此滿足 $2^m\geqslant n$ ，即 $m\geqslant log_2{n}$

SICP習題1.45 為什麼做average damp的次數需要大於等於log2n

簡單粗暴的標題，這個推導是我看了知乎的兩篇文章之後，結合自己的理解寫的，多畫了幾個影象，力求通俗易懂知乎的參考文章如下 https://zhuanlan.zhihu.com/p/25601871 https://www.zhihu.com/question/28838814/answe

SICP 習題（1.41）解題總結

結果調用 -m name 方法定義 height double plain SICP 習題1.41 看似和周邊的題目沒有關系，突然叫我們去定義一個叫double的過程，事實上這道題的核心還是高階函數。題目要求我們定義一個過程double，它以一個過程作為參數，這個作為

SICP習題解答2.1-2.6

ex2.1-2.4 #lang racket ; exercise 2.1 (define (make-rat n d) (define g (gcd n d)) (let ((g1 (if (< d 0) (- g) g))) (cons (/ n

SICP習題解答1.9-1.19

ex1.9-1.10 #lang racket ; exercise 1.9 ;; 第一個是遞迴計算過程，第二個是迭代計算過程 ; exercise 1.10 (define (A x y) (cond ((= y 0) 0) ((= x 0)

SICP 習題（1.43）解題總結

ews define 通過嵌套總結 com plain 例如 pan SICP 習題 1.43 是前面兩道題的延續，習題要求我們定義一個過程(repeat f n

C/C++算法競賽入門經典Page16 習題1-4 正弦和余弦

pan 輸出定義 c++ col return 算法競賽入門經典算法 printf 題目:輸入正整數 n（n<360）,輸出n度的正弦余弦數值。提示：用數學函數需要用到的數學函數: acos,sin,cos 函數頭文件:math.h #include <

C/C++算法競賽入門經典Page16 習題1-3 連續和

n) 練習 clu amp ... blog 變量入門經典 c++ 題目:輸入正整數 n，輸出1+2+3+...+n的值。提示：目標是解決問題，而不是練習編程為了輸出1+2+3+...+n,可以使用首項加末項乘以項數除以2的公式 (1)聲明1+2+3+...+n的和m,

C/C++算法競賽入門經典Page16 習題1-5 打折

span 博客保留 double int 保留兩位小數 c++ 小數 pan 題目:一件衣服 95元，若消費滿300元，可打八五折。輸入購買衣服件數，輸出需要支付的金額（單位元），保留兩位小數. 通過分支結構可以判斷是否要打八五折. (1)聲明衣服件數a,需要支付的金額b

C/C++算法競賽入門經典Page16 習題1-6 年份

blog std ear code http ret 判斷 ++ class 題目:輸入年份，判斷是否為潤年，如果是，輸出yes，否，輸出no。提示：簡單的除以4是不夠的分析閏年的條件: (1)能被4整除. (2)整百的年份一定要能被400整除. 步驟: (1)聲明一個i

C/C++算法競賽入門經典Page16 習題1-6 三角形

include 是否 angle not 習題 bsp 如果著作權能夠題目:輸入三角形 3 條邊長度值（均為正整數），判斷是否能為直角三角形的3個邊長。如果可以，則輸出yes，不能輸出no，如果無法構成三角形，則輸出 not a triangle 分析能構成三角形的條

習題1

由於整除循環之間累乘 src color 技術分享 break 題目如下編程實現：設n為1~1500間的所有整數，按照如下要求計算和輸出。 1、判斷並輸出n中能被3和5和7同時整除的數的個數； 2、判斷並輸出n中能被3或5或7整除的數的個數；

2014-10-30NOIP復習題1

tarjan oid incr 排序 int() 高精有向圖詳細 one Problem 1 Graph (graph.cpp/c/pas) 【題目描述】給出 N 個點，M 條邊的有向圖，對於每個點 v，求 A(v) 表示從點 v 出發，能到達的編號最大的點。

C 語言習題 1-12

ret blog 1-1 main stat pan tdi 編寫一個程序 std 練習 1-12 編寫一個程序，以每行一個單詞的形式打印其輸入。 1 #include <stdio.h> 2 3 #define IN 1 4 #define OUT

SCIP習題 1.4

運算 fin sci def + - 語句習題 define 1.4 (define (a-plus-abs-b a b) (if (> b 0) (+ a b) (- a b)) )(define (a-plus-abs-b-book a b) ((if (

SCIP習題 1.1

href def con else 代碼 http adt nts 後來 (define a 3) (define b (+ a 1)) b (+ a b(* a b)) (= a b) (if (and (> b a) (< b (* a b))) b

SCIP習題 1.3

習慣 define 函數感覺 sum 形參 fin 1.3 and (define (SumOfTwoLarge x y z) (if (and (<= x y) (<= x z)) (+ (* y y) (* z z)) (SumOfT

SCIP習題 1.10

意思 con 習題 else 什麽 def cond sci efi (define (A x y) (cond ((= y 0) 0) ((= x 0) (* 2 y)) ((= y 1) 2) (else (A (- x 1)

sicp 習題

lan 版本 define sca fine cal ack asc row 1.11 求f(n)=f(n-1)+2*f(n-2)+3*f(n-3) 1 #lang racket 2 (define (fff n) 3 (define (fff-iter a

生成器習題1

eve pla 根據 open display pycharm expr 輸出結果次循環預熱 li = [x * x for x in range(10)] print(li) li = (x * x for x in range(10)) print(

在CentoOS 7上使用Nginx 1.12做TCP反向代理

toc repo pac nginx ace yum源 linu 配置yum _for 1、關閉selinux、firewall2、修改/etc/security/limits.conf添加以下內容：* soft nofile 65536* hard nofi

SICP習題1.45 為什麼做average damp的次數需要大於等於log2n

相關推薦