ПОИСК Статьи Рисунки Таблицы Распределение молекул идеального газа по скоростям (распределение Максвелла) из "Компьютеры Применение в химии" Прежде чем составлять программу, надо твердо установить для себя, в каких единицах измерения будут выражены величины, входящие в уравнение, и в какой системе единиц предполагается вводить исходные данные и выводить результат вычислений. В программе есть подсказчик , который напоминает, в каких единицах следует вводить исходные параметры (строки 2(Ю, 4(Ю, 500), и указывает размерность выходных данных (строки 4600—5500). Для ввода и вывода параметров используют общепринятые или удобные единицы измерения, а для вычислений все величины непосредственно в программе переводят в систему СГС. Так, скорость в м/с (переменная U) вводится оператором INPUT в строке 500 и в первой части строки 3000 переводится в см/с (переменная UH). Во второй части строки 3000 размерность универсальной газовой постоянной пересчитывается из Дж/(К моль) (переменная R) в эрг/(К моль) (переменная RH). Кроме того, приведенное выще выражение разделено на число Авогадро (L), чтобы получить ответ в мольных долях. Поскольку М = mL (М — молекулярная масса) и R = kL, отношение т/к ь приведенной выше формуле заменено равным ему отношением M/R. [c.33] Вывод текста программы по частям. Программа МАКСВЕЛЛ занимает довольно много места и целиком не помещается на экране. Для работы с программой на экране ее вызывают по частям, дополнив команду LIST номерами соответствующих строк. Для тех же целей можно использовать замедленный построчный вывод программы на экран. [c.34] БЕЙСИК позволяет записывать в одной строке друг за другом несколько операторов, разделяя их двоеточиями ( ) (см. разд. 5.3.1). Однако следует иметь в виду, что при выполнении оператора СОТО имеет значение только номер той строки, первому оператору которой передается управление, и нет никакой возможности передать управление оператору, стоящему во второй половине строки. [c.35] Для вычисления по довольно длинной формуле в программе МАКСВЕЛЛ используется новый, ранее не встречавщийся прием. В строке 4000 вычисляется первая часть формулы и промежуточный результат присваивается переменной ОН в строке 4100 вычисляется вторая часть формулы, умножается на ВЫ, и произведение запоминается как новое значение переменной ОЫ. В строке 4200 результат, вычисленный в предыдущей строке, для получения ответа в метрах умножается на 100. [c.35] Задание 8. Используйте программу МАКСВЕЛЛ для расчета распределения частиц по скоростям для аргона при 1000 °С для различных значений скорости, постройте график функции распределения. [c.35] Задание 9. Используя профамму МАКСВЕЛЛ , найдите при заданной скорости и молекулярной массе такую скорость, для которой функция распределения максимальна. Решите задачу методом проб и ошибок, не составляя новой программы (программа для нахождения максимума функций будет рассмотрена ниже). Как смещается максимум при изменении температуры Ответы на эти вопросы можно получить в аналитической форме, однако очень простые численные методы в некоторых случаях дают ответ более коротким путем и позволяют избежать ошибок. [c.35] Вернуться к основной статье