Не отобразилась форма расчета стоимости? Переходи по ссылке

Не отобразилась форма расчета стоимости? Переходи по ссылке

Магистерская диссертация на тему «Методы для решения краевых задач, в том числе ‘жестких’ краевых задач»

На примере системы дифференциальных уравнений цилиндрической оболочки ракеты – системы обыкновенных дифференциальных уравнений 8-го порядка (после разделения частных производных).

Написание диссера за 10 дней

1. Введение.

На примере системы дифференциальных уравнений цилиндрической оболочки ракеты – системы обыкновенных дифференциальных уравнений 8-го порядка (после разделения частных производных).

Система линейных обыкновенных дифференциальных уравнений имеет вид:

Y(x) = A(x) ∙ Y(x) + F(x),

где Y(x) – искомая вектор-функция задачи размерности 8х1, Y(x) – производная искомой вектор-функции размерности 8х1, A(x) – квадратная матрица коэффициентов дифференциального уравнения размерности 8х8, F(x) – вектор-функция внешнего воздействия на систему размерности 8х1.

Здесь и далее вектора обозначаем жирным шрифтом вместо черточек над буквами

Краевые условия имеют вид:

U∙Y(0) = u,

V∙Y(1) = v,

где

Y(0) – значение искомой вектор-функции на левом крае х=0 размерности 8х1, U – прямоугольная горизонтальная матрица коэффициентов краевых условий левого края размерности 4х8, u – вектор внешних воздействий на левый край размерности 4х1,

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Y(1) – значение искомой вектор-функции на правом крае х=1 размерности 8х1, V – прямоугольная горизонтальная матрица коэффициентов краевых условий правого края размерности 4х8, v – вектор внешних воздействий на правый край размерности 4х1.

В случае, когда система дифференциальных уравнений имеет матрицу с постоянными коэффициентами A=const, решение задачи Коши имеет вид [Гантмахер]:

Y(x) = e∙ Y(x)  +  e∙ e∙ F(t) dt,

где

e= E + A(x-x) + A (x-x)/2! + A (x-x)/3! + …,

где E это единичная матрица.

Матричная экспонента ещё может называться матрицей Коши или матрициантом и может обозначаться в виде:

K(x←x) = K(x — x) = e.

Тогда решение задачи Коши может быть записано в виде:

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)  ,

где Y*(x←x) = e∙ e∙ F(t) dt   это вектор частного решения неоднородной системы дифференциальных уравнений.

2. Случай переменных коэффициентов.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Этот вариант рассмотрения переменных коэффициентов проверялся в кандидатской диссертации.

Из теории матриц [Гантмахер] известно свойство перемножаемости матричных экспонент (матриц Коши):

e= e∙ e ∙ … ∙ e ∙ e,

K(x←x) = K(x←x) ∙ K(x←x) ∙ … ∙ K(x←x) ∙ K(x←x).

В случае, когда система дифференциальных уравнений имеет матрицу с переменными коэффициентами A=A(x), решение задачи Коши предлагается искать при помощи свойства перемножаемости матриц Коши. То есть интервал интегрирования разбивается на малые участки и на малых участках матрицы Коши приближенно вычисляются по формуле для постоянной матрицы в экспоненте. А затем матрицы Коши, вычисленные на малых участках, перемножаются:

K(x←x) = K(x←x) ∙ K(x←x) ∙ … ∙ K(x←x) ∙ K(x←x),

где матрицы Коши приближенно вычисляются по формуле:

K(x←x) = e,      где ∆x= x- x.

3. Формула для вычисления вектора частного решения неоднородной системы дифференциальных уравнений.

Эта очень простая формула еще не обсчитана на компьютерах. Вместо неё обсчитывалась значительно ранее выведенная и гораздо более сложная формула, приведенная в:

Численный метод переноса краевых условий для жестких дифференциальных уравнений строительной механики
Журнал «ММ», Том: 14 (2002), Номер: 9, 3 стр. 1409-003r.pdf

Вместо формулы для вычисления вектора частного решения неоднородной системы дифференциальных уравнений в виде [Гантмахер]:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Y*(x←x) = e∙ e∙ F(t) dt

предлагается использовать следующую формулу для каждого отдельного участка интервала интегрирования:

Y*(x←x) = Y*(x- x) = K(x- x) ∙K(x- t) ∙ F(t) dt .

Правильность приведенной формулы подтверждается следующим:

Y*(x- x) = e∙e∙ F(t) dt ,

Y*(x- x) = e∙e∙ F(t) dt ,

Y*(x- x) = e∙ F(t) dt ,

Y*(x- x) = e∙ F(t) dt ,

Y*(x- x) = e∙ e∙ F(t) dt ,

Y*(x←x) = e∙ e∙ F(t) dt,

что и требовалось подтвердить.

Вычисление вектора частного решения системы дифференциальных уравнений производиться при помощи представления матрицы Коши под знаком интеграла в виде ряда и интегрирования этого ряда поэлементно:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

Y*(x←x) = Y*(x- x) = K(x- x) ∙K(x- t) ∙ F(t) dt =

= K(x- x) ∙ (E + A(x- t) + A (x- t)/2! + … ) ∙ F(t) dt =

= K(x- x) ∙ (EF(t) dt  + A∙(x- t) ∙ F(t) dt  + A/2! ∙(x- t) ∙ F(t) dt  + … ) .

Эта формула справедлива для случая системы дифференциальных уравнений с постоянной матрицей коэффициентов A=const.

Вектор F(t) может рассматриваться на участке (x- x) приближенно в виде постоянной величины F(х)=constant, что позволят вынести  его из под знака интеграла, что приводит к совсем простому ряду для вычислений на рассматриваемом участке.

Для случая дифференциальных уравнений с переменными коэффициентами в приведенной выше формуле для каждого участка может использоваться осредненная матрица А: A=А(х)  коэффициентов системы дифференциальных уравнений.

Приведем (итерационные или рекуррентные) формулы вычисления вектора частного решения, например, Y*(x←x)  на рассматриваемом участке (x←x) через вектора частного решения Y*(x←x), Y*(x←x), Y*(x←x) соответствующих подучастков (x←x), (x←x), (x←x).

Имеем:

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)  ,

Также имеем формулу для отдельного подучасточка:

