VPF::Нахождение корня методом деления пополам

Elfet

Дата 3.5.2005, 20:08 (ссылка)

(нет голосов)

Загрузка ...

Белый и Пушистый

Профиль
Группа: Awaiting Authorisation
Сообщений: 3776
Регистрация: 2.4.2003

Репутация: нет
Всего: 16

Вот сказали сделать smile

(процедуру нахождение корня методом деления пополам)
Как это выглядит и работает на листке я понимаю smile

а вот запрограмить - никак smile

Прошу помощи! smile

Вот полное задание:

Цитата

Построить процедуру нахождения в интервале [a,b] единственного корня уровнения вида f(x)=0, с точностью E (для значения функции в окресности корня) или N (для значения в корне) , применив метод деления пополам.

Вот то что у меня не работает smile

Код


#include <iostream>
#include <cmath>
using namespace std;

float found(float a, float b, float eps);
float f(float x);
float mod(float x){return sqrt(x*x);};

void main()
{
    int e;
    float a,b,eps;
    
    cout << "Found X." << __DATE__ << endl;
    cout << "Write 'a' and push Enter:" << endl;
    cin >> a;
    cout << "Write 'b' and push Enter:" << endl;
    cin >> b;
    cout << "Write 'eps' and push Enter:" << endl;
    cin >> eps;
    
    cout << endl << "f( " << found(a,b,eps) << " )=0";
    
    cin >> e;
}

float found(float a, float b, float eps)
{

    float x,y,z;
    z=f(a);
    do
    {
     x=(a+b)/2; y=f(x);
     if(y*z>0)
     { 
         a=x;
     }
     else 
     {
         b=x;
     }

    }while( (mod(y) < eps) && (mod(a-b) < eps) );

    return x;
}

float f(float x)
{
    return 3*x-2;
}

--------------------

PHP Чат Скрипт

Elfet

Дата 4.5.2005, 19:33 (ссылка)

(нет голосов)

Загрузка ...

Белый и Пушистый

Профиль
Группа: Awaiting Authorisation
Сообщений: 3776
Регистрация: 2.4.2003

Репутация: нет
Всего: 16

Нужна помощь! smile

Зачёт в четверг! А у меня всё не выходит smile

--------------------

PHP Чат Скрипт

Pakshin A. S.

Дата 4.5.2005, 20:58 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Участник Клуба
Сообщений: 5056
Регистрация: 16.2.2003

Репутация: 5
Всего: 61

А можно объяснить сам метод нахождения?

Elfet

Дата 5.5.2005, 21:19 (ссылка)

(нет голосов)

Загрузка ...

Белый и Пушистый

Профиль
Группа: Awaiting Authorisation
Сообщений: 3776
Регистрация: 2.4.2003

Репутация: нет
Всего: 16

Всё разобрался!

Цитата(Elfet @ 3.5.2005, 21:08)

}while( (mod(y) < eps) && (mod(a-b) < eps) );

на

Цитата(Elfet @ 3.5.2005, 21:08)

}while( (mod(y) > eps) && (mod(a-b) > eps) );

Сам метод очень прост:
Мы делем участок пополам и смотрим который ближе у ОХ, затем берём этот участок и также анализируем smile

Пока не дойдём до точности вещественного нуля smile

ну или до одной из границ smile

--------------------

PHP Чат Скрипт

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »