Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
FloraWinx
4 / 4 / 1
Регистрация: 05.10.2013
Сообщений: 123
#1

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

20.02.2014, 02:12. Просмотров 620. Ответов 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).
Конец.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.02.2014, 02:12
Я подобрал для вас темы с готовыми решениями и ответами на вопрос Метод бисекции для нахождения корней в нелинейном уравнении (C++):

Нахождения корней уравнения: метод половинного деления (бисекции) или метод хорд - C++
Разработать программу нахождения корней уравнения f(x) =0 на интервале с точностью e = 0,001 (интервал или подобрать самостоятельно). При...

Мой код - метод бисекции, метод секущих (метод хорд) - C++
Всем привет!!! Изучаем в институте С++. Сделал код, и там, и там одна и та же проблема - при любых вбиваемых значениях программа делает...

Определить функцию, позволяющую распознавать наличие вещественных корней в квадратном уравнении - C++
Даны три квадратных уравнения ax^2+bx+с, bx^2+ax+с, cx^2+ax+b. Сколько из них имеют вещественные корни. (определить функцию, позволяющую...

Определить функцию, позволяющую распознавать наличие вещественных корней в квадратном уравнении - C++
Даны три квадратных уравнения ах^2 + bх + с = 0, bх^2 + ах + с = 0, сх^2 + ах + b = 0. Сколько из них имеют вещественные корни?...

Для нахождения корней квадратного уравнения использовать функцию - C++
Даны вещественные числа a,b,c,d,e,f.Переменноый S присвоить значение 1, если оба уравнения ах2 + bх + с = 0 и dx2 + ex + f = 0.для...

Написать перегруженные шаблоны функций для нахождения корней линейного (a*x + b = 0) уравнения - C++
Кто знает как решить задачу на С++, при помощи шаблонных функций? Написать перегруженные шаблоны функций для нахождения корней линейного ...

1
UnsKneD
алкокодер
155 / 151 / 12
Регистрация: 27.12.2012
Сообщений: 550
20.02.2014, 07:06 #2
FloraWinx,
У вас в цикле ни одно из значений fabs(f_midle)>eps не изменяется.
Обратите внимание на
Пока |F(xi)| > epsy повторять: и xi := xn + dx
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.02.2014, 07:06
Привет! Вот еще темы с ответами:

Метод бисекции - C++
Доброго времени суток! Очень нуждаюсь в вашей помощи. Поиск корней F(x) c помощью метода бисекции(деление отрезка пополам) + скользящее...

Метод бисекции - C++
Код: #include &lt;iostream.h&gt; #include &lt;float.h&gt; #include &lt;cmath&gt; #include &lt;float.h&gt; #include &lt;iomanip.h&gt; float f(float x) { ...

C++ метод бисекции - C++
Добрый день, помогите пожалуйста в решении задачи. fmin(x), f(x) = x^3 – 3 sinx,

Метод итераций для нахождения обратной матрицы - C++
Люди, подскажите как написать программу нахождения обратной матрицы методом итераций, или хотя бы решения системы уравнений


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

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru