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

Добрый вечер! Помогите пожалуйста численно решить систему дифференциальных уравнений с запаздыванием в пакете МАTHCAD.



Заранее большое спасибо за помощь!

0

Сообщение от anna_gorina численно решить систему дифференциальных уравнений с запаздыванием

Вы забыли указать в каком интервале надо решать задачу и каковы начальные условия.
Я выбрал свои. Для больших Т получаются очень большие числа.

1

Уважаемый Symon, огромное Вам спасибо! Не могли бы Вы прислать мне результат в форме прикрепленного файла, так как все очень мелко и неразборчиво. Пыталась увеличить масштаб, однако изображение стало ненамного лучше. Заранее благодарю!

0

Сообщение от anna_gorina Не могли бы Вы прислать мне результат

Пжлста
Система ДУ с запаздыванием.zip

0

Добрый день Symon! К сожалению, у меня не открывается посланный Вами файл с расширением .xmсd. Не могли бы Вы его переформатировать в .mсd и прислать? Надеюсь на Вашу помощь. Заранее большое спасибо.

0

Сообщение от anna_gorina не открывается посланный Вами файл с

Пжлста
Система ДУ с запаздыванием1.zip

0

Здравствуйте, Симон! К сожалению, то, что Вы мне прислали, тоже не открывается. Вероятно, потому что у меня старая версия - Mathcad 2001i. Какую версию Вы мне посоветовали бы? Кроме этого хотелось бы получить нормальный листинг решенной Вами системы, то есть две (или более) страниц полностью, так как в присланном Вами листинге в векторном уравнении F(t,w,f) никак не определена функция f, которая в Вашем текстовом комментарии равна u(t-T). К тому же, в листинге отсутствует функция решения системы дифференциальных уравнений rkfixed, или подобная ей, и, по-видимому, много еще чего.
Еще, некорректно и произвольным образом задана функция u0(t), которую, по моему мнению, надо находить, решая сначала аналогичную систему ДУ, но без запаздывания, а, затем подставляя ее в u0(t).
Надеюсь, вы мне ответите, посоветуете и пришлете. Заранее благодарю!

0

Сообщение от anna_gorina Какую версию Вы мне посоветовали бы?

Маткад 15 m020 или последняя m045

Сообщение от anna_gorina в присланном Вами листинге в векторном уравнении F(t,w,f) никак не определена функция f,

При определении функции, в частности, векторной F(t,w,f), в качестве аргументов можно использовать любые буквы (хоть русские). Что они из себя представлют и как они используются, становится понятным из разбора тела функции в правой части знака присвоения (по контексту). А когда наступает пора использовать эту функцию, вместо формальных параметров подставляем настоящие, которые должны быть определены выше или непосредственно в F. Например, вместо f подставляем u0.

Сообщение от anna_gorina некорректно и произвольным образом задана функция u0(t), которую, по моему мнению, надо находить, решая сначала аналогичную систему ДУ, но без запаздывания, а, затем подставляя ее в u0(t).

При постановке задачи, кроме ДУ обязательно указываются начальные и\или краевые условия. У вас их нет.
Об этом я уже писал. Нет и интервала, где нужно искать решение. Поэтому я выбрал свои. Ищу решение на промежутке [0,T]. Для уравнения с запаздыванием надо еще знать начальную функцию u0(t) на промежутке [-T,0]. Когда t меняется в промежутке [o,T], t-T меняется в [-T,0] Так как в уравнении фигурирует u(t-T), то не зная искомую функцию на [-T,0], не сможем определить и на [0,T].
Внимательно посмотрите оригинал постановки задачи, там должны быть начальные данные и подставляйте вместо моих. Если их там нет, скажите об этом преподу.

Сообщение от anna_gorina в листинге отсутствует функция решения системы дифференциальных уравнений rkfixed, или подобная ей, и, по-видимому, много еще чего.

Все, что надо, тут есть! Встроенные функции типа rkfixed для ДУ с запаздыванием не используются. По крайней мере я не знаю, как это сделать. Я решаю простым методом Эйлера с помощью рекуррентной фаормулы.

Сообщение от anna_gorina все очень мелко и неразборчиво.

Присылаю увеличенный вариант и с други начальным условием

0