- •Понятие о моделировании.
- •2.1. Виды моделирования.
- •3.1. Математическое моделирование. Источники ошибок. Необходимость тестирования.
- •4.1. Сравнительные характеристики пакетов mathcad, matlab, maple, mathematica. Основные приемы работы в mathcad.
- •5.1. Понятие о системах с сосредоточенными и распределенными параметрами.
- •6.1. Сведение системы оду произвольного порядка к системе оду первого порядка в нормализованной форме Коши.
- •7.1. Процедуры решения ду в среде mathcad. Примеры соответствующих документов mathcad.
- •8.1. Последовательность получения математической модели колебательной механической системы с сосредоточенными параметрами.
- •9.1. Численные методы решения оду.
- •10.1. Метод Эйлера и его модификации для решения оду
- •11.1.Методы типа Рунге-Кутта для решения оду
- •13.1. Понятие о граничных задачах для оду. Метод стрельбы для решения граничной задачи.
- •14.1. Построение уравнения изогнутой оси балки при различных условиях на ее концах на основе метода стрельб.
- •15.1.Оценка погрешности решения оду. Способ Рунге для оценки такой погрешности.
- •16.1.Понятие о жестких дифференциальных уравнениях. Процедуры для решения таких уравнений в среде mathcad.
- •17.1. Классификация методов решения оду.
- •19.1. Задача Бюффона как пример использования случайных величин при решении детерминированной задачи.
- •20.1.Вычисление площадей и объемов с использованием случайных величин.
- •21.1.О выборе количества экспериментов для получения заданной степени точности при использовании случайных величин.
- •22.1Метод Судзуо-Какутани для решения граничных задач теории потенциалов
- •23.1. Понятие о конкурирующих стратегиях
- •24.1. Моделирование смо.
- •25.1. Приближение инженерных данных.Виды приближения.
- •Поточечное среднеквадратическое приближение.
- •Непрерывное приближение в среднеквадратичном смысле.
- •Равномерное приближение.
- •1.2.Узловой метод получения математической модели системы.
- •2.2. Метод получения топологических уравнений с использованием м-матрицы.
- •3.2. Получение эквивалентной системы технических объектов.
- •4.2. Компонентные уравнения для различного типа подсистем.
- •Механическая поступательная система.
- •Механическая вращательная подсистема
- •5.2. Основные положения получения математической модели технического объекта на макроуровне.
- •6.2. Особенности выбора узлов сетки при интерполяции различными сплайнами.
- •7.2. Параметрический рациональный сплайн.
- •8.2. Обобщенные кубические сплайны.
- •9.2. Параметрический Эрмитов кубический сплайн.
- •10.2. Интерполяция кривых локальными сплайнами.
- •11.2. Вычисление интеграла по таблице значений функции с использованием интерполирующего полинома Эрмита.
- •12.2. Кубический сплайн дефекта 2 (s3,2(X)).
- •13.2. Определение сплайна. Сплайн первой степени.
- •14.2. Алгоритм «прогонка» для решения системы линейных уравнений с диагональным преобладанием.
- •15. 2. Сплайны. История возникновения. Понятие о звене сплайна.
- •16.2. Использование обратного интерполирования для решения уравнений.
- •17.2. Метод Мюллера для решения трансцендентных уравнений.
- •18.2 Области использования интерполирования.
- •19.2. Тригонометрическая интерполяция.
- •20.2. О наилучшем выборе узлов интерполирования.
- •21.2. Остаточный член интерполяционной формулы Лагранжа.
- •24.2. Разложение аппроксиматора по системе базисных функций.
7.2. Параметрический рациональный сплайн.
Среди всех сплайнов используемых для приближения кривых, рациональный сплайн обладает наиболее универсальными свойствами. Параметрический рациональный сплайн:
Где - суммарная длина хорд, т.о. между точками .
Где задаются одинаковыми для пары ,
Для вычисленных параметров сплайна воспользуемся формулами для рационального сплайна, подставив в их вместо .Для построения рац сплайна нужно задать краевые условия, далее выражают A B через C D, используя условие непрерывности сплайна слева и справа от точки, далее чтобы найти C D используют условие непрерывности производной первой и второй в точках слева и справа, найдя ABCD строим сплайн.
8.2. Обобщенные кубические сплайны.
- это такие сплайны, которые при определенных значениях входящих в них параметров, переходят в обычный кубический сплайн, при этом многие задачи интерполирования решаются с их помощью с лучшими результатами. Одним из достоинств является простота реализации на ЭВМ.
Рациональные сплайны.
- это одна из разновидностей обобщенных кубических сплайнов. Они позволяют более эффективно учитывать точки излома функции, отслеживать большие градиенты, хорошо учитывать выпуклости и вогнутости кривой. Кроме того, позволяет исключить осцилляции на прямолинейных участках.
Пусть на интервале [a,b] задана система узлов ∆ так, что a= x0<x1<…<xn=b.
Рациональный сплайн – функция SR(x), которая
На каждом из подинтервалов [xi, xi+1] имеет вид:
SR(x) = Ai t +Bi (1- t) + + (1)
t =
hi= xi+1 – xi
pi ,qi – параметры сплайна (заданные числа) -1< pi ,qi< ∞;
SR(x) C2 [a,b], т.е.непрерывен вместе со 2-ой производной на [a,b].
Для определения сплайна необходимо найти значения Ai ,Bi ,Ci ,Di .
9.2. Параметрический Эрмитов кубический сплайн.
Этот сплайн представляет собой совокупность двух Эрмитовых сплайнов.
(7)
- производная по S от x(S) в точке Si.
Для возможности вычисления по формуле (7) необходимо определить .
Поскольку в реальных задачах информация о наклонах обычно отсутствует, то, как и в случае обычного Эрмитова сплайна , заменим их приближенными значениями.
Поскольку точное значение параметра вычислить невозможно, то будем строить Эрмитов сплайн близкий к сплайну 7 в некотором смысле.
Во-первых, для описания сплайна введем параметризацию по суммарной длине хорд.
!!!Там, где сверху квадратики, должно быть наверху.
где
Во-вторых, точные значения производных заменим по приближенным разностным формам.
(8)
(9)
(10)
где
Эти формулы используются в том случае, если кривая не замкнута. Если кривая замкнута, то вместо формул 8 и 10 используем:
(11)
Рекомендации по выбору узлов:
1. Следует выбирать узлы так, чтобы (то есть, чтобы длины звеньев были практически одинаковы);
2. В точках излома кривой следует вводить по два близких узла. В этом случае будет снижена асцилляция кривой, заключающаяся в том, что сплайн сильно уклониться от истинной кривой.