Главная » Просмотр файлов » Гурский Д., Турбина Е. - Вычисления в MathCad 12

Гурский Д., Турбина Е. - Вычисления в MathCad 12 (1077322), страница 93

Файл №1077322 Гурский Д., Турбина Е. - Вычисления в MathCad 12 (Гурский Д., Турбина Е. - Вычисления в MathCad 12) 93 страницаГурский Д., Турбина Е. - Вычисления в MathCad 12 (1077322) страница 932018-01-10СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

12.2). Посмотрите, на каком широком промежутке полином ведет себя точно так же, как синус. Удивительно, не правда ли?Лsm(x)Рис. 12.2. Кривая синуса и график аппроксимирующего ее полиномаПомимо приблизительного нахождения значений функций в неаналитических точках,ряды Тейлора активно используются при вычислении пределов. Данный подход хорош12.6. Разложение функций в ряды Фурье • 4 1 1своей универсальностью и простотой, поэтому именно его использует символьныйпроцессор Mathcad при подсчете пределов. Приведем и мы пример такого применениярядов Тейлора.Пример 1 2 . 1 1 . Нахождение пределов с использованием разложения в рядТейлораПусть стоит задача без применения оператора Km найти следующий предел:-х,.е- cos(x)hmх-» °x3-sin(x)Чтобы найти данный предел, нужно раскрыть неопределенность -1/0, что довольно непростосделать посредством стандартных средств.

Однако если разложить функцию в ряд Тейлора, задача будет решена довольно просто:-хе "- cos(x).1series,х,9 ->1212х +294х3 .' .18030240х •sin(x)Полагая х равным 0, находим, что предел будет равен 1/12. Этот же ответ будет выдан и оператором lim Mathcad:-х2е- cos(x)I— -> —hmхx3-sin(x)~> °12Приведем примеры расчета с использованием ряда Тейлора еще двух пределов:11yc-(sin(x)-x)-1/x-(sin(x)-x)..Xseries,x,54cos(x) + —limcos(x) нх-> 0•i —1Xhm v + x /x-> 01,X-> e22(е2}2fу +x /,' "2series, x, 3 -> e12.6.

Разложение функций в ряды ФурьеРазложение в ряд Тейлора периодических функций — это не очень эффективный ход,так как при этом хороший уровень аппроксимации будет достигнут лишь в окрестности точки разложения. При отдалении же от нее приближение будет плохим, хотя участки412 •Глава 12. Ряды и пределыкривой будут идентичны участку, на котором проводилось разложение. Очевидно, чтодля того, чтобы хорошо приближать периодическую функцию, аппроксимирующийполином сам должен быть периодической функцией с таким же периодом, как у приближаемой функции.

Чтобы удовлетворить этому условию, нужно разложить функцию не по степеням х, а по тригонометрическим функциям. Полученный при этом рядназывается рядом Фурье. Он имеет следующий общий вид (для функций с периодом 2л):°°а0f(x) = — + \*(an-cos(n-x) + bn-sin(n-x))n=0•В математическом анализе доказано, что погрешность приближения функции рядомФурье минимальна тогда и только тогда, когда коэффициенты а0, b n и ап определяютсяследующим образом:f(x)dxan = - -1Гf(x)-cos(n-x)dxn JnJ-n- nВ Mathcad нет оператора наподобие series, с помощью которого можно было бы производить разложение функций в ряд Фурье. Одиако это не значит, что данная проблемане решаема.

Используя вычислительные возможности программы, с ней вполне можно справиться. Как это сделать, показано в примере 12.12.Пример 12.12. Разложение функции в ряд ФурьеПусть стоит задача — разложить на промежутке от -л до л в ряд Фурье функцию f(x)=x. Для этогосоздадим функцию Fourier(x, k), основываясь на приведенных выше формулах:fi»:=xгпFourier(x,k) : = —fl[x)dxH}Лf(x)cos(nx)dxcos(n-x)(n-x) ++ fljx)-sin(nx) dx-sin(n-x))Чтобы получить ряд из к членов в символьном виде, функцию Fourier следует вычислить аналитически.

Найдем девять членов разложения f(x)=x в ряд Фурье:z(x) := Fourier(x, 9) factor —> 2-sin(x) — sin(2x) нsin(3x)sin(4-x) ...+ -sin(5x)sin(6x) + --sin(7-x)sin(8-x) + --sin(9-x)53749Тригонометрические ряды сходятся медленно, а полученный — вообще на уровне знакопеременного гармонического. Поэтому точность приближения им функции при суммировании девятичленов разложения будет невысока:z ( l ) = 0.988z ( 2 ) = 2.044То, что девяти членов разложения явно недостаточно для качественной аппроксимации, можноувидеть, построив график z(x) (рис. 12.3).Что же делать, если функцию нужно аппроксимировать с высокой точностью? Идею получитьнеобходимый для этого ряд в явном виде отбросим сразу (оперировать выражением из десятков,а то и сотен элементов очень сложно).

Можно попробовать просчитать функцию Fourier численно. Посмотрим, какую точность обеспечит суммирование 20, 50 и 100 членов разложения:Fourier(-l,20) = -0.9447022Fourier(-1,50) = -0.9946389Fourier(-1,100) = -1.000320221412.6. Разложение функций в ряды Фурье-5* 413хРис. 12.3. Из-за медленной сходимости ряда девять членов — это слишком мало для хорошегоприближения функции (обратите внимание на колебания, которые совершает кривая)В общем, если функцию нужно приблизить с точностью до 0.001-0.0001, численный расчет вполне приемлем.

