
- •Предисловие
- •Введение
- •Зачем нужен такой практикум?
- •О чем сказано далее
- •О системе 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
ПРИЛОЖЕНИЕ
A. Дополнительные задачи
Здесь приведены примеры задач, предлагаемых студентам для самостоятельной работы.
Свободные колебания
1.Изобразите на экране монитора нормальные колебания двойного маятника (см. [9, задача 6.3]).
2.Колебания молекулы CO2 можно представлять как наложение гармонических колебаний. Изобразите каждое из этих колебаний и их суперпозицию (см. [7, §24, задача 1]).
Электрические и магнитные поля
3.Изобразите силовые линии магнитного поля, создаваемого парой прямых токов и однородным полем, перпендикулярным этим токам (магнитная ловушка), для следующих случаев:
а) токи текут в направлении, перпендикулярном плоскости ХУ; б) токи текут по прямым, скрещивающимся в пространстве.
Частица в магнитной ловушке
4.Исследуйте движение заряженной частицы в магнитных полях, заданных в задачах 3а, 3б.
E |
= |
e |
E |
Сила, действующая на частицу в магнитном поле, - сила Лоренца F |
c |
[EvB]. |
Так как она зависит от скорости, схема расчета, используемая в работе “МАЯТНИК” и “ПЛАНЕТА”, оказывается здесь недостаточно точной. Можно рекомендовать схему расчета, основанную на равенстве вида
x(t + h) = x(t) + hx˙ (t) + |
h2 |
h3 |
˙ |
4 |
|
|
|
x¨(t) + |
|
x¨(t) + O(h |
), |
||
2 |
6 |
˙
причем x˙ (t) = v(t), x¨(t) = f(t) = F (t)/m, а x¨(t) с той же точностью
выражается как f(t)−f(t−h).
h
80

5.Исследуйте движение заряженной частицы в поле магнитного диполя (моделируя тем самым движение частиц в радиационных поясах).
Фокусировка пучков частиц
6.Пучок частиц, летевших параллельно оси X, фокусируется в поле U(x, r). Изобразите движение группы частиц, стартовавших с большого расстояния с разными прицельными параметрами одновременно:
а) U = |
Ay2 |
|
|
Ay2 |
|
|
, |
б) U = |
|
. |
|
r2+B2 |
(r2 +B2)2 |
7.Пучку частиц, фокусируемых в поле, можно сопоставить картину фокусировки волн в неоднородной среде. Можно показать (но мы не будем этого делать), что длина волны должна быть обратно пропорциональна скорости
частицы. Иначе говоря, фронт волны должен смещаться в направлении Ev, но
на расстояние
∆Er Ev∆t/v2.
Получите картину волновых поверхностей, соответствующих фокусировке пучка в полях, заданных в задаче 9.
8. Пучок частиц летит параллельно оси X и рассеивается в поле
√ |
|
|
|
2 |
2 |
|
|
U(r) = −α/ r |
+ a . |
Определить распределение рассеянных частиц по углам. Считать, что падающий поток частиц равномерно распределен по площади своего поперечного сечения.
9.Движение частиц в стоячей волне на поверхности жидкости
x = x0 + aekz0 cos kx0 cos ωt, z = z0 + aekz0 sin kx0 cos ωt.
Получите картину движения частиц. Дуги траекторий, прорисованные частицами за время τ 1/ω, дадут представление о линиях тока жидкости.
Концентрация частиц
10.Рассматривается диффузия (случайные блуждания) частиц в поле U = k2 (x2+ y2 + z2). Определите наблюдаемую концентрацию частиц в зависимости от r = √x2 + y2 + z2.
81