3854

Метод хорд і дотичних

Лабораторная работа

Информатика, кибернетика и программирование

Метод хорд і дотичних Мета роботи: вивчити і засвоїти ітераційні методи розв’язування алгебраїчних і трансцендентних рівнянь. Короткі теоретичні відомості. Метод хорд. Метод хорд – метод лінійної інтерполяції (метод пропорційних частин, ме...

Украинкский

2012-11-09

135 KB

50 чел.

Метод хорд і дотичних

Мета роботи: вивчити і засвоїти ітераційні методи розв’язування алгебраїчних і трансцендентних рівнянь.

Короткі теоретичні відомості.

Метод хорд.

Метод хорд – метод лінійної інтерполяції (метод пропорційних частин, метод січних). Ідея методу полягає в тому, що на досить малому проміжку  дугу кривої  заміняють хордою, а за наближене значення кореня приймають точку перетину хорди з віссю Ох.

Значення функції визначаються в точках, які розташовані на осі Ох через рівні інтервали. Це робиться доти, поки кінці інтервалів ,  не будуть мати різні знаки. Тоді пряма, що проведена через ці дві точки, перетинає вісь абсцис у точці

                                               (1 )

Формула (1) відображає метод хорд, який є ітераційним методом першого порядку. Для досить малих проміжків  можна використати оцінку похибки

                                                 (2)

Рис. 1. Алгоритм методу половинного ділення.

Після цього знаходять значення f () і порівнюють його з f (). Надалі користуються  замість того значення, з яким воно збіглося за знаком. Якщо   ≤ ε , то вся процедура повторюється спочатку (рисунок 2). Алгоритм методу хорд подібний до попереднього, за винятком процедури оцінки .

Треба також враховувати, що в алгоритмі обчислень за цим методом контроль похибки проводиться за тим кінцем інтервалу, що рухається. В випадку, що показаний на рисунку1,  аналізуються послідовні наближення: на першому кроці x1 x2 ≤ ε , на другому – x1 x3 ≤ ε , на третьому – x3 x4 ≤ ε і т. д.

Похибка розв’язку оцінюється за формулою:

                                                         (3),

де M1 ,m1 – відповідно, найбільше та найменше значення модуля першої похідної на відрізку.

Рис. 2. Метод хорд.

Метод дотичних.

У методі дотичних здійснюється екстраполяція за допомогою дотичної до кривої в даній точці (рисунок 3):

В основі цього методу лежить розвинення функції в ряд Тейлора

Доданки, що містять h у другому і вищих степенях, відкидаються, внаслідок чого одержується наведена вище наближена формула для оцінки .

Швидкість збіжності цього алгоритму значною мірою залежить від вірного вибору початкової точки. Коли в процесі обчислень кут нахилу дотичної f′(x) перетворюється на нуль, застосування цього методу ускладнюється. Можна також показати, що у випадку дуже великих значень f′′(x) (опуклість функції) чи кратних коренів метод Ньютона стає неефективним.

Початкове наближення слід вибирати з умови

Похибка методу оцінюється як:

де M2 – найбільше за модулем значення другої похідної на інтервалі [ , ] .

Рис. 3. Метод дотичних.

Завдання лабораторної роботи

Знайти один із коренів рівняння f(x) = 0, використовуючи методи дотичних, хорд і простої ітерації з точністю . Вивести кількість ітерацій, необхідних для досягнення заданої точності для кожного з методів.

f(x)=x2-cos(5x)

Лістинг програми

#include <conio.h>

#include <iostream.h>

#include <math.h>

double f(double x)

{

return (x*x-cos(5*x));

}

double f1(double x)

{

return (2*x+5*sin(5*x));

}

double f2(double x)

{

return (2+25*cos(5*x));

}

int main()

{

  double a,b,c,exp=0.00001;

  cout<<"Vvedit znachennia"<<endl;

  cout<<"a="; cin>>a;

  cout<<"b="; cin>>b;

  {

do

  {

  c=a-f(a)/(f(b)-f(a))*(a-b);

  if (f(c)*f(a)>0) a=c;

  else b=c;

  }

while (fabs(f(c))>=exp);

  cout<<"Metod xord->"<<c<<"\n";

   }

   {

 if (f(a)*f2(a)>0) c=a;

  else c=b;

do

    {

      c=c-(f(c)/f1(c));

    }

while (fabs(f(c))>=exp);

   cout<<"Metod dotuchnux->"<<c<<"\n";

   }

getch();

return 0;

}

Вигляд програми

Висновок: на дані лабораторній роботі я вивчив два методи знаходження розвязку функції з значеннями на проміжку – це метод Хорд і метод Дотичних(метод Ньютона).


 

А также другие работы, которые могут Вас заинтересовать

39957. Газодинамика как раздел механики сплошных сред 907.5 KB
  Краткий очерк развития механики жидкости и газа. Математический аппарат используемый в механике жидкости и газа [1. Газодинамика как раздел механики сплошных сред Многие машины и аппараты созданные к настоящему времени характеризуются перемещением газа или жидкости внутри их или перемещением самого аппарата в среде газа или жидкости. Целью курса Газодинамика является изучение явлений протекающих в газе и жидкости и закономерностей которым эти явления подчиняются.
39958. УРАВНЕНИЯ ГАЗОВОЙ ДИНАМИКИ ДЛЯ ЕДИНИЧНОЙ СТРУЙКИ 401.5 KB
  Предельная скорость движения газа. Уравнение неразрывности Выведем основные уравнения газовой динамики для элементарной струйки газа поперечные размеры которой настолько малы что в каждом ее сечении можно считать постоянными все основные параметры потока: скорость давление температуру и плотность газа. Чтобы получить уравнение неразрывности рассмотрим стационарное установившееся движение элементарной струйки газа рис. Элементарная струйка Рассмотрим некоторый участок струйки между двумя нормальными к поверхности тока сечениями 1 и...
39959. Элементы гидродинамики 441 KB
  Cилы действующие в жидкости 3.1 Элементарный параллелепипед в потоке жидкости Грани бесконечно малой частицы жидкости имеющей в начале движения форму прямого параллелепипеда с ребрами dx dy dz с течением времени могут скашиваться и растягиваться рис.8 представляет собой уравнение неразрывности жидкости.9 Здесь под плотностью жидкости понимается предел отношения массы частицы к ее объему 3.
39960. ГАЗОДИНАМИЧЕСКИЕ ФУНКЦИИ 81 KB
  ГАЗОДИНАМИЧЕСКИЕ ФУНКЦИИ План лекции. Зависимость параметров потока в функции числа M. Зависимость параметров потока в функции скоростного коэффициента. Зависимость параметров потока в функции числа M.
39961. ДЕТАЛИ МАШИН И ОСНОВЫ КОНСТРУИРОВАНИЯ 10.06 MB
  1 а е: Ft Н окружная сила на барабане ленточного или на звездочке цепного конвейера; V м с скорость движения ленты или цепи; Dб мм диаметр барабана; Zзв число зубьев тяговой звездочки; Рзв мм шаг тяговой цепи.2 Вид передачи Твердость зубьев Передаточное число Uрек Uпред Зубчатая цилиндрическая: тихоходная ступень во всех редукторах uт 350 НВ 40. Термообработка зубчатых колес редуктора улучшение твердость зубьев 350НВ. Первая группа колеса с твердостью поверхностей зубьев Н  350 НВ Применяются в слабо и...
39962. Специализированный вычислитель (СВ) 194 KB
  При обращении ВчУ в режиме Чтение к ОЗУ по адресу 034320 обращение происходит в ячейке ДЗУ с адресом 134320. Специализированный вычислитель СВ относится к классу специализированных ЭВМ и предназначен для решения специфических задач обработки информации: 1. Отображение информации на рабочих местах РМ лиц боевого расчета; 3. Вычислительное устройство ВчУ является основным операционным устройством СВ предназначенным для обработки цифровой и логической информации реагирования на сигналы прерывания внешних устройстви управления...
39963. Методы локализации неисправностей в аппаратуре СВ и РМ 47 KB
  Наиболее склонными к поломке элементами являются транзисторы. Основные же мероприятия по устранению неисправности на принципиальном уровне сводятся к выпаиванию неисправного элемента и впаиванию на его место нового в случае необходимости замены элемента резисторы транзисторы диоды и другие. На принципиальном уровне неисправными элементами могут быть транзисторы на платах: ВУ2: Т1 Т2 Т3 либо Т4. Более полная информация о неисправных транзисторах находится в перечне элементов схемы.
39964. Отчет по учебной геологической практике 69 KB
  Целью проведения полевой практики по инженерной геологии является закрепление теоретического материала и ознакомление с природными условиями залегания различных типов горных пород а также с формами проявления геологических и инженерногеологических процессов. Ее учебными задачами являются: Приобретение навыка визуального определения геологических особенностей горных пород. В течении практики в полевых условиях изучаются: Вещественный состав и строение пород. Условия формы залегания пород.
39965. Учебная геологическая практика 865 KB
  4 Порядок проведения практики. Оценка практики. Цели и задачи практики Учебная геологическая практика проводится в летнее время после изучения студентами курса Инженерная геология.