Кейс: нужно численно решить уравнение x^3 - 2x + 2 = 0 — сравните применимость методов Ньютона, секущих и бисекции, назовите критерии выбора метода и условия сходимости для каждого
Функция и интервал: f(x)=x3−2x+2f(x)=x^3-2x+2f(x)=x3−2x+2. Имеется знакопеременация на [−2,−1][-2,-1][−2,−1]: f(−2)=−2, f(−1)=3f(-2)=-2,\; f(-1)=3f(−2)=−2,f(−1)=3, значит есть корень в этом интервале (единственный вещественный корень). Про методы — кратко и по пунктам. 1) Бисекция - Условие применимости: fff непрерывна на [a,b][a,b][a,b] и f(a)f(b)<0f(a)f(b)<0f(a)f(b)<0. Для нашего случая можно взять [a,b]=[−2,−1][a,b]=[-2,-1][a,b]=[−2,−1]. - Условие сходимости: гарантированная (глобальная) сходимость к корню при указанных условиях. - Скорость: линейная, погрешность уменьшается вдвое за итерацию: после nnn итераций погрешность ≤b−a2n\le \frac{b-a}{2^n}≤2nb−a. - Преимущества/недостатки: очень надёжна, не требует производной; медленная. - Рекомендация для задачи: надёжный «фолбэк»; даст корень с любой точностью, но потребуется много итераций. 2) Метод Ньютона - Формула: xk+1=xk−f(xk)f′(xk)x_{k+1}=x_k-\dfrac{f(x_k)}{f'(x_k)}xk+1=xk−f′(xk)f(xk), где f′(x)=3x2−2f'(x)=3x^2-2f′(x)=3x2−2. - Условие применимости: fff дифференцируема, нужно вычислять f′f'f′. - Условие сходимости (локальное): если корень простой (f′(α)≠0f'(\alpha)\ne0f′(α)=0) и начальное x0x_0x0 достаточно близко к α\alphaα, то сходимость квадратичная: ∣xk+1−α∣≤C∣xk−α∣2|x_{k+1}-\alpha|\le C|x_k-\alpha|^2∣xk+1−α∣≤C∣xk−α∣2. Более формально — при f∈C2f\in C^2f∈C2 в окрестности корня и f′(α)≠0f'(\alpha)\ne0f′(α)=0 существует окрестность, в которой любой старт сходится квадратично. - Риски: может расходиться или прыгать к другому корню при плохом старте; проблемы если f′(x)f'(x)f′(x) близко к нулю. - Для данной функции: в корне α≈−1.7693 \alpha\approx -1.7693α≈−1.7693 имеем f′(α)=3α2−2≈7.39f'(\alpha)=3\alpha^2-2\approx 7.39f′(α)=3α2−2≈7.39 (не мало), значит метод должен быстро сходиться при разумном старте (например x0=−1.8x_0=-1.8x0=−1.8 или −1.7-1.7−1.7). Быстрая (квадратичная) сходимость и малое число итераций. 3) Метод секущих - Формула: xk+1=xk−f(xk)(xk−xk−1)f(xk)−f(xk−1)x_{k+1}=x_k-\dfrac{f(x_k)(x_k-x_{k-1})}{f(x_k)-f(x_{k-1})}xk+1=xk−f(xk)−f(xk−1)f(xk)(xk−xk−1). - Условие применимости: нужен только доступ к значениям fff (две начальные точки x0,x1x_0,x_1x0,x1). - Условие сходимости: локальная сходимость порядка φ=1+52≈1.618\varphi=\dfrac{1+\sqrt5}{2}\approx1.618φ=21+5≈1.618 для простого корня при f∈C2f\in C^2f∈C2 и достаточно близких начальных приближениях. Может потерпеть неудачу, если знаменатель близок к нулю или начальные точки плохи; сохраняет брутальность меньше чем бисекция. - Преимущества/недостатки: быстрее бисекции, не требует производной; медленнее Ньютона по скорости, но дешевле, если вычисление производной дорого. - Для задачи: хорошие стартовые точки, например x0=−2, x1=−1.5x_0=-2,\; x_1=-1.5x0=−2,x1=−1.5, обычно приведут к быстрой сходимости к тому же корню. Критерии выбора в практике - Нужна гарантия и простота → бисекция (например при наличии только бракетирования). - Доступен/недорог в вычислении производной и есть разумный начальный приближение → метод Ньютона (быстрее). - Производная недоступна или дорого вычисляется, но нужен быстрый метод → секущие (компромисс). - Для надёжности используют защищённые схемы: комбинируют Ньютон/секущие с бисекцией (например, если шаг Ньютона выходит из бракетирования — выполнить бисекцию). Короткое практическое заключение по уравнению x3−2x+2=0x^3-2x+2=0x3−2x+2=0: - Надёжно: бисекция на [−2,−1][-2,-1][−2,−1]. - Быстро и эффективно: Ньютон с x0=−1.8x_0=-1.8x0=−1.8 (или −1.7-1.7−1.7) — квадратичная сходимость, корень α≈−1.7693 \alpha\approx -1.7693α≈−1.7693. - Без производной: секущие с x0=−2, x1=−1.5x_0=-2,\;x_1=-1.5x0=−2,x1=−1.5 — сверхлинейная сходимость порядка ≈1.618\approx1.618≈1.618.
f(x)=x3−2x+2f(x)=x^3-2x+2f(x)=x3−2x+2. Имеется знакопеременация на [−2,−1][-2,-1][−2,−1]: f(−2)=−2, f(−1)=3f(-2)=-2,\; f(-1)=3f(−2)=−2,f(−1)=3, значит есть корень в этом интервале (единственный вещественный корень).
Про методы — кратко и по пунктам.
1) Бисекция
- Условие применимости: fff непрерывна на [a,b][a,b][a,b] и f(a)f(b)<0f(a)f(b)<0f(a)f(b)<0. Для нашего случая можно взять [a,b]=[−2,−1][a,b]=[-2,-1][a,b]=[−2,−1].
- Условие сходимости: гарантированная (глобальная) сходимость к корню при указанных условиях.
- Скорость: линейная, погрешность уменьшается вдвое за итерацию: после nnn итераций погрешность ≤b−a2n\le \frac{b-a}{2^n}≤2nb−a .
- Преимущества/недостатки: очень надёжна, не требует производной; медленная.
- Рекомендация для задачи: надёжный «фолбэк»; даст корень с любой точностью, но потребуется много итераций.
2) Метод Ньютона
- Формула: xk+1=xk−f(xk)f′(xk)x_{k+1}=x_k-\dfrac{f(x_k)}{f'(x_k)}xk+1 =xk −f′(xk )f(xk ) , где f′(x)=3x2−2f'(x)=3x^2-2f′(x)=3x2−2.
- Условие применимости: fff дифференцируема, нужно вычислять f′f'f′.
- Условие сходимости (локальное): если корень простой (f′(α)≠0f'(\alpha)\ne0f′(α)=0) и начальное x0x_0x0 достаточно близко к α\alphaα, то сходимость квадратичная: ∣xk+1−α∣≤C∣xk−α∣2|x_{k+1}-\alpha|\le C|x_k-\alpha|^2∣xk+1 −α∣≤C∣xk −α∣2.
Более формально — при f∈C2f\in C^2f∈C2 в окрестности корня и f′(α)≠0f'(\alpha)\ne0f′(α)=0 существует окрестность, в которой любой старт сходится квадратично.
- Риски: может расходиться или прыгать к другому корню при плохом старте; проблемы если f′(x)f'(x)f′(x) близко к нулю.
- Для данной функции: в корне α≈−1.7693 \alpha\approx -1.7693α≈−1.7693 имеем f′(α)=3α2−2≈7.39f'(\alpha)=3\alpha^2-2\approx 7.39f′(α)=3α2−2≈7.39 (не мало), значит метод должен быстро сходиться при разумном старте (например x0=−1.8x_0=-1.8x0 =−1.8 или −1.7-1.7−1.7). Быстрая (квадратичная) сходимость и малое число итераций.
3) Метод секущих
- Формула: xk+1=xk−f(xk)(xk−xk−1)f(xk)−f(xk−1)x_{k+1}=x_k-\dfrac{f(x_k)(x_k-x_{k-1})}{f(x_k)-f(x_{k-1})}xk+1 =xk −f(xk )−f(xk−1 )f(xk )(xk −xk−1 ) .
- Условие применимости: нужен только доступ к значениям fff (две начальные точки x0,x1x_0,x_1x0 ,x1 ).
- Условие сходимости: локальная сходимость порядка φ=1+52≈1.618\varphi=\dfrac{1+\sqrt5}{2}\approx1.618φ=21+5 ≈1.618 для простого корня при f∈C2f\in C^2f∈C2 и достаточно близких начальных приближениях. Может потерпеть неудачу, если знаменатель близок к нулю или начальные точки плохи; сохраняет брутальность меньше чем бисекция.
- Преимущества/недостатки: быстрее бисекции, не требует производной; медленнее Ньютона по скорости, но дешевле, если вычисление производной дорого.
- Для задачи: хорошие стартовые точки, например x0=−2, x1=−1.5x_0=-2,\; x_1=-1.5x0 =−2,x1 =−1.5, обычно приведут к быстрой сходимости к тому же корню.
Критерии выбора в практике
- Нужна гарантия и простота → бисекция (например при наличии только бракетирования).
- Доступен/недорог в вычислении производной и есть разумный начальный приближение → метод Ньютона (быстрее).
- Производная недоступна или дорого вычисляется, но нужен быстрый метод → секущие (компромисс).
- Для надёжности используют защищённые схемы: комбинируют Ньютон/секущие с бисекцией (например, если шаг Ньютона выходит из бракетирования — выполнить бисекцию).
Короткое практическое заключение по уравнению x3−2x+2=0x^3-2x+2=0x3−2x+2=0:
- Надёжно: бисекция на [−2,−1][-2,-1][−2,−1].
- Быстро и эффективно: Ньютон с x0=−1.8x_0=-1.8x0 =−1.8 (или −1.7-1.7−1.7) — квадратичная сходимость, корень α≈−1.7693 \alpha\approx -1.7693α≈−1.7693.
- Без производной: секущие с x0=−2, x1=−1.5x_0=-2,\;x_1=-1.5x0 =−2,x1 =−1.5 — сверхлинейная сходимость порядка ≈1.618\approx1.618≈1.618.