12469

Розв’язання функціональних рівнянь з однією змінною

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

Математика и математический анализ

Лабораторна робота №3 Розвязання функціональних рівнянь з однією змінною Мета роботи: ознайомитися з методами розв‘язання рівнянь з однією змінною розглянути реалізацію цих методів у середовищі MatLab. Задачі лабораторної роботи: реалізувати один з методів у ві

Украинкский

2013-04-27

372.82 KB

48 чел.

Лабораторна робота №3

Розв’язання функціональних рівнянь з однією змінною

Мета роботи: ознайомитися з методами розв‘язання рівнянь з однією змінною, розглянути реалізацію цих методів у середовищі MatLab.

Задачі лабораторної роботи: реалізувати один з методів у відповідності до варіанту в середовищі MatLab, випробувати його роботу на прикладі розв‘язання конкретного рівняння.

Теоретичні відомості

Постановка задачі

Важливу задачу прикладного аналізу представляє розв’язок функціональних рівнянь з однією змінною (фізика, механіка, економіка та ін.)

У загальному вигляді функціональне рівняння записується у вигляді:

     (3.1)

де F(х) визначена і неперервна на скінченому чи нескінченному інтервалі [a,b].

Число називається коренем r кратності, якщо воно задовольняє рівнянням:

  (3.2)

Однократний корінь називається простим.

Два рівняння F(х) та G(х) називаються рівносильними, якщо будь-який розв’язок кожного з них є розв’язком і для другого.

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

Шляхом алгебраїчних перетворень із будь-якого алгебраїчного рівняння можна отримати рівняння в канонічній формі:

   (3.3)

де n – степінь алгебраїчного рівняння .

Кожне алгебраїчне рівняння має хоча б один дійсний корінь чи пару комплексних.

При зведенні алгебраїчного рівняння до канонічної форми отримаємо ті ж корені, що і для початкового рівняння, але при цьому можуть з’явитися і сторонні корені. Наприклад:

 (3.4)

Якщо функція F(х) не є алгебраїчною, то рівняння (3.1) називається трансцендентним. Наприклад:

    (3.5)

В деяких випадках розв’язування трансцендентних рівнянь зводиться до розв’язання алгебраїчних.

Більшість функціональних рівнянь з однією змінною не розв’язуються точними методами шляхом аналітичних перетворень. На практиці для їх розв’язання застосовують чисельні методи. Розв’язати таке рівняння – значить встановити, чи має воно корені на вказаному проміжку, кількість коренів, і відшукати значення коренів із заданою точністю (-окіл).

Задача чисельного знаходження дійсних і комплексних коренів рівняння зазвичай складається з двох етапів:

  1.  відокремлення коренів – знаходження достатньо малих околів у заданій області, в яких знаходиться один корінь.
  2.  уточнення коренів – обчислення коренів із заданим ступенем точності в деякому околі.

Для знаходження дійсних коренів рівняння (3.1) застосовуються наступні чисельні методи:

  1.  метод половинного ділення;
  2.  метод хорд;
  3.  метод дотичних (Ньютона);
  4.  комбінований метод;
  5.  метод Рибакова, знаходження всіх дійсних коренів;
  6.  метод простої ітерації.

Відокремлення коренів

Відокремлення коренів – встановлення „тісних” проміжків, що вміщують кожен проміжок – один корінь.

В багатьох випадках цей етап проводять графічно. Враховуючи, що дійсні корені рівняння (3.1) – це точки перетину графіка функції F(х) з віссю абсцис, досить побудувати графік F(х) і відмітити на осі OX відрізки, кожен з яких вміщує по одному кореню.

Рис. 3.1. Графічний спосіб відокремлення коренів.

Побудову графіка вдається сильно спростити, якщо замінити рівняння (3.1) рівносильним йому рівнянням F1(х)= F2(х). В цьому випадку будують графіки функцій F1(х) та F2(х), потім на вісі OX відмічають відрізки, що вміщують абсциси точок перетину цих графіків.

Рис. 3.2. Варіанти розв’язання рівняння при заміні функції рівносильними функціями F1(х) та F2(х) на проміжку [ 1; 3 ].

При написанні програми, яка виконує відокремлення коренів графічний спосіб не є зручним, тому застосовують найпростіший алгоритм відокремлення коренів: проміжок, на якому відшукуються розв’язки, ділять рівномірно на невеликі відрізки, на кінцях кожного відрізку розраховують значення функції, якщо функція на кінцях відрізку змінює знак – відрізок містить корінь.

