Основные задачи математического моделирования - Экономико-математическое моделирование курсовая работа

Главная

Экономико-математическое моделирование
Основные задачи математического моделирования

Общая постановка задачи линейного программирования (ЛП). Приведение задачи ЛП к стандартной форме. Теоремы двойственности и их использование в задачах ЛП. Транспортная задача и её решение методом потенциалов. Интерполирование табличных функций.


посмотреть текст работы


скачать работу можно здесь


полная информация о работе


весь список подобных работ


Нужна помощь с учёбой? Наши эксперты готовы помочь!
Нажимая на кнопку, вы соглашаетесь с
политикой обработки персональных данных

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ЖЕЛЕЗНОДОРОЖНОГО ТРАНСПОРТА ФИЛИАЛ ГОУ ВПО «УРАЛЬСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ПУТЕЙ СООБЩЕНИЯ»
Математические модели вагонов и процессов
Тема: Основные задачи математического моделирования
1. Общая постановка задачи линейного программирования (ЛП)
2. Приведение задачи линейного программирования к стандартной форме
3. Примеры экономических задач, приводящихся к задачам ЛП
3.1 Графический метод решение задач ЛП
3.2 Теоремы двойственности и их использование в задачах ЛП
3.3 Транспортная задача и её решение методом потенциалов
4. Метод аналитического представления экспериментальных данных
4.1 Интерполирование табличных функций
4.2 Эмпирические формулы табличных функций
Математическое моделирование - процесс разработки, отладки, оперирования математическими моделями с целью получения данных о свойствах объекта проектирования. В свою очередь математические модели - совокупность математических выражений определяющих свойства объекта проектирования: универсальность, экономичность, адекватность, точность.
Наиболее эффективное применение вычислительная техника находит при проведении трудоемких расчетов в научных исследованиях. При решении задачи основная роль принадлежит человеку, а машина выполняет роль по разработанной программе. Основные этапы математического проектирования: постановка задачи, которая определяет условия и конечную цель решения; построение математической модели - формулировка задачи с помощью формул и алгебраических (физических) законов; разработка численного метода - позволяет свести задачу к некоторому вычислительному алгоритму; разработка алгоритма и построение блок-схемы; анализ результатов. Процесс решения задачи записывается в виде последовательности элементарных арифметических и логических операций, приводящих к конечному результату.
Основное требование применяемое к математической модели - она должна достаточно точно отражать характерные черты явления, вместе с тем должна обладать сравнительно простой доступностью исследования.
Решения поставленной задачи производятся с помощью численных методов. Для решения математических задач используются следующие основные группы методов: графические - позволяют оценить порядок искомой величины - решение находится путем геометрических построений; аналитические - решение задачи удается выразить с помощью формул; численные - основной инструмент для решения сложных математических задач, позволяют свести решение задачи к выполнению конечного числа арифметических действий над числами.
1. Общая постановка задачи линейного программирования (ЛП)
Задача линейного программирования (ЛП) состоит в нахождении минимума (или максимума) линейной функции при линейных ограничениях.
Общая форма задачи имеет вид: найти при условиях
Здесь и далее нам удобнее считать с и а і вектор - строками, а x и b= (b 1 ,...,b m ) T - вектор столбцами.
Наряду с общей формой широко используются также каноническая и стандартная формы. Как в канонической, так и в стандартной форме
т.е. все переменные в любом допустимом решении задачи должны принимать неотрицательные значения (такие переменные принято называть неотрицательные в отличие от так называемых свободных переменных, на область значений которых подобное ограничение не накладывается). Отличие же между этими формами состоит в том, что в одном случае I 2 = 0, а в другом - I 1 = 0.

