LEC-18 (1014362), страница 7
Текст из файла (страница 7)
МЕТОД НАИМЕНЬШИХ КВАДРАТОВ ПРИ В=0: a=-0.08371.
Анализируя полученные значения коэффициентов для разных методов видим явно несовпадающие значения коэффициента А.
В0
МЕТОД ВЫБРАННЫХ ТОЧЕК: a= -2.91832; b= 3.88663.
МЕТОД СРЕДНИХ: a= -1.41004; b= 3.08540.
МЕТОД НАИМЕНЬШИХ КВАДРАТОВ: a= -2.98023; b= 3.95859.
Анализируя полученные значения коэффициентов для разных методов, видим, что метод средних дает для коэффициентов А и В выпадающие значения. Не принимая их во внимание получим: А=-2.95; В = 3.92.
7.2.3 Нелинейная парная регрессия
ВЫБОР ЭМПИРИЧЕСКИХ ФОРМУЛ ДЛЯ АНАЛИЗА НЕЛИНЕЙНЫХ ЗАВИСИМОСТЕЙ
ТЕОРЕТИЧЕСКАЯ ЧАСТЬ
Выше была рассмотрена линейная зависимость вида у=Ах+В для случаев, когда А0; В=0 и А0, В0. Но, к сожалению, построение этой зависимости не дает ответа на вопрос о том, какая аналитическая зависимость наилучшим образом подходит (описывает) имеющееся распределение. Наиболее популярные на практике эмпирические зависимости имеют вид:
1) линейная функция: у = Ах + В;
2) показательная функция: у = АВх ;
3) дробно-рациональная функция: у = (Ах + В)-1;
4) логарифмическая функция: у = А . ln(х) + В;
5) смешанная функция: у = АхВ.
В зависимости от параметра В она определяет параболическую зависимость (В>0), гиперболическую зависимость (В < 0) и линейную зависимость (В= 0);
6) гиперболическая функция: у = А + В/х;
7) дробно-рациональная функция: у = х/(Ах + В).
Для того, чтобы выбрать теперь вид аналитической зависимости, которая наилучшим образом соответствует исходным экспериментальным данным, поступим следующим образом. Выполним промежуточные вычисления. На области определения независимой переменной (мы ранее условились, что это будет хi ) выберем две точки, достаточно надежные и по возможности как можно дальше отстоящие друг от друга. Обозначим их Х1 и Х2. Этим точкам соответствуют значения Y1 и Y2. Найдем теперь среднее арифметическое, среднее геометрическое и среднее гармоническое для выбранных точек:
Построим график, который, по нашему мнению, наиучшим образом будет соответствовать имеющимся экспериментальным данным. И, зная XAP , XГЕОМ , и ХГАРМ , найдем из графика приближенные Y*АР, Y*ГЕОМ и Y*ГАРМ . При построении графика можно использовать метод построения интерполяционной кривой по выбранным точкам (Дорот и др., 1977; Крылов и др., 1972; Троицкий, Иванова 1975) или методы п.1.
Теперь найдем погрешности результатов сравнений:
| Y*АР- YАР| =1; | Y*АР - YГЕОМ| = 2 ; | Y*АР - YГАРМ| = 3 ;
| Y*ГEOM - YАР | = 4 ; | Y*ГЕОМ - YГЕОМ | = 5 ;
| Y*ГАРМ - YАР | = 6 ; | Y*ГАРМ - YГАРМ | = 7
и выберем = min { 1, 2, ..., 7 }.
1. Если наименьшим среди всех абсолютных значений окажется 1 , то в качестве аналитической зависимости для данных точек будет служить линейная функция вида у = Ах + В.
2. Если наименьшей абсолютной ошибкой является 2, то в качестве эмпирической зависимости следует выбрать показательную функцию у = АВx.
3. В том случае, если наименьшая из абсолютных ошибок есть 3, то искомая эмпирическая зависимость определяется дробно - рациональной функцией вида у = (Ах + В) -1 .
4. Если наименьшая из абсолютных ошибок есть 4 , то хорошим приближением будет служить логарифмическая функция у=А . ln(х) + В.
5. Для случая, когда наименьшей абсолютной ошибкой окажется 5 , то в качестве эмпирической зависимости рекомендуется выбрать смешанную функцию у = АхB .
6. Если наименьшей из абсолютных ошибок окажется 6 , то за искомую зависимость следует выбрать гиперболическую функцию у= А + В/х.
7. И, наконец, в том случае, если наименьшая из всех абсолютных ошибок есть 7, то в качестве зависимости следует выбрать дробно - рациональную функцию вида у = х/(Ах + В).
Для уточнения коэффициентов выбранной аналитической зависимости у = f(х, А, В) воспользуемся, как и в п.1, тремя методами.
Метод выбранных точек
На кривой, которую предварительно построим по множеству экспериментальных точек, выберем две произвольные S1 (х1*, у1*) и S2 (х2*, у2*). Зная вид зависимости f(х,А,В), составим систему
разрешая которую относительно параметров А и В, находим их числовые значения.
Метод средних
В эмпирическую формулу у = f(х, А, В) подставляем последовательно хi и получаем уi , которые будут отклоняться от табличных на e i= уi - f(хi, А, В). Согласно методу средних, надо определить так А и В, чтобы e = 0. Для этого вся совокупность значений pазбивается на две группы так, чтобы алгебраическая сумма уклонений в каждой группе равнялась нулю. Таким образом, для определения параметров А и В имеем
откуда получаем из совместного решения системы значение двух параметров А и В.
Метод наименьших квадратов
Согласно этому методу А и В должны быть определены так, чтобы выполнялось условие минимума функции
В силу необходимого условия экстремума функции находим частные производные функции F по неизвестным коэффициентам А и В и приравниваем их к нулю, откуда получаем систему
из решения которой находим А и В.
Ниже в таблице 1 приводятся явный вид коэффициентов А и В для всех рассматриваемых здесь видов зависимостей.
Таблица 1
Таблица 1
(продолжение)