МОСКОВСКИЙ ИНЖЕНЕРНО-ФИЗИЧЕСКИЙ ИНСТИТУТ технический университет Курсовой проект на тему Моделирование значений случайных векторов Студент Барбот А. В. Группа А9-01 Руководитель Власов В. А. Москва 2002 г. Содержание. стр.1. Аннотация. 32. Введение. 43. Необходимые сведения. 64. Исходные данные и обозначения 75. Вывод неизвестных коэффициентов системы уравнений. 86. Реализация программы в среде Matlab. 117. Примеры работы программы. 138. Заключение. 199. Список литературы. 1. Аннотация. Решение многих прикладных задач требует моделирования случайных векторов.
В работе приводится метод моделирования случайных векторов с одинаковым для всех координат одномерным законом распределения, заданной матрицей ковариации и математическим ожиданием составляющих. Для решения этой задачи используется система алгебраических уравнений с неизвестными коэффициентами.
По соответствующему алгоритму разработана программа имитации значений векторов по заданной ковариационной матрице и математическим ожиданиям составляющих с треугольной матрицей преобразования. Изучена возможность покоординатных преобразований. Проведена проверка датчика псевдослучайных чисел системы MATLAB. 2.
Введение. Решение многих прикладных задач, таких как проведение модельных машинных экспериментов с помощью математического моделирования требует моделирования случайных векторов.
Предполагая определенные свойства объекта исследования и характеристики измерительной аппаратуры, исследователь имитирует результаты измерений, обрабатывает их тем или иным способом и сравнивает результат с заложенными ранее характеристиками объекта. Особенно необходимы такие эксперименты при решении некорректных обратных задач. При этом необходимо моделировать не только закономерное влияние на результат измерения свойств объекта исследования и аппаратные искажения, но и случайные погрешности измерений, т.е. случайные величины вектора с заданным законом распределения.
Результат эксперимента, как правило, представляет собой массив отсчетов вольтамперная характеристика, спектр излучения источника света, пространственное распределение яркости в изображении и т.п Если отсчеты считать независимыми случайными величинами их средние значения отражают какие-то закономерности, но к средним прибавлена случайная погрешность, то задача сводится к генерации значений независимых случайных величин погрешностей с нулевым средним и заданным законом распределения.
В общем случае эту задачу легко решить с помощью генератора случайных чисел, равномерно распределенных в интервале 0,1, который встроен практически во все языки программирования высокого уровня. Однако в реальных экспериментах, особенно если они выполняются быстро с помощью автоматизированных измерительных систем, погрешности измерения в различных экспериментальных точках могут быть коррелированны.
Ниже описывается метод моделирования случайных векторов с одинаковым для всех координат одномерным законом распределения, заданной матрицей ковариации и математическим ожиданием составляющих. Для решения этой задачи предлагается использовать систему алгебраических уравнений с неизвестными коэффициентами. Алгоритм получения очередного случайного вектора заключается в следующем по заданным ковариационным матрицам и математическим ожиданиям составляющих случайных векторов вычисляются значения неизвестных коэффициентов системы линейных алгебраических уравнений моделируется случайный вектор, координаты которого независимы и имеют заданное одномерное распределение с помощью указанной системы алгебраических уравнений получается случайный вектор. Доказано, что при выполнении условий реализуемости системы линейных алгебраических уравнений закон распределения координат совпадает с одномерным законом распределения координат, а значения коэффициентов ковариации любой пары равны соответствующим элементам заданной матрицы коэффициентов ковариации.
Моделирующая программа, использующая предложенный метод, определяет значения коэффициентов системы линейных алгебраических уравнений и проверяет выполнение условий реализуемости этой системы.
В случае невыполнения условий реализуемости программа указывает на необходимость корректировки задаваемой матрицы коэффициентов ковариации. Если указанные условия реализуемости выполнены, то программа позволяет выбрать количество объем выборки и размерность моделируемых векторов.
По окончании моделирования программа проверяет соответствие параметров закона распределения координат исходным требованиям, а также находит оценки для полученных в результате моделирования коэффициентов ковариации координат. Программа реализована в вычислительной среде MATLAB. 3. Необходимые сведения. Ниже приводятся необходимые сведения и определения из линейной алгебры и теории вероятности.
Математическое ожидание случайной величины обладает следующими свойствами. 1. Математическое ожидание постоянной равно этой постоянной, т.е 2. Постоянный множитель можно вынести за знак математического ожидания, т.е 3. Математическое ожидание суммы двух случайных величин равно сумме их математических ожиданий, т.е 4. Математическое ожидание произведения двух независимых случайных величин равно произведению их математических ожиданий, т.е Можно доказать, что для случайных величин и для независимых случайных величин. Дисперсия случайной величины обладает следующими свойствами. 1. Дисперсия постоянной величины равно нулю, т.е 2. Постоянный множитель можно вынести за знак дисперсии возведя его в квадрат, т.е 3. Дисперсия суммы двух независимых случайных величин равно сумме их дисперсий, т.е Можно доказать, что для случайных независимых величин Линейное преобразование случайных векторов.
Предположим, что - случайный - мерный вектор с математическим ожиданием и корреляционной матрицей. Введем матрицу преобразования размером и сформируем мерный вектор Можно показать справедливость следующих выражений вектор математического ожидания Если - случайный - мерный вектор, координаты которого являются центрированными случайными величинами, то для выражения справедливо . 4.
Исходные данные и обозначения. Исходными данными для поставленной зада... Вектор задается с помощью генератора случайных чисел, встроенного в си... Так как, аналогично, используя приведенные выше свойства математическо... Так как нижнетреугольная матрица и, то. Эти рекуррентные соотношения позволяют найти элементы матрицы по элеме...
Входные данные размерность вектора, объем выборки Гистограмма нормальн... ков. Входные данные размерность вектора, объем выборки Выходные данные , Ка... Гистограмма распределения координат для выборки Гистограмма распределе... матрица covksi covksi построение выборочной ковариационной матрицы dis...
Заключение. Итогом курсовой работы по теме Моделирование значений случайных векторов явилась разработка алгоритма моделирования значений случайных векторов, с одинаковым для всех координат одномерным законом распределения, заданной матрицей ковариаций и математическим ожиданием составляющих.
Для достижения поставленной цели была изучена система MatLab, с помощью которой была создана моделирующая программа, позволяющая моделировать вектора по заданным параметрам.
В программе присутствуют средства проверки соответствия полученных данных условиям поставленной задачи.
Результаты работы программы представлены в виде графических данных и в виде сохранения массива полученных векторов в файл.
Список литературы . 1. Гусак А.А. Высшая математика.
В 2-х т. Т. 2.Учеб. Пособие для студентов вузов.
Мн. ТетраСистемс, 1998. 448 с. 2. Лазарев Ю.Ф. MatLAB 5.x К. Издательская группа BHV, 2000. 384 с. 3. Потемкин В.Г. Система инженерных и научных расчетов MATLAB 5.x. В 2-х т. Т. 2 М. ДИАЛОГ-МИФИ, 1999. 303 с. 4. Вентцель Е.С. Теория вероятностия.
М. Государственное издательство физико- математической литературы, 1962. 564 с. Моделирование значений случайных векторов Система линейных уравнений, или в матричной форме ,Примеры работы программы. Входные данные размерность вектора, объем выборки Выходные данные , Объем выборки недостаточен для оценки полученных распределений координат, так как ковариационная матрица измерений значительно отличается от заданной.
Поэтому увеличим объем выборки.где Рассматривая ковариацию случайных величин, получаем выражение для элементов ковариационной матрицы, причем справедливо выражение Далее выражаем неизвестные коэффициенты Эти рекуррентные соотношения позволяют найти элементы матрицы по элементам ковариационных матриц Выборочная ковариационная матрица определятся выражением Гистограмма распределения координат для Гистограмма распределения координат для Входные данные размерность вектора, объем выборки Выходные данные , C увеличением объема выборок до элементы выборочной ковариационной матрицы отличаются от заданных не более чем на 5. Гистограмма распределения координат для Гистограмма распределения координат для Документация по программе MonteCarlo.m Входные данные N количество моделируемых случайных векторов. covksi ковариационная матрица элементов моделируемых случайных векторов.
Mksi матрица математического ожидания выходных векторов n размерность случайных векторов вычисляется в программе по заданной ковариационной матрице.
Выходные данные ksi матрица размерности, столбцами которой являются моделируемые случайные векторы. out.txt файл распечатка случайных векторов ksi сохраняется в той же директории, что и программа.