12979

Математичне моделювання та диференціальні рівняння

Лекция

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

Лекція 1 Математичне моделювання та диференціальні рівняння. 1.1. Поняття математичного моделювання. Поняття математичного моделювання трактується різними авторами по своєму. Ми будемо його повязувати з нашою спеціалізацією прикладна математика. Під ма

Украинкский

2013-05-07

300.5 KB

0 чел.

Лекція 1

Математичне моделювання та диференціальні рівняння.

1.1. Поняття математичного моделювання.

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

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

  а) постановка задачі;

  б) побудова математичної моделі;

  в) перевірка її адекватності;   

  г) узагальення та теоретичне дослідження данного класу задач;

  д) створення програмного забезпечення;

  е) проведення обчислювального експеременту;

  ж) впровадження цих результатів в виробнитство.

  Розглянемо питання використання диіеренціальних рівнянь в деяких предметних областях.

1.2. Диференціальні рівняння в екології.

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

  Опишемо математично процес розмноження чи вмирання популяцій.

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

    (1.1)

В (1.1) і  можуть залежити від . Наприклад:

    (1.2)

Де  – коефіцієнт народжуваності,  – смертності. Маємо з (1.2)

    (1.3)

Розвязок диференціального рівняння запишемо в вигляді

З розв’язку (1.4) видно, що при  популяція вижчваюча, а при  – вмираюча.

    (1.4)

  Рівняння (1.3) в деяких випадках береться нелінійне

    (1.5)

  Це рівняння Беруллі при  і його розв’язок запишеться в такому вигляді

    (1.6)

З формули (1.6) видно, що при . При цьому можливі випадки

, та

Рівняння (1.5) описує.

  Можна говорити і про більш складні рівняння, системи рівнянь.

  Розглянемо більш детально двух видову модель «хижак-жертва», яка була побудована для виявлення коливань рибних уловів в Адріатичному морі.

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

    (1.7)

де додатні константи.

  В (1.7) доданок виражає залежність прирісту великих риб від числа малих,  – зменьшення числа малих риб від великих.

1.3. Закони Кеплера руху планет.

  Згідно закону всесвітнього тяжіння два тіла, які знаходятся на віддалі  друг від друга і які мають маси  і  притягаються з силою

    (1.8)

де - константа тяжіння.

  Опишемо рух планети з масою навколо Сонця маси . Вплив других планет на них не будемо враховувати. (Мал 1.1).

Сонце знаходиться в початку координат, а планета має положення  в момент часу . Використавши другий закон Ньютона маємо:

    (1.9)

Враховуючи, що

Позначимо , прийдемо до системи

    (1.10)

Без обмеження загальності візьмемо початкові умови:

при     (1.11)

Перейдемо до полярних координат:

Позначивши отримані вирази в (1.10) будемо мати

Помножимо перше рівняння на ,друге на  і складемо:

    (1.12)

Домножимо перше рівняння на ,друге на  і складемо:

    (1.13)

Перепишемо в нових змінних умови (1.11):

Рівняння (1.13) перепишемо у вигляді

    (1.14)

    (1.15)

Звідки маємо

    (1.6)

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

Звідки

    (1.17)

,або

Останній вираз означає секторну швидкість. З (1.16) випливає, що вона являється постійною. Це означає, що радіус-вектор “замітає” за рівні проміжки часу рівні площі.

  1-ій закон Кеплера: кожна із планет рухається по плоскій кривій відносно Сонця так, що радіус-вектор, який зв’язує Сонце і кожну з планет, “замітає” рівні площі за рівні проміжки часу.

  Задачу Кощі (1.12)-(1.14) можна розвязати. Розвзок має еліпсоідальну форму, на основі цього робиться наступний висновок:

  2-ій закон Кеплера: траєкторії планет рухаються по еліпсам, в одному з фокусів яких знаходиться Сонце.

  З аналізу траєкторій випливає таке твердження:

  3-ій закон Кеплера: квадрати періодів обертання планет пропорційні кубам великих осей їх орбіт.

1.4. Диференціальні рівняння закону пропиту і пропозиції в економічних дослідженнях.

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

    Нехай  – ціна, наприклад, на фрукти, – тенденція формування ціни. Тоді, як попит так і пропозиція будуть функціями введених величин. Як показує практика, ці функції можуть бути різними. Часто попит  і пропозиція задаються лінійними

    (1.17)

залежностями. Наприклад:

Для того, щоб попит відповідав пропозиції необхідно:

Звідки

    (1.8)

Припустимо, що в момент 1кг фруктів коштував 1крб. Тоді , , отже

   (1.19)

Це закон зміни цін, щоб між попитом і пропозицією була рівновага.

1.5 Найпростіші рівняння руху частинок в електромагнитних поясах.

  Швидкість зміни імпульсу частинки

дорівнює силі Лоренса, яка діє на неї

    (1.20)

де зарядове число, заряд частинки,  – вектор напруженності прискорюючого поля,  – вектор магнітної індукції,  – вектор швидкості частинки.

де  маса спокою, -приведена енергія частинки.

- векторний добуток двох змінних.

З (1.20) маємо:

    (1.21)

  Рівняння (1.21) не враховує власного поля пучка(кулонівських сил).

  Систему (1.21) перепишемо в скалярній формі:

    (1.22)

Визначимо

тобто

так як , то визначимо:

Тому

    (1.23)

Підставляючи (1.23) в (1.22) отримаємо рівняння руху.

  Але в ці складні рівняння ще входять компоненти електромагнітного поля, які визначаються рівняннями максвела:

    (1.24)

Тут  – електрична і магнітна сталі,  – обємна густина заряду,  – вектор густини струму, - знак транспонування.

  А (1.24) – це рівняння в частинних похідних з складними граничними умовами. Задача заключається не тільки в моделюванні рівнянь руху, а й в розрахунках оптимальних систем.

1.6. Використання диференціальних рівнянь в біології і математичних обчисленнях.

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

 де     (1.25)

 const, , коефіцієнт пропорційності; розвязуючи рівняння (1.25) ми отримаємо таку залежність:

    (1.26)

  Математика. Обчислити невласний інтеграл

    (1.27)

залежний від параметра .

Знайдемо похідну:

Отримали диференціальне рівняння

    (1.28)

При цьому відомо:

    (1.29)

Розв’язуючи задачу Коші (1.28),(1.29), отримаємо:

    (1.30)

1.7. Побудова диференціальнихрівнянь з заданими параметричними сімействами кривих.

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

   (1.31)

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

    (1.32)

Тоді з (1.31) та (1.32) як з системи рівнянь, вилучаємо сталу  і отримаємо шукане диференціальне рівняння першого порядку.

  Якщо ж задано - параметричне сімейство кривих:

    (1.33)

то до (1.33) додаються дані співвідношення:

    (1.34)

з(1.33) та (1.34), як з системи рівнянь, кількість яких , вилучаються сталі  і отримане таким чином співвідношення між

    (1.35)

і буде шуканим диференціальним рівняння  -го порядку.

  В (1.32) та (1.34) означають частинні похідні відповідних порядків за вказаними змінними. При цьому припускаємо, що похідні існують, тобто функції (1.32) та (1.34) являються диференційовними відповідну кількість разів.

  Аналогічно поступають і при складанні систем рівнянь.

Приклад 1.1. Знайти диференціальне рівняння першого порядку, розв’язками якого буде однопараметричне сімейство

    (1.36)

Розвязання. Продиференйіюємо за  праву частину нашого співвідношення в припущенні, що .

    (1.37)

  Враховуючи (1.36) рівність (1.37) перепишемо таким чином:

    (1.38)

З (1.38) знаходимо  

і підставивши в (1.36) отримаємо шукане диференціальне рівняння

    (1.39)

Приклад 1.2. Знайти диференціальне рівняння другого порядку, розв’язками якого буде двопараметричне сімейство

   (1.40)

Розвязання. Згідно описаного вище складаємо систему рівнянь:

    (1.41)

З якої вилучивши  і  знаходимо шукане диференціальне рівняння:

    (1.42).


 

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

