Модераторы: Poseidon
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> [c++]Метод Дихотомии, нелинейные уравнения 
:(
    Опции темы
Dukalis
  Дата 25.5.2008, 14:20 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

Репутация: нет
Всего: нет
Здравствуйте, кто поможет и скинет у кого может есть код (С, С++), решение нелинейных уравненией методом дихотомии (он же Метод деления отрезка пополам или метод бисекции)
.Спасибо
Извините , я не там тему создал, модераторы можете пожалуйста перекинуть тему в раздел Помощи, спасибо.

Это сообщение отредактировал(а) Dukalis - 25.5.2008, 14:22
PM MAIL   Вверх
maxdiver
Дата 25.5.2008, 21:15 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 381
Регистрация: 29.1.2008
Где: Саратов

Репутация: нет
Всего: 18
Ну а что тут особого?
Насколько я себе представляю метод деления пополам, предполагается, что функция f(x) - монотонна на отрезке [a;b], на котором и ищут её корень (ну и непрерывна, естественно).
Пусть для определённости функция f(x) возрастает, тогда имеем такой код:
Код
double l=a, r=b; // границы отрезка
const double EPS = 1E-9; // точность
while (r-l > EPS) {
   double m = (l + r) / 2;
   if (f (m) > 0)
      r = m;
   else
      l = m;
}
// результат (X корня) содержится в l (впрочем, и в r тоже, они же отличаются только на EPS :) )
PM MAIL WWW ICQ   Вверх
_shef_
Дата 26.5.2008, 16:21 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


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

Репутация: нет
Всего: 1
 Простите, ноя позволю себе уточнить код

Код

double l=a, r=b; // границы отрезка
const double EPS = 1E-9; // точность
while (r-l > EPS) {
   double m = (l + r) / 2;
   if (f (m)*f (r) > 0)
      r = m;
   else
      l = m;
}


Иначе возможно что решение не будет найдено.
PM MAIL WWW   Вверх
maxdiver
Дата 26.5.2008, 16:43 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 381
Регистрация: 29.1.2008
Где: Саратов

Репутация: нет
Всего: 18
_shef_
Почему?
Если функция, как я написал, возрастает, то можно просто f(m) > 0.
Вот если функция монотонна, но неизвестно, как именно - тогда да, надо f(m) f ®.
Разве я не прав?

Добавлено через 2 минуты и 14 секунд
Да, кстати. После цикла есть смысл проверить, действительно ли найденное l является корнем (просто посчитать f(l)). А то, если корня всё-таки не существует, то получится палево.
PM MAIL WWW ICQ   Вверх
_shef_
Дата 28.5.2008, 10:37 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


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

Репутация: нет
Всего: 1
maxdiver
Да Вы конечно правы. Просто у нас учат студентов решать в общем случае, когда не известно возрастает или убывает функция. Согласитесь что это более правильно, а то придется сначала искать производные. А по поводу отсутствия решения - у нас обычно в задании дают интервал на котором оно гарантированно существует. Иначе придется сначала изолировать корень.
PM MAIL WWW   Вверх
  
 Ответ в темуСоздание новой темы Создание опроса
Правила форума "Центр помощи"

ВНИМАНИЕ! Прежде чем создавать темы, или писать сообщения в данный раздел, ознакомьтесь, пожалуйста, с Правилами форума и конкретно этого раздела.
Несоблюдение правил может повлечь за собой самые строгие меры от закрытия/удаления темы до бана пользователя!


  • Название темы должно отражать её суть! (Не следует добавлять туда слова "помогите", "срочно" и т.п.)
  • При создании темы, первым делом в квадратных скобках укажите область, из которой исходит вопрос (язык, дисциплина, диплом). Пример: [C++].
  • В названии темы не нужно указывать происхождение задачи (например "школьная задача", "задача из учебника" и т.п.), не нужно указывать ее сложность ("простая задача", "легкий вопрос" и т.п.). Все это можно писать в тексте самой задачи.
  • Если Вы ошиблись при вводе названия темы, отправьте письмо любому из модераторов раздела (через личные сообщения или report).
  • Для подсветки кода пользуйтесь тегами [code][/code] (выделяйте код и нажимаете на кнопку "Код"). Не забывайте выбирать при этом соответствующий язык.
  • Помните: один топик - один вопрос!
  • В данном разделе запрещено поднимать темы, т.е. при отсутствии ответов на Ваш вопрос добавлять новые ответы к теме, тем самым поднимая тему на верх списка.
  • Если вы хотите, чтобы вашу проблему решили при помощи определенного алгоритма, то не забудьте описать его!
  • Если вопрос решён, то воспользуйтесь ссылкой "Пометить как решённый", которая находится под кнопками создания темы или специальным флажком при ответе.

Более подробно с правилами данного раздела Вы можете ознакомится в этой теме.

Если Вам помогли и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, Poseidon, Rodman
 
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Центр помощи | Следующая тема »

Подписаться на тему | Подписка на этот форум | Скачать/Распечатать тему


 




[ Время генерации скрипта: 0.0524 ]   [ Использовано запросов: 22 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
 По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.