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

Все доброго дня

Возникла следующая проблема:
Пытаюсь решить систему уравнений, состоящую из 6 дифференциальных, и 15 алгебраических уравнений через встроенную функции Odesolve.
Известны все начальные условия.
После ввода всех необходимых данных в функцию Odesolve возникает ошибка : "This value must be a finction, but has form: any1"

Прошерстил кучу источников, почитал книжки по Маткаду
В итоге не нашел у себя причин к появлению данной ошибки.
Операторы дифференцирования введены с вкладки "Calculus", не в ручную, что обычно приводит к появлению данной ошибки.

Также пробовал решить "укороченную" систему уравнений состоящую из 4 дифф. уравнений и 8 алгебраических с помощью встроенной функции Radau (для жестких систем) и rkfixed

В первом случае (Radau) маткад ругается на то что неверная размерность у матрицы М - "These array dimensions do not match" , а в случае использования функции rkfixed возникает ошибка "This function cannot be used here"

Уже который день бьюсь с данными проблемами и не могу понять, в чем проблема

В архиве лежат 2 файла, соответственно для решения через Odesolve и другой с решением через Rkfixed и Radau

Миниатюры

     

 

Вложения

SDE_problem.ZIP (161.3 Кб, 23 просмотров)

0

У Вас интегрирование внутри блока Odesolve идет, но срывается, видимо из-за переполнения (деление на 0). Надо проверить все правые части системы (у Вас действительно деление на 0).

0

А разве при этом программа не выдала бы сообщение: "Devided by zero"?
Там где возможно деление на ноль в правой части уравнения, я задал ненулевые начальные условия.
Мое предположение о причине появления данной ошибки состоит в том, что Маткад не распознает одно из написанных уравнений как функцию от аргумента "t".
Но в каком месте ?Не могу понять

0

Сообщение Divided by zero не выдается при работе программного блока Odesolve.
У Вас ещё похоже проблема с последним уравнением wg(t)=w*psi(t), потому что оно не задает уравнение как таковое (w - константа). Уберите эту функцию.

0

Не помогло.
Какие по вашему мнению могут быть причины появления данной ошибки кроме тех, что вы уже написали ?

0

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

0