
- •I. Система обчислень mathcad
- •1. Основи роботи з MathCad
- •1.1. Математичні вирази
- •1.2. Типи даних
- •1.3. Функції
- •1.4. Текстові фрагменти
- •1.5. Графічні області
- •1.6. Побудова фігур, що перетинаються
- •1.7. Приклад
- •1.8. Створення анімаційного кліпу
- •1.9. Вправи для самостійної роботи.
- •1.10. Контрольні запитання
- •2. Розв’язок рівнянь засобами mathcad
- •2.9. Вправи для самостійної роботи.
- •2.10. Контрольні запитання
- •2.1. Чисельний розв’язок нелінійного рівняння
- •2.2. Відсутність збіжності функції root
- •2.3. Рекомендації з використання функції root
- •2.4. Знаходження коренів полінома
- •2.5. Розв’язок систем рівнянь
- •2.6. Розв’язок матричних рівнянь
- •2.7. Наближені розв’язки
- •2.8. Символьний розв’язок рівнянь
- •2.9. Вправи для самостійної роботи.
- •2.10. Контрольні запитання
- •3. Символьні обчислення
- •3.7. Вправи для самостійної роботи
- •3.8. Контрольні запитання
- •3.1. Виділення виразів для символьних обчислень
- •3.2. Символьні операції
- •3.2.1. Операції з виділеними виразами
- •3.2.2. Операції з виділеними змінними
- •3.2.3. Операції з виділеними матрицями
- •3.2.4. Операції перетворення
- •3.3. Стиль представлення результатів обчислень
- •3.4. Приклади символьних операцій у командному режимі
- •3.5. Оператори обчислення границь функцій
- •3.6. Завдання операторів користувача
- •3.7. Вправи для самостійної роботи
- •3.8. Контрольні запитання
- •II. Програмування в системі mathcad
- •1. Програмування з використанням програм-функційMathCad
- •1.1. Опис програми - функції і локального оператору присвоєння
- •1.2. Звернення до програми-функції MathCad
- •1.3. Програмування в програмі-функції лінійних алгоритмів
- •1.4. Програмування в програмі-функції алгоритмів, що розгалужуються
- •1.5. Програмування в програмі-функції циклічних алгоритмів
- •1.5.1. Програмування циклу типу арифметичної прогресії
- •1.5.2. Програмування ітераційних циклів
- •1.6. Можливі використання умовного оператора if
- •1.7. Додаткові оператори програмування циклів у пакеті Mathcad 7 Professional
- •2. Модульне програмування в Mathcad
- •2.1. Модульне програмування одного документу
- •2.2. Модульне програмування у декількох документах MathCad
- •2.1. Модульне програмування одного документу
- •2.2. Модульне програмування у декількох документах MathCad
- •III. Автоматизація статистичних обчислень.
- •1. Вибірки і їхнє представлення
- •1.1. Основні поняття
- •1.1. Основні поняття
- •2.1 Властивості оцінок
- •2.2. Теоретичне порівняння оцінок
- •2.3. Статистичне порівняння оцінок
- •2.4. Завдання для самостійної роботи
- •3. Довірчі межі та інтервали.
- •3.1. Основні положення
- •3.1.1. Побудова довірчих границь і інтервалів
- •3.2. Рівень довіри
- •3.3. Завдання для самостійної роботи
- •4. Методи побудови оцінок
- •4.1. Метод моментів
- •4.2. Метод найбільшої правдоподібності
- •4.3. Властивості оцінок найбільшої правдоподібності
- •4.1. Метод моментів
- •4.2. Метод найбільшої правдоподібності
- •4.3. Властивості оцінок найбільшої правдоподібності
- •5. Граничні теореми
- •5.1. Теорема Бернуллі
- •5.2. Закон великих чисел у формі Чебишева
- •5.3. Реалізація практично достовірної події
- •5.4. Стиск розподілу з ростом числа доданків
- •5.5. Посилений закон великих чисел
- •5.6. Теорема Гливенко основна теорема статистики
- •5.7. Центральна гранична теорема
- •6. Однаково розподілені доданки.
- •6.1. Різнорозподілені доданки
- •6.1. Різнорозподілені доданки
- •Додаток 1 Системні змінні
- •Додаток 2
- •Оператори
- •Показові і логарифмічні функції
- •Функції роботи з частиною числа (округлення й ін.)
2.2. Теоретичне порівняння оцінок
Усі три оцінки незсунені, що можна перевірити методами теорії імовірностей. Обчислимо дисперсії оцінок :
Dâ1 =
D(
)
=
,
Dâ2 =
D(
max xi )
=
,
Dâ3 =
D(x(k) +
x(k+1))
,
Звідси випливає, що â2 — найбільш точна оцінка, а â3 — найменш.
Пояснимо наведені формули для дисперсій .
Перша :
Dâ1 =
=
=
=
.
Друга. Визначимо функцію розподілу статистики max xi :
F(z) P{ max xi <
z} =
P{x1 <
z, ..., xn <
z} =
=
;
щільність розподілу
p(z)
= F(z)
=
,
z[0, a].
Далі
Mâ 2 =
M(
max xi )
=
=
,
Mâ22 =
M
=
,
Dâ2 =
Mâ22 — (Mâ2)2=
.
Третя. Використовуємо
теорему Крамера, відповідно до якої
вибіркова p
-квантиль
має дисперсію, рівну приблизно
,
де xp — p-квантиль, f(x) -
щільність розподілу спостережень
вибірки. У нашому випадку (при n
= 2k)
статистика 0.5 (x(k) +x (k+1) ) m є
вибірковою медіаною (p
= 0.5)
, f(x0.5)
= 1/a , â3=
2m,
і тому
Dâ3=Dm =
=
.
2.3. Статистичне порівняння оцінок
Далеко не завжди вдається аналітично обчислити дисперсію оцінки. Як експериментально визначити, яку з оцінок використовувати? По одній вибірці не можна судити про розкид значень оцінки, оскільки значення усього одне; необхідно мати декілька вибірок, наприклад, k = 20, (чи хоча б 5 10), оцінити розкид значень для кожної оцінки і віддати перевагу тій оцінці (той спосіб оцінювання), для якої розкид менше. Якщо ж вибірка усього одна, то випливає (якщо n досить велике) розбити її випадковим чином на декілька вибірок, і по них порівнювати якість оцінок.
Сформуємо k =20 вибірок з розподілу R[0, a=10] обсягу n для різних n=10, 40, 160 і визначимо розкид оцінок. Характеристиками розкиду значень а1,...,аk оцінки âбудемо вважати розмах w = max ai - min ai і середньоквадратичне відхилення (СКВ)
Sa=
,
.
Як приклад у табл.1 і на мал.1 наведено результати порівняння трьох оцінок.
Таблиця 1. Розкид значень оцінок.
|
|
â1 |
â2 |
â3 |
|
amin |
7.98 |
9.21 |
6.04 |
n = 10 |
amax |
13.80 |
10.98 |
15.69 |
|
w |
5.82 |
1.77 |
9.65 |
|
Sa |
1.51 |
0.53 |
2.35 |
|
amin |
8.59 |
9.77 |
7.02 |
n = 40 |
amax |
11.35 |
10.24 |
12.89 |
|
w |
2.76 |
0.47 |
5.86 |
|
Sa |
0.84 |
0.14 |
1.56 |
|
amin |
9.12 |
9.85 |
8.67 |
n = 160 |
amax |
11.26 |
10.06 |
12.24 |
|
w |
2.14 |
0.21 |
3.57 |
|
Sa |
0.50 |
0.05 |
0.94 |
Порівняння значень розмахів w і СКО Sа для 3 оцінок показує, що оцінка â2(х1, ... , хn) є найбільш точною, а оцінка â3(х1, ... , хn) - найменш.
Наведені результати експериментального порівняння 3 способів обробки спостережень показують наступне.
1. Значення оцінок концентруються в околиці оцінюваного параметра (прояв властивості незсуненості оцінок).
2. З ростом числа спостережень точність (величина розкиду) оцінок поліпшується (прояв властивості спроможності).
3. Різні оцінки розрізняються по величині середньої помилки, звідки ясно, що різні способи обробки спостережень потрібно порівнювати по величині середнього значення деякого критерію якості, наприклад, середнього значення квадрату помилки.