- •Предисловие
- •Введение
- •Зачем нужен такой практикум?
- •О чем сказано далее
- •О системе MATLAB
- •Первые задачи
- •Фигуры Лиссажу
- •Биения
- •Волны
- •Основные графические объекты и их использование
- •Отрисовка движущихся кривых
- •Об использовании безразмерных переменных
- •Маятник
- •Свободные колебания
- •Вынужденные колебания
- •Переходные колебания
- •Резонанс
- •О случайном движении
- •Движение частиц в центральном поле
- •Траектория финитного движения
- •Влияние малого возмущения
- •Движение двух частиц
- •Случайные блуждания и диффузия
- •Закономерности случайных блужданий
- •Оценка параметров движения броуновской частицы в жидкости
- •Программа, изображающая случайные блуждания
- •Броуновские частицы в поле тяжести
- •Броуновское движение
- •Случайные силы
- •Корреляционные функции
- •Шары
- •Расчет движения шаров
- •Алгоритм расчета
- •Процедура Balls
- •Динамический хаос
- •Почему движение шаров становится непредсказуемым?
- •Как убедиться в появлении хаоса?
- •Функции распределения
- •Стохастический нагрев и стохастическое охлаждение
- •Потери пучка при прохождении через вещество
- •Эффективные сечения
- •Потери частиц пучка при прохождении слоя
- •Потери энергии
- •Распределение по углам и энергиям
- •Работа с сигналами и модель диодного выпрямителя
- •Работа с сигналами
- •Расчет простейших цепей
- •Статическая модель диода. Решение нелинейных уравнений
- •ПРИЛОЖЕНИЕ
- •Свободные колебания
- •Электрические и магнитные поля
- •Частица в магнитной ловушке
- •Фокусировка пучков частиц
- •Концентрация частиц
- •Приближенные методы решения систем дифференциальных уравнений
- •Моделирование распределения случайных величин
- •Компиляция файлов MATLAB
- •Список литературы
- •Работа в командном окне
- •Вход в систему MATLAB
- •Интерактивный доступ к справочной информации и документации
- •Команда hеlр
- •Команда lookfor
- •Меню Help
- •Редактирование и перевызов командной строки
- •Формат вывода
- •Копия протокола сессии
- •Введение матриц
- •Явное определение матриц
- •Функции построения матриц
- •Операции, выражения и переменные
- •Правила записи операторов
- •Матричные операции
- •Операции с массивами
- •Сохранение данных из рабочей области
- •Операторы for, while, if, case и операторы отношения
- •Цикл for
- •Цикл while
- •Условный оператор if
- •Оператор переключения case
- •Условия (операторы отношения)
- •Функция find
- •Функции MATLAB
- •Скалярные функции
- •Векторные функции
- •Матричные функции
- •M-файлы
- •Файлы-программы, или сценарии
- •Файлы-функции
- •Текстовые строки, сообщения об ошибках, ввод
- •Работа с m-файлами
- •Список путей доступа
- •Работа со списком путей доступа
- •Текущий каталог
- •Средство просмотра и редактирования путей доступа Path Browser
- •Использование редактора/отладчика
- •Отладка m-файлов
- •Сравнение алгоритмов: flops и etime
- •Графика
- •Плоские графики
- •Команда plot
- •Разметка графика и надписи
- •Управление осями при выводе графиков
- •Несколько графиков на листе
- •Специальные виды графиков
- •Столбиковые диаграммы
- •Ступенчатые кривые
- •Гистограммы
- •Изображение функций
- •Трехмерные изображения
- •Одномерная кривая
- •Сеточные поверхности
- •Изолинии
- •Дескрипторная графика (графика низкого уровня)
- •Графические объекты и их иерархия
- •Дескрипторы и работа с ними
- •Свойства графических объектов и работа с ними. Функции get и set
- •Движущиеся графики (анимация)
- •Разработка графического интерфейса пользователя
- •Создание внешнего вида интерфейса
- •Способы взаимодействия графического интерфейса с функциями пользователя
- •Общая структура функции NAME.M графического интерфейса
- •Функционирование графического интерфейса
- •Разработка функции Run
- •Разработка функции Exit
- •Разработка функции Edit
Лучше обойти эти трудности, исследуя движение в пространстве и под действием возмущающей силы f , не параллельной плоскости траектории. В этом случае орбита искажается и поворачивается, но в отрезок не вырождается.
Для наблюдения орбиты, если она оказывается пространственной кривой, удобно изображать сразу три ее проекции на плоскости (X,Y), (Y,Z) и (X,Z), разместив их в трех расположенных рядом графических окнах. (Необходимо согласовать направление осей координат в разных окнах друг с другом.) Можно также с помощью свойства кривой ’ZData’ строить пространственное изображение траектории (подробнее о построении трехмерной кривой см. Дополнение, п. 8.3).
В этом случае становится ясно, что анализ результатов – задача не менее сложная, чем построение модели. (Аналогично обстоит дело и в обычном, «железном» эксперименте.)
Задание 4. Повторите исследование влияния возмущения δU = −fr на финитное движение, если исходное поле сильно отличалось от кулоновского. Можно положить, например, U0 = −α/r + β/r2.
5.3.Движение двух частиц
Движение нескольких взаимодействующих друг с другом тел лишь в исключительных случаях можно рассчитать аналитически. В то же время численный расчет оказывается совсем немногим более сложен, чем для одного тела, т.е. вполне доступен. При исследовании возможных движений в этом случае наиболее сложной частью задачи оказывается представление результатов – перебор и описание различных видов движения, понимание, какие из движений нужно считать качественно различными. Разумеется, такая постановка вопроса не предполагает однозначного ответа, да и исследование в рамках данного практикума не может быть полным. Это скорее иллюстрация наличия практически неограниченного числа возможностей.
Задание 5. Изобразите на экране движение одновременно двух точек в поле
U = −α/r.
Получите траекторию одной точки с точки зрения другой (т.е. траекторию Венеры, например, с точки зрения земного наблюдателя).
Введите взаимодействие частиц друг с другом: U12 = β/|r1 − r2|.
"Выключив"поле U(r) и сохраняя лишь взаимодействие двух тел, получите их траектории.
40