Главная » Просмотр файлов » Дьяконов В.П. Maple 9.5 и 10 в математике, физике и образовании

Дьяконов В.П. Maple 9.5 и 10 в математике, физике и образовании (1185901), страница 67

Файл №1185901 Дьяконов В.П. Maple 9.5 и 10 в математике, физике и образовании (Дьяконов В.П. Maple 9.5 и 10 в математике, физике и образовании.djvu) 67 страницаДьяконов В.П. Maple 9.5 и 10 в математике, физике и образовании (1185901) страница 672020-08-25СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 67)

Реп(ение дифференциальных уравнений специального вида 467 хождению вычислительного процесса, в ходе которого решение грубо отлично от существующего. Мар!е в большинстве случаев дает верное решение даже без указания метода решения. Это связано с тем, что система дифференциальных уравнений при решении его функцией Жо)че анализируется и в зависимости от результатов анализа выбирается наиболее подходящий метод решения. Кроме того, большинство методов (например, самый распространенный г[(Г45) реализуют алгоритмы контроля погрешности вычислений и дробления шага решения, если погрешность оказывается больше заданной. В связи с указанным решение жестких систем дифференциальных уравнении средствами системы Мар!е не вызывает особых трудностей и может быть осуществлено даже при выборе не вполне удачного метода. Однако при этом возможны следующие ситуации: ° может резко возрасти время вычислений из за чрезмерно сильного уменьшения шага решения; ° может оказаться превышенным число итерации в ходе дробления шага; лля «особо жестких» систем адаптивныи выбор шага может не помочь и погрешность решения будет большой.

Во избежание этого рекомендуется при решении жестких систем дифференциальных уравнений все же пользоваться специально для них созданными методами, например методом Розенброка (опция п)е(Ьо([ = гоаепЬгос)( для функции озо)че). 7.7.2. Примеры решения жестких систем дифференциальных уравнений В качестве первого примера исследуем и решим следующую систему дифференциальных уравнений (файл Ыез): > йед2:= и1ЕГ (и (Е), Е) = -11*и (1) +9*ч (е), с)1ЕГ (ч(1), Е) 9*и (Е) -11*ч(Ф.) ( Нед2:= — и(() =-1! и(()+9ч((), — ч(() = 9 и(()-11чЯ (! д с(( 4( Загрузив пакет!гпа!й вычислим собственные значения матрицы данной системы дифференциальных уравнений: > еьсЬ (11па1д): И: =~паьгьх (2, 2, [-11, 9, 9, -11) ); -11 9 > де."=еьдеача1иеа(И) ее:=-2, -20 Они оказались отрицательными.

Кроме того. очевидно, что значение жесткости данной системы з = ! О. Его трудно назвать очень большим, но в целом условия жесткости лля данной системы выполняются. Теперь решим зту систему методом Розенберга. Решение представлено на рис. 7.22. Обратите внимание на то, что представлены две точки и график решения. К достоинствам реализации примененного метода относится отсутствие необходимости в составлении матрицы Якоби. которую приходится задавать при использовании ряда функций системы Ма(- Ьса([, имеющихся для решения жестких систем дифференциальных уравнений (9!.

Еше один пример задания и решения жесткой системы дифференциальных уравнений представлен на рис. 7.23. Собственные значения матрицы этой системы 7. 7. Решение дифференциальных уравнений еиециальнага вида 4б9 равны — 2 и -1000, а жесткость системы з = 500 (проверьте сами по аналогии с ранее приведенным примером). Таким образом, зта система намного жестче, чем система из первого примера. Обратите внимание на то, что она решается без задания метода решения, но с опцией мИТ= ггве, вынуждаюшей Мар!е выбирать метод для решения жесткю систем дифференциальных уравнений.

7.7.3 Пример решения системы жестких дифференциальных уравнений химической кинетики Жесткие системы дифференциальных уравнений, часто описывают кинетику химических процессов, например, растворение веществ в растворах нли смешивание газов. На рис. 7.24 показано решение жесткой сисгемы из трех дифференциальных уравнений, описывающих один из типовых химических процессов — какой именно в данном случае не важно.

"„;,:"р: „,:.т'...,,. ° .„.„';..а ,"';:",! с л.-', "" °:,.',~,'; ":.'::,':.,д.::;:: 47() Глава 7. Решение ди(7)ференциальных уравнений описываюшие динамику развития колебаний в различных колебательных системах, например, автогенераторах на электронных лал(пах, полевых и биполярных транзисторах. Пример задания и решения дифференциального уравнения Ван-Дер Поля при сравнительно малом п)ц = ! (и при выборе метода решения по умолчанию) представлен на рис. 7.25. Нетрудно заметить, что выбор Мар!е пал на метод г)((45 и что этот метод не очень удачен даже для этого метода с (пц = 1.

Хотя обшая форма колебаний (близкая к синусоидальной, но все же заметно искаженная) в интервале г от 0 до 20 просматривается, уже в данном случае видна нестабилыюсть колебаний. При увеличен максимального значения т до !00 и более, нестабильность колебаний становится весьма заметна (проверьте это сами). Гн ен 1 ) с г с о.- ь, )нь я ь [ Рец)ение дифференциального уравения Ван-Дер Поля при малом параметре п)Ц .г! > теататт) аечз: — а(гг(у(ь),ь,м-ни*()-у(М "2) а)тт(у(ь),т)~у(т) 7.7. Решение дифференциа[ьньп)равнеуюий енециальнага вида 471 МС вес*у с и [ Решение дифференциального уравения Ван-Дер Поля при большом параметре гпо ~ ! > се>С>гС; Оеоа:- О111(у(С),с,С)-о >[1 у[С) "У) *О)[1[у(С),С) «у[С) ЬуЗ =~ „у(г)/ — р([ — у(г) )[--у(г))слг)=о а > ео:-аооо:[оз у- (у[о) -д,о[у) (о)-о); к) = ( у(0) = 1.

[)() ХО) = 0 ) > О>013: о>о1не( (с[етз) 00100 1оз, ош ес[о, гн1ое 0..6000, >СЗГЕ Сгые) г [ свои = р и оно~Мосс) ееа ргос > р)оСе[ое рзоС) [«[еозЗ, [С,у(С) ),оо)о — Ьз .Ш: с; 'э"~,' зн; ;",.;"д",0 сг;; н.>.Г' .л ." '....:. '.

', '": . ', "',::...":,"'у."., .,)рй 472 Глава 7. Решение дифференциа[ьных уравнений е [а [ог Г И С КОРаВ 1 СГ .С ~4 ОСЬ Решение дифференциального уравнения с двумя ираевь[ми условиями > вес] с- 4 О[[1[О[11[к(Г],Г],Х]+х[Х] - Г[ 449 = 4[ х(г) а «(г) = г 2 > бхо1 с — Охо1ое ((Оес[, х (О) — 5, х (5] 151, ооеех1о] с с(ао].= ргос(х аср) гед ргос > с[хо1(0.);с1хо1(5.] г с( Г = 0 . х(Г) = 5 ОООООО КИКККК Х88, — х(Л = 12 701?27872980]914 ' с]г [- с] г = 5, х(с) =! 5, — х(г) = -9 87054432828710482 ~ се > ю[Ю(р1охе] годер1ос[дхо1, [(С,х[Х]], [1,15]],О..б,оо1ог=-ыао[с] еггга о, г огс ь .са г г ха о ссог о сг ог сг о* ш«, с г ь ао)о*с а ;, сг ',, ка 74 *Ф;5:;,сс"„:„; 7.8.

Региение дифференииальных уравнений с частными производнъгми 473 ~ 1цпсь — (опция) множество или список с неопределенными функциями или именами; ° Н!ЙТ вЂ” (опция) равенство в форме Н1)чТ = агяцгпепц где аргумент может быть символом '+', ", любым алгебраическим выражением или строкой 'ыг)р'; ° ПЧТЕОВАТŠ— (опция) задает автоматическое интегрирование для множества ООЕь (если РОЕ решается при разделении переменных; е Ьцй6 — опция, задающая попытку построения явного выражения для неопределенной функции, независимо от общности найденного решения; ° пцгпепс — ключевое слова, задающее решение в численном виде; ° о1пег орйопь — другие опции.

7.8.2. Инструментальный пакет расширения РОЕтоо1 Для решения дифференциальных уравнений с частными производными и его визуализации в Мар!е 9.5 служит специальный инструментальный пакет РОЕгоо1: > пзев(РРВгоо1о1; (РОЕргог, Ьийа', сагезрйл сйатггр, дсйапее, дсоеЯЪ, дес1аге, с(Яогдег. аро!у~опп, дзиЬз, тарсге, зерагаЬгйгу, зр)йвгг1р, зрнувуз, ипаесlаге) Ввиду небольшого числа функций этого пакета приведем их определения: Ьцйб(ьо1) — конструирует улучшенную форму решения, полученного функцией рбьо!че; саьеьрй1(ьуь, о1, 02,...) — преобразует форму дифференциального уравнения; сйагь1пр(РОЕ 1) — находит характеристическую последовательность, дающую дифференциальное уравнение первого порядка; 6спапде(1г,ехрг,о1,02,...) — выполняет замену переменных в математических выражениях или функциях; 6соеП(ехрг,у(х)) — возвращает коэффициенты полиномиала дифференциального уравнения; 6ес(аге(ехрг) и др. — задает функцию для компактного ее отображения; 61Пог6ег(а,х) — возвращает порядок дифференциала в алгебраическом выражении а; 6ро!у(опп(ьуь,по Еп,ор1ь) — возвращает пслиномиальную форму для заданной системы ьуь не полиномиальных дифференциальных уравнений; 6ьцЬь(6епч1=а,...,ехрг) — выполняет дифференциальные подстановки в выражение ехрг; барбе(РОЕ,1п1о,() — создает карту РОЕ в различных форматах 1пго с опциональным заданием имени неизвестной функции Г; ьерагаЬППу(РОЕ, Г(х,у,...), "') — определяет условия разделения для сумм или произведений РОЕ; ьр!йь1пр(РОЕ,() — разделяет характеристическую последовательность на несоединенные поднаборы; ьрй1ьуь(ьуь,(цпсь) — разделяет наборы уравнений (алгебраические и дифференциальные) на несоединенные поднаборы; цп6ес1аге(1(х),...) и др.

— отменяет задание функции для компактного ее отображения. 474 Глава 7. Реигение дифференциальнвгх уравнений 7.8.3. Примеры решения дифференциальных уравнений с частными производными Примеры решения лиффереициальных уравнений и систем с частными производными представлены ниже (файл рг)е): > гевеагГ! чЙГЬ(РОЕГоо1в] ! > РРЕ у= х*сИЙЙ(Й(х, у), у) -сИЙЙ(Й(х, у),х) =Й(х, у) "2*5 [х] IЬ(у) ! РОЕ:= л1 — Г(х,у)~ -< — Г(х,у)~ = ( д ) ( д ) Г(х„у) а(х) [,ду ',) !,ду ' 3 )](у) > аов г= рг)во1че(РПЕ) анв:= Г(х.

у)— < а х )! -=+у+— 2 2 > РРЕ:= 5 [х,у) *гИЙЙ(5[х, у), у,х) + г]гЙЙ (5 (х, у), х) *гИЙЙ (5 (х, у), у] = 1г РОЕ:= Я(х,у' 3(х,у) + — о(х,у) — Б(х,у) = ! [ дудх ~ [,ду Лду > вегас г= рг]во1че(РВЕ,Н1НТ Й(х]*9(у]) в(гас: = (Б(х, у) = Г(х) а(у)) д[ю(]еге с г( с! г( ! ! ( — Г(х) = =', — 5(у) =— с(х Г(х) 4у 2 6(у) с, ~ > Ьог1с](вггчс)! ,]2, + гу,~,у гу $(х, у)— > раво1че (РОЕ, Н1НТ=Р (х, у) " (1/2) ] Б(х,у) = > РОЕ г= с]гЙЙ(Й(хуУуг),х) + гИЙЙ(Й(х,У,г),У) "2 Й(х,У,г)+г; РОЕ:= — Г(л,у,г) + — Г(х,у,Й) =Г(х,у,2)+2 ' <,дх ' ' ~ <.дх > рг]во1че (РОЕ, Нтнт=веглр] ! р( г)= СЗе-*, р,( ю)= С4е-', г( в)= С5, х( в)= в+ Сб, (2 у) у( г)=2 С4е-'+ С2,Г( в)= СЗе-' — С5+е СГ) ), д д агап ( р = — Г(х,у,г), р = — Г(х,у,г)) 7.8.

Характеристики

Тип файла
DJVU-файл
Размер
14,75 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Как Вы думаете, сколько людей до Вас делали точно такое же задание? 99% студентов выполняют точно такие же задания, как и их предшественники год назад. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6390
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее