- •ПРЕДИСЛОВИЕ
- •ОПРЕДЕЛЕНИЕ И НАЗНАЧЕНИЕ МОДЕЛИРОВАНИЯ
- •1.1. ЧТО ТАКОЕ МОДЕЛЬ?
- •Место моделирования среди методов познания
- •Определение модели
- •Определение модели
- •Цели моделирования
- •1.2. КЛАССИФИКАЦИЯ МОДЕЛЕЙ
- •Материальное моделирование
- •Идеальное моделирование
- •Когнитивные, концептуальные и формальные модели
- •1.3. Классификация математических моделей
- •Классификационные признаки
- •Классификация математических моделей в зависимости от оператора модели
- •Классификация математических моделей в зависимости от целей моделирования (рис. 1.11)
- •Классификация математических моделей в зависимости от методов реализации (рис. 1.12)
- •Задания для самостоятельного выполнения
- •Глава 2
- •ЭТАПЫ ПОСТРОЕНИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ
- •2.1. ОБСЛЕДОВАНИЕ ОБЪЕКТА МОДЕЛИРОВАНИЯ
- •2.2. КОНЦЕПТУАЛЬНАЯ ПОСТАНОВКА ЗАДАЧИ МОДЕЛИРОВАНИЯ
- •2.3. МАТЕМАТИЧЕСКАЯ ПОСТАНОВКА ЗАДАЧИ МОДЕЛИРОВАНИЯ
- •2.4. ВЫБОР И ОБОСНОВАНИЕ ВЫБОРА МЕТОДА РЕШЕНИЯ ЗАДАЧИ
- •stop
- •2.5. РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ В ВИДЕ ПРОГРАММЫ ДЛЯ ЭВМ
- •2.6. ПРОВЕРКА АДЕКВАТНОСТИ МОДЕЛИ
- •2.7. ПРАКТИЧЕСКОЕ ИСПОЛЬЗОВАНИЕ ПОСТРОЕННОЙ МОДЕЛИ И АНАЛИЗ РЕЗУЛЬТАТОВ МОДЕЛИРОВАНИЯ
- •Задания для самостоятельного выполнения
- •ПРИМЕРЫ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
- •3.1. СТАТИЧЕСКИЙ АНАЛИЗ КОНСТРУКЦИЙ
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Математическая постановка задачи
- •Решение задачи
- •Анализ результатов
- •Концептуальная постановка
- •Математическая постановка задачи
- •Методика решения задачи
- •Анализ результатов
- •3.3. ДИНАМИКА ПОПУЛЯЦИИ
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Решение задачи
- •Анализ результатов
- •Математическая постановка задачи для модели Ферхюльста
- •Решение задачи
- •Анализ результатов
- •Численное исследование модели Ферхюльста
- •3.4. МОДЕЛЬ КОНКУРЕНЦИИ ДВУХ ПОПУЛЯЦИЙ
- •Математическая постановка задачи
- •Качественный анализ задачи
- •Численное исследование модели конкуренции популяций
- •3.5. ГАРМОНИЧЕСКИЙ ОСЦИЛЛЯТОР
- •Содержательная постановка задачи
- •Концептуальная постановка задачи
- •Математическая постановка задачи
- •Решение задачи
- •Качественный анализ задачи
- •Численное исследование модели
- •Качественный анализ задачи
- •Решение задачи
- •Задания для самостоятельного выполнения
- •Глава 4
- •СТРУКТУРНЫЕ МОДЕЛИ
- •4.1. ЧТО ТАКОЕ СТРУКТУРНАЯ МОДЕЛЬ?
- •4.2. СПОСОБЫ ПОСТРОЕНИЯ СТРУКТУРНЫХ МОДЕЛЕЙ
- •4.3. ПРИМЕРЫ СТРУКТУРНЫХ МОДЕЛЕЙ
- •Задания для самостоятельного выполнения
- •МОДЕЛИРОВАНИЕ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ
- •5.1. ПРИЧИНЫ ПОЯВЛЕНИЯ НЕОПРЕДЕЛЕННОСТЕЙ И ИХ ВИДЫ
- •5.2. МОДЕЛИРОВАНИЕ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ, ОПИСЫВАЕМОЙ С ПОЗИЦИЙ ТЕОРИИ НЕЧЕТКИХ МНОЖЕСТВ
- •5.3. МОДЕЛИРОВАНИЕ В УСЛОВИЯХ СТОХАСТИЧЕСКОЙ НЕОПРЕДЕЛЕННОСТИ
- •5.4. МОДЕЛИРОВАНИЕ МАРКОВСКИХ СЛУЧАЙНЫХ ПРОЦЕССОВ
- •Задания для самостоятельного выполнения
- •ЛИНЕЙНЫЕ И НЕЛИНЕЙНЫЕ ЯВЛЕНИЯ
- •6.1. О ЗАКОНЕ ГУКА И ГРАНИЦАХ ЛИНЕЙНОСТИ
- •6.3. О ПОСТРОЕНИИ СПЛОШНЫХ МОДЕЛЕЙ. ВЫВОД ВОЛНОВОГО УРАВНЕНИЯ ИЗ ЗАКОНОВ МЕХАНИКИ
- •6.4. РЕШЕНИЕ ВОЛНОВОГО УРАВНЕНИЯ МЕТОДОМ ФУРЬЕ
- •6.6. УРАВНЕНИЯ МАКСВЕЛЛА
- •6.7. О КЛАССИФИКАЦИИ КВАЗИЛИНЕЙНЫХ СИСТЕМ
- •6.8. СВЯЗЬ НЕПРЕРЫВНОГО И ДИСКРЕТНОГО НА ПРИМЕРАХ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ И УРАВНЕНИЯ ШРЕДИНГЕРА
- •6.9. О ПОЛЬЗЕ ФЕНОМЕНОЛОГИИ ПРИ ПОСТРОЕНИИ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
- •6.10. АНАЛИЗ ПОДОБИЯ И РАЗМЕРНОСТИ
- •6.11. АВТОМОДЕЛЬНОСТЬ
- •6.12. САМООРГАНИЗАЦИЯ И СТРУКТУРЫ В НЕЛИНЕЙНЫХ СРЕДАХ
- •6.13. О НЕЛИНЕЙНЫХ ВОЛНАХ В СПЛОШНЫХ СРЕДАХ
- •6.14. ИЕРАРХИЧЕСКИЕ МОДЕЛИ ТУРБУЛЕНТНОСТИ И МНОГОМАСШТАБНЫЕ ФУНКЦИОНАЛЬНЫЕ БАЗИСЫ
- •Иерархический базис для турбулентных полей
- •Одномерный иерархический базис
- •Двумерный базис
- •6.15. ВЕЙВЛЕТЫ
- •Непрерывное вейвлет-преобразование
- •6.16. ВЕЙВЛЕТ-АНАЛИЗ ВРЕМЕННЫХ КОЛЕБАНИЙ
- •6.17. О ФРАКТАЛАХ И ИХ ПРИМЕНЕНИИ
- •Примеры фракталов
- •Подобие и скейлинг
- •Множества Мандельброта и Жюлиа
- •Фрактальная размерность кластеров
- •Экспериментальные методы определения фрактальной размерности
- •Модель случайных фракталов для описания растущих дендритных структур
- •Результаты применения модели случайных фракталов
- •6.18. НЕЛИНЕЙНЫЕ МОДЕЛИ ДНК
- •Структура и физические свойства ДНК
- •Модель Пейрара-Бишопа
- •Задания для самостоятельного выполнения
- •Глава 7
- •7.2. ИМИТАТОР СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ
- •7.3. КЛЕТОЧНЫЕ АВТОМАТЫ
- •7.4. МОДЕЛИРОВАНИЕ ДИСЛОКАЦИЙ В МЕТАЛЛЕ
- •Самоорганизация дислокаций в модели клеточных автоматов
- •Задания для самостоятельного выполнения
- •ЯЗЫК ФОРМАЛЬНОГО ОПИСАНИЯ АЛГОРИТМОВ
- •====== Приложение 2
- •П2.1. Решение уравнений высоких степеней и трансцендентных уравнений с одним неизвестным
- •П2.2. Решение систем линейных уравнений
- •БИБЛИОГРАФИЧЕСКИЙ СПИСОК
- •ПРЕДМЕТНЫЙ УКАЗАТЕЛЬ
- •ОГЛАВЛЕНИЕ
- •Введение в математическое моделирование
По образцу, данному выше для квазилинейных систем (добав ляя к независимым переменным v), или элементарными сложени ем и вычитанием уравнений получаем:
£ (P + v )+ (v + n )|-(P + v ) = О,
Э/ |
дх |
| ( P |
_ v) + (v_0)|-(i» _ v )= о. |
dt |
дх |
Очевидно (см. (6.31),(6.34)), характеристиками системы будут
кривые, dx/dt = v ± а, на которщх постоянны соответствующие ве
личины P±v. Здесь в отличие от системы, описывающей малые возмущения, а представляет собой поле местной скорости звука,
±а —относительные скорости распространения волн по газу. Вдоль первого семейства характеристик распространяется волна разрежения, вдоль второго —волна сжатия [63].
6.8. СВЯЗЬ НЕПРЕРЫВНОГО И ДИСКРЕТНОГО НА ПРИМЕРАХ УРАВНЕНИЯ КОЛЕБАНИЙ СТРУНЫ И УРАВНЕНИЯ ШРЕДИНГЕРА
Лишь первый исчез, как сейчас же в ином положении
Новый родится за ним, а нам кажет ся, - двинулся первый.
Лукреций
Вопросы о движении, переходе постепенных количественных изменений в качественные, появлении у целого свойств, которы ми не обладает ни одна из его частей, являются одними из клю чевых вопросов современного фундаментального естествознания. Эти вопросы имеют глубокие философские корни. Развитие есте ственных наук заставляет еще и еще раз возвращаться к ним. Уче ных XIX в. поразило наличие волновых свойств у света, который они представляли как поток дискретных частиц. К глубокому пе ресмотру фундаментальных понятий привело в XX в. создание квантовой механики. Оказалось, что дискретные и непрерывные свойства материи нельзя противопоставлять друг другу, что они неразрывно связаны между собой. Выяснилось и другое важное об
стоятельство. Анализ многих явлений требует сочетания дискрет ного и непрерывного подходов. И вопрос о соотношении тех и других свойств при построении теории оказывается далеко не про стым. От его успешного решения часто зависит, насколько глубо ко нам удается разобраться в изучаемом объекте.
Рассмотрим задачу, при решении которой был развит ряд клю чевых идей. Проведем анализ поведения упругой струны, по ко торой ударили в начальный момент времени. Остановимся преж де на системе, состоящей из точечного груза массой т, к которо му прикреплены две одинаковые упругие горизонтальные нити длиной IQ / 2 , натянутые силой F0 (сила тяжести отсутствует). При отклонении груза от положения равновесия появляется возвраща
ющая сила, пропорциональная отклонению F = -2 F0sinа - 4F Q U / I Q ,
тогда согласно второму закону Ньютона имеем:
(6.46)
Уравнение (6.46) описывает колебания с круговой частотой со:
и - A sinart+В costo/,
где константы А и В определяются начальными положением и скоростью груза. Мы решили задачу о колебании струны, вся масса которой сосредоточена в центре. Для однородной струны с массой т и длиной /0 это слишком грубое приближение; разумнее заме нить ее набором из N шариков с массой ц = m/N, расположенных на расстоянии h = /0/JVи соединенных нитями, натянутыми силой F0 (так рассуждали, в частности, при выводе уравнений колебаний струны Иоганн и Даниил Бернулли). Если ик —отклонение к-то шарика от положения равновесия, то при условии, что разница в отклонениях соседних шариков мала,
|
|
|
ик^ -Ъ 1 к +икЛ |
l0 dt2 |
|
~ 0 |
h |
^ ЭЧ _ Ы |
u k +1 - 4 + u k -1 |
||
Эt2 |
m |
|
h2 |
В пределе при h -»0 |
получаем уже известное нам волновое |
уравнение
Э 2 и
c2 = FQlQ/m.
Этот вывод уравнения впервые был сделан Даламбером, кото рый не только записал указанное уравнение, но и нашел его об щее решение (6.39) в виде суперпозиции двух волн
и = f(x+ ct) + g(x-ct).
Однако использовать данное решение для струны конечных размеров непросто. Действительно, после удара по струне вправо и влево идут волны. Они доходят до концов струны, отражаются, идут в обратную сторону —устанавливается некий режим, описы вать который с помощью полученной формулы неудобно. Возмо жен другой путь, предложенный Фурье. Так как струна соверша ет колебательные движения, решение задачи ищут в виде
и = {Asin сЩ+ В cos a>t)z(x),
что дает в итоге (см. разд. 6.4)
и = X (At s i n |
+вп cos« у )sin ~Т~- |
п=0 |
Ч> |
Решение нашлось в виде суперпозиции стоячих волн, для ко торых на длине струны укладывается целое число п полуволн. Соб ственные значения Хп определяют, с какой частотой может коле баться струна в конфигурации л-й стоячей волны, форму которой описывает собственная функция z„ = sin (nnx/l0).
Найденные решения по виду сильно различаются. Равноправ ность их не очевидна. Этот вопрос стал причиной дискуссии в середине XVIII в. (спор о струне) между Эйлером, Даламбером и Лагранжем. Дискуссия позволила убедиться не только в эквивален тности двух решений (ведь исходная задача имеет решение и оно единственное!), но и лучше разобраться в уравнениях. Решение, полученное Фурье, дает возможность выяснить соотношение не прерывного и дискретного в этой задаче. Если мы имеем автомо дельное решение волнового уравнения в виде стоячей волны
w = cn(Osin(n«x//0), то окажется, что для разных cn(t) получаются
уравнения
d 2c (t)
— ^ -- + co2cn(/) = 0, n =0, 1, 2,...
Но это опять уравнения колебаний. Значит, струна оказывается эквивалентной бесконечному множеству независимых колеблю щихся грузов.
Интересно и другое: в непрерывной задаче, описывающей ко лебания струны, есть дискретный набор собственных частот. Не прерывное и дискретное вновь оказываются тесно связанными.
Глубокая связь дискретного и непрерывного отмечена и в физике микромира, где в одних случаях материю удобно рассмат ривать как электромагнитную волну, а в других —как поток час тиц (квантов). Микрочастицы в некоторых опытах ведут себя как волны, например, испытывая дифракцию и интерференцию. В то же время было установлено, что свет (волна) квантуется, регист рируются дискретные порции света —фотоны. Оказалось, что и здесь довольно гармоничный дуализм может быть описан с помо щью линейного уравнения математической физики —уравнения Шредингера. Выражение простейшей (плоской) волны, описыва ющей колебания в пространстве с частотой со и волновым векто ром к, записывается как exp (/foe —/со/).
Сопоставим этой волне частицу с энергией Е= цсо и импуль
сом p = h k , где h —постоянная Планка. Выразив со и £ из после дних соотношений и подставив их в формулу для волны, получим волновую функцию для частицы с энергией Е и импульсом р:
■Р .Е .
у(х,/) = ехр l-X -1 —t
Эта комплексная функция определяет плотность вероятности нахождения частицы во времени и пространстве. Нетрудно убе диться в справедливости тождеств
и - к
/ Э/
Но это означает, что энергия Е является собственным значением
оператора |
, а квадрат импульса р2 —собственным значени- |
/ |
dt |
ем оператора -й2Д\у, при этом в обоих случаях выступает в ка честве собственной функции.
Если частица массы т движется в потенциальном поле V(x),
то в силу закона сохранения энергии Е = р2/2т+У Чтобы полу
чить волновой аналог этого соотношения, мы должны заменить Е, р2 и V соответствующими дифференциальными операторами:
А |
Эц/ |
А2 |
т. |
- 7 |
¥ - 2 ^ 4 v + ^ ’ |
откуда и следует фундаментальное уравнение волновой и кванто вой механики
= + (6 -4 7)
открытое австрийским физиком Эрвином Шредингером в 1926 г. Оно описывает движение частицы в заданном потенциальном поле.
Квадрат амплитуды волновой функции J\|i(x,t)^*(x,t)dx = P(G,t)
G
определяет вероятность P{G,t) нахождения частицы в момент вре мени / в области G. Для оператора (6.47) можно поставить задачу на собственные значения:
tP" |
°° |
—— Aw + K\|f = Е у , |
f \|A|/*fltc=l, |
2m |
J |
|
—oo |
в результате чего получим спектр Еп значений энергии частицы, которые она может принимать, двигаясь в заданном потенциаль ном поле. Соответствующие им собственные функции \|/„ показы вают, с какой вероятностью можно обнаружить частицу в разных точках пространства.
Решив уравнение Шредингера для кулоновского потенциала V= —е2//* (в —заряд электрона, г —расстояние до протона), мож но найти энергетические уровни атома водорода, с огромной точ ностью согласующиеся с экспериментом.
До появления уравнения Шредингера существовала матричная квантовая механика Гейзенберга, использующая в качестве аппа рата пространство бесконечномерных векторов с ограниченным скалярным квадратом. Теория же Шредингера использует вектор ное пространство функций с ограниченным скалярным квадратом