Поиск:

Ответ в темуСоздание новой темы Создание опроса
> численные методы (Builder c++) Часть 2, модифицированный метод Ньютона 
:(
    Опции темы
leahov
Дата 4.12.2007, 18:10 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


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

Репутация: нет
Всего: нет



имеется код решения нелинейного уравнения методом Ньютона (он же метод касательных)

Код

//---------------------------------------------------------------------------

#include <vcl.h>
#pragma hdrstop
#include<math.h>

#include "urav.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
    : TForm(Owner)
{
}
//---------------------------------------------------------------------------
float f(float x)
{
//    return x*x*x + 3*x - 2;
 return x*x*x + x - 5;
}
float fdx(float x)
{//    return 3*x*x +3;
  return 3*x*x +1;
}
float fd2x(float x)
{
    return 4*x;
}
float absolute(float x)
{    return x < 0 ? -x : x;
}
void __fastcall TForm1::kasat(float a, float b, float eps)
{
 if(f(a)*f(b) < 0)
 {
  ShowMessage("Касательные: a и b неправильные");
  return;
 }
 float x0 = a;
 if(f(a)*fd2x(a) < 0)
 {
  x0 = b;
 }
 int step = 1;
 float x1 = x0 - f(x0)/fdx(x0);
 while(absolute(x1-x0) > eps)
 {
  x0 = x1;
  x1 = x0 - f(x0)/fdx(x0);
  step = step + 1;
 }
 Label17->Caption = x1;
 Label18->Caption = f(x1);
 Label25->Caption = step;
}

void __fastcall TForm1::Button1Click(TObject *Sender)
{
   float a, b, eps;
 //  ShowMessage("Vvedite a, b, epsilon: ");
   a = StrToFloat(Edit1->Text);
   b = StrToFloat(Edit2->Text);
   eps = StrToFloat(Edit3->Text);
   kasat(a, b, eps);

}
//---------------------------------------------------------------------------



как сделать из него модифицированный метод Ньютона, если кто знает подскажите

P.S. суть методов такова, метод ньютона (он же метод касательных) это Метод основан на замене функции f(x) на каждом шаге итерационного процесса поиска касательной, пересечение которой с осью абсцисс дает приближение корня. его формула Xn+1 = X n - (f (Xn / f ' (Xn), модифицированный метод ньютона это Этот метод заключается в том, что вместо вычисления производной на каждом шаге поиска находится ее приближенное значение с помощью конечных разностей и формула его Xn+1 = X n - (f (Xn / f ' (X0) )
Короче в формуле f` вместо (xn) становится (x0)
PM MAIL   Вверх
Ionitte
Дата 7.12.2022, 01:55 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

Репутация: нет
Всего: нет




Модератор: Сообщение скрыто.

PM MAIL   Вверх
scalosale
Дата 8.2.2023, 14:12 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



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

Репутация: нет
Всего: нет




Модератор: Сообщение скрыто.

PM MAIL   Вверх
  
Ответ в темуСоздание новой темы Создание опроса
Правила форума "С++ Builder"
Rrader

Запрещается!

1. Публиковать ссылки на вскрытые компоненты

2. Обсуждать взлом компонентов и делиться вскрытыми компонентами

  • Литературу по С++ Builder обсуждаем здесь
  • Действия модераторов можно обсудить здесь
  • С просьбами о написании курсовой, реферата и т.п. обращаться сюда
  • Настоятельно рекомендуем заглянуть в DRKB (Delphi Russian Knowledge Base) - крупнейший в рунете сборник материалов по Дельфи


Если Вам понравилась атмосфера форума, заходите к нам чаще! С уважением, Rrader.

 
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | C++ Builder | Следующая тема »


 




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


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

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