Форум программистов, компьютерный форум CyberForum.ru

Метод бисекции для нахождения корней в нелинейном уравнении - C++

Восстановить пароль Регистрация
 
FloraWinx
 Аватар для FloraWinx
4 / 4 / 1
Регистрация: 05.10.2013
Сообщений: 122
20.02.2014, 02:12     Метод бисекции для нахождения корней в нелинейном уравнении #1
Решил реализовать метод бисекции для нахождения корней в нелинейном уравнении.
Открыл псевдо код на википедии, но почему то оно у меня зацикливаеться, помогите найти ошибку.
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
#include <iostream>
#include <conio.h>
#include <math.h>
 
 
 
double sign(double mak){
    double x=mak;
    if (x > 0) return 1;
if (x < 0) return -1;
return 0;}
 
using namespace std;
 
double bisection_count(double end , double start , double eps ){
    double res;
    double temp;
    double temp_midle;
    double f_st , f_end;
    f_st=start; 
    f_end= end;
    double f_mid_i;
 
    f_st = (pow(start,4) +3*pow(start,3)-2*start+4);
    f_end = ( pow(end,4) +3*pow(end,3)-2*end+4);
    if(f_st==0) {
        cout<<"Root equels"<<endl;
        return f_st;}
 
    else if (f_end==0) {
        cout<< "Root equels 0"<<endl; 
         return f_end;
    }
 
 
    temp=end-start;
    temp_midle=(end-start)/2;
 
    double f_midle =    ( pow(temp_midle,4) +3*pow(temp_midle,3)-2*temp_midle+4);
 
    while (fabs(f_midle)>eps){
 
        temp=temp/2;
        temp_midle=end+temp;
 
        f_mid_i= ( pow(temp_midle,4) +3*pow(temp_midle,3)-2*temp_midle+4);
 
      if(sign(f_st)!=sign(f_mid_i))
      {end=temp_midle;}
      cout<<"Root is found:";
    }
 
    return temp_midle;
}
 
 
int main(){
double start;
double end;
double eps;
cin>>eps>>start>>end;
double z;
z=bisection_count(end,start,eps);
cout<<"res:"<<z;
getch();
 
    return 0;}
вот псевдо код который я использовал
Начало.
Ввод xn, xk, epsy.
Если F(xn) = 0, то Вывод (корень уравнения — xn).
Если F(xk) = 0, то Вывод (корень уравнения — xk).
dx := xk — xn.
Пока |F(xi)| > epsy повторять:
dx := dx / 2;
xi := xn + dx;
если sign(F(xn)) ≠ sign(F(xi)), то xk := xi;
иначе xn := xi.
конец повторять
Вывод (Найден корень уравнения — xi с точностью по y — epsy).
Конец.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.02.2014, 02:12     Метод бисекции для нахождения корней в нелинейном уравнении
Посмотрите здесь:

для нахождения корней квадратного уравнения использовать функцию C++
C++ метод бисекции C++
Метод итераций для нахождения обратной матрицы C++
Метод бисекции C++
Нахождения корней уравнения: метод половинного деления (бисекции) или метод хорд C++
Определить функцию, позволяющую распознавать наличие вещественных корней в квадратном уравнении C++
Определить функцию, позволяющую распознавать наличие вещественных корней в квадратном уравнении C++
C++ Метод бисекции

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
UnsKneD
алкокодер
 Аватар для UnsKneD
153 / 149 / 11
Регистрация: 27.12.2012
Сообщений: 548
20.02.2014, 07:06     Метод бисекции для нахождения корней в нелинейном уравнении #2
FloraWinx,
У вас в цикле ни одно из значений fabs(f_midle)>eps не изменяется.
Обратите внимание на
Пока |F(xi)| > epsy повторять: и xi := xn + dx
Yandex
Объявления
20.02.2014, 07:06     Метод бисекции для нахождения корней в нелинейном уравнении
Ответ Создать тему
Опции темы

Текущее время: 20:42. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru