Файл: Гинзбург Г.Л. Руководство по использованию метода наименьших квадратов в экономико-статистических расчетах учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 27.07.2024
Просмотров: 84
Скачиваний: 0
матичеокой школе, основоположником и создателем которой является Чебышев.
Теория нанлучшего приближения функций в настоящее время нашла широкое практическое применение в экономике и технике. Круг конкретных задач, для разрешения которых Чебышев ее создал, значительно расширился, а эффективные практические методы неустанно разрабатываются.
В связи с появлением ЭВМ ученые, работающие в этой области, начинают все чаще и чаще обращаться к работам Чебышева. На праздновании 200-летнего юбилея Академии наук СССР, а также на заседании, посвященном 50-летию со дня смерти Чебышева, многие ученые говорили о настоятель ной необходимости самого полного и тщательного изучения научного .наследия этого выдающегося математика (акаде мики Н. Г. Бруевич, И. И. Артоболевский и др.).
Основная установка Чебышева состоит в следующем: «Практическая деятельность человека представляет чрезвы чайное разнообразие и для удовлетворения ее требований, разумеется, недостает науке многих и различных метод. Но из них особенную важность имеют те, которые необходимы для решения различных видоизменений одной и той же за дачи, общей для всей практической деятельности человека: как располагать средствами своими для достижения по воз можности большей выгоды». [15]
Из этой оснозной установки вытекает лейтмотив всех ра бот Чебышева в области теории наилучшего приближения функций — заменить точное решение задачи приближенным так, чтобы погрешность, получающаяся при этой замене, бы ла наименьшей.
В 1946 году академик В. С. Немчинов изложил систему статистическо-математических вычислений в полиномах Че бышева, позволяющую производить все расчеты на ЭВМ. [10]
Необходимо упомянуть и работу В. И. Хотимского [13], в которой даются весьма удобные таблицы для нахождения уравнений параболических кривых.
Та же задача — упрощения вычислений при использовании метода наименьших квадратов — рассматривалась в работе «О связном анализе статистических рядов», написанной груп пой авторов — бывших студентов статистического отделения МГУ: Боярским, Писаревым, Старовским под руководством профессора Б. С. Ястремского. [19].
Весьма значительными представляются работы советских ученых Н. С. Четверикова — представителя «чупровской» школы и В. Д. Пирятина из Харьковского университета.
М. В. Игнатьев и Н. С. Четвериков составили и опубликова
10
ли очень удобные таблицы множителей для вычисления пара метров параболических уравнений. 1
В 1963 году вышла книга Н. С. Четверикова [17], в кото рой помещена небольшая по объему, но весьма содержатель ная статья «О технике вычисления параболических кривых». В этой статье приводятся те же таблицы множителей.
В. Д. Пирятин дает подробную разработку формул, уп рощающих технику вычислений при использовании метода наименьших квадратов, и приводит общее решение полиномов степени т. [11]
Д. Р. Каустон предлагает таблицы и способ их использо вания для построения полиномами кривых регрессий .по ме тоду наименьших квадратов. Таблицы дают возможность подбирать многочлен до 4-го порядка только при двенадца ти равноотстоящих значениях абсцисс, но тем не менее его работа имеет общий характер. 2
В последнее время в НИИ ЦСУ СССР группа научных сотрудников под руководством Е. А. Александрова вела ра боту эвристического характера с целью испытания способа выравнивания динамических рядов (экспериментальная пси хология). Группе испытуемых студентов были предложены задачи на выравнивание с целью объединить прошлый опыт или привычки в законе гомогенности (опыт в психологическом поле).
Был разработан метод трех параметров, метод ветвящихся процессов и т. д. Полученные кривые математического урав нения не имеют, упрощения в вычислениях при ручном счете
пока нет, но работа на обозримом отрезке |
времени может |
дать положительный результат. |
|
С появлением ЭВМ открылись новые возможности приме |
|
нения методов математической статистики |
в планировании |
и управлении народным хозяйством. |
|
Оказалась возможной значительная формализация вычис лительного процесса. Это обстоятельство намного расширило круг задач, допускающих практическое решение. Однако из этого вовсе не следует вывод об эффективности ЭВМ во всех случаях. Известно, что надежды, возлагавшиеся на внедрение ЭВМ, оправдывались не всегда.
Поэтому к помощи ЭВМ следует прибегать только тогда, когда это экономически обосновано.
§ 3. Постановка задачи
Пусть дана система, состоящая из т функций
___________ |
<P o(M )q?l(*), |
. . . . ф т -1 ( х ) , |
1 М. В. |
Игнатьев и Н. С. |
Четвериков. «Вопросы конъюнк |
туры», М., 1926, т., 2, вып. 1. |
|
2 «Биометрике», Нью-Йорк, 1968, № 2.
11
определенных на отрезке [а, в] и линейно независимых, т. е. таких, что ни одна из этих функций не является линейной комбинацией остальных. Выражение вида
|
|
m—1 |
|
называется обобщенным |
полиномом порядка |
т — 1. Здесь |
|
коэффициенты |
а0) а{. .., |
ai ...am— произвольные числа. |
|
В конкретных экономических исследованиях часто быва |
|||
ет необходимо аппроксимировать функцию f(x), |
заданную на |
||
дискретном множестве точек Х = [х1: х г.......х п\ |
обобщенны |
||
ми полиномами |
(1,1) (очевидно должно выполняться т ^.п). |
||
Наиболее распространенными аналитическими способами |
аппроксимации таких функций являются: 1) наилучшее сред нее приближение (метод наименьших квадратов) и 2) наи
лучшее равномерное приближение |
(метод Чебышева). |
|
§ 4. Наилучшее среднее приближение |
|
|
метод наименьших квадратов |
|
|
Пусть на множестве точек |
Х = {jcx, х 2, ..., х п] |
задана |
функция y = f(x). |
|
от дан- |
За меру отклонения обобщенного полинома (1,1) |
||
ной функции принимается величина |
|
|
|
2 |
|
f{xt)— Pm- 1 (Xi)
;= 1
называемая квадратичным отклонением.
Метод наименьших квадратов заключается в минимиза ции величины е.
В качестве системы { <р,-(л')} на практике используются системы функций двух видов
а) фо(х) = 1, (pi (х) =х, ..., <pm-i (х) = х т~1 ;
б) ортогональные полиномы Чебышева.
Каждая из этих систем применяется в методе наимень ших квадратов.
В случае а) обобщенным полиномом является обычный
полином |
|
|
т ~1 |
(1,2) |
Pm-i (х) = а0 + агх + |
х |
|
||
Для того, чтобы полином (1, 2) |
стал |
оптимальным, |
необхо |
|
димо выбрать такие коэффициенты а0, яь ..., |
ат - п р |
и ко |
||
торых квадратичное отклонение s достигает |
минимума. |
12
Для достижения минимума квадратичного отклонения приравняем к нулю частные производные от величины
e = S | |
— (а0 + а1х + |
|
|
|
|
|
||||
взятые по всем переменным at |
( /= 0,1, ..., |
т — 1). |
|
|
||||||
Получим систему т уравнений с п неизвестными |
|
|
||||||||
|
|
а0, #1> •' •1 &т—1 .• |
|
|
|
|
||||
да„ |
— 2 |
I?(-Ч)—(<r0+ |
a |
i |
X |
') |
= |
0 |
||
— |
=2 V |
^(a ,-)—K + ^ iA + ... + a m_i a" |
') U = 0 |
|||||||
^■1 |
^TTli |
|
|
|
|
|
|
|
|
(1.3) |
|
г=1 |
|
|
|
|
|
|
|
|
|
^am-1 |
= 2 У |
©(a ,) — ( a 0 + |
a aA + |
... + |
a m- i A ' n- |
|) U m~ 1--=0 |
||||
' |
|
|
|
|
|
|
|
|
|
|
Обозначим |
x\ _)_ Ag-j- ... -j- Xln —S; |
|
|
|
|
|||||
|
|
|
|
|
|
|||||
|
x / /(-Ч) 4“ -*2/ (-*2) H~ |
+ Xaf {xn) = |
6;. |
|
|
|||||
|
|
(/ = |
0,1,... , m— l ) . |
|
|
|
|
|||
Тогда система |
(1,3) |
примет вид |
|
|
|
|
|
|||
S qCLq "Т $1а 1”Ь •••~"ЬS m~l О-т- 1= Ьо |
|
|
|
|
||||||
До 4" S,at -}-... -f- Smam-i — b} |
|
|
|
|
(1,4) |
|||||
|
|
|
|
|
|
|
|
|
|
|
S m —1 |
|
-f- . . . |
- ( - S 2 m -2 |
|
----b m |
|
|
|||
\ |
|
|
|
|
|
|
|
|
|
|
Здесь S0= n |
|
|
|
|
|
|
|
|
|
|
В результате |
решения |
этой системы |
получим |
искомые |
||||||
коэффициенты а0 = ао, ai=ai, •••, flm-i = a OT-i |
|
|
{©Да)} |
|||||||
Случай б), когда за |
основную систему функций |
приняты ортогональные полиномы Чебышева, принципиаль но не отличается от разобранного выше случая а), однако, имеет преимущества при машинном и ручном счете.
13