Однако более высокой точности он обеспечить не может, что связано с погрешностью численных методов интегрирования. Ее, конечно, можно уменьшить, присвоив TOL минимальное значение. Однако при этом численные алгоритмы интегрирования перестанут сходитьсяуже при п порядка 50.Если стоит задача аппроксимировать функцию рядом Фурье максимально точно, нужно вывести формулу общего члена и использовать ее в сочетании с оператором суммы. При этом можнобудет суммировать тысячи и даже миллионы членов разложения буквально за секунды. В нашемслучае найти формулу общего члена будет очень просто, если вынести за скобку 2:sin(x)-^p + ^ M - ^+^. . .

. + ( — 1)sin(n-x)На основании данной формулы создадим функцию, которая позволит суммировать произвольное количество членов ряда:(-1)-sin(nx)п=0Проверим, какую точность аппроксимации обеспечит суммирование 1000,10 000 и 10 000 000 (!)членов ряда:F(-l, 1000)= -0.998866428398192F(-l, 10000)= -1.00008257593928F(-l, 10000000)= -1.00000000750972Суммирование 1000 членов ряда Фурье для f(x)=x дает достаточную точность, чтобы можнобыло построить «правильный» график (рис.

12.4).F(x.lOOO)^Рис. 12.4. «Периодическая линейная» функцияГлава 13. Исследование функцийи оптимизацияИсследование функций различной степени сложности и построение их графиков является важнейшей задачей математического анализа, с которой приходится сталкиваться любому старшекласснику и студенту естественной специальности вуза.В основе методов исследования поведения функций как одной, так и многих переменных лежит вычисление производных (первого и второго порядков), пределов, а также(для функций нескольких переменных) решение систем уравнений.

Изучив матершшпредыдущих глав, вы наверняка убедились в том, что с подобными операциями Mathcadпод вашим чутким руководством справляется весьма эффективно, поэтому длительную и, в общем-то, не требующую аналитического подхода процедуру исследованиянекоторой функции вы сможете провести с его помощью, затратив всего пару минут.Благодаря удивительным графическим возможностям системы вы без труда построите кривую или поверхность любой сложности, описывающую анализируемую функцию — график наглядно отражает особенности ее поведения, что может значительноупростить проведение исследования.13.1. Исследование функций однойпеременнойВ данном разделе мы рассмотрим наиболее часто встречающиеся в практикуме по высшей математике задачи, связанные с исследованием функций одной переменной: определением точек экстремума и перегиба, промежутков монотонности, выпуклостии вогнутости, нахождением асимптот. К сожалению, в Mathcad имеются встроенныесредства лишь для отыскания максимумов и минимумов функций (особенностям ихиспользования мы уделим внимание в соответствующем разделе), поэтому большинство этапов исследования функции в Mathcad придется проводить, используя классические формулы и определения.Пример 1 3 .

1 . Исследовать функцию3и построить ее график13.1. Исследование функций одной переменной * 4 1 5Поскольку алгоритм исследования функции включает в себя большое количество пунктов, разберемкаждый из них по отдельности, прибегая, в случае необходимости, к теоретическим пояснениям.Для облегчения поставленной задачи построим в первую очередь график функции (рис.

13.1). Заметьте, при проведении исследования на бумаге такой возможности у вас нет.„-2хРис. 13.1. График исследуемой функцииВо-первых, найдем область определения и область значений функции. X может принимать любые значения (в выражении функции отсутствуют корни четной степени и знаменатель), поэтому функция определена на всей числовой оси. Очевидно, область значений функции принадлежит промежутку (-<»; <*>).Во-вторых, определим, является ли функция четной, нечетной или периодической.Видно, что график несимметричен относительно оси ординат и начала координат, значит этофункция общего вида.

Приведем аналитическое подтверждение.1Г2lf(a) -• |_(а + 3)-а J31Г2lf(-a) -» [(-а + 3)-а J13Г2l3-fl» -• -|_(а + 3)-а JКак видно из полученных выражений, условия четности (f(-a)=f(a)) и нечетности (f(-a)=-f(a))не соблюдаются, значит, это функция общего вида.Докажем непериодичность исследуемой функции. Функция является периодической, если длялюбого х существует такое действительное число Т, отличное от нуля, что f(x+T)=f(x). Доказательством периодичности/непериодичности функции является нахождение величины Т из уравнения f(x+T)"f(x). Если Т является действительным числом, функция периодична. В случае непериодичности Т будет являться функцией от х.123)) x2= V (x+ T + 3)(х+ Т) 2 solve,T -•22+22416* Глава 13. Исследование функций и оптимизацияПри решении уравнения f(x+T)=f(x) относительно Т мы получили три значения. По определениюTV0, остальные же выражения зависят от х, следовательно, функция является непериодической.В-третьих, точек разрыва функция не имеет, поскольку непрерывна на промежутке (-«>; то).В-четвертых, найдем точки пересечения графика функции с осями координат.|(х + 3)х solve,х ->ОГрафик функции пересекает ось Ох в точках (-3; 0), (0; 0).цо) = оГрафик функции пересекает ось Оу в точке (0; 0).В-пятых, определим точки экстремума, значения функции в этих точках, а также интервалы возрастания и убывания функции.Необходимым условием существования экстремума в точке х(1 является равенство нулю или несуществование (равенство бесконечности) производной в данной точке.

Достаточным условиемсуществования экстремума в точке х() является смена знака производной при переходе через точку х0.х+ 2— f(x) simplify -» хdxх+ 2[(x+3)-x2]solve,x —> - 2[(x+3)-x2]Производная равна нулю при х=-2, при х=-3 и х=0 она обращается в бесконечность. Значит,данные критические точки являются точками перегиба или экстремума.

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

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

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

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