Metodichka_lab_Ekonometria
.pdfРаздел 1. Элементы линейной алгебры и математической статистики – математическая основа курса «Эконометрия». Функции программы MS Excel выполнения операций с матрицами и категории
«Статистические»
1.1 Теоретические замечания
При решении ряда задач в курсе «Эконометрия» используются операции с матрицами, которые выполняются согласно правилам линейной алгебры. Реализация этих операций в программной среде MS Excel имеет свои особенности. Рассмотрим функции MS Excel и порядок их применения для выполнения операций с матрицами.
1.2 Выполнение операций с матрицами с помощью функций MS
Excel
Порядок выполнения операций с матрицами в табличном редакторе
MS Excel
1. Выделить диапазон ячеек под результат операции. Размерность диапазона определяется по правилам выполнения операций с матрицами
− если матрица A размерности n × m , то матрица AT будет иметь размерность m × n ;
− если матрица A размерности n × m , матрица B размерности m ×l , то произведение матриц A B будет иметь размерность n ×l ;
− если матрица A размерности n × n то матрица A−1 (если она существует) также будет иметь размерность n × n .
2.Вызвать необходимую функцию (табл. 1.1) или набрать ее в командной строке.
3.Ввести диапазоны ячеек (или сделать ссылки на ячейки), в которых хранятся преобразуемые матрицы.
4.Удерживая клавиши Cntrl + Shift, нажать клавишу Enter.
|
|
|
|
Таблица 1.1 |
|
Функции программы MS Excel для работы с матрицами |
|||||
|
|
|
|
|
|
Операция |
|
Функции в MS Excel |
|
|
|
Транспонирование |
ТРАНСП (массив) (категория «Ссылки и мас- |
|
|||
матрицы |
сивы») |
|
|
|
|
Умножение матриц |
МУМНОЖ |
(массив 1, |
массив 2) |
(категория |
|
«Математические») |
|
|
|
||
|
|
|
|
||
Вычисление |
МОПРЕД |
(массив) |
(категория |
«Матема- |
|
определителя |
|
||||
тические») |
|
|
|
|
|
матрицы |
|
|
|
|
|
|
|
|
|
|
|
Нахождение обратной |
МОБР (массив) (категория «Математические») |
|
|||
матрицы |
|
|
|
|
|
9
Рассмотрим пример выполнения лабораторной работы 1 для варианта 40.
Исходная матрица A размерности 3 × 4 при m = 40 примет вид
41 42 43 44 A = 44 45 46 47 .
47 48 49 49
Для выполнения лабораторной работы нужно создать новую рабочую книгу MS Excel, в ячейки которой сначала вводятся исходные данные задачи.
1. Введем элементы матрицы A в ячейки листа MS Excel (ячейки B6 : E8 рис. 1.1). Отметим, что ее можно ввести непосредственно или получить для своего варианта, используя вспомогательную матрицу, которая хранится в ячейках C2 : F 4 (рис. 1.1).
Рис. 1.1 – Организация данных и выполнение операции транспонирования средствами программы MS Excel
Для выполнения операции транспонирования выделим диапазон
4 ×3 в ячейках G5: I8 (рис. 1.1) для транспонированной матрицы AT . Вызовем надстройку «Мастер функций» кликнув на иконку fx или
выбрав в меню «Вставка» подменю «Функция». В результате появится окно как на рис. 1.2 слева. В поле «Категория» нужно выбрать «Ссылки и массивы» или «10 недавно использовавшихся», а затем в списке функций выбрать функцию «ТРАНСП» и нажать «ОК». В результате появится новое окно как на рис. 1.2 справа. В поле «Массив» нужно ввести ссылку на ячейки, в которых хранится матрица A . Сделать это можно двумя способами: первый – непосредственно с клавиатуры набрать адреса левой верхней и правой нижней ячейки, разделенные двоеточием; второй – выделив мышью ячейки, в которых хранится матрица A, курсор при этом должен находиться в поле «Массив».
10
Рис. 1.2 – Последовательность выполнения операции транспонирования
После нажатия клавиш Cntrl+Shift+Enter в ячейках G5: I8 (рис. 1.1) появится результат транспонирования матрицы A
|
41 |
44 |
47 |
|
AT = |
42 |
45 |
48 |
. |
43 |
46 |
49 |
|
|
|
44 |
47 |
49 |
|
Пункт 1 задания 1 выполнен полностью.
2. Для выполнения второго задания – нахождения матрицы B = AAT необходимо воспользоваться функцией «МУМНОЖ» (массив 1, массив 2) (категория «Математические») (табл. 1.1), порядок работы с которой такой же, как и с функцией «ТРАНСП» (см. задание 1, лаб. раб. 1). Чтобы задать массив 1 нужно сделать ссылку на ячейки B6 : E8, в которых хранится матрица A, массив 2 – на ячейки G5 : I8 , в которых хранится
матрица AT . После выполнения последовательности действий как показано на рис. 1.3, найдем матрицу B
|
7230 |
7740 |
8206 |
|
B = |
7740 |
8286 |
8785 |
. |
|
8206 |
8785 |
9315 |
|
Пункт 2 задания 1 выполнен полностью.
11
Рис. 1.3 – Организация данных и выполнение операции умножения матриц средствами программы MS Excel
3. Для выполнения третьего задания – вычисления определителя матрицы B воспользуемся функцией «МОПРЕД» (массив) (табл. 1.1). Поставив курсор в ячейку B17 (рис. 1.4), в которой будет храниться значение определителя, необходимо вызвать надстройку «Мастер функций» а затем в списке функций выбрать функцию «МОПРЕД» и в поле «Массив» ввести ссылку на ячейки B11: D13, в которых хранится матрица B , нажать «ОК». После выполнения этих действий в ячейке B17 получим результат: det B = 54. Пункт 3 задания 1 выполнен полностью.
Рис. 1.4 – Организация данных и вычисление определителя матрицы средствами программы MS Excel
4. В четвертом задании для нахождения матрицы обратной к матрице
Bвоспользуемся функцией «МОБР» (массив) (табл. 1.1), порядок работы
скоторой такой же, как и с функцией «ТРАНСП» (см. задание 1, лаб. раб. 1). В поле «Массив» нужно ввести ссылку на ячейки B11: D13, в которых хранится матрица B .
12
Рис. 1.5 – Организация данных и нахождение обратной матрицы средствами программы MS Excel
После выполнения этих действий в ячейках G11: I13 получим результат
145,65 –155,37 18,22
B −1 = –155,37 166,93 –20,56 . 18,22 –20,56 3,33
Пункт 4 задания 1 выполнен полностью.
5. В пятом задании, чтобы найти результат умножения матрицы B−1 , которая хранится в ячейках G11: I13 (рис. 1.4) на число m = 40 , которое хранится в ячейке B3, по правилам умножения матрицы на число, нужно
каждый элемент матрицы B−1 умножить на это число. Для этого в ячейки
G16 : I18, в которых будут хранится элементы матрицы m B−1 , нужно ввести формулы B3 G11, … B3 I18. Как видно эти формулы отличаются только вторым множителем. Поэтому формулу можно ввести один раз, а затем скопировать ее, сделав абсолютную ссылку на ячейку B3. Для этого нужно поставить курсор в ячейку G16 , содержащую формулу $B$3 G11 и протянуть ее за правый нижний угол вправо до ячейки I16 , а затем вниз до ячейки I18 . После этих действий экранная форма решения пункта 5 задания 1 примет вид, такой как на рис. 1.6.
Задание 1 выполнено полностью.
Рис. 1.6 – Организация данных и выполнение операции умножения матрицы на число средствами программы MS Excel
13
1.3 Обработка наблюдений статистических признаков
Рассмотрим пример выполнения задания 2 для варианта 40. Исходные данные при m = 40 следующие (табл. 1.2):
Таблица 1.2
Наблюдения статистических признаков X и Y
X |
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
49 |
50 |
Y |
8 |
12 |
12 |
12 |
16 |
16 |
16 |
20 |
20 |
24 |
Для выполнения задания 2 нужно создать или открыть новый рабочий лист в текущей рабочей книге MS Excel, в ячейки которого сначала нужно ввести названия строк и столбцов из табл. 1.2, а затем числовые значения статистических признаков X и Y . В результате получим экранную форму как на рисунке 1.7.
1. Сначала вычислим числовые характеристики статистических признаков X и Y в ячейках B13: B18, C13: C16 с помощью встроенных функций MS Excel (табл. 1.3), которые хранятся в надстройке «Мастер функций» в категории «Статистические».
10
Для ячейки B13 формула имеет вид ∑xi , а ее реализация в ячейке:
i=1
СУММ(B3:B12). В ячейку C13 нужно ввести формулу СУММ(С3:С12). Это можно сделать, протягивая формулу из ячейки B13 за правый нижний угол вправо на одну ячейку.
1 10
Для ячейки B14 формула имеет вид x = 10 i∑=1xi , а ее реализация в
ячейке: СРЗНАЧ(B3:B12) (в ячейках B3:B12 хранятся значения признака
X ).
|
|
|
1 |
10 |
(x − x)2 |
|
|
Для ячейки |
B15 |
формула имеет вид DX = |
∑ |
, а ее |
|||
|
|||||||
|
|
|
10 i =1 |
i |
|
||
реализация в ячейке: ДИСПР(B3:B12).
Для ячейки B16 формула имеет вид σ(X ) =
DX , а ее реализация в ячейке: СТАНДОТКЛОНП(B3:B12).
|
|
|
1 |
10 |
|
|
|
|
|
|||
Для ячейки B17 формула имеет вид cov(X ,Y ) = |
|
|
|
|
∑(x |
i |
− x)( y |
j |
− y) , |
|||
|
|
|
|
|||||||||
|
|
|
10 i =1 |
|
|
|
||||||
а ее реализация в ячейке: КОВАР(B3:B12;C3:C12). |
|
|
|
|
|
|||||||
|
|
|
|
cov(X ,Y ) |
|
|
|
|||||
Для ячейки |
B18 |
формула имеет вид r |
= |
|
|
, |
|
а ее |
||||
σ(X ) σ(Y ) |
|
|||||||||||
|
|
XY |
|
|
|
|
|
|||||
реализация в ячейке: КОРРЕЛ(B3:B12;C3:C12).
14
Рис. 1.7 – Организация данных задания 2 на рабочем листе MS Excel
Таблица 1.3
Функции программы MS Excel для расчета числовых характеристик статистических признаков
Характери- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Функции в MS |
стики стати- |
|
|
Расчетная формула |
|
|||||||||||
|
|
|
Excel (категория |
||||||||||||
стических |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
«Статистические») |
|
признаков |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
n |
|
|
|
|
СРЗНАЧ (число1; |
|
Среднее |
|
|
|
x = |
|
|
|
x |
|
|
|
число2; ...) |
|||
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
n i∑=1 |
|
i |
|
||||
|
|
|
n |
(xi − x)2 = |
|
n |
2 −(x)2 |
|
|||||||
Дисперсия |
DX = |
1 |
∑ |
1 ∑xi |
ДИСПР (число1; |
||||||||||
|
|
n i=1 |
|
|
|
|
|
|
|
n i=1 |
|
число2; ...) |
|||
Стандартное |
|
|
|
σ(X ) = |
DX |
|
СТАНДОТКЛОНП |
||||||||
отклонение |
|
|
|
|
(число1; число2; ...) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Ковариация |
cov(X ,Y ) = |
1 |
|
n |
|
− x)(y j − y) |
КОВАР (массив1; |
||||||||
|
|
∑(xi |
массив2) |
||||||||||||
|
|
|
|
|
n i =1 |
|
|
|
|
|
|||||
Коэффициент |
|
|
r XY = |
|
|
cov(X ,Y ) |
|
|
КОРРЕЛ (массив1; |
||||||
корреляции |
|
|
σ(X ) σ(Y ) |
|
массив2) |
||||||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
15 |
В ячейках C13: C16 для расчета числовых характеристик статистического признака Y формулы такие же, как и в соответствующих ячейках B13: B16, только вместо ссылки на массив B13: B12 должна быть ссылка на массив C3: C12, в котором хранятся значения признака Y . Поэтому эти формулы можно скопировать в соответствующих ячейках B13: B16, протягивая за правый нижний угол вправо на одну ячейку.
Далее выполним расчеты в ячейках D3 : G12 . Сначала введем формулы в ячейки D3 : G3 как указано в табл. 1.4. Чтобы вычислить значения в ячейках D4 : G12 достаточно формулы в 3–й строке скопировать в эти ячейки, протягивая каждую формулу за правый нижний угол вниз до 12–й строки включительно.
|
|
|
|
Таблица 1.4 |
|
Реализация в MS Excel формул задания 2 |
|||||
|
|
|
|
|
|
Адрес ячейки |
Формула |
Реализация в MS Excel |
|||
D3 |
X 2 |
= B3^2 |
|
||
E3 |
Y 2 |
= C3^2 |
|
||
F3 |
XY |
= B3 * C3 |
|
||
G3 |
X − |
|
|
= B3 − $B$14 |
|
X |
|
||||
Осталось заполнить ячейки D13: G13 . Для этого можно скопировать формулу СУММ(С3:С12) из ячейки C13 , протягивая ее за правый нижний угол вправо до ячейки G13.
В результате получим экранную форму как на рисунке 1.8.
2. Для построения корреляционного поля показателей X и Y воспользуемся надстройкой «Мастер диаграмм», которую можно вызвать, кликнув на соответствующую иконку. После этого появится окно как показано слева на рисунке 1.9. Далее в категории «Стандартные» нужно выбрать тип диаграммы «Точечная» и «Далее». В новом окне (правое окно на рис. 1.9), в поле «Диапазон» нужно ввести ссылку на массив $B$3:$C$12, в котором хранятся значения признаков X и Y , а также указать, что ряды данных в столбцах. После нажатия кнопки «ОК» будет построено корреляционное поле.
Для подписи осей диаграммы, на поле диаграммы нужно кликнуть правой кнопкой мыши и в появившемся окне выбрать функцию «Параметры диаграммы». После этого появится окно как на рисунке 1.10, в поля которого нужно ввести название диаграммы и осей. Для редактирования осей нужно кликнуть на оси правой кнопкой мыши и выбрать функцию «Формат оси» и в появившемся окне задать вид, цену деления начало и конец оси. В результате этих операций корреляционное поле признаков X и Y примет такой вид как на рисунке 1.11. Задание 2 выполнено полностью.
16
Рис. 1.8 – Результат выполнения задания 2 на рабочем листе MS Excel
Рис. 1.9 – Построение корреляционного поля средствами MS Excel
17
Рис. 1.10 – Задание параметров корреляционного поля средствами MS |
|||||||
|
|
|
|
Excel |
|
|
|
|
30 |
|
|
|
|
|
|
|
25 |
|
|
|
|
|
|
|
20 |
|
|
|
|
|
|
Y |
15 |
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
40 |
42 |
44 |
46 |
48 |
50 |
52 |
|
|
|
|
X |
|
|
|
|
|
Рис.1.Корреляционное поле |
|
|
|||
Рис. 1.11 – Результат построения корреляционного поля средствами MS
Excel
1.4 Нахождение значений статистических критериев
Для выполнения задания 3 нужно создать или открыть новый рабочий лист в текущей рабочей книге MS Excel, в ячейки которого сначала нужно ввести уровни значимости α1 = 0,05 и α1 = 0,01 (ячейки B3
и E3, рис. 1.12) и номер варианта m = 40 (ячейка H 3). В ячейки B5 и E5 ввести формулы H3+1 и 50–H3 для расчета степеней свободы по формулам k1 = m +1 и k2 =50 − m , соответственно. Для нахождения значений
18
