
- •Лекция 1.
- •1.1. Основные понятия и определения.
- •1.2. Требования к модели. Функции модели
- •1.3. Классификация моделей
- •1.3.1. Статические и динамические модели
- •1.3.2. Детерминированные и стохастические модели
- •1.3.3. Классификация математических моделей.
- •Лекция 2. Разновидности математических задач, возникающих при моделировании эмс.
- •2.1. Приближение функций. Интерполяция, экстраполяция, аппроксимация. Приближение периодических функций.
- •2.2. Алгебра комплексных чисел.
- •2.3. Решение систем линейных алгебраических уравнений (слау), матричная алгебра.
- •2.4. Решение систем нелинейных алгебраических уравнений.
- •2.5. Решение систем обыкновенных дифференциальных уравнений.
- •2.6. Решение систем дифференциальных уравненийв частных производных.
- •Лекция 3. Ошибки вычислений.
- •3.1. Общие характеристики вычислительных процессов.
- •3.2. Классификация погрешностей.
- •3.3. Абсолютная и относительная погрешности. Точные десятичные знаки.
- •Лекция 4. Приближение функций
- •4.1. Каноническая форма интерполяционного полинома.
- •4.2. Интерполяционный полином Лагранжа.
- •4.3. Интерполяция сплайнами.
- •4.3.1. Линейный сплайн
- •4.3.2. Кубический сплайн
- •Лекция 5. Аппроксимация функций.
- •5.1. Степенной базис
- •5.2. Базис в виде классических ортогональных полиномов
- •5.3. Малая помехоустойчивость метода наименьших квадратов при решении задач идентификации
- •5.3.1. Теория множественности моделей
- •Лекция 6. Приближение периодических функций.
- •6.1. Общие сведения.
- •6.2. Ряды Фурье.
- •6.3. Функции Уолша.
- •Лекция 7. Решение систем линейных алгебраических уравнений.
- •7.1. Область применения слау в задачах математического моделирования эмс.
- •7.2. Прямые методы решения слау.
- •7.3. Итерационные методы.
- •Лекция 8. Решение нелинейных уравнений.
- •8.1. Отделение корней уравнения.
- •8.1.1. Графический метод отделения корней.
- •8.1.2. Аналитический метод отделения корней.
- •8.2. Метод половинного деления (метод дихотомии).
- •8.3. Метод хорд.
- •8.4. Метод касательных (метод Ньютона-Рафсона).
- •Лекция 9. Численное интегрирование и дифференцирование.
- •9.1. Метод прямоугольников.
- •9.2. Метод трапеций.
- •9.3. Метод Симпсона.
- •9.4. Численное дифференцирование.
- •Лекция 10. Решение систем обычных дифференциальных уравнений (оду).
- •10.1. Метод Эйлера.
- •10.2. Методы Рунге-Кутты.
- •10.2.1. Метод Рунге-Кутты-Мерсона
- •10.3. Метод Адамса.
- •10.4. Визуализация решений оду.
- •Лекция 11. Визуальное моделирование динамических систем.
- •Лекция 12. Численное решение систем дифференциальных уравнений в частных производных.
- •12.1. Уравнения математической физики.
- •12.1.1. Уравнения параболического типа.
- •12.1.2. Уравнения гиперболического типа.
- •12.1.3. Уравнения эллиптического типа
- •12.2. Основные понятия метода сеток.
- •Лекция 13. Решение оптимизационных задач.
- •13.1. Методы безусловной одномерной оптимизации
- •13.1.1. Постановка задачи.
- •13.1.2 Метод обратного переменного шага.
- •13.1.3. Метод половинного деления
- •13.1.4. Метод квадратичной аппроксимации (метод Пауэлла).
- •13.2 Методы оптимизации многомерных функций.
- •13.2.1. Метод покоординатного спуска.
- •13.2.2. Метод наискорейшего спуска (метод градиентов)
- •13.2.3. Метод Нелдера-Мида.
- •13.2.3. Метод пчелиного роя.
- •Лекция 14. Идентификация параметров эмс.
- •14.1. Аппроксимация переходных характеристик элементарными динамическими звеньями
- •14.1.1. Апериодическая переходная характеристика
- •14.1.2.Колебательная переходная характеристика.
5.3. Малая помехоустойчивость метода наименьших квадратов при решении задач идентификации
Изменение
среднеквадратичной ошибки будет
следующим: для уравнения первой степени
δ2
= = 25%, второй— δ2
= 15%, третьей — δ2
= 0%, четвертой δ2
= 0% и т. д. (показано на рис. 5.1 пунктиром).
Таким образом, при полном отсутствии
помех метод позволяет «открыть» истинное
уравнение. Легко убедиться, что при
наложении на исходные данные очень
небольшого шума (Δy
= ±0,01) получается кривая, показанная
сплошной линией (рис. 5.1). Кривая
подтверждает правило «чем сложнее
модель, тем она точнее».
При наличии даже очень малых неточностей измерений (помех) с помощью критерия Δ2(N) нельзя отличить слишком сложные (переусложненные) модели от модели оптимальной сложности, отвечающей действительному процессу. Задачу открытия единственной модели оптимальной сложности решают методы самоорганизации моделей, основанные на помехоустойчивых, «внешних» критериях. Как увидим ниже, для экстраполяции (прогнозирования) метод наименьших квадратов тоже не пригоден.
Это не исключает эффективности метода наименьших квадратов при решении задачи приближенной аппроксимации функций внутри области интерполяции, где достаточно густо заданы опытные точки. Собственно говоря, метод и был задуман Гауссом для решения задачи интерполяции.
5.3.1. Теория множественности моделей
Теория множественности моделей утверждает, что по экспериментальным данным принципиально нельзя найти единственную модель. Например, в качестве полинома регрессии можно взять полином любого вида и любой степени, и для каждого из них регрессионный анализ укажет значения коэффициентов. В любом достаточно сложном уравнении подбираются оценки коэффициентов так, чтобы ошибка на интервале наблюдения (интерполяции) была мала или даже равна нулю. Например, если задано пять экспериментальных точек, то уравнение, имеющее пять коэффициентов (или параметров), дает кривую, точно проходящую через все заданные точки. В этом и состоит теория множественности моделей. Отсюда следует, что для каждого объекта, рассматриваемого как некоторый «черный ящик», можно составить не одну единственную, а бесконечное множество моделей, имеющих одинаковые или почти одинаковые внешние проявления. Так, ясно, что не раскрыв электрический выключатель, мы никаким способом не можем установить, сколько разрывов электрической цепи он выполняет: один, два или три.
По этому поводу У. Р. Эшби пишет: «Время простых моделей прошло... С тех пор, как было доказано (Шенноном в 1938 г. для релейных сетей и Мак-Каллоком и Питтсом в 1943 г. для сетей из простых нейроноподобных элементов), что любое поведение машины, поддающейся описанию, может быть реализовано не одной машиной, а бесконечным рядом их, разработка новой модели показывает только, что ее автор не усвоил значения работ Шеннона, Мак-Каллойа и Питтса» [89, 90].
Во многих научных работах еще до сих пор авторы указывают, что их теория или формула дала кривую, которая весьма близко прошла через экспериментальные точки. Такого рода «доказательства» правильности теорий нельзя принимать серьезно.
Решение вопроса о выборе единственного уравнения регрессии оптимальной сложности дает принцип внешнего дополнения.
Только внешнее дополнение может служить арбитром при сравнении противоречащих теорий. Ошибка, измеренная на всех экспериментальных точках, не является внешним дополнением.