Метод половинного ділення

Припустимо, що рівняння має на відрізку [a, b] єдиний корінь. Функція F(x) на цьому відрізку неперервна. Поділимо відрізок [a, b] точкою навпіл (рис. 3.3). Якщо , то є коренем рівняння. Якщо F(c)  0, то вибираємо відрізок [a, с] чи [с, b], на якому функція змінює знак, тобто F(а) F(с) < 0 чи F(с) F(b) < 0. Новий зменшений відрізок знову поділимо навпіл і проводимо повторно ті ж самі дії. В результаті здобудемо на деякому етапі або точний корінь або послідовність вкладених один в одного відрізків , , … , таких, що  .

Рис. 3.3. Знаходження кореня рівняння методом половинного ділення.

Процес поділу продовжимо до тих пір, поки довжина відрізка не стане менше заданої точності :

|a – b| <      (3.6)

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

Метод половинного поділу практично зручно застосовувати для грубого знаходження кореня даного рівняння, так як при збільшені точності значно збільшується об’єм обчислювальної роботи.

Метод хорд

Даний метод ґрунтується на лінійній інтерполяції функції F(x)=0 по двом значенням, що мають протилежні по знаку значення функції F(a) та F(b). Метод хорд швидше збігається до розв‘язку навіть при досить малих значеннях . Потрібно знайти корінь рівняння F(x)=0 на проміжку [a; b], і якщо відомо, що F(x) неперервна на [a; b] і F(a)∙F(b) < 0. Крім того,  перша F'(x) і  друга F''(x) похідні функції F(x) зберігають на проміжку [a; b] свій знак. Замінимо функцію F(x) лінійною функцією, яка проходить через вузлові точки (a , F(a)) і (b , F(b)) :

     (3.7)

Лінійна функція P(x) на кінцях відрізку [a; b] приймає такі ж значення, як і функція F(x)=0.

Рис. 3.4. Знаходження кореня за методом хорд.

В якості першого наближення при знаходженні кореня функції F(x)=0 візьмемо точне значення кореня функції P(x)=0, тобто х1, яке розрахуємо з рівняння:

                (3.8)

При подальшому дослідженні відрізків [a; х1] та 1; b], виберемо той, на якому функція змінює знак. З рис. 3.4 бачимо, що таким відрізком є [a; х1] . Для вибраного відрізка побудуємо лінійне наближення функції за формулою (3.7), виконаємо розрахунки кореня для лінійного наближення за формулою (3.8). Розрахунки припинимо, коли .

Метод січних

Метод січних подібний до методу хорд, тільки точки (х0 , F(х0)) і (х1 , F(х1)) взяті з одного боку від кореня рівняння F(x)=0. Геометрична інтерпретація методу представлена на рис.3.5.

В якості початкового наближення обираємо точки (х0 , f(x0)) та (х1, f(x1)).

Рис. 3.5. Знаходження кореня за методом січних.

Через точки (х0 , f(x0)) та (х1 , f(x1)) проводимо січну до графіку функції, яка перетинає вісь ох в точці 2, 0). Перевіряємо виконання умови , якщо вона не виконується, проводимо січну через точки (х1 , f(x1)) та (х2 , f(x2)), знаходимо точку перетину січної з віссю ох (точка 3, 0)) і перевіряємо виконання чергової умови , і так до виконання умови виходу з ітераційного процесу.

Для математичного опису методу січних отримаємо формулу прямої, що проходить через дві точки (х0 , f(x0)) та (х1 , f(x1)) :

           (3.9)

Враховуючи, що f(xn+1) = 0 отримаємо загальну формулу для методу січних:

               (3.10)

Метод січних має досить високу збіжність до розв‘язку у випадку, коли F(x)=0 - гладка функція, але в процесі розв’язання деяких рівнянь швидкість збіжності може знижатись (наприклад, на ділянках функції, близьких до функції х = const).

Метод дотичних (Ньютона)

Метод дотичних базується на заміні функції F(x)=0 у точці початкового наближення х0 дотичною, яка при перетині з віссю ох в точці 1, 0) дає перше наближення. У цьому методі Ньютон замість інтерполяції використав екстраполяцію, що знаходиться за допомогою дотичної у визначеній точці. Геометрична інтерпретація методу дотичних (Ньютона) показана на рис. 3.6. В якості початкового наближення обираємо абсцису х0 = b. У точці (х0 f(x0)) проводимо дотичну до графіку функції, яка перетинає вісь ох в точці 1, 0).

