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

Форум пользователей пакета Maple

Модератор: Admin

fenix85
Сообщения: 141
Зарегистрирован: Вт мар 01, 2011 12:06 pm

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

Сообщение fenix85 » Вт июл 17, 2018 7:13 pm

Подскажите пожалуйста, как получить интеграл от численного решения диф. уравнений. У меня есть решение - функция Z6(x), нужно ее проинтегрировать на интервале и получить график этого интеграла. Я не знаю как это сделать...(((
Вложения
eq.rar
(34.02 КБ) 21 скачивание

Markiyan Hirnyk
Сообщения: 1221
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение Markiyan Hirnyk » Чт июл 19, 2018 9:21 pm

Добавьте в конце вашего кода команду

Код: Выделить всё

plot(x->int(Z6(s), s = 0 .. x, numeric, epsilon = 0.1e-1) , 0 .. tau, numpoints = 5, colour = [black], thickness = [4, 4], labeldirections = [horizontal, vertical], labels = [s*M, `dθ`/ds]);

fenix85
Сообщения: 141
Зарегистрирован: Вт мар 01, 2011 12:06 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение fenix85 » Пт июл 20, 2018 3:03 pm

Зделал и получил вот такое сообщение Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct
Смотрел на странице документации по этой ошибке но так ничего и не понял(
У меня Maple 13.
Последний раз редактировалось fenix85 Пт июл 20, 2018 3:21 pm, всего редактировалось 1 раз.

Markiyan Hirnyk
Сообщения: 1221
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение Markiyan Hirnyk » Пт июл 20, 2018 3:21 pm

Пожалуйста, представьте через Dropbox или иной файлообменник mw file с выполненным кодом. Возможно также, что вы не выполнили все ваши предыдущие команды. Конечно же, программа должна начинаться с команды restart.

Markiyan Hirnyk
Сообщения: 1221
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение Markiyan Hirnyk » Сб июл 21, 2018 7:43 am

fenix85, спасибо за файл, переданный через Гугл. Действительно, при выполнении команды в Мэйпле 16.02 ( версии 13 не имею) получается та же ошибка. Работает в Мэйпле 2018.1.Возможно, изменился синтаксис Мэйпла. Подумаю над этим.

Markiyan Hirnyk
Сообщения: 1221
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение Markiyan Hirnyk » Сб июл 21, 2018 8:11 am

Команда

Код: Выделить всё

plot(evalf(Int('Z6(s)', s = 0 .. x)), x = 0 .. tau, numpoints = 5);
работает в Мэйпле 16.02, но медленно: выполняется примерно 6 минут на моем компе.

fenix85
Сообщения: 141
Зарегистрирован: Вт мар 01, 2011 12:06 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение fenix85 » Сб июл 21, 2018 7:36 pm

Благодарю за помощь! У меня тоже заработало, на моем компьютере выполняется примерно 7 минут. Правда меня удивляет что интеграл получается прямая линия, тогда как сама функция осциллирует.

Markiyan Hirnyk
Сообщения: 1221
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение Markiyan Hirnyk » Сб июл 21, 2018 8:25 pm

Амплитуда колебаний порядка 10^(-7) около среднего значения 0.14e-4, поэтому график первообразной и близок к прямой линии.

Markiyan Hirnyk
Сообщения: 1221
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение Markiyan Hirnyk » Сб июл 21, 2018 10:28 pm

Еще один подход изложен здесь. Идея состоитв следующем: добавить

Код: Выделить всё

diff(z9(x), x) = z6(x)
в , в и в c понятными дальнейшими изменениями. Проверял, работает.

fenix85
Сообщения: 141
Зарегистрирован: Вт мар 01, 2011 12:06 pm

Re: интеграл от решения дифференциальных уравнений

Сообщение fenix85 » Пн июл 23, 2018 8:23 pm

спасибо большое!