Метод бисекции

@noobchik · Регистрация: 15.06.2012

Студворк — интернет-сервис помощи студентам

Код:

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
#include <iostream.h>
#include <float.h>
#include <cmath>
#include <float.h>
#include <iomanip.h>
float f(float x)
{
    return x*x-4;
}
 
 
 
float fsgn(float t)
{
    if (t<0)
        t=-1.0;
    if (t>0)
        t=1.0;
    if (t=0)
        t=0.0;
 
    
    return t;
}
 
float q(float a,float b)
{
 
 
 
  float abseps = (FLT_EPSILON + 2.0 * FLT_EPSILON) * (1.0 + fmax(fabs(a), fabs(b)));
  float fa = f(a);
    float m;
  while( b - a > abseps ){
    float m = a + (b - a)/2.0; 
    float fm = f(m);
    if( fsgn(fa) * fsgn(fm) <= 0.0 )
      b = m;
    else{
      a = m;
      fa = fm;
    }
 
   
  }
  return m;
}
 
 
int main()
{
    cout << fixed;
    cout <<setprecision(10);
    cout << q(0, 10) << endl;
 
}

Выводит 0. В чем дело? Заранее спасибо.

@Catstail · 18.06.2012, 18:15

C++
1
2
3
4
5
6
7
8
9
10
11
12
float fsgn(float t)
{
    if (t<0)
        t=-1.0;
    if (t>0)
        t=1.0;
    if (t==0)   // Ошибка была здесь! 
        t=0.0;
 
    
    return t;
}

@noobchik · 18.06.2012, 18:31 **[ТС]**

Да я это условие, от отчаяния написал, если оно не больше и не меньше нуля то это и есть ноль.

Добавлено через 43 секунды
Проблема остается открытой.

@Catstail · 18.06.2012, 18:35

Но получается-то верный ответ... Хотя проверять число с плав. точкой на равенство - не вполне корректно.

Добавлено через 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
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
#include "iostream.h"
 
float f(float x)
{
    return x*x-4;
}
 
float fsgn(float t)
{
    if (t<0)
        t=-1.0;
    if (t>0)
        t=1.0;
    //if (t==0)
    //    t=0.0;
 
    
    return t;
}
 
float q(float a,float b)
{
  
  float abseps = 0.000001; //(FLT_EPSILON + 2.0 * FLT_EPSILON) * (1.0 + fmax(fabs(a), fabs(b)));
  float fa = f(a);
  float m,fm;
  
  while((b - a) > abseps )
 
  {
        //cout << m << endl;
        m = 0.5*(a+b); 
        fm = f(m);
        if( fsgn(fa) * fsgn(fm) <= 0.0 )
          b = m;
        else
        {
          a = m;
          fa = fm;
        }
   }
  return m;
}
 
void main()
{
    cout << q(0, 10) << endl;
 
}

Взгляни на картинку (ответ = 2):

@noobchik · 18.06.2012, 19:40 **[ТС]**

Спасибо.

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11680&d=1772460536 Одним из. . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .

@noobchik 2 / 2 / 0 Регистрация: 15.06.2012 Сообщений: 23
	18.06.2012, 18:31 [ТС]
	Да я это условие, от отчаяния написал, если оно не больше и не меньше нуля то это и есть ноль. Добавлено через 43 секунды Проблема остается открытой. 0

@noobchik 2 / 2 / 0 Регистрация: 15.06.2012 Сообщений: 23
	18.06.2012, 19:40 [ТС]
	Спасибо. 0