Рис. 3.6. Знаходження кореня за методом дотичних.

Перевіряємо виконання умови , якщо вона не виконується, проводимо дотичну в точці (х1 , f(x1)), знаходимо точку перетину дотичної з віссю ох (точка 2, 0)) і перевіряємо виконання чергової умови , і так до виконання умови виходу з ітераційного процесу.

Для математичного опису методу дотичних отримаємо формулу дотичної, що проведена в точці (х0 , f(x0)) і має кутовий коефіцієнт f'(x0):

            (3.11)

У точці перетину цієї дотичної з віссю ох (точка 1, 0)) y = 0, тому формулу 6 можемо записати у вигляді:

             (3.12)

Для загального випадку формулу Ньютона можемо записати у вигляді:

          (3.13)

Метод дотичних (Ньютона) має високу збіжність до розв‘язку, але час виконання ітерації дещо збільшується за рахунок необхідності обчислення похідної f'(xk).

Метод простої ітерації

Замінимо рівняння F(x)=0 рівносильним рівнянням:

x=f(x)       (3.14)

Припустимо w – корінь рівняння, а x0 – отримане будь-яким способом нульове наближення до кореня w. Підставимо x0 в праву частину рівняння (3.14), отримаємо x1=f(x0); x2=f(x1); … ; xn=f(xn-1). Цю числову послідовність називають послідовністю наближень чи ітераційною послідовністю. Послідовність наближень може бути збіжною чи розбіжною.

 Теорема: ітераційна послідовність буде збіжною при будь-якому початковому значенні x0[a; b], якщо для рівняння x=f(x), що має єдиний розв’язок на проміжку [a; b], виконуються умови:

  1.  f(x) визначена і диференційована на [a; b];
  2.  f(x) визначена для всіх x[a; b];
  3.  існує таке дійсне q, що для всіх x[a; b].

Швидкість збіжності характеризується нерівністю (умова Ліпшиця):

    (3.15)

Умовою припинення ітерацій є .

При розв’язанні поставленого в лабораторній роботі завдання потрібно ознайомитися з особливостями програмування на внутрішній мові середовища MatLab.

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

Створити програму на внутрішній мові середовища МatLAB, що реалізує метод, заданий по варіанту (таблиця 3.1), провести тестування створеної програми на прикладі.

Варіанти завдань до лабораторної роботи №3.

Таблиця 3.1

Варіант

Метод

Тестовий приклад

1

половинного ділення

2

січних

3

хорд

4

дотичних

5

простої ітерації

6

половинного ділення

7

січних

8

хорд

9

дотичних

10

простої ітерації

Створені програми перевірити на прикладі алгебраїчних рівнянь, що вводяться користувачем з клавіатури.

Додаток А

Лабораторна робота № _.

Тема

Мета роботи: ознайомитися з методами розв‘язання рівнянь з однією змінною, розглянути реалізацію цих методів у середовищі MatLab.

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

<У розділі наводяться ті теоретичні відомості, що використовуються при розв‘язанні задачі по варіанту лабораторної роботи>*

2. Хід виконання лабораторної роботи.

<У розділі пояснюються особливості алгоритму реалізації задачі у середовищі MATLAB, детально розписуються етапи розв‘язання поставленого завдання>

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

<У розділі наводиться лістінг реалізації задачі у середовищі MATLAB, в якому повинні бути основні коментарі по етапам виконання поставленого завдання>

4. Результати виконання програми.