34109. Пять базовых техник психотерапии:суггестия, абреакция, манипуляция, разъяснение, интерпретация 95 KB
  Абреакция как основа эмоционального состояния и способности пациента справляться с ним. кратко говоря различных ментальных процессов терапевтом индивидуум в авторитетном положении у пациента индивидуума в зависимом положении независимо от или с исключением рационального или критического реалистического мышления последнего. Техническое использование суггестии может быть главным образом формальным к примеру чтобы в общем индуцировать пациента к фантазии или сну какой бы не была эта фантазия или сон или главным образом...
34110. Понятие регрессии. Роль регрессии в развитии психоаналитической терапии 48 KB
  Понятие регрессии. Роль регрессии в развитии психоаналитической терапии. Процесс регрессии как временный постоянный защитный топический ситуационный. Патологическая и нормальная регрессии их формирование в процессе развития и их значение в функционировании психического аппарата и формировании различных уровней психопатологии.
34111. Принцип психоаналитической нейтральности. Реакции аналитика на пациента: рациональные аффективные, комплиментарные, эмпатические, контрпереносные 48 KB
  Принцип психоаналитической нейтральности. В данной теме особое внимание следует уделить пониманию центрального базового значения психоаналитического понимания нейтральности. Слово НЕЙТРАЛЬНОСТЬ neutrlity и концепция ПСИХОАНАЛИТИЧЕСКОЙ НЕЙТРАЛЬНОСТИ были амбивалентными с самого момента рождения психоанализа. Приветствуемая одно время как настолько фундаментальная что принимается как данность к нейтральности тут же стали тихо относиться как мифу.
34112. Психоаналитическое понятие тревоги и ее типы 85.5 KB
  Тревога и процесс регрессии в психоаналитической ситуации. Тревога рассматривается как архаичный аффект оторвавшийся от первоначального смыслового контекста. Объективная тревога это тревога вызванная известной опасностью. Невротическая тревога вызвана неизвестной опасностью.
34113. Неспецифические аспекты психоаналитической терапии 77 KB
  В данной теме необходимо сформировать четкое представление о неспецифических формах взаимодействия аналитика и пациента. Данный раздел дает четкое представление о вспомогательных формах и методах во взаимодействии аналитика и пациента в рамках психоаналитической терапии. Если он будет это делать с неохотой аналитик может сказать что его интересуют факты. Пациента увязнувшего в неискренней похвале своих родителей можно спросить: Ваши родители действительно замечательные люди Расспрашивание для прояснения очевидности: Вместо того чтобы...
34114. Роль сновидений в психоаналитической терапии и техника работы с ними 73 KB
  Работа сновидения. Роль сновидения в работе психического аппарата. Развитие понимание сновидения и его роли в терапевтическом процессе от З. Классические подходы к пониманию сновидения его роль в общей структуре психики.
34115. Психоанализ и психоаналитическая терапия, основные принципы 67.5 KB
  Основные принципы классического психоанализа разработанного в наследие З. Основные отличия внешние организационные и методологические основы клинического психоанализа психоаналитической терапии. Обратить особое внимание на основные принципы классического психоанализа разработанного З. Предлагаю обсудить вопрос который постоянно в той или иной форме возникает в ходе как профессиональных так и студенческих обсуждений отголоски этой дискуссии звучат и в раздающихся все чаще и чаще утверждениях о том что под брендом психоанализа скрывается...
34116. Показание и противопоказания психоаналитической терапии 62 KB
  Некоторые особенности российского пациента. Так же следует обратить особое внимание на особенности российского пациента и особенности построения терапии в зависимости от психологической конституции. Фрейд полагал что последние две силы связаны между собой и что существует некоторое соответствие внешней реальности и психологической предрасположенности самого пациента Тем самым предполагалось наличие патогенных компонентов в прошлом которые должны предопределять повышенную чувствительность по отношению к определенным обстоятельствам в...
34117. Сеттинг. Определение, взаимозависимость терапевтической задачи и сеттинга 46.5 KB
  Роль сеттинга в построение переходного пространства в рамках котрого происходит развертывание фантазий пациента и осуществляется работа с переносом и сопротивлением. Следует разобраться в ключевой роли сеттинга для формирования у пациента способности восприимать и продуцировать символическую организацию мира. Пациент лежит на кушетке или софе а психоаналитик сидит позади него оставаясь большей частью вне поля зрения пациента стараясь вмешиваться в процесс мышления пациента настолько мало насколько это возможно и не иначе как посредством...