Задача.Изобразить на графике приблизительную форму электронных облаков в атомах.
Анализ. По современным представлениям электронные уровни в атоме определяются четырьмя квантовыми числами. Форма электронного облака определяется двумя из этих чисел:
• число l определяет тип орбитали (значения 0-3 соответствуют s-, p-, d- и f- орбиталям);
• число m определяет магнитный момент электрона и может изменяться в диапазоне от -l до l.
При m=0 форма электронного облака определяется на основе многочленов Лежандра первого рода:
, где l – степень многочлена.
В этом случае .
Параметрическое задание соответствующей поверхности имеет следующий вид:
Углы изменяются в диапазоне от 0 до 2л.
1. Запустите программу MathCad.
2. Определите переменную l, которая укажет тип орбитали.
l:=3
3. Построение поверхности будем производить по точкам. Задайте два диапазона, которые будут определять изменение параметров , задающих поверхность. Удобно определить границы диапазона в целых числах (через точку с запятой, на экране изображаются две точки), а затем произвести перемасштабирование на отрезок [0; 2л].
i:=0..100 j:=0..100
4. Определите двумерные матрицы, определяющие значения координат х, у и z в зависимости от значения параметров. Используйте названия переменных Х0, Y0 и Z0.
Рис. 6. Трехмерное изображение электронной f-орбитали
5. Дайте команду Insert > Graph > Surface Plot (Вставка > График > Поверхность) или воспользуйтесь кнопкой Surface Plot (Поверхность) на панели инструментов Graph (График).
6. В появившейся области графика вместо заполнителя укажите имена отображаемых матриц через запятую, заключив все их в скобки: (X0,Y0,Z0).
7. Чтобы изменить формат построенного графика, дважды щелкните на его области. Откроется диалоговое окно 3-D Plot Format (Формат трехмерного графика).
8. На вкладке General (Общие) установите флажок Equal Scales (Равный масштаб), чтобы обеспечить одинаковый масштаб по осям координат.
9. На вкладке Appearance (Оформление) установите переключатель Fill Surface (Заливка поверхности), чтобы обеспечить заливку построенной поверхности.
10. На вкладке Lighting (Подсветка) включите режим освещения поверхности. Установите флажок Enable
Lighting (Включить подсветку), отключите все источники света, кроме первого.
11. На панели Light Location (Размещение источника)задайте координаты источника света. Используйте кнопку Применить, чтобы сразу видеть последствия сделанных настроек. По окончании настройки закройте диалоговое окно щелчком на кнопке ОК.
12. Путем протягивания мыши в области графика измените направление осей координат, чтобы изображение было видно наиболее отчетливо.
13. Изменяя значение l, можно увидеть форму электронных облаков для разных орбиталей, в том числе и не встречающихся в природе.
Мы научились строить трехмерные графики с изображением объемных поверхностей, заданных параметрически. Это фактически означает умение изображать любые фигуры, которые могут потребоваться в ходе практической работы.
Упражнение 8. Решение дифференциальных уравнений (30 мин)
Задача. Найти функцию у(х), удовлетворяющую дифференциальному уравнению
и имеющую значение 0 при х = 0 .
Анализ. Это простое дифференциальное уравнение допускает точное аналитическое решение, однако в данном упражнении предполагается использование стандартной функции программы MathCad, осуществляющей численное решение данного уравнения. Результат вычислений можно после этого сравнить с точным решением.
1. Запустите программу MathCad.
2. Задайте начальное значение функции как элемент вектора y, размерность которого соответствует числу решаемых уравнений (в данном случае единице): yo := 0.
3. Создайте функцию Т(х, у), которая вычисляет значение производной при заданных значениях независимой переменной и неизвестной функции:
4. Определите начальное (точка 0) и конечное значение отрезка интегрирования. а := 0, Ь := 12-p.
5. Укажите число шагов интегрирования.
K:=20.
6. Вычислите численное решение уравнения при помощи функции rkfixed.
Рис. 7. Графики численного и точного решения дифференциального уравнения
Результат вычислений — матрица Z с двумя столбцами, первый из которых содержит значения независимой переменной, а второй — соответствующие значения функции.
7. Постройте график полученного решения.
8. Определите аналитическое решение данного уравнения при тех же начальных условиях.
9. Нанесите аналитическую кривую на тот же график и сравните поведение численного и точного решения.
10. Измените число шагов, на которые делится отрезок интегрирования, и исследуйте, как изменяется результат расчета при уменьшении и увеличении этого параметра.
Мы научились численно решать дифференциальные уравнения первого порядка с помощью программы MathCad. Использованный метод без изменений переносится на системы, содержащие два или большее число дифференциальных уравнений. Увеличение величины шага интегрирования ускоряет получение результата, но снижает его точность. При слишком большой величине шага результат расчетов может вообще не соответствовать реальному решению.
Практикум по математическойсистеме Mathcad: Метод, указ. по дисциплине «Информатика» для студентов специальностей 1105, 1106, 1703 и 1204 очно-заочной формы обучения
Составитель В.Г. Рыжков
Издается полностью в авторской редакции
Тем.план 2001, поз. 368
Лицензия ЛР № 020377 от 22.01.97
Подписано в печать 27.11.2001. Формат 60х84 1/16. Бумага тип.№ 1
, где l – степень многочлена.
.
изменяются в диапазоне от 0 до 2л.
и имеющую значение 0 при х = 0 .
Рис. 7. Графики численного и точного решения дифференциального уравнения