Y*(x←x) = Y*(x- x) = K(x- x) ∙K(x- t) ∙ F(t) dt.

Можем записать:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)  ,

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)  .

Подставим Y(x) в Y(x) и получим:

Y(x) = K(x←x) ∙ [K(x←x) ∙ Y(x)  +  Y*(x←x)]  +  Y*(x←x) =

= K(x←x) ∙ K(x←x) ∙ Y(x)  +  K(x←x) ∙ Y*(x←x) +  Y*(x←x).

Сравним полученное выражение с формулой:

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)

и получим, очевидно, что K(x←x) = K(x←x) ∙ K(x←x) и, самое главное здесь — для частного вектора получаем формулу:

Y*(x←x)  = K(x←x) ∙ Y*(x←x) +  Y*(x←x).

То есть вектора подучастков Y*(x←x) и  Y*(x←x) не просто складываются друг с другом, а с участием матриц Коши подучастков.

Аналогично запишем:

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

И подставим сюда формулу для Y(x) и получим:

Y(x) = K(x←x) ∙ [K(x←x) ∙ K(x←x) ∙ Y(x)  +  K(x←x) ∙ Y*(x←x) +  Y*(x←x)]  +  Y*(x←x) =

= K(x←x) ∙ K(x←x) ∙ K(x←x) ∙ Y(x)  +  K(x←x) ∙ K(x←x) ∙ Y*(x←x) +  K(x←x) ∙ Y*(x←x)  +  Y*(x←x).

Сравнив полученное выражение с формулой:

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)

очевидно, получаем, что K(x←x) = K(x←x) ∙ K(x←x) ∙ K(x←x) и вместе с этим получаем формулу для частного вектора:

Y*(x←x)  = K(x←x) ∙ K(x←x) ∙ Y*(x←x) +  K(x←x) ∙ Y*(x←x)  +  Y*(x←x).

То есть именно так (по своеобразным рекуррентным формулам) и вычисляется частный вектор – вектор частного решения неоднородной системы дифференциальных уравнений, то есть так вычисляется, например, частный вектор Y*(x←x)  всего участка (x←x)  на основе вычисленных векторов Y*(x←x), Y*(x←x), Y*(x←x) подучастков (x←x), (x←x), (x←x).

4. Метод «переноса краевых условий» в произвольную точку интервала интегрирования.

Метод обсчитан на компьютерах. По нему уже сделано 3 кандидатских физ-мат диссертации.

Метод подходит для любых краевых задач. А для «жестких» краевых задач показано, что метод считает быстрее, чем метод С.К.Годунова до 2-х порядков (в 100 раз), а для некоторых «жестких» краевых задач не требует ортонормирования вовсе. Смотри:

Численный метод переноса краевых условий для жестких дифференциальных уравнений строительной механики
Журнал «ММ», Том: 14 (2002), Номер: 9, 3 стр. 1409-003r.pdf

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Полное решение системы дифференциальных уравнений имеет вид

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)  .

Или можно записать:

Y(0) = K(0←x) ∙ Y(x)  +  Y*(0←x)  .

Подставляем это выражение для Y(0) в краевые условия левого края и получаем:

U∙Y(0) = u,

U∙[ K(0←x) ∙ Y(x)  +  Y*(0←x) ] = u,

[ U∙ K(0←x) ] ∙ Y(x)  = u U∙Y*(0←x)  .

Или получаем краевые условия, перенесенные в точку x:

U∙ Y(x)  = u  ,

где U= [ U∙ K(0←x) ] и u = u U∙Y*(0←x) .

Далее запишем аналогично

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)

И подставим это выражение для Y(x) в перенесенные краевые условия точки x

U∙ Y(x)  = u,

U∙ [ K(x←x) ∙ Y(x)  +  Y*(x←x) ]  = u  ,

[ U∙ K(x←x) ] ∙ Y(x)  = u — U∙ Y*(x←x)   ,

Или получаем краевые условия, перенесенные в точку x:

U∙ Y(x)  = u  ,

где U= [ U∙ K(x←x) ] и u = u — U∙ Y*(x←x)   .

И так в точку x переносим матричное краевое условие с левого края и таким же образом переносим матричное краевое условие с правого края и получаем:

U∙ Y(x)  = u  ,

V∙ Y(x)  = v  .

Из этих двух матричных уравнений с прямоугольными горизонтальными матрицами коэффициентов очевидно получаем одну систему линейных алгебраических уравнений с квадратной матрицей коэффициентов:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Y(x) =  .

А в случае «жестких» дифференциальных уравнений предлагается применять построчное ортонормирование матричных краевых условий в процессе их переноса в рассматриваемую точку. Для этого формулы ортонормирования систем линейных алгебраических уравнений можно взять в [Березин, Жидков].

То есть, получив

U∙ Y(x)  = u,

применяем к этой группе линейных алгебраических уравнений построчное ортонормирование и получаем эквивалентное матричное краевое условие:

U∙ Y(x)  = u.

И теперь уже в это проортонормированное построчно уравнение подставляем

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x)  .

И получаем

U∙ [ K(x←x) ∙ Y(x)  +  Y*(x←x) ]  = u  ,

[ U∙ K(x←x) ] ∙ Y(x)  = u — U∙ Y*(x←x)   ,

Или получаем краевые условия, перенесенные в точку x:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

U∙ Y(x)  = u  ,

где U= [ U∙ K(x←x) ] и u = u — U∙ Y*(x←x)   .

Теперь уже к этой группе линейных алгебраических уравнений применяем построчное ортонормирование и получаем эквивалентное матричное краевое условие:

U∙ Y(x)  = u.

И так далее.

И аналогично поступаем с промежуточными матричными краевыми условиями, переносимыми с правого края в рассматриваемую точку.

В итоге получаем систему линейных алгебраических уравнений с квадратной матрицей коэффициентов, состоящую из двух независимо друг от друга поэтапно проортонормированных матричных краевых условий, которая решается любым известным методом для получения решения Y(x) в рассматриваемой точке x:

Y(x) =  .

5. Второй вариант метода «переноса краевых условий» в произвольную точку интервала интегрирования.

Этот вариант метода еще не обсчитан на компьютерах.

Предложено выполнять интегрирование по формулам теории матриц [Гантмахер] сразу от некоторой внутренней точки интервала интегрирования к краям:

Y(0) = K(0←x) ∙ Y(x) +  Y*(0←x)  ,

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Y(1) = K(1←x) ∙ Y(x) +  Y*(1←x)  .

Подставим эти формулы в краевые условия и получим:

U∙Y(0) = u,

U∙[ K(0←x) ∙ Y(x)  +  Y*(0←x) ] = u,

[ U∙ K(0←x) ] ∙ Y(x)  = u U∙Y*(0←x)  .

и

V∙Y(1) = v,

V∙[ K(1←x) ∙ Y(x)  +  Y*(1←x) ] = v,

[ V∙ K(1←x) ] ∙ Y(x)  = v V∙Y*(1←x)  .

То есть получаем два матричных уравнения краевых условий, перенесенные в рассматриваемую точку x:

[ U∙ K(0←x) ] ∙ Y(x)  = u U∙Y*(0←x)  ,

[ V∙ K(1←x) ] ∙ Y(x)  = v V∙Y*(1←x)  .

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

Эти уравнения аналогично объединяются в одну систему линейных алгебраических уравнений с квадратной матрицей коэффициентов для нахождения решения Y(x)  в любой рассматриваемой точке x:

Y(x)  =  .

В случае «жестких» дифференциальных уравнений предлагается следующий алгоритм.

Используем свойство перемножаемости матриц Коши:

K(x←x) = K(x←x) ∙ K(x←x) ∙ … ∙ K(x←x) ∙ K(x←x)

и запишем выражения для матриц Коши, например, в виде:

K(0←x) = K(0←x) ∙ K(x←x) ∙ K(x←x),

K(1←x) = K(1←x) ∙ K(x←x) ∙ K(x←x) ∙ K(x←x),

Тогда перенесенные краевые условия можно записать в виде:

[ U∙ K(0←x) ∙ K(x←x) ∙ K(x←x) ] ∙ Y(x)  = u U∙Y*(0←x)  ,

[ V∙ K(1←x) ∙ K(x←x) ∙ K(x←x) ∙ K(x←x) ] ∙ Y(x)  = v V∙Y*(1←x)

или в виде:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

[ U∙ K(0←x) ∙ K(x←x) ∙ K(x←x) ] ∙ Y(x)  = u*  ,

[ V∙ K(1←x) ∙ K(x←x) ∙ K(x←x) ∙ K(x←x) ] ∙ Y(x)  = v*  .

Тогда рассмотрим левое перенесенное краевое условие:

[ U∙ K(0←x) ∙ K(x←x) ∙ K(x←x) ] ∙ Y(x)  = u*  ,

[ U∙ K(0←x) ] ∙ { K(x←x) ∙ K(x←x) ∙ Y(x) }  = u*  ,

[     матрица      ] ∙ {                     вектор                        }  = вектор  .

Эту группу линейных алгебраических уравнений можно подвергнуть построчному ортонормированию, которое сделает строчки [матрицы] ортонормированными, {вектор} затронут не будет, а вектор получит преобразование. То есть получим:

[ U∙ K(0←x) ] ∙ { K(x←x) ∙ K(x←x) ∙ Y(x) }  = u*  .

Далее последовательно можно записать:

[[ U∙ K(0←x) ] ∙  K(x←x) ] ∙ { K(x←x) ∙ Y(x) }  = u*  ,

[                      матрица                          ] ∙ {         вектор          }  = вектор  .

Аналогично и эту группу линейных алгебраических уравнений можно подвергнуть построчному ортонормированию, которое сделает строчки [матрицы] ортонормированными, {вектор} затронут не будет, а вектор получит преобразование. То есть получим:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

[[ U∙ K(0←x) ] ∙  K(x←x) ]  ∙ { K(x←x) ∙ Y(x) }  = u*  ,

Далее аналогично можно записать:

[[[ U∙ K(0←x) ] ∙  K(x←x) ]  ∙ K(x←x) ] ∙ {  Y(x)   }  = u*  ,

[                                        матрица                                         ] ∙ { вектор}  = вектор  .

Аналогично и эту группу линейных алгебраических уравнений можно подвергнуть построчному ортонормированию, которое сделает строчки [матрицы] ортонормированными, {вектор} затронут не будет, а вектор получит преобразование. То есть получим:

[[[ U∙ K(0←x) ] ∙  K(x←x) ]  ∙ K(x←x) ]  ∙ Y(x)  = u*  .

Аналогично можно проортонормировать матричное уравнение краевых условий и для правого края независимо от левого края.

Далее проортонормированные уравнения краевых условий:

[ U∙ K(0←x) ] ∙ Y(x)  = u*  ,

[ V∙ K(1←x) ] ∙ Y(x)   =  v*

как и ранее объединяются в одну обычную систему линейных алгебраических уравнений с квадратной матрицей коэффициентов для нахождения искомого вектора Y(x) :

Y(x)  = .

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

6. Метод дополнительных краевых условий.

Этот метод еще не обсчитан на компьютерах.

Запишем на левом крае ещё одно уравнение краевых условий:

M ∙ Y(0) = m .

В качестве строк матрицы M можно взять те краевые условия, то есть выражения тех физических параметров, которые не входят в параметры краевых условий левого края L или линейно независимы с ними. Это вполне возможно, так как у краевых задач столько независимых физических параметров какова размерность задачи, а в параметры краевых условий входит только половина физических параметров задачи. То есть, например, если рассматривается задача об оболочке ракеты, то на левом крае могут быть заданы 4 перемещения. Тогда для матрицы М можно взять параметры сил и моментов, которых тоже 4, так как полная размерность такой задачи – 8. Вектор m правой части неизвестен и его надо найти и тогда можно считать, что краевая задача решена, то есть сведена к задаче Коши, то есть найден вектор Y(0)  из выражения:

Y(0) = ,

то есть вектор Y(0)  находится из решения системы линейных алгебраических уравнений с квадратной невырожденной матрицей коэффициентов, состоящей из блоков U и M.

Аналогично запишем на правом крае ещё одно уравнение краевых условий:

N ∙ Y(0) = n ,

где матрица N записывается из тех же соображений дополнительных линейно независимых параметров на правом крае, а вектор n неизвестен.

Для правого края тоже справедлива соответствующая система уравнений:

Y(1) = .

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

Запишем Y(1) = K(1←0) ∙Y(0) + Y*(1←0)  и подставим в последнюю систему линейных алгебраических уравнений:

∙ [ K(1←0) ∙Y(0) + Y*(1←0) ]  = ,

∙ K(1←0) ∙Y(0)  =  —  Y*(1←0),

∙ K(1←0) ∙Y(0)  =  ,

∙ K(1←0) ∙Y(0)  =  .

Запишем вектор Y(0)  через обратную матрицу:

Y(0) = ∙

и подставим в предыдущую формулу:

 

∙ K(1←0) ∙ ∙  = .

Таким образом, мы получили систему уравнений вида:

В ∙  = ,

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

где матрица В известна, векторы u и s известны, а векторы m и t неизвестны.

Разобьем матрицу В на естественные для нашего случая 4 блока и получим:

∙  = ,

откуда можем записать, что

В11 ∙ u + B12 ∙ m = s,

B21 ∙ u + B22 ∙ m = t.

Следовательно, искомый вектор m вычисляется по формуле:

m = B12 ∙ (s – B11∙ u).

А искомый вектор n вычисляется через вектор t:

t  = B21 ∙ u + B22 ∙ m,

n = t + N ∙ Y*(1←0).

В случае «жестких» дифференциальных уравнений предлагается выполнять поочередное построчное ортонормирование.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Запишем приведенную выше формулу

∙ K(1←0) ∙ ∙  =

в виде:

∙ K(1←x2) ∙ K(x2←x1) ∙ K(x1←0) ∙ ∙  = .

Эту формулу можно записать в виде разделения левой части на произведение матрицы на вектор:

[ ∙ K(1←x2) ]     ∙   { K(x2←x1) ∙ K(x1←0) ∙ ∙  }    =

[     матрица       ]      ∙   {                        вектор                        }     =    вектор

Эту группу линейных алгебраических уравнений можно подвергнуть построчному ортонормированию, которое сделает строчки [матрицы] ортонормированными, {вектор} затронут не будет, а вектор получит преобразование. То есть получим:

[ ∙ K(1←x2) ]     ∙   { K(x2←x1) ∙ K(x1←0) ∙ ∙  }    =

Здесь следует сказать, что подвектор t подвергать преобразованию не нужно, так как невозможно, так как его первоначальное значение не известно. Но подвектор t нам оказывается и не нужен для решения задачи.

Далее запишем:

[[ ∙ K(1←x2) ] ∙ K(x2←x1)]     ∙      { K(x1←0) ∙ ∙  }    =

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

[                       матрица                    ]      ∙     {                вектор            }     =    вектор

Аналогично и эту группу линейных алгебраических уравнений можно подвергнуть построчному ортонормированию, которое сделает строчки [матрицы] ортонормированными, {вектор} затронут не будет, а вектор получит преобразование. То есть получим:

[[ ∙ K(1←x2) ] ∙ K(x2←x1)]      ∙   { K(x1←0) ∙ ∙  }    =    .

И так далее.

В результате поочередного ортонормирования получим:

В ∙  = ,

∙  =  .

Следовательно, искомый вектор m вычисляется по формуле:

m = B12 ∙ (s – B11∙ u).

7. Формула для начала счета методом прогонки С.К.Годунова.

Эта формула обсчитана на компьютерах в кандидатской диссертации.

Рассмотрим проблему метода прогонки С.К.Годунова.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Предположим, что рассматривается оболочка ракеты. Это тонкостенная труба. Тогда система линейных обыкновенных дифференциальных уравнений будет 8-го порядка, матрица A(x) коэффициентов будет иметь размерность 8х8, искомая вектор-функция Y(x)  будет иметь размерность 8х1, а матрицы краевых условий будут прямоугольными горизонтальными размерности 4х8.

Тогда в методе прогонки С.К.Годунова для такой задачи решение ищется в следующем виде:

Y(x)  = Y(x) c +  Y(x) c + Y(x) c + Y(x) c + Y*(x),

или можно записать в матричном виде:

Y(x)  = Y(x) ∙ c  +  Y*(x),

где векторы  Y(x), Y(x), Y(x), Y(x) – это линейно независимые вектора-решения однородной системы дифференциальных уравнений, а вектор Y*(x) – это вектор частного решения неоднородной системы дифференциальных уравнений.

Здесь Y(x)=|| Y(x), Y(x), Y(x), Y(x)  || это матрица размерности 8х4, а c  это соответствующий вектор размерности 4х1из искомых констант c,c,c,c.

Но вообще то решение для такой краевой задачи с размерностью 8 (вне рамок метода прогонки С.К.Годунова) может состоять не из 4 линейно независимых векторов Y(x), а полностью из всех 8 линейно независимых векторов-решений однородной системы дифференциальных уравнений:

Y(x)=Y(x)c+Y(x)c+Y(x)c+Y(x)c+

+Y(x)c+Y(x)c+Y(x)c+Y(x)c+Y*(x),

И как раз трудность и проблема метода прогонки С.К.Годунова и состоит в том, что решение ищется только с половиной возможных векторов и констант и проблема в том, что такое решение с половиной констант должно удовлетворять условиям на левом крае (стартовом для прогонки) при всех возможных значениях констант, чтобы потом найти эти константы из условий на правом крае.

То есть в методе прогонки С.К.Годунова есть проблема нахождения таких начальных значений Y(0), Y(0), Y(0), Y(0), Y*(0) векторов Y(x), Y(x), Y(x), Y(x), Y*(x), чтобы можно было начать прогонку с левого края x=0, то есть чтобы удовлетворялись условия U∙Y(0) = u на левом крае при любых значениях констант c,c,c,c.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Обычно эта трудность «преодолевается» тем, что дифференциальные уравнения записываются не через функционалы, а через физические параметры и рассматриваются самые простейшие условия на простейшие физические параметры, чтобы начальные значения Y(0), Y(0), Y(0), Y(0), Y*(0) можно было угадать. То есть задачи со сложными краевыми условиями так решать нельзя: например, задачи с упругими условиями на краях.

Ниже предлагается формула для начала вычислений методом прогонки С.К.Годунова.

Выполним построчное ортонормирование матричного уравнения краевых условий на левом крае:

U∙Y(0) = u,

где матрица U прямоугольная и горизонтальная размерности 4х8.

В результате получим эквивалентное уравнение краевых условий на левом крае, но уже с прямоугольной горизонтальной матрицей U размерности 4х8, у которой будут 4 ортонормированные строки:

U∙Y(0) = u,

где в результате ортонормирования вектор u преобразован в вектор u.

Как выполнять построчное ортонормирование систем линейных алгебраических уравнений можно посмотреть в [Березин, Жидков].

Дополним прямоугольную горизонтальную матрицу U до квадратной невырожденной матрицы W:

W = ,

где матрица М размерности 4х8 должна достраивать матрицу U до невырожденной квадратной матрицы W размерности 8х8.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

В качестве строк матрицы М можно взять те краевые условия, то есть выражения тех физических параметров, которые не входят в параметры левого края или линейно независимы с ними. Это вполне возможно, так как у краевых задач столько независимых физических параметров какова размерность задачи, то есть в данном случае их 8 штук и если 4 заданы на левом крае, то ещё 4 можно взять с правого края.

Завершим ортонормирование построенной матрицы W, то есть выполним построчное ортонормирование и получим матрицу W размерности 8х8 с ортонормированными строками:

W = .

Можем записать, что

Y(0) = (М)транспонированная = М.

Тогда, подставив в формулу метода прогонки С.К.Годунова, получим:

Y(0) = Y(0) ∙с + Y*(0)

или

Y(0) = М∙с + Y*(0).

Подставим эту последнюю формулу в краевые условия U∙Y(0) = u и получим:

U∙ [ М∙с + Y*(0) ]= u.

Отсюда получаем, что на левом крае константы c уже не на что не влияют, так как

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

U∙ М = 0 и остается только найти Y*(0) из выражения:

U∙ Y*(0) = u.

Но матрица U имеет размерность 4х8 и её надо дополнить до квадратной невырожденной, чтобы найти вектор Y*(0) из решения соответствующей системы линейных алгебраических уравнений:

Y*(0) = ,

где 0 – любой вектор, в том числе вектор из нулей.

Отсюда получаем при помощи обратной матрицы:

Y*(0) = ∙ ,

Тогда итоговая формула для начала вычислений методом прогонки С.К.Годунова имеет вид:

Y(0) = М∙с + ∙ .

8. Второй алгоритм для начала счета методом прогонки С.К.Годунова.

Этот алгоритм обсчитан на компьютерах в кандидатской диссертации.

Этот алгоритм требует дополнения матрицы краевых условий U до квадратной невырожденной:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

 

Начальные значения Y(0), Y(0), Y(0), Y(0), Y*(0) находятся из решения следующих систем линейных алгебраических уравнений:

Y*(0) = ,

Y(0) = , где i = , , , ,

где 0 – вектор из нулей размерности 4х1.

9. Замена метода численного интегрирования Рунге-Кутта в методе прогонки С.К.Годунова.

Эта замена формул Рунге-Кутта на формулу теории матриц обсчитана на компьютерах в кандидатской диссертации.

В методе С.К.Годунова как показано выше решение ищется в виде:

Y(x)  = Y(x) ∙ c  +  Y*(x).

На каждом конкретном участке метода прогонки С.К.Годунова между точками ортогонализации можно вместо метода Рунге-Кутта пользоваться теорией матриц и выполнять расчет через матрицу Коши:

Y(x) = K(x- x) ∙Y(x).

Так выполнять вычисления быстрее, особенно для дифференциальных уравнений с постоянными коэффициентами.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

И аналогично через теорию матриц можно вычислять и вектор Y*(x) частного решения неоднородной системы дифференциальных уравнений. Или для этого вектора отдельно можно использовать метод Рунге-Кутта, то есть можно комбинировать теорию матриц и метод Рунге-Кутта.

10. Метод половины констант.

Этот метод пока не обсчитан на компьютерах.

Выше было показано, что решение системы линейных обыкновенных дифференциальных уравнений можно искать в виде только с половиной возможных векторов и констант. Была приведена формула для начала вычислений:

Y(0) = М∙с + ∙ .

Из теории матриц известно, что если матрица ортонормирована, то её обратная матрица есть её транспонированная матрица. Тогда последняя формула приобретает вид:

Y(0) = М∙с + U∙u

или

Y(0) = U∙u +  М∙с

или

Y(0) =  ∙ ,

Таким образом записана в матричном виде формула для начала счета с левого края, когда на левом крае удовлетворены краевые условия.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Далее запишем     V∙Y(1) = v      и    Y(1) = K(1←0) ∙Y(0) + Y*(1←0)   совместно:

V∙ [ K(1←0) ∙Y(0) + Y*(1←0) ]  = v

V∙ K(1←0) ∙Y(0)   = v — V∙Y*(1←0)

и подставим в эту формулу выражение для Y(0):

V∙ K(1←0) ∙ ∙ = v — V∙Y*(1←0).

V∙ K(1←0) ∙ ∙ = p.

Таким образом мы получили выражение вида:

D ∙ = p,

где матрица D имеет размерность 4х8 и может быть естественно представлена в виде двух квадратных блоков размерности 4х4:

∙ = p.

Тогда можем записать:

D1∙ u + D2 ∙ c = p.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Отсюда получаем, что:

c =  D2 ∙ ( p — D1∙ u )

Таким образом, искомые константы найдены.

Далее показано как применять этот метод для решения «жестких» краевых задач.

Запишем

V∙ K(1←0) ∙ ∙ = p.

совместно с K(1←0)  = K(1←x2) ∙ K(x2←x1) ∙ K(x1←0) и получим:

V∙ K(1←x2) ∙ K(x2←x1) ∙ K(x1←0) ∙ ∙ = p.

Эту систему линейных алгебраических уравнений можно представить в виде:

[ V∙ K(1←x2)  ] ∙  {  K(x2←x1) ∙ K(x1←0) ∙ ∙  } =  p.

[     матрица    ] ∙  {                                   вектор                                     } = вектор

Эту группу линейных алгебраических уравнений можно подвергнуть построчному ортонормированию, которое сделает строчки [матрицы] ортонормированными, {вектор} затронут не будет, а вектор получит преобразование. То есть получим:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

