Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Справочный материал к контрольным, домашним работам.doc
Скачиваний:
6
Добавлен:
20.05.2014
Размер:
987.14 Кб
Скачать

Домашняя работа по теме «числовые ряды» (2 семестр)

1. Необходимый признак сходимости. Если ряд ап сходится, то ап0. Обратное неверно – если ап0 , то о сходимости ряда ничего сказать нельзя. Следствие (достаточный признак расходимости). Если ап не стремится к нулю, то ряд расходится. Поэтому исследование ряда на сходимость в первой задаче производится по программе: «Если ап0 , то нет ответа, в прочих случаях – расходимость».

Итак, для данного ряда задаем общий член ряда

> a[n]:=(2*n+1)/(3*n+4)*sin(n);

Находим предел этой числовой последовательности

> L:=limit(a[n],n=infinity);

Применяем программу

> if L=0 then print ('не*знаю') else print ('pacx') end if;

2. Предуведомление. В условии второй задачи имеется требование «Применяя предельный признак сравнения, упростить ряды», т.е. сохранить в них главные слагаемые. Это требование относится к докомпьютерной эпохе, и выполнять его не обязательно.

Признак Коши для рядов с положительными членами. Пусть существует предел . Тогда рядап сходится, если L<1, и расходится, если L>1. При L=1 признак ответа не дает. Признак удобно применять к рядам, общий член которых представляет собой выражение в некоторой степени.

Пример. Пусть общий член ряда имеет вид

> a[n]:=((2*n+1)/(3*n+2))^sqrt(n^2+2);

Находим предел выражения

> L:=limit(a[n]^(1/n),n=infinity);

Составляем программу, соответствующую признаку Коши

> if L<1 then print ('cx') elif L>1 then print ('pacx') else print ('не*знаю') end if;

Признак Даламбера для рядов с положительными членами. Пусть существует предел . Тогда рядап сходится, если L<1, и расходится, если L>1. При L=1 признак ответа не дает. Признак применяют к рядам, общий член которых содержит экспоненты.

Пример. Пусть общий член ряда имеет вид

> a:=n->n*3^n/sqrt(2^n+5^n);

Вычисляем предел выражения и переводим его в десятичный вид

> L:=evalf(limit(a(n+1)/a(n),n=infinity));

Составляем программу, соответствующую признаку Даламбера

> if L<1 then print ('cx') elif L>1 then print ('pacx') else print ('не*знаю') end if;

Примечание. MAPLE не «умеет» сравнивать числа вида и 1. Он может срав­нивать только рациональные или десятичные числа. Поэтому вычисление предела следует применять совместно с командой> evalf(…);

3. Первые три примера этой задачи сводятся при помощи предельного признака сравнения (сохранение главных слагаемых, соотношения эквивалентности) к ряду Дирихле , который сходится приp>1 и расходится в остальных случаях. Исходный ряд ведет себя в смысле сходимости так же, как эквивалентный ему ряд Дирихле.

Последний пример при помощи предельного признака сравнения сводится к ряду , который должен быть исследован при помощи интегрального при­знака Коши: Ряд ведет себя в смысле сходимости так же, как.

4. Исследовать оба знакочередующихся ряда на абсолютную сходимость. Для этого рассматриваем ряды, составленные из модулей. Это ряды с положительными членами, и к ним можно применять признаки, рассмотренные в предыдущих задачах. Если сходится ряд, составленный из модулей, то сам ряд также сходится. Он называется абсолютно сходящимся. Если же ряд из модулей расходится, то про сам ряд ничего сказать нельзя. К нему нужно применять признак Лейбница: «Если |an| убывает, начиная с некоторого номера (его производная по п начиная с некоторого номера отрицательна), и стремится к нулю, то ряд сходится».

Для вычисления суммы S ряда берут некоторое, достаточно большое, число слагаемых и вычисляют частичную суммуSn. При этом, конечно же, допускается ошибка|S – Sn| , которая не превосходит абсолютной величины первого отброшенного слагаемого. В данной задаче следует для каждого из рядов вычислить сумму, взяв 1000 слагаемых, оценить ошибку (она|a1001| ). По этой ошибке следует выяснить, сколько верных знаков содержит результат.

Пример. Вычислить сумму ряда .

Записываем частичную сумму S1000

> S[1000]:=Sum(n*(-1)^n/(n^2+1),n=1..1000);

Вычисляем её значение

> evalf(S[1000]);

Оцениваем ошибку

> delta<evalf(abs(1001*(-1)^1001/(1001^2+1)));

Ответ: Сумма ряда S  – 0, 27 . Все знаки верные.