Автор: Пользователь скрыл имя, 14 Декабря 2010 в 22:43, курсовая работа
Метод Монте-Карло использующий математическое моделирование, реализуют, как правило, на ЭВМ. Отечественный и зарубежный опыт показывают, применение ЭВМ для целей оценки точности выходных параметров способствует повышению технического уровня изделий, так как на ЭВМ можно смоделировать и просчитать большее число альтернативных вариантов и выбрать из множества вариантов наиболее лучший (оптимальный или близкий к оптимальному).
Вероятностный метод, а иначе расчетно-аналитический метод с учетом вероятностного рассеивания первичных параметров, решается с помощью аналитических методов.
МИНИСТЕРСТВО ОБРАЗОВАНИЯ РЕСПУБЛИКИ БЕЛАРУСЬ
ИНФОРМАТИКИ И РАДИОЭЛЕКТРОНИКИ”
К КУРСОВОМУ ПРОЕКТУ ПО УЧЕБНОЙ ДИСЦИПЛИНЕ “ТОКТиН” НА ТЕМУ:
“Сравнительная
оценка точности выходного
параметра, полученная
по методу Монте–Карло
и рассчитанная вероятностным
методом”
Научный руководитедь:
В настоящее время, задачи
Точность выходных параметров характеризует степень приближения его истинного значения к номинальному, при отклонении первичных параметров, соответствующих производственным погрешностям. Из всего многообразия методов решения этой задачи большой интерес представляют следующие методы: метод Монте-Карло (метод статических испытаний) и вероятностный метод.
Метод Монте-Карло использующий математическое моделирование, реализуют, как правило, на ЭВМ. Отечественный и зарубежный опыт показывают, применение ЭВМ для целей оценки точности выходных параметров способствует повышению технического уровня изделий, так как на ЭВМ можно смоделировать и просчитать большее число альтернативных вариантов и выбрать из множества вариантов наиболее лучший (оптимальный или близкий к оптимальному).
Вероятностный
метод, а иначе расчетно-
1.Постановка задачи
Исходными данными для
1.РЭС – мультивибратор.
2.Электрическая принципиальная схема.
3.Математическая
модель для выходного
Uвых=1+С1/С2 | (1.1) |
4.Сведения о первичных параметрах:
С1=0.47 пФ±10%;
С2=0.68 пФ±20%;
Исполнение резисторов - дискретное.
Микросхема DА1к140УД8.
U1= 1.2B; U2= 1.2B.
Из задания на курсовое
Для конденсаторов с допускоми на первичный параметр в десять и в двадцать процентов воспользуемся нормальным законом распределения.
В результате решения нашей
задачи двумя методами мы получим
значения выходного параметра (Uвых)-математическое
ожидание и среднее квадратическое отклонение(СКО).
Сравнив полученные результаты двумя
методами, мы сможем сделать вывод о точности
методов.
2.Краткое
описание используемых методов
2.1.Метод Монте-Карло
Этот
метод иначе называется методом
статистических испытаний. С помощью
этого метода можно оценить М(y)-
Основу
метода Монте-Карло составляет процесс
получения случайных реализаций
устройства или процесса. Каждая реализация
описывается значением
Статистическая обработка этого ряда позволяет определить характеристики M(y) и .
При практической реализации метода Монте-Карло используют математическое или физическое моделирование устройств или процессов. В данном курсовом проекте будем использовать математическое моделирование.
Получив требуемые исходные данные можно приступать к определению точности выходного параметра.
Будем действовать следующим образом:
1. Для начала задаемся числом реализации процесса N.
2. Используя генераторы случайных чисел, получаем случайную комбинацию первичных параметров для первой реализации процесса.
3. Подставляем полученную
4. Действия, описанные в пунктах 2…3, повторяем N раз. В итоге получим ряд
5. Выполняем статистическую
6. Проверяем условие
где
- заданная до проведения моделирования
допустимая погрешность (ошибка) в определении
характеристики М(у).
=e-допустимая
ошибка в определении среднего значения
параметра , т.е. разница между оценкой
(y) и истинным значением математического
ожидания М(y) , которая еще допускается.
Допустимая ошибка назначается из условия
Если условие выполняется, то устанавливаем
значение допуска на выходной параметр.
Если условие не выполняется, увеличиваем
число реализаций процесса, корректируем
значение
и вновь по условию (2.1) проверяем, достигнута
ли заданная точность.
2.2.Вероятностный метод
Запишем
уравнение относительной
Коэффициент влияния i-того первичного
параметра
определяется как
При вероятностном методе записанным уравнением воспользоваться сразу не представляется возможным, так как неясно, какие конкретно численно значения требуется подставлять в записанное уравнение в силу случайности этих величин.
Поэтому также оказывается случайной, и для количественного ее описания используют две характеристики:
-математическое ожидание
d(Δy/y)пр – половина поля допуска выходного параметра(Δy/y)пр.
Указанные характеристики могут использоваться для оценки точности выходного параметра. В промышленности в качестве комплексной оценки точности выходных параметров используется производственный допуск на выходной параметр, который устанавливается на основе двух выше записанных характеристик.
Интересующие
нас расчетные соотношения
где -математическое ожидание относительной производственной погрешности i-того первичного параметра;
-коэффициент влияния i-того первичного параметра.
Для d(Δy/y)пр воспользуемся выражением :
, (2.6)
где d(Δx/x)пр –половина поля допуска первичного параметра;
-коэффициент корреляции между i-м и
j-м первичными параметрами.
3.Решение задачи вероятностным
методом
Для решения задачи нам нужно
определить коэффициент
BС1 = С1/(С1+С2) = 0.47/(0.47+0.68) = 0.4 |
BС2 = -С1/(С1+С2) = - 0.47/(0.47+0.68) = - 0.4 |
Определяем математическое ожидание относительной производственной погрешности выходного параметра по формуле (2.6)
Определяем
половину поля допуска выходного
параметра по формуле (2.6).
В оканчательном
виде производственный допуск устанавливается
как
4. Метод Монте-Карло
Вначале мы вводим значения первичных параметров М(xi) и s(xi), после организуем цикл n=50 реализаций. Для каждой из реализаций выходного параметра генерируются значения первичных параметров
с учетом нормального закона распределения по выражению
где x - нормально распределенная случайная величина;
- математическое ожидание и СКО x;
- равномерно
распределенное число в диапазоне (0…1).
Далее генерируются значения
первичных параметров с учетом
нормального закона