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

Метод половиного деления. - C++

Восстановить пароль Регистрация
 
Luntic
0 / 0 / 0
Регистрация: 11.05.2010
Сообщений: 7
05.06.2010, 00:32     Метод половиного деления. #1
Задача: найти нуль функции использовав метод половинного деление. Вооопщем с заданием справился. Не понятно вот это:-надо проверить имеет ли корень(т.е. нуль функции на заданом промежутке), и если неимеет, то вывести "корней нет." Пока не могу придумать(ну кроме производной), как сделать до цикла.

Добавлено через 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
#include "stdafx.h"
 
 
 
int _tmain(int argc, _TCHAR* argv[])
{
 
    double a,b,c,eps;
    cout<<"input a=";cin>>a;
    cout<<"input b=";cin>>b;
    cout<<"input eps=";cin>>eps;
 
c=poldel(a,b,eps);
cout<<"c="<<c<<endl;
 
 
 
 
_getch();
    return 0;
}
 
double fun(double x)
 
{
return x*x-4;
    }
 
double poldel(double a,double b,double eps)
        {
double c;
//bool k=0;
c=(a+b)/2;
/*if ((fun(a)>=0) && (fun(c)>=0))   
{
    cout<<"net kornej"<<endl;
    return 0;
}*/
else 
do
 
            {
     c=(a+b)/2; 
      if (fun(c)*fun(a)<=0) 
      {
          b=c;
          //k=1;
      }
       else a=c;
            }
    while (fabs(a-b)>=eps);
    /*if ( k==0)
    {
        cout<<"net kornej"<<endl;
        return 0;
    }*/
    return c;
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
лендер
46 / 46 / 2
Регистрация: 12.01.2010
Сообщений: 183
05.06.2010, 01:11     Метод половиного деления. #2
можна сделать так:

Цитата Сообщение от Luntic Посмотреть сообщение
double poldel(double a,double b,double eps)
{
double c;
int k = 0;
c=(a+b)/2;

do

{
k++;
c=(a+b)/2;
if (fun(c)*fun(a)<=0)
{
b=c;
}
else a=c;
}
while (fabs(fun(c))>=eps && k < 100 );

if( k == 100 )
cout << "корней нет";

return c;
}[/CPP]
mexx-nikitin
78 / 78 / 3
Регистрация: 28.03.2010
Сообщений: 302
05.06.2010, 12:18     Метод половиного деления. #3
У меня была подобная РГР, может поможет
метод.срр
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
#include "stdafx.h"
#include "conio.h"
#include "locale.h"
#include "iostream"
#include "math.h"
#include "function.h"
using namespace std;
 
int _tmain(int argc,_TCHAR* argv[])
{
    setlocale (LC_CTYPE,"Russian_Rusia.1251");
        double a=0;
        double b=2.0;
        double e=0.1;
        double x;
        for(;;)
        {
            double x1=(a+b)/2-e/2;
            double x2=(a+b)/2+e/2;
        if(f(x2)>f(x1))
            b=x2;
        else
            a=x1;
        double l=b-a;
        if(l<0.2)
            break;
    }
    x=(a+b)/2; 
    double min=f(x); 
    cout<<"\t"<<"Значение x при котором функция имеет минимум "<<x<<endl;
    cout<<"\t"<<"Значение функции в этой точке Y(x)="<<min<<endl;
    _getch();
    return 0;
}
function.cpp
C++
1
2
3
4
5
6
7
8
9
10
#include "stdafx.h"
#include "cmath"
 
double f (double x)
{   
    double y;
    y=x*x*x*x-4/3*x*x*x+20/3*x*x-112/27*x+10;
 
    return(y);
}
functin.h
C++
1
double f (double x);
Luntic
0 / 0 / 0
Регистрация: 11.05.2010
Сообщений: 7
05.06.2010, 15:49  [ТС]     Метод половиного деления. #4
спасибо=)) но всё таки это не то что хочет препод(а что он хочет я не могу понять, а он объяснить=) )

Добавлено через 51 минуту
оказывается надо найти корень ур-е первой степени. Жалко что уточнили сейчас, но не в этом суть. Тогда просто, вот:

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
#include "stdafx.h"
 
 
 
int _tmain(int argc, _TCHAR* argv[])
{
 
    double a,b,c,eps;
    cout<<"input a=";cin>>a;
    cout<<"input b=";cin>>b;
    cout<<"input eps=";cin>>eps;
 
c=poldel(a,b,eps);
cout<<"c="<<c<<endl;
 
 
 
 
_getch();
    return 0;
}
 
double fun(double x)
 
{
return x*3-12;
    }
 
double poldel(double a,double b,double eps)
        {
double c;
c=(a+b)/2;
if ((fun(a)*fun(b)>0))   
{
    cout<<"net kornej"<<endl;
    return 0;
}
else 
do
 
            {
     c=(a+b)/2; 
      if (fun(c)*fun(a)<=0) 
      {
          b=c;
          //k=1;
      }
       else a=c;
            }
    while (fabs(a-b)>=eps);
    return c;
}
Добавлено через 2 минуты
а это с проверкой в цикле :
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
#include "stdafx.h"
 
 
 
int _tmain(int argc, _TCHAR* argv[])
{
 
    double a,b,c,eps;
    cout<<"input a=";cin>>a;
    cout<<"input b=";cin>>b;
    cout<<"input eps=";cin>>eps;
 
c=poldel(a,b,eps);
cout<<"c="<<c<<endl;
 
 
 
 
_getch();
    return 0;
}
 
double fun(double x)
 
{
return x*3-12;
    }
 
double poldel(double a,double b,double eps)
        {
double c;
bool k=0;
c=(a+b)/2;
do
 
            {
     c=(a+b)/2; 
      if (fun(c)*fun(a)<=0) 
      {
          b=c;
          k=1;
      }
       else a=c;
            }
    while (fabs(a-b)>=eps);
    if ( k==0)
    {
        cout<<"net kornej"<<endl;
        return 0;
    }
    return c;
}
Yandex
Объявления
05.06.2010, 15:49     Метод половиного деления.
Ответ Создать тему
Опции темы

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