ЗАДАЧА 14.2. Динамическое определение энтропии
а. Возьмите N = 10 и с помощью программы entropy вычислите час-тоту совпадений Rn и энтропию Sn для каждого макросостояния мо-дели частиц в ящике.
б. Сравните полученные приближенные результаты для с точными,приведенными в табл. 14.1.
в. Оцените погрешность своих вычислений Sn, пропустив программус более длинными последовательностями.
г. Если возможно, повторите все указанные выше вычисления длябольших значений N. Какое наибольшее значение TV можно рассматри-вать, используя величину micro, как оиа определена в (14.4) и впрограмме entropy? Существует ли практическое ограничение сверхуиа значения N, которое можно рассматривать?
В программе entropy формирование конфигураций осуществляется гораздо быстрее, чем подсчет совпадений. Первое пропорционально числу обменов, а второе пропорционально квадрату числа обменов. Поэтому, вероятно, было бы удобно запомнить значения micro, порождаемые подпрограммой exchange, в отдельном файле. Таким образом мы можем отделить анализ данных от их формирования и при необходимости, использовать различные методы анализа. В языке True BASIC файл открывается с помощью инструкции:
OPEN #1: name "datal", access output, create new
Запись значений micro в файл #1 производится инструкцией:
PRINT #1: micro(iexch) Файл #1 закрывается инструкцией:
CLOSE *1
Чтобы прочитать данные, используем инструкции:
OPEN #1: name "datal", access input INPUT #1: micro(iexch) CLOSE #1