Метод наименьших квадратов
Сущность метода наименьших квадратов состоит в том, что наивероятнейшими значениями аргументов искомой аналитической зависимости будут такие, при которых сумма квадратов отклонений экспериментальных значений функции от значений самой функции , будет наименьшей: . Обычно функция является функцией нескольких аргументов , где - неизвестные коэффициенты многочлена. Тогда на основании экспериментальных пар следует определить искомых аргументов аналитической зависимости, которая наилучшим образом описывает массив , т.е. в этом случае метод наименьших квадратов требует выполнения условия:
.
Применение метода наименьших квадратов при статистической обработке результатов измерений требует учета ряда условий:
§ Значения аргументов известны точно;
§ Результаты измерений независимы и содержат лишь случайные погрешности с одинаковыми дисперсиями;
§ Погрешности измерения имеют нормальное распределение.
Первое условие приближенно выполняется за счет измерения значения с меньшей погрешностью, чем . Наличие только случайных погрешностей обеспечивается исключением из результатов измерений возможных систематических погрешностей.
На основе метода наименьших квадратов можно выполнять аппроксимацию различных аналитических зависимостей, например, выражаемых такими полиномами:
, где a, b, c, d, …, g – константы.
Рассмотрим случай, когда искомая зависимость имеет линейный характер вида .
При использовании метода наименьших квадратов необходимо по набору из экспериментальных координат ( ) найти такие оценки неизвестных постоянных a и b, при которых получится прямая линия, наилучшим образом отражающая истинную анализируемую линию.
График функции – прямая линия с коэффициентом b = tgα, пересекающая ось ординат в точке а.
В соответствии с методом наименьших квадратов наилучшим оценкам a и b соответствуетминимальное значение выражения , где - отклонение измеренных значений от вычисленных при .
Эта величина минимальна, если ее частные производные равны нулю:
и = 0.
Решая систему этих двух уравнений, находим формулы для оценок значений a и b:
При этом:
; ; ; ; .
Степень приближения найденных значений a и b к истинным значениям
Этих величин оценивается с помощью их СКО и :
, ,
где - СКО погрешности измерения величины y, значение которой можно получить из паспортных данных на средство измерения или вычислить по формуле:
.
В качестве примера практического применения метода наименьших квадратов рассмотрим аппроксимацию нагрузочной характеристики одного из устройств преобразовательной техники. Для построения нагрузочной характеристики измеряют 5 ... 10 пар значений выходного напряжения и тока нагрузки . Индекс соответствует текущему измерению (n – число измерений).
В данном примере было снято десять пар (n=10) экспериментальных точек и напряжения и тока соответственно.
i | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
,B | 5,35 | 5,20 | 5,11 | 4,92 | 4,76 | 4.76 | 4,50 | 4,33 | 4,30 | 4,04 |
,mA | 105 | 110 | 129 | 148 | 154 | 181 | 190 | 206 | 225 | 241 |
По полученным экспериментальным результатам построим графическое отображение полученных пар и аппроксимируем их.
Рис. 1 Аппроксимация исследуемой зависимости методом наименьших квадратов
Из расположения экспериментальных точек видно, что аппроксимирующим уравнением может быть полином первой степени типа . Таким образом, на основании массива экспериментальных данных по уравнениям вычисляем коэффициенты a и b:
;
Следовательно, исследуемое устройство имеет нагрузочную характеристику, аналитически описываемую как
ПРИМЕР 2
Требуется установить реальную зависимость сопротивления металлического проводника от температуры по результатам совместных измерений температуры и сопротивления. При этом теоретическая зависимость определена как , где - сопротивление проводника при температуре 0 0С; - температурный коэффициент сопротивления проводника; - температура в 0С.
,0С | 10 | 15 | 20 | 25 |
, Ом | 10,3 | 10,9 | 11,3 | 11,6 |
Преобразуем заданную зависимость в вид , где , .
Расчеты по формулам
; ; ; ;
при п=4, = ti,, и дают следующие результаты: а = 9б52 Ом, b = 0,09 Ом/градус.
Пусть средство измерения имеет СКО Ом. Тогда, проведя вычисления по формулам
, ,
получим: Ом, Ом/ градус.
Окончательно получаем:
Дата добавления: 2015-10-22; просмотров: 36; Опубликованный материал нарушает авторские права?.