[ V∙ K(1←x2)  ]  ∙  {  K(x2←x1) ∙ K(x1←0) ∙ ∙  } =  p.

И так далее.

В итоге поочередного вычленений матриц слева из вектора и ортонормирования получим систему:

D ∙ = p,

Отсюда получаем, что:

c =  D2 ∙ (pD1∙ u)

Таким образом, искомые константы найдены.

11. Применяемые формулы ортонормирования.

Эти формулы обсчитаны в кандидатской диссертации.

Взято из: Березин И.С., Жидков Н.П. Методы вычислений, том II, Государственное издательство физико-математической литературы, Москва, 1962 г. 635 стр.

 

Пусть дана система линейных алгебраических уравнений порядка n:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

Здесь над векторами поставим черточки вместо их обозначения жирным шрифтом.

Будем рассматривать строки матрицы А системы как векторы:

=(,,…,).

Ортонормируем эту систему векторов.

Первое уравнение системы А= делим на .

При этом получим:

++…+=,    =(,,…,),

где =, =, =1.

Второе уравнение системы заменяется на:

++…+=,    =(,,…,),

где =, =,

=-(,),  =-(,).

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

Аналогично поступаем дальше. Уравнение с номером i примет вид:

++…+=,    =(,,…,),

где =, =,

=-(,)-(,)-…-(,),

=-(,)-(,)-…-(,).

Процесс будет осуществим, если система линейных алгебраических уравнений линейно независима.

В результате мы придем к новой системе С=, где матрица С будет с ортонормированными строками, то есть обладает свойством С*С= E, где Е – это единичная матрица.

(Таким образом, решение системы можно записать в виде = С.)

12. Вывод формул, позаимствованный из «Теории матриц» Гантмахера.

Система линейных обыкновенных дифференциальных уравнений с постоянными коэффициентами имеет вид:

Y(x) = A Y(x) + F(x).          (1)

Разложим Y(x) в ряд Маклорена по степеням x:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Y(x)=Y + Yx + Yx/2!  +  …,  где Y=Y(0), Y= Y(0), …         (2).

Из (1) почленным дифференцированием при А=const и F(x)=0 получим:

Y= AY= AY,    Y= A Y = AY,             (3)

Положив в (3) x=0 и подставив в (2) получим:

Y(x) = Y + Ax Y + A x/2!  Y + … = e Y,      (4)

где e = E + Ax + A x/2!  + …, где Е – единичная матрица. (5)

Если принять x=x, то (4) заменится на

Y(x) = e Y(x),       (6)

Рассмотрим случай A=const и F≠0.

Введем в рассмотрение вектор-функцию Ya(x) в виде: Y(x)= eYa(x).         (7)

Продиффренцируем (7) и подставим в (1). Получим:

eYa(x) = F(x).   (8)

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

При получении (8) учитывалось, что:

=  = A + A x + A x/2!  + … = A e.

Из (8) следует, что:

Ya(x) = c + .            (9)

Подставим в (7) и получаем:

Y(x) = ec + e.          (10)

Положив x=x  в (10) получим:

c = e Y(x).     (11)

Окончательно получаем:

Y(x) = e Y(x) + e.               (12)

ЛИТЕРАТУРА

  1. Гантмахер Ф.Р. Теория матриц. – М.: Наука, 1988. – 548 с.
  2. Березин И.С., Жидков Н.П. Методы вычислений, том II, Государственное издательство физико-математической литературы, Москва, 1962 г., 635 с.

© к.ф.-м.н. Алексей Юрьевич Виноградов. Февраль 2010.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Смотрите мои сайтики:

www.AlexeiVinogradov.narod.ru

www.VinogradovAlexei.narod.ru

www.Vinogradov-Alexei.narod.ru

www.Vinogradov-math.narod.ru

Пишите мне:

[email protected]

P.S. 28 февраля 2010:

Мой отец (доктор физико-математических наук профессор МГТУ им. Баумана Юрий Иванович Виноградов) предложил использовать и другую (гораздо более эффективную по времени счета) матричную формулу вместо матричной экспоненты – что-то на основе Вольтерра. Это есть в статье в журнале «Математическое моделирование»:

Численный метод переноса краевых условий для жестких дифференциальных уравнений строительной механики
Журнал «ММ», Том: 14 (2002), Номер: 9, 3 стр. 1409-003r.pdf

P.P.S. 22 апреля 2010:

Я тут как-то на днях поговорил с отцом и вот что любопытное насчет метода выяснилось.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Когда я сам решал своим методом «переноса краевых условий» «жесткие» краевые задачи, то я в каждой рассматриваемой точке x=x* решал соответствующую систему линейных алгебраических уравнений для нахождения решения Y(x*) в каждой рассматриваемой точке x=x*.

А мой отец утверждает, что «жесткими» бывают только краевые задачи, а начальные задачи «жесткими» не бывают. И поэтому он находил решение Y(x*) в какой-то одной точке x=x* по моему методу, а дальше решал (влево и вправо от рассматриваемой точки x=x*) как задачу Коши (от найденных начальных условий Y(x*) в этой одной точке): Y(x)=K(х←x*)·Y(x*)+Y*(xx*). И если это так, то так решать, естественно, гораздо быстрее, так как надо только домножать на матрицу Коши (матрициант, матричную экспоненту) вместо решения систем линейных алгебраических уравнений.

Я думаю, что это очень существенное уточнение по скорости счета.

P.P.P.S. Метод решения «жестких» начальных задач, то есть «жестких» задач Коши. Придуман вечером 16 января 2008 года.

Метод переноса краевых условий работает успешно. Следовательно, по аналогии можно пытаться решать «жесткие» начальные задачи, то есть «жесткие» задачи Коши.

Пусть дана начальная задача:

Y(x) = A·Y(x),  Y(0) = Yнач,

Можем записать:

Y(0) = K(0←x) · Y(x),

K(0←x) · Y(x) = Yнач

K(0←x3) · K(x3←x2) · K(x2←x) · Y(x) = Yнач

[ K(0←x3) ] · { K(x3←x2) · K(x2←x) · Y(x) } = Yнач

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

[ матрица  ]    {                   вектор                  }     вектор

Выполняем построчное ортонормирование этой системы линейных алгебраических уравнений с квадратной матрицей коэффициентов и получаем систему с ортонормированными строчками в квадратной матрице:

[ K(0←x3) ]орт · {K(x3←x2) · K(x2←x) · Y(x)} = Yнач_орт

Аналогично записываем

[[ K(0←x3) ]орт · K(x3←x2) ] ·{ K(x2←x) · Y(x) } = Yнач

[                  матрица               ]   {         вектор        }     вектор

Далее выполняем построчное ортонормирование и получаем:

[[ K(0←x3) ]орт · K(x3←x2)]орт · {K(x2←x) · Y(x)} = Yнач_2орт

Аналогично получаем

[[ K(0←x3) ]орт · K(x3←x2)]орт · K(x2←x)]орт · Y(x) = Yнач_3орт

То есть получили итоговую систему линейных алгебраических уравнений для нахождения вектора Y(x). Для этого вектор Yнач_3орт  надо домножить слева на матрицу транспонированную по отношению к левой ортонормированной матрице, так как если матрица ортонормированна, то обратная ей есть транспонированная матрица.

P.P.P.P.S. 11 сентября 2009:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диссер

Долго было лень записывать, но вдруг кто-то сам не сразу догадается – ещё один метод решения «жестких» начальных задач, то есть «жестких» задач Коши. Хотя у меня есть такое подозрение, что «жесткими» бывают только краевые задачи, а начальные задачи «жесткими» НЕ бывают. Но на всякий случай приведу метод решения «жестких» начальных задач.

Начальные условия имеют вид:

Y(0) = Yнач.

Полное решение системы обыкновенных дифференциальных уравнений

Y(x) = A · Y(x) + F(x)

имеет вид:

Y(x) = K(x←0)  · Y(0) + Y*(x←0).

Или можно записать:

Y(0) = K(0←x1)  · Y(x1) + Y*(0←x1).

Подставим это выражение в краевые условия и получим:

K(0←x1)  · Y(x1) + Y*(0←x1) = Yнач

или

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

K(0←x1)  · Y(x1)  = Yнач — Y*(0←x1)

или

K1 · Y(x1) = Y1.

Проортонормируем это выражение построчно и получим эквивалентное выражение:

K1орто · Y(x1) = Y1орто.

Тогда

Y(x1) = (K1орто)транспонир · Y1орто.

Подставим вместо Y(x1) выражение через Y(x2) и получим:

K(x1←x2)  · Y(x2) + Y*(x1←x2) = (K1орто)транспонир · Y1орто

или

K(x1←x2)  · Y(x2) = (K1орто)транспонир · Y1орто Y*(x1←x2)

или

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

K2 · Y(x2) = Y2.

Проортонормируем построчно и получим эквивалентное выражение:

K2орто · Y(x2) = Y2орто.

Тогда:

Y(x2) = (K2орто)транспонир · Y2орто.

И так далее.

P.P.P.P.P.S. Метод для численного интегрирования дифференциальных уравнений.

Читали нам как-то в бауманке численные методы решения дифференциальных уравнений. И, кажется, приводили аналитический вывод формул одного из авторов. Или это просто мелькнуло в учебнике (я имею в виду вывод формул). Уже не очень помню. Запомнилась только собственная мысль, что людям вообще-то проще всего даются геометрические аналогии и выводы, сделанные на основе понятных геометрических картинок. Ну, вот тогда я и нарисовал один из вариантов численного решения дифференциальных уравнений и помню даже перевёл геометрические картинки в буквенные формулы приближённых вычислений. Сейчас повторно выводить буквенные формулы для численного интегрирования дифференциальных уравнений мне не кажется интересным. А вот привести картинки тех студенческих мыслей вполне можно для обсуждения.

Далее идёт картинка с текстом и с рисунками численного интегрирования:

17 сентября 2010: Забыл ранее кое-что добавить – добавляю сегодня насчет обратной матрицы.

Для однородной системы дифференциальных уравнений имеем:

Y(x) = K(x←x) ∙ Y(x).

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

Можем записать:

Y(x) = K(x←x) ∙ Y(x) и

Y(x) = K(x←x) ∙ Y(x).

Подставляем одну формулу в другую и получаем:

Y(x) = K(x←x) ∙ Y(x) = K(x←x) ∙ K(x←x) ∙ Y(x),

то есть получаем:

Y(x) = K(x←x) ∙ K(x←x) ∙ Y(x),

но последнее возможно только когда

K(x←x) ∙ K(x←x) = Е – единичная матрица,

то есть матрицы K(x←x) и K(x←x) взаимообратны.

То есть доказано, что

K(x←x) = K(x←x),

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

то есть

K(x-x) = K(x-x).

18 сентября 2010: Еще кой о чем вспомнил – вычисление матрицы Коши методами типа Рунге-Кутта.

Я сейчас ниже выскажу мысли, которые мне кажутся очевидными, но они приводятся без доказательства и они не проверялись вычислительными экспериментами – так что все в Ваших руках (можно проверить численно на компьютере и можно написать и опубликовать соответствующие статьи).

Итак. Матрица Коши K(x-x) вычисляется как матричная экспонента или по формуле Вольтерра (смотри выше). Я же считаю очевидным, что матрицу Коши можно вычислять и методами типа Рунге-Кутта и другими аналогичными численными методами, включая мой «геометрический» (неопробованный) численный метод, который приводится на картинке выше.

Матрица Коши K(x-x) обладает тем свойством, что ее значение при нулевом аргументе (x-x)=∆х=0 это есть единичная матрица K(0)=Е.

Тогда мне кажется очевидным, что вектора, составляющие матрицу Коши K(x-x) при НЕнулевом аргументе ∆х=(x-x)≠0 можно получить следующим путем: берем из единичной матрицы Е вертикальный вектор (столбец) с номером i и методом Рунге-Кутта от этого взятого начального значения (i-го столбца из единичной матрицы) вычисляем некий вектор-столбец на выбранном шаге интегрирования ∆х=(x-x)≠0 и записываем этот вектор-столбец в результирующую матрицу Коши на свое же i-ое место, то есть полученный вектор-столбец записываем в качестве i-го столбца результирующей матрицы Коши. И так формируем все столбцы вычисляемой (результирующей) матрицы Коши.

