Реферат: Решение одного нелинейного уравнения
.
Иначе можно получить уравнение (4) следующим способом: левую и правую часть уравнения (1) умножить на произвольную константу и прибавить к левой и правой части х, т.е. получаем уравнение вида: (6), где .
На заданном отрезке [a; b] выберем точку х0 - нулевое приближение - и найдем: х1 = f (x0 ), потом найдем: х2 = f (x1 ), и т.д. Таким образом, процесс нахождения корня уравнения сводится к последовательному вычислению чисел: хn = f (xn-1 ) n = 1,2,3… Если на отрезке [a; b] выполнено условие: |f ' (x) |<=q<1 то процесс итераций сходится, т.е. . Процесс итераций продолжается до тех пор, пока |xn - xn-1 |<=, где - заданная абсолютная погрешность корня х. При этом будет выполняться: .
Применение метода к конкретной задаче (анализ).
Дано уравнение вида x² - ln (1+x) - 3 = 0 при x [2,3]. Задача состоит в том, чтобы решить это нелинейное уравнение 4 известными способами: метод половинного деления, метод касательных, метод секущих и метод простой итерации.
Изучив методы и применив их к данному уравнению приходим к такому выводу: при решении данного уравнения 4 известными способами результат одинаков во всех случаях. Но количество итераций при прохождении метода значительно отличается. Зададим приближенную точность = . Если в случае половинного деления количество итераций составляют 20, при методе простых итераций равно 6, при методе секущих они составляют 5, а при методе касательных их количество равно 4. Из полученного результата видно, что более эффективным методом является метод касательных. В свою очередь метод половинного деления является более неэффективным, затрачивающий больше времени на выполнение, но являющийся самым простым из всех перечисленных методов при исполнении. Но не всегда результат будет таковым. Подставляя другие нелинейные уравнения в программу, в результате получается, что при методе простой итерации при разных видах уравнений количество итераций колеблется. Количество итераций может быть значительно больше, чем в методе половинного деления и меньше, чем в методе касательных.
Листинг программы:
1. Метод половинного деления
#include <stdio. h>
#include <math. h>
#include <conio. h>
#define e 0.000001
double func (double x)
{
return ( ( ( (x*x) - (log (1+x))) - 3));
}
void main ()
{
FILE *res;
res=fopen ("bisekciy. txt","w");
int x1,x2;
double x;
int k;
k=0;
x1=2;
x2=3;
x=0;
double a,b,c;
a=x1;