В обоих случаях А есть матрица размерности m x n, i-я строка которой совпадает с вектором а i .
Задача ЛП в общей форме сводится (в определенном смысле) к задаче ЛП в канонической (стандартной) форме. Под этим понимается существование общего способа построения по исходной задаче (в общей форме) новой задачи ЛП (в нужной нам форме), любое оптимальное решение которой "легко" преобразуется в оптимальное решение исходной задачи и наоборот. (Фактически, связь между этими задачами оказывается еще более тесной). Тем самым мы получаем возможность, не теряя общности, заниматься изучением задач ЛП, представленных либо в канонической, либо в стандартной форме. Ввиду этого наши дальнейшие рассмотрения задач ЛП будут посвящены, главным образом, задачам в канонической форме.
2. Приведение задачи линейного программирования к стандартной форме
Любая задача линейного программирования приводится к стандартной (канонической) форме основной задачи линейного программирования, которая формулируется следующим образом: найти неотрицательные значения переменных X 1 , X 2 , X n , удовлетворяющих ограничениям в виде равенств:
A 11 X 1 + A 12 X 2 + … + A 1n X n = B 1 ;
A 21 X 1 + A 22 X 2 + … + A 2n X n = B 2 ;
A m1 X 1 + A m2 X 2 + … + A mn X n = B m ;
и обращающих в максимум линейную функцию этих переменных:
E = C 1 X 1 + C 2 X 2 + … + C n X n Ю max
При этом также требуется, чтобы правые части равенств были неотрицательны, т.е. должны соблюдаться условия:
Приведение к стандартной форме необходимо, так как большинство методов решения задач линейного программирования разработано именно для стандартной формы. Для приведения к стандартной форме задачи линейного программирования может потребоваться выполнить следующие действия:
перейти от минимизации целевой функции к ее максимизации;
изменить знаки правых частей ограничений;
перейти от ограничений-неравенств к равенствам;
избавиться от переменных, не имеющих ограничений на знак.
3. Примеры экономических задач, приводящихся к задачам ЛП
Предприятие выпускает два вида продукции А 1 и А 2 , используя при этом три вида сырья S 1 , S 2 , S 3 . Известны запасы сырья: S 1 = 40ab +12 a 2 (У.Е.) , S 2 = 56ab (У.Е.), S 3 = 46ab +20 b 2 (У.Е.). Расход сырья вида S 1 на производство единицы продукции А 1 составляет 2b+a (У.Е.); на А 2 составляет 2а (У.Е.). Расход сырья вида S 2 на производство единицы продукции А 1 составляет 2b (У.Е.); на А 2 составляет 4а (У.Е.). Расход сырья вида S 3 на производство единицы продукции А 1 составляет 2b (У.Е.); на А 2 составляет 2b+3a (У.Е.). Доход от реализации единицы продукции А 1 составляет 3b (У.Е.), единицы продукции А 2 составляет 2b+a (У.Е.). При a=3; b=t+4; t= № варианта.
Составить такой план производства продукции, при котором доход будет максимальным.
3.1 Графический метод решение задач ЛП
Обозначим x 1 - количество продукции А 1 , необходимой изготовить, x 2 - количество продукции А 2 , необходимой изготовить, a=3, b=4+4=8.
На основе данной таблицы составим систему и целевую функцию вида:
Решим задачу графическим методом. Для этого построим график
(приложение 1), три прямые из системы:
Пересечение прямых образует четырехугольник: OABC который является решением системы ограничений.
Целевая функция f x =24х 1 +19x 2 достигает наибольшего значения (max) в одной из вершин этого четырехугольника (B), поэтому для ответа на поставленный вопрос необходимо найти координаты пересечения прямых y1 и y2.
Координаты пересечения прямых y1 и y3 найдем, решая систему уравнений:
х 1 = -4752 / -132 х 2 = (1068- 19х 1 ) / 6
Таким образом, функция f x =24х 1 +19x 2 принимает максимальное значение в точке (36; 64), отсюда можно сделать вывод, что доход будет максимальным при выпуске продукции А 1 =36 у.е и А 2 =64 у.е.
При этом максимально возможная прибыль составит:
f(X max) = f (36; 64) =24 · 36+19 · 64 = 2080 ( У.Е.)
При этом на выполнение данного плана потребуется ресурсов:
S 1 - 19·36+6·64= 1068 - израсходуется полностью
S 2 - 16·36+12·64= 1344 - израсходуется полностью
S 2 - 16·36+25·64= 2176 - останется 208 ( У.Е.)
3.2 Теоремы двойственности и их использование в задачах ЛП
Каждой задаче линейного программирования можно определенным образом сопоставить некоторую другую задачу (линейного программирования), называемую двойственной или сопряженной по отношению к исходной или прямой задаче. Дадим определение двойственной задачи по отношению к общей задаче линейного программирования, состоящей, как мы уже знаем, в нахождении максимального значения функции.
Задача: В условиях предыдущей задачи, найти двойственные оценки цен на сырьё, из решения двойственной задачи по теореме двойственности.
Чтобы составить двойственную задачу составим расширенную матрицу:
Транспонируем матрицу и сформулируем двойственную задачу:
g(y)= 1068y 1 +1344y 2 +2384y 3 min
19·36+6·64 = 1068 - следовательно, y 1 >0
16·36+12·64 = 1344 - следовательно, y 2 >0
16·36+25·64 = 2176 - меньше 2384, следовательно y 3 =0
19 y 1 +16y 2 +0 = 24 / х 6 = 114 y 1 +96 y 2 = 144
6 y 1 +12y 2 +0 = 19 / х 19 = 114 y 1 +228 y 2 = 361
y 2 = -217 / -132 = 1,6439393939393939393939393939394;
y 1 = -0.121212121212121212121212121213
Проверим правильность найденного решения, найденные значения подставим в целевую функцию.
g(y)= 1068y 1 +1344y 2 +2384y 3 = 1068 · (-0.12121212121212121212121212121213) + 1344· 1.6439393939393939393939393939394 + 2384· 0 = 2080 (у.е.)
Ответ: y 1 = -0.1212121212121213(у.е); y 2 =1.643939393939394(у.е); y 3 = 0(у.е.) , увеличение запаса сырья вида S 3 , прибыли не увеличит.
3.3 Транспортная задача и её решение методом потенциалов
Транспортная задача ставится следующим образом: имеется m пунктов отправ- ления А1, А2, ..., Am. B которых сосредоточены запасы каких-то однородных грузов в количестве соответственно a1, а2, ..., аm единиц. Имеется n пунктов назначения B1, В2, …, Вn подавшие заявки соответственно на b1, b2, ..., bn единиц груза. Известны стоимости Сi,j перевозки единицы груза от каждого пункта отправления Ai до каждого пункта назначения Bj. Все числа Сi,j, образующие прямоугольную таблицу заданы. Требуется составить такой план перевозок (откуда, куда и сколько единиц поставить), чтобы все заявки были выполнены, а общая стоимость всех перевозок была минимальна.
Рассмотрим решение «закрытой транспортной задачи», т.е. когда сумма всех заявок равна сумме всех запасов.
Решение транспортной задачи методом потенциалов.
Для того, чтобы некоторый допустимый план X={xi,j}m.n транспортной зада- чи был оптимальным необходимо и достаточно, чтобы ему соответствовала система из m+n чисел U1, U2, ..., Um; V1, V2, ..., Vn, удовлетворяющих условиям Vj-Ui<=Cij (i=1,m; j=1,n) (1), а для всех Xij>0 имело бы место строгое равенство Vj-Ui=Cij. Числа Ui, Vi называются потенциалами соответственно пунктов отправления, а условия (1) (2) называются условиями потенциальности системы.
Для оптимальности плана транспортной задачи необходимо и достаточно, чтобы он был потенциальным. Алгоритм метода потенциалов состоит из предварительного и повторяющегося общего шага.
Предварительный план состоит из следующих операций: составление первоначального ациклического плана перевозок; построение для полученного плана системы m+n чисел U1 ,U2, ..., Um; V1, V2, …, Vn таких, чтобы выполнялись условия Vj-Ui=Cij для всех базисных клеток; проверка построенной системы на потенциальность.
Если система нe потенциальна, т.е. план Х не оптимален, переходим к об- щему шагу.
Общий шаг повторяется до тех пор, пока система не станет потенциальной.
Он состоит из следующих операций: улучшение плана, т.е. замена плана Х новым планом X' со стоимостью перевозок, не превышающей стоимость плана X; построение для X' новой системы потенциалов U'i, V'j путем перестроения старой; проверка системы U'i, V'j на потенциальность.
Предложенный алгоритм сходится за конечное число шагов.
Решение транспортной задача начинается с нахождения опорного плана. Для этого существуют различные способы. Например, способ "северо-западного угла", способ минимальной стоимости по строке, способ минимальной стоимости по столбцу и способ минимальной стоимости таблицы.
Будем заполнять таблицу перевозками постепенно начиная с левой верхней ячейки (“ceвеpo-западного” угла таблицы). Будем рассуждать при этом следующим образом. пункт B1 подал заявку на 18 единиц груза. Удовлетворим эту заявку за счёт запаса 48, имеющегося в пункте A1 , и запишем перевозку 18 в клетке (1:1). После этого заявка пункта B1 удовлетворена, а в пункте A1 осталось ещё 30 единиц груза. Удовлетворим засчёт них заявку пункта В2, (27 единиц), запишем 27 в клетке (1: 2); оставшиеся 3 единицы пункта а1 назначим пункту В3. В составе заявки пункта В3 остались неудовлетворенными 39 единиц. Из них 30 покроем за счет пункта А2, чем его запас будет исчерпан, и еще 9 возьмём из пункта Аз. Из оставшихся 18 единиц пункта Аз 12 выделим пункту В4; оставшиеся 6 единиц назначим пункту В5, что вместе со всеми 20 единицами пункта А4 покроет его заявку. За этом распределение запасов закончено; каждый пункт назначения получил груз согласно своей заявки. Это выражается в том, что сумма перевозок в каждой строке равна соответствующему запасу, а в столбце -- заявке. Таким образом, нами сразу же составлен план перевозок, удовлетворяющий балансовым условиям. Полученное решение является опорным решением транспортной задачи:
На станциях сосредоточен однородный груз в количестве единиц груза, который требуется перевезти на станции назначения в соответствии с заявками каждой станции единиц груза. Известны затраты на перевозку единицы груза с любой станции на любую станцию .
Требуется составить такой план перевозок, чтобы весь груз был вывезен, все заявки были удовлетворены, а суммарные затраты были бы минимальны.
Таблица 1. Запасы и заявки на станциях - участников процесса перевозок.
Таблица 2. Стоимость перевозки единицы груза
Стоимость доставки единицы груза из каждого пункта отправления в соответствующие пункты назначения задана матрицей тарифов
Проверим необходимое и достаточное условие разрешимости задачи: , .
Используя метод наименьшей стоимости, построим первый опорный план транспортной задачи.
В результате получен первый опорный план , который является допустимым, так как все грузы из баз вывезены, потребность магазинов удовлетворена, а план соответствует системе ограничений транспортной задачи. Подсчитаем число занятых клеток таблицы, их 7, а должно быть . Следовательно, опорный план является невырожденным.
Проверим оптимальность опорного плана. Найдем предварительные потенциалы по занятым клеткам таблицы, в которых , полагая, что .
Опорный план не является оптимальным, так как существуют оценки свободных клеток, для которых . Выбираем максимальную оценку свободной клетки (1;5): 6. Для этого в перспективную клетку (1;5) поставим знак «+», а в остальных вершинах многоугольника чередующиеся знаки «-», «+», «-».
Из грузов , стоящих в минусовых клетках, выбираем наименьшее, то есть .
Прибавляем 14 к объемам грузов , стоящих в плюсовых клетках и вычитаем 14 из , стоящих в минусовых клетках. В результате получим новый опорный план.
Проверим оптимальность опорного плана. Найдем предварительные потенциалы по занятым клеткам таблицы, в которых , полагая, что .
линейный программирование транспортный задача
Опорный план является оптимальным, так все оценки свободных клеток удовлетворяют условию .
Оптимальный план предписывает из цеха А 1 перевозить на склад B 2 4 единицы груза, на склад B 4 12 единиц груза, на склад B 5 14. Из цеха А 2 перевозить на склад B 1 18 единиц груза, на склад B 2 10 единиц груза. Из цеха А 3 перевозить на склад B 3 19 единиц груза, на склад B 5 3 единицы груза.
4. Метод аналитического представления экспериментальных данных
Первоначально, исследование любого объекта начинается с накопления информации, как правило экспериментальным путем. Этот этап математического моделирования называется математической обработкой экспериментальных данных. Изучим описание этих данных в аналитическом виде ( с помощью математических формул).
Допустим, исследуется экспериментальными методами, зависимость между двумя величинами x и y. Любой опыт может быть поставлен конечное число раз, это означает что из эксперимента связь между указанными величинами может быть известна в виде конечного количества пар их числовых значений.
Пусть в ходе эксперимента изменили величину Х и измерили зависимость от этих изменений величину Y. Эту зависимость можно рассматривать как функцию y=y(x), значение которой известны только при конечном количестве изолированных значений аргумента X, т.е. как функцию заданную в виде таблицы:
Для построения математической модели объекта иследователю неоходимо научиться определять значения функций, заданной таблично при значения аргумента не входящих в эту таблицу. Наилучший способ это сделать - найти аналитическое выражение для зависимости y=y(x). Геометрически, это означает, найти кривую y=y(x) проходящую через заданную систему точек (X i ; Y i ) i=0, n.
В этом случае можно получить значения функции для любых значений аргумента Х, взятых из области определения функции.
Функцию, в аналитическом виде удобно исследовать, т.е. можно использовать обширный аппарат математического анализа. (дифференцировать, интегрировать)
В инженерной практике существует два подхода к этой задаче:
1.Построение интерполяционного многочлена.
2. Построение эмпирической формулы, по методу наименьших квадратов.
4.1 Интерполирование табличных функций
Пусть функция y=y(x) заданная в виде таблицы (1), требуется найти значения функции при значениях аргумента не входящих в таблицу, с определенной степенью точности. Будем подбирать функцию Ц(x), приближающую табличную функцию y(x) так , чтобы в точках Х i она принимала те же значения что и y i , (как исходная функция).
Этот процесс называется интерполяцией. Функция Ц(x), называется интерполирующей функцией.
Табличные значения аргументов x 0 , x 1, x 2, …. x n , (n+1)называются узлами интерполяции.
Чаще всего интерполирующую функцию ищут в виде интерполяционного многочлена. Многочленом (полиномом) n-ой степени, от неизвестной x называется функция вида:
P n (x) = a 0 + a 1 x+ a 2 x 2 +….. + a n · x n (2)
Действительные числа a 0 , a 1 , a 2 ….. a n - называются коэффициентами многочлена, их количество на единицу больше, чем его порядок (степень).
P 1 (x) = a 0 + a 1 x - линейная функция
P 2 (x) = a 0 + a 1 x+ a 2 x 2 - квадратичная функция
Интерполяционным многочленом называется многочлен, коэффициент которого выбрали таким образом, чтобы значение многочлена в узлах таблицы, в точности совпадали со значениями табличной функции.
Интерполяционный многочлен в форме Лагранжа.
Применяется для интерполирования табличных функций, с небольшим количеством узлов, этот многочлен составляют на основе таблицы (1), по формуле:
Многочлен Лагранжа может быть составлен если в таблице (1) нет одинаковых аргументов, в противном случае, в одном из знаменателей окажется нулевой сомножитель. Каждое слагаемое в правой части формулы (3), представляет собой многочлен в степени n, следовательно и сам многочлен Лагранжа является многочленом n - ой степени.
Для функции, заданной таблицей, составить интерполяционный многочлен Лагранжа. С его помощью найти приближенное значение функции в точке .
с помощью полученного многочлена, можно получить приближенное значение функции y=y(x), при любом значении x, из промежутка для которого составлена таблица [2.5; 3.2].
Интерполяционный многочлен в форме Ньютона.
Существует несколько разновидностей многочленов Ньютона.
P n (x) = С 0 + С 1 (x - x 0 ) + C 2 (x - x 0 ) (x - x 1 ) + ...+ C n (x - x 0 )(x - x 1 ) ... (x - x n -1) (4)
где С i - коэффициенты которые определяются из условий совпадения значений многочлена в узлах интерполяции с соответствующими значениями таблицы функции.
Формула (4) называется интерполяционным многочленом Ньютона для интерполирования вперед. Его применяют когда требуется вычислить значения от аргументов близких к x 0 .
Если нужно вычислить значения функции для аргумента находящегося вблизи x n , то используют многочлен Ньютона для интерполирования назад.
P n (x) = С 0 + С 1 (x - x n ) + C 2 (x - x n ) (x - x n -1) + ...+ C n (x - x n )(x - x n -1)· ... ·(x - x 2 )( x - x 1 ) (5)
Интерполяционный многочлен Ньютона при равноотстоящих узлах.
Пусть функция y=y(x) задана таблицей (1), предположим что разность между любыми соседними значениями аргумента есть постоянное число, называемое шагом таблицы, с обозначением:
Разность между значениями функции в соседних узлах интерполяции, называется разностями первого порядка:
Разностями второго порядка, называются разности полученные из разностей первого порядка, в соседних узлах интерполяции:
Д 2 y i =y i+1 -y i , где i=0,…n-1.
Разностями порядка m называются разности полученные из разностей порядка m -1, в соседних узлах интерполяции.
Д m y i =Д m-1 ·y i-1 - Д m+1 ·y i , где I =0 , …, n-m.
С учетом этого многочлен Ньютона, для интерполирования вперед (4) примет вид:
Многочлен Ньютона для интерполирования назад (5) примет вид:
Для ее записи используются данные записанные из таблицы (7)
Интерполирование с заданной точностью.
Пусть задана функция в виде таблицы (1), с большим количеством равноотстоящих узлов. Требуется найти значение функции y(x), при промежуточном значении аргумента с заданной точностью е нужно:
1) В каждое слагаемое найденного P n (x), в место x подставить , вычислить.
2) При вычислении суммы, отбросить те слагаемые, которые меньше заданной точности.
В условиях предыдущей задачи, построить многочлен Ньютона и чертеж (приложение 2).
Функция задана таблицей. Требуется составить многочлены Ньютона для интерполирования вперед и назад. С их помощью вычислить значения функции в точках и с погрешностью не более чем 0.005.
4.2 Эмпирические формулы табличных функций
Требуется подобрать формулу приближенно описывающую функциональную зависимость заданную этой таблицей. При этом не требуется чтобы приближающая функция в узлах интерполяции совпадала с табличными значениями y i , но не далеко отклонялась от них, кроме этого ее аналитическая формула не должна одержать большое количество параметров, а их количество не должно зависеть от количества табличных точек. Всем этим требованиям удовлетворяет эмпирическая формула.
Задача отыскания эмпирической формулы для табличной функции, состоит из следующих этапов:
1. Определение вида эмпирической формулы.
Реализация этого этапа начинается с построения в прямоугольно-декартовой системе координат, точечного графика функции, по исходным данным таблицы (1). Затем, построенный точечный график, сравнивается с различными кривыми, графики которых известны (линейных, квадратичных, показательных, логарифмических и д.р.) Это дает основание для выбора возможного типа формул.
Например, по результатам таблицы (эксперимента) получен точечный график.
Эту зависимость можно попытаться описать:
а) y=a 0 + a 1 x + a 2 x 2 - уравнение параболы
в) y=a 0 + e bx , b<0 - показательная функция
В формуле (а) три неизвестных параметра (a 0 a 1 a 2 ), в остальных по два. Таким образом, у поставленной задачи может быть несколько решений, начинать исследование нужно всегда с более простой формулы. При рассмотрении графика всегда следует иметь в виду, что при использовании эмпирической формулы используется лишь часть кривой, которая соответствует интервалу изменения аргумента х.
2) Определение параметров эмпирической формулы, методом наименьших квадратов.
Метод наименьших квадратов позволяет подобрать оптимальные значения параметров эмпирической формулы таким образом, чтобы сумма квадратов, всех уклонений значений эмпирической формулы, от опытных данных была бы наименьшей.
Пусть эмпирическая формула является линейной функцией y=a 0 + a 1 x (1) Для определения неизвестных параметров, составляется система линейных алгебраических уравнений (СЛАУ).
Где n - количество узлов в исходной таблице данных. Для составления системы вида (2), на практике обычно составляют таблицу вида:
В табл.3 сначала заполняют первые 5 столбцов, последняя строка которых содержит коэффициент системы (2), по этим данным записывается система (2), затем по формулам Крамера из системы определяется значение параметров a 0 , a 1 . Далее заполняется 6-ой столбец таблицы: при заданных значениях х i , при найденных параметрах a 0 a 1 , вычисляется значение линейной функции. В 7-ом столбце находятся разности между табличными значениями функции (3 столб.) и значениями эмпирической формулы (6 столб.) Последний столбец таблицы - все значения е возведенные в квадрат.
Три последних столбца таблицы (3), заполняются для вычисления среднеквадратичного уклонения (СКУ), которое определяется по формуле:
Оценка точности эмпирической формулы.
Для оценки точности эмпирической формулы используется понятие (СКУ) среднеквадратичного уклонения. СКУ показывает среднюю величину уклонения опытных значений исследуемой зависимости, от расчетных, полученных по эмпирической формуле. Каждое измерение в эксперименте, всегда проводится с некоторой погрешностью и табличные значения функции отличаются от истинных. Одной из задач построения эмпирической формулы, является сглаживание случайных погрешностей измерения, путем подбора такой зависимости между (x) и (y), которая наилучшим образом приближает не только исходные но и промежуточные значения.
Величину СКУ используют для определения пригодности эмпирической формулы. Если ее значение приблизительно равно погрешности экспери- ментальных данных и число параметров много меньше точек в таблице, то эмпирическая формула может использоваться.
Способы сведения сложных эмпирических формул к линейным.
Если зависимость между (x) и (y) не является линейной или квадратичной, то для подбора оптимальных параметров эмпирической формулы, переходят с помощью соответствующей замены переменных, к линейной или квадратичной функции. Неизвестные параметры определяются по методу наименьших квадратов. После этого возвращаются к исходным переменным. Оценка точности полученной эмпирической формулы, определяется по значению (СКУ), которое вычисляется для исходного вида эмпирической формулы.
Пример сложных эмпирических формул и их преобразования.
1) y(x) = a·x b - степенная функция, неизвестные параметры (a) и (b). Для их нахождения прологарифмируем обе части выражения.
Y=A+b·X - преобразовали к линейной функции, для нее СЛАУ будет иметь вид:
где неизвестные (А и b), столбцы для таблицы (3) - i, x i , y i , X i =ln x i , Y i =ln y i , , X i · Y i . Систему решим по формулам Крамера. Находим (А и b). Зная (А), учитывая формулу A = ln a , находим a = e A .
При проведении опыта получена таблица значений двух величин. Задана эмпирическая формула, которая в пределах опыта достаточно точно определяет зависимость между ними.
Требуется, используя метод наименьших квадратов, найти параметры эмпирической формулы, вычислить среднеквадратическое уклонение и построить на одном чертеже графики эмпирической и табличной зависимости.
В таблице приведены результаты наблюдений за n- числом групп вагонов, поступающих на данный путь сортировочного парка станции, от m - числа вагонов в группе:
1. Прологарифмируем заданную нелинейной зависимость . Проведем замену и перейдем к линейной зависимости . Используя исходную таблицу с данными m i , n i составим для неё таблицу значений m i , N i для линейной зависимости.
2. Возьмем прямоугольную систему координат, построить точки (m i , N i ) из таблицы линейной зависимости (приложение 3).
Как видно точки (m i , N i ) располагаются вдоль некоторой прямой.
3. Найдем методом наименьших квадратов коэффициенты в линейной зависимости.
Продифференцировав функцию S по b и A , получим систему уравнений:
4. Найдем уклонения эмпирической зависимости от значений исходной таблицы е i = f(m) - n i , i= 1,2 .., 6.
Вычислим среднеквадратичное уклонение:
Вывод: найденное значение СКУ определяет пригодность нашей эмпирической формулы, т.к. ее значение приблизительно равно погрешности экспериментальных данных.
Построим график эмпирической формулы и на этом же чертеже нанесем точки (m i ; n i ) исходной таблицы (приложение 4).
По результатам наблюдений, проведенным на железнодорожной станции, составлена таблица зависимости времени расформирования составов на сортировочной горке от числа вагонов в составе:
Требуется найти зависимость времени расформирования t (мин) от числа вагонов m, в виде трех формул:
Вычислить среднеквадратичные уклонения (СКУ) и выбрать наиболее подходящую эмпирическую формулу. Построить графики эмпирических зависимостей вместе с точками исходной таблицы. Написать вывод о пригодности эмпирической формулы по значению СКУ.
t(мин) - время расформирования состава
1. Т.к. зависимость линейная преобразования не требуются.
2. Найдем методом наименьших квадратов коэффициенты в линейной зависимости.
Продифференцировав функцию S по b и a , получим систему уравнений:
3. Найдем уклонения эмпирической зависимости от значений исходной таблицы
е i = f(m) - y i , i= 1,2 .., 10.
4. Вычислим среднеквадратичное уклонение:
5. Построим график эмпирической формулы и на этом же чертеже нанесем точки (m i ; t i ) исходной таблицы (приложение 5).
1. Перейдем с помощью преобразования от заданной нелинейной зависимости к линейной . Используя исходную таблицу с данными m i , t i составим для неё таблицу значений M i , t i для линейной зависимости.
2. Возьмем прямоугольную систему координат, построим точки (M i ; t i ) из таблицы линейной зависимости (приложение 6).
Как видно точки (M i , t i ) располагаются вдоль некоторой прямой.
3. Найдем методом наименьших квадратов коэффициенты в линейной зависимости.
Продифференцировав функцию S по b и a , получим систему уравнений:
4. Найдем уклонения эмпирической зависимости от значений исходной таблицы е i = f(m) - y i , i= 1,2 .., 10.
Вычислим среднеквадратичное уклонение:
Построим график эмпирической формулы и на этом же чертеже нанесем точки (m i ; t i ) исходной таблицы (приложение 7).
2. Найдем методом наименьших квадратов коэффициенты в полиномиальной зависимости.
Продифференцировав функцию S по с, b и a , получим систему уравнений:
3. Найдем уклонения эмпирической зависимости от значений исходной таблицы е i = f(m) - y i , i= 1,2 .., 10.
4. Вычислим среднеквадратичное уклонение:
5. Построим график эмпирической формулы и на этом же чертеже нанесем то
Основные задачи математического моделирования курсовая работа. Экономико-математическое моделирование.
Реферат: Об американской прессе. Скачать бесплатно и без регистрации
Как Писать Вступление К Сочинению Егэ
Пособие по теме Учёт, финансирование и кредитование
Реферат по теме Политические воззрения К.П. Победоносцева
Дипломная работа по теме Разработка эффективной стратегии бренда на основе социальных коммуникаций
Курсовая работа: Становление регулярной полиции в России
Контроль Реферат
Сочинения На Тему Мильон Проблем
Доклад: Рунеберг, Йохан Людвиг
Расторжение Трудового Договора Курсовая Работа
Дипломная работа по теме Устаревшая лексика в художественном тексте (на материале поэм А.С. Пушкина)
Итоговое Сочинение 2022 Гроза Островский
Реферат Цели Задачи Методы
Социально Политические Взгляды Аристотеля Реферат
Мотивация Туралы Эссе
Реферат по теме Фонетика эмоциональной речи в ее устной и письменной реализации
Реферат: Фразеология русского языка. Скачать бесплатно и без регистрации
Доклад по теме Конституция Российской Федерации и ее развитие
Моя Любимая Кукла Сочинение
Лирический Герой В Поэзии Есенина Реферат
Концепція спонтанного характеру ринкового порядку - Экономика и экономическая теория реферат
Развитие коммуникативных умений младших школьников с нарушением слуха в процессе внеурочной деятельности - Педагогика курсовая работа
Политология как наука - Политология презентация