Например, для дифференциальных уравнений цилиндрической оболочки ракеты матрица Коши имеет размерность 8х8, то есть состоит из 8 нужных нам столбцов размерности 8х1. И методами типа Рунге-Кутта мы вычисляем соответственно 8 столбцов матрицы Коши K(x-x) на выбранном шаге интегрирования ∆х=(x-x)≠0, беря в качестве начальных значений векторов-столбцов для метода Рунге-Кутта столбцы размерности 8х1 из единичной матрицы размерности 8х8.

Кстати, для этого случая круговой цилиндрической оболочки мы имеем систему дифференциальных уравнений с постоянными коэффициентами и поэтому можем вычислить матрицу Коши однажды на одном маленьком интервале ∆х=(x-x)≠0, а на всем интервале задачи (1-0) получаем полную матрицу Коши перемножением самой на себя единожды вычисленной матрицы Коши малого участка: K(1-0)= K(∆х)∙ K(∆х)∙…∙ K(∆х), что очень сильно ускоряет вычисления по сравнению со случаями переменных коэффициентов (конус, сфера), так как при случае переменных коэффициентов приходится вычислять матрицы Коши K(x-x) независимо для каждого маленького отдельного шага (x-x) полного интервала интегрирования всей задачи с их аналогичным последующим перемножением.

То есть для случая постоянных коэффициентов системы дифференциальных уравнений не особенно важна скорость вычисления матрицы Коши K(∆х) одного отдельного малого участка ∆х, так как на малом участке матрица Коши K(∆х) вычисляется единожды, так как потом полученная матрица Коши малого участка перемножается сама на себя.

Тогда получается, что для вычисления матрицы Коши можно применять все известные численные методы, начиная с метода Рунге-Кутта и заканчивая моим предложенным выше и пока неопробованным «геометрическим» численным методом.

Скажу дополнительно, что касательно вектора частного решения неоднородной системы дифференциальных уравнений у меня нет таких очевидных догадок, какие высказаны выше касательно матрицы Коши.

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

То есть я не догадываюсь как вычислять методами типа Рунге-Кутта вектор Y*(x←x):

Y*(x←x) = e∙ e∙ F(t) dt   —

вектор частного решения неоднородной системы дифференциальных уравнений.

Хотя вполне может оказаться, что все на самом деле просто и надо применять метод Рунге-Кутта к начальному полностью нулевому вектору, то есть вектору, состоящему полностью из нулей. Но это не факт. То есть мы так очевидно получим вариант частного решения неоднородной системы дифференциальных уравнений, но не факт, что этот частный вектор будет удовлетворять формуле:

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x).

То есть, интегрируя Рунге-Куттом от полностью нулевого вектора неоднородную систему дифференциальных уравнений, мы получим таки некий вектор частного решения, но не факт, что этот вектор уложиться без поправок в формулу:

Y(x) = K(x←x) ∙ Y(x)  +  Y*(x←x).

То есть касательно частного вектора Y*(x←x) тут тоже все в Ваших руках – можете проверить численно разные численные методы типа Рунге-Кутта и может быть можете придумать коррекцию (если она потребуется) и можно публиковать статьи.

19 сентября 2010: еще об ускорении вычислений – применение «параллельных» вычислений.

В современной математике для современных компьютеров разрабатываются и применяются так называемые методы «параллельных вычислений». Особенно для промышленных суперкомпьютеров на основе многопроцессорной технологии. Эти методы сильно ускоряют вычисления для практических задач в НИИ и КБ.

Приведенные здесь методы как будто специально предназначены для «параллельных вычислений» даже без дополнительной модификации.

В приведенных методах используется матрица Коши и ее свойство перемножаемости:

Нужна помощь в написании магистерской?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена магистерской

K(x←x) = K(x←x) ∙ K(x←x) ∙ … ∙ K(x←x) ∙ K(x←x).

Так вот, очевидно, что множители K(x←x), K(x←x), … , K(x←x), K(x←x), входящие в формулу для вычисления матрицы Коши K(x←x), вычисляются полностью независимо друг от друга и поэтому могут вычисляться «параллельно» без привлечения дополнительных математических приемов.

То есть отдельные матрицы Коши на отдельных участках (для случая системы дифференциальных уравнений с переменными коэффициентами) могут вычисляться «параллельно» и например, просто на разных процессорах многопроцессорного компьютера какого-нибудь КБ (конструкторского бюро).

Аналогично «параллельно» может вычисляться и вектор частного решения неоднородной системы дифференциальных уравнений, так как подвектора Y*(x←x), из которых складывается частный вектор, вычисляются полностью независимо, то есть «параллельно»:

Y*(x←x) = Y*(x- x) = K(x- x) ∙K(x- t) ∙ F(t) dt =

= K(x- x) ∙ (E + A (x- t) + A (x- t)/2! + … ) ∙ F(t) dt =

= K(x- x) ∙ (EF(t) dt  + A∙(x- t) ∙ F(t) dt  + A/2! ∙(x- t) ∙ F(t) dt  + … ) ,

где A=А(х) – матрица коэффициентов системы дифференциальных уравнений, вычисленная в приближении ее постоянства на малом подучастке (x- x) и вычисленная на малом подучастке (x- x), например, в начальной точке x рассматриваемого подучастка (x- x). Хотя матрица А подучастка может быть вычислена и, например, в точке x подучастка (x- x): A=А(x) или в точке (x- x)/2: A=А((x- x)/2).

Кстати, в приведенной выше формуле осреднению может подвергаться не только матрица А: A=А(х)  коэффициентов системы дифференциальных уравнений, но и вектор F(t) может рассматриваться на подучастке (x- x) приближенно в виде постоянной величины F(х)=constant, что позволят вынести  его из под знака интеграла, что приводит к совсем простому ряду для вычислений на рассматриваемом подучастке.

То есть все приведенные выше методы Алексея Юрьевича Виноградова очень удачно встраиваются в современную тенденцию высокоскоростных «параллельных вычислений».

Средняя оценка 0 / 5. Количество оценок: 0

Поставьте оценку первым.

Сожалеем, что вы поставили низкую оценку!

Позвольте нам стать лучше!

Расскажите, как нам стать лучше?

929

Закажите такую же работу

Не отобразилась форма расчета стоимости? Переходи по ссылке

Не отобразилась форма расчета стоимости? Переходи по ссылке