Читайте также: |
|
Данные экспериментальных дифрактограмм хранятся в виде таблицы в текстовом файле. В начале происходит выбор файла с помощью графического окна (функция uigetfile), затем данные из файла считываются в два массива: в первом хранится угол 2θ, во втором интенсивность I. Экспериментальную кривую интенсивности I (2 θ) разделили на поляризационный множитель, который вычислили по формуле: , где α =13.30– угол отражения монохроматора. Получили кривую I' (2 θ) = I (2 θ)/ P (2 θ) (рис. 1).
Далее необходимо «сгладить» полученную экспериментальную кривую, т.е. избавится от шумов, сохранив при этом все особенности функции. Для этой операции было рассмотрено несколько методов.
Аппроксимация различными сплайнами
Рис.2. Сплайн шестой степени.
На большей части кривой данный метод даёт приемлемый результат, однако в некоторых местах возникают точки разрыва, кроме того в самом начале возникает ложный локальный экстремум, а начальный участок представляет собой наиболее ценную для получения дальнейшей ФРРА информацию. Данное явление заметно снижается при повышении степени сплайна, однако эффект сглаживания так же ослабляется.
Рис.3. Сплайн двадцатой степени.
Рис.4. Виден ложный пик при малых углах.
Дата добавления: 2015-07-08; просмотров: 129 | Нарушение авторских прав
<== предыдущая страница | | | следующая страница ==> |
Пакет Mathematica | | | Взвешенная локальная регрессия |