< У розділі наводиться результат реалізації задачі у середовищі MATLAB на тестовому прикладі у відповідності з варіантом (у числовому та графічному вігляді.>

5. Висновок.

<У розділі робиться висновок, у якому проводиться аналіз результатів і оцінка точності та сфери застосування методу.>

_____________________

* Примітки: при захисті лабораторної роботи студент повинен знати суть всіх методів, що вивчаються в рамках лабораторної роботи.


 

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

22504. ЧТО ТАКОЕ СТРАТЕГИЧЕСКОЕ УПРАВЛЕНИЕ 343.32 KB
  Сущность стратегического управления Было бы наверно считать что еще 2030 лет назад фирмы совершенно не имели стратегического поведения. Разработка идей стратегического управления нашла отражение в работах таких авторов как Fannhofs и Granger 1971; Schendel и Hatten 1972; Irwin 1974; и др. В качестве ведущей идеи отражающей сущность перехода к стратегическому управлению от оперативного управления явилась идея необходимости переноса центра внимания высшего руководства на окружение для того чтобы соответствующим образом и своевременно...
22505. АНАЛИЗ СРЕДЫ 354.03 KB
  Внутренняя среда организации является источником ее жизненной силы. Она заключает в себе тот потенциал который дает возможность организации функционировать а следовательно существовать и выживать в определенном промежутке времени. Но внутренняя среда может также быть и источником проблем и даже гибели организации в том случае если она не обеспечивает необходимого функционирования организации. И на них претендуют многие другие организации находящиеся в этой же среде.
22506. ОПРЕДЕЛЕНИЕ МИССИИ И ЦЕЛЕЙ ОРГАНИЗАЦИИ 219.47 KB
  Но целевое начало в деятельности организации возникает отнюдь не только потому что ей нужно иметь ориентиры чтобы не погибнуть в изменяющемся окружении. В первую очередь целевое начало в деятельности организации возникает потому что организация это объединение людей преследующих определенные цели. Люди создают организации для того чтобы с их помощью решать свои проблемы.
22507. ЧТО ТАКОЕ СТРАТЕГИЯ ФИРМЫ И КАК ОНА ВЫРАБАТЫВАЕТСЯ 494.62 KB
  Сущность стратегии организации . Выбор стратегии и ее реализация составляют основную часть содержания стратегического управления. Существуют два противоположных взгляда на понимание стратегии. Первое понимание стратегии базируется на следующем процессе.
22508. КОНЦЕПЦИЯ ПРОДУКТА В СТРАТЕГИЧЕСКОМ УПРАВЛЕНИИ 466.76 KB
  Понятие продукта Понятие продукта очень сильно варьируется в зависимости от того кто вкладывает смысл в это понятие в контексте каких процессов оно рассматривается относительно к каким видам продукции применяется и т. Каждое данное понимание продукта отвечает на вопрос о том что делает фирма и принципиально отличается одно от другого в зависимости от того кто и в каком контексте дает ответ. Свой определенный смысл в понимание продукта вкладывает разработчик продукта посвоему смотрит на продукт плановик и наконец свой...
22509. ДИНАМИКА ПРОДУКТА 311.55 KB
  От того как фирма смотрит на производимый ею продукт как она подходит к выработке своей стратегии продукта во многом зависит сумеет ли фирма найти свое место на рынке сможет ли обеспечить себе устойчивые связи с окружением осуществить своевременные и адекватные запросам среды изменения. Но это было статическое рассмотрение продукта. Другим очень важным направлением рассмотрения продукта в стратегическом управлении является его динамическое рассмотрение в соответствии с которым продукт предстает как изменяющееся во времени...
22510. СТРАТЕГИЯ ПРОДУКТА 404.19 KB
  Ответ на эти вопросы относящиеся к стратегии продукта формируется под влиянием нескольких факторов. В первую очередь стратегия продукта зависит от стратегии фирмы. Однако хотя стратегия продукта и является подстратегией общей стратегии фирмы и соответственно определяется характером ее целей и ее стратегией тем не менее можно указать на отдельные факторы оказывающие непосредственное влияние на выработку стратегии продукта.
22511. КАК СТРОИТСЯ ВЗАИМОДЕЙСТВИЕ ЧЕЛОВЕКА И ОРГАНИЗАЦИИ 368.1 KB
  Для того чтобы понять как строится взаимодействие человека с организацией необходимо не только уяснить в чем суть проблемы этого взаимодействия а также то что в личности человека определяет его поведение в организации и какие характеристики организационного окружения влияют на то как происходит включение человека в деятельность организации. Подходы к построению взаимодействия человека и организационного окружения Взгляд на поведение человека в организации может быть осуществлен с двух позиций: с позиции человека взаимодействующего с...
22512. ВХОЖДЕНИЕ ЧЕЛОВЕКА В ОРГАНИЗАЦИЮ 134.02 KB
  Быть членом организации совсем не одно и то же что входить в организацию становиться ее членом. Втретьих это проблема изменений и модификаций в организации которые происходят даже тогда когда организация уже имеет свободное место для человека и сама принимает человека на это место в соответствии с ее потребностями и критериями отбора. От решения данных проблем зависит не только то сможет ли человек войти в организацию но и то как человек будет функционировать в организации как будет строиться его взаимодействие с организационным...