2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 Вопрос по системе ДУ (maple)
Сообщение29.04.2012, 09:09 


08/04/09
43
Магнитогорск
Здравствуйте!

Создал новую тему (хотя похожие есть), т.к. к имеющимся уж очень давно в последний раз обращались! Решение о правильности моего поступка, по всей видимости, остается за модераторами. Очень нужна помощь в этом вопросе, поскольку обратиться больше не к кому. Вопрос в общем-то не сложный, но странный :-) Есть система ДУ с условиями:
Код:
sys:={14.52*a[3](t)+18.60*a[1](t)+2.35*(diff(a[3](t), t)), 18.60*a[3](t)+6.28*(diff(a[1](t), t))+33.61*a[1](t)};
ics:={a[1](0) = .63, a[3](1) = .21}

Решая ее в Maple 13 и 15 получаю
Код:
{a[1](t) = Float(undefined)*exp(-.9051675582*t)+Float(undefined)*exp(-10.60750253*t), a[3](t)=Float(undefined)*exp(-.9051675582*t)+Float(undefined)*exp(-10.60750253*t)}

где Float(undefined) -- насколько я понимаю, не определенное вещественное число. Возникает такое, вроде как, когда пытаешься найти значение функции, которое не принадлежит ее области определения. Однако, решая систему без условий получаем
Код:
dsold := evalf(dsolve(sys);assign(dsold):
{a[1](t) = _C1*exp(-.9051675582*t)+_C2*exp(-10.60750253*t), a[3](t) = -1.501489565*_C1*exp(-.9051675582*t)+1.776014119*_C2*exp(-10.60750253*t)}

далее подставляем $t=0,t=1$ соответственно в $a_1(t),a_3(t)$
Код:
raz := subs(t = 0, a[1](t)):dva := subs(t = 1, a[3](t)):
solve({dva = .21, raz = .63}, {_C1, _C2});
{_C1 = -.35, _C2 = .99}

находим $C_1,C_2$. Как это объяснить? И почему Maple не может решить эту систему с условиями сразу, найдя эти константы?
Задача у меня именно такая: одни условия задаются в начале временного промежутка, другие -- в конце. Система решается "в лоб" мейплом, если задавать $a_1(0)=0.63,a_3(0)=0.21$, но только в этом случае. Если оба условия задавать в $t=1$ -- тоже не решается. С чем это связано и почему так происходит? Может я все-таки что-то делаю не правильно?

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ 1 сообщение ] 

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group