- 10•x•e 2x = 10*x*exp(2*x)
- x•e -x +cos(3x) = x*exp(-x)+cos(3*x)
- x 3 -x 2 +3 = x^3-x^2+3
- Выражение 0.9*x=sin(x)+1 необходимо преобразовать к виду: sin(x)+1-0.9*x . Аналогично, x^2-7=5-3x к виду x^2+3x-12 .
Пусть дано уравнение f(x)=0 , где f(x) определено и непрерывно в некотором конечном или бесконечном интервале a ≤ x ≤ b . Всякое значение ξ, обращающее функцию f(x) в нуль, то есть такое, что f(ξ)=0 называется корнем уравнения или нулем функции f(x) . Число ξ называется корнем k -ой кратности, если при x = ξ вместе с функцией f(x) обращаются в нуль ее производные до (k-1) порядка включительно: f(ξ)=f’(ξ)= … =f k-1 (ξ) = 0 . Однократный корень называется простым.
Приближенное нахождение корней уравнения складывается из двух этапов:
- Отделение корней, то есть установление интервалов [αi,βi] , в которых содержится один корень уравнения.
- f(a)•f(b) , т.е. значения функции на его концах имеют противоположные знаки.
- f’(x) сохраняет постоянный знак, т.е. функция монотонна (эти два условия достаточны, но НЕ необходимы) для единственности корня на искомом отрезке).
- f”(x) сохраняет постоянный знак, т.е. функция выпукла вверх, либо – вниз.
- Уточнение приближенных корней, то есть доведение их до заданной точности.
Геометрическая интерпретация метода Ньютона (метод касательных)
Критерий завершения итерационного процесса имеет вид
Локализировать наименьший положительный корень уравнения
Вопросы для самоподготовки
1. Что значит найти корень уравнения с точностью ε ?
2. Каковы этапы приближенного решения нелинейных уравнений? Какова цель каждого этапа?
3. Теорема о существовании и единственности корня на отрезке. Аналитическое и графическое отделение корней.
4. Метод половинного деления (алгоритм, геометрическая иллюстрация, условие окончания вычислений).
5. Метод хорд (алгоритм, геометрическая иллюстрация, условие окончания вычислений).
6. Метод касательных (условия применимости, алгоритм, геометрическая иллюстрация, условие окончания вычислений).
7. Комбинированный метод (условия применимости, алгоритм, геометрическая иллюстрация, условие окончания вычислений).
8. Метод итераций (алгоритм, геометрическая иллюстрация, условие окончания вычислений, достаточное условие сходимости итерационного процесса).
9. Сравнительная оценка методов уточнения корней.
Пример выполнения работы. Найти наименьший положительный корень уравнения
Найти наименьший положительный корень уравнения
.
1. Область определения функции ее производные
2. Строим графики функций: находим точки пересечения графиков. Из рис. 2.6 видно, что наименьший положительный корень данного уравнения лежит внутри отрезка . Проверим аналитически, что корень отделен на этом отрезке. Вычисляем: Поскольку и функция непрерывна, то в силу теоремы 1 внутри отрезка имеются корни. Поскольку и для всех , следовательно для всех , т. е. функция возрастающая на . Поэтому на основании теоремы 2 внутри этого отрезка имеется один корень уравнения и он может быть взят в качестве начального.
Рис. 2.6. Графики функций и
3. С помощью микрокалькулятора делаем 3 шага методом половинного деления; результаты заносим в табл. 2.2.
Таблица 2.2
Уточнение начального отрезка методом половинного деления
N | |||||||
0,5 | 1,0 | 1,5 | 1,0 | – 0,974 | – 0,386 | 1,766 | |
1,0 | 1,25 | 1,5 | 0,5 | – 0,386 | 0,449 | 1,766 | |
1,0 | 1,125 | 1,25 | 0,25 | – 0,386 | – 0,023 | 0,449 | |
1,125 | 1,1875 | 1,25 | 0,125 |
В результате получаем: уточненный отрезок [1,125; 1,250]; приближенное значение корня ; погрешность корня равной . При этом были введены следующие дополнительные обозначения: .
Дальнейшее уточнение корня проводим комбинированным методом. Так как , то левый конец отрезка [1,125; 1,250] уточняем методом хорд, а правый – методом Ньютона. Поэтому используем формулы (2.13). Результаты вычислений заносим в табл. 2.3.
Таблица 2.3
Уточнение корня комбинированным методом
N | ||||||
1,125 | 1,250 | 0,125 | – 0,023092 | 0,444045 | 4,243056 | |
1,131114 | 1,144170 | 0,013056 | – 0,002622 | 0,041961 | 3,460994 | |
1,131882 | 1,132046 | 0,000164 |
Так как , вычисления прекращаем на втором шаге. Находим корень уравнения
4. Продолжаем выполнение работы в компьютерном классе. Запускаем программу Mathcad. Открываем файл Lab2.mcd. Вводим функцию
Строим график функции на найденном начальном интервале [0,5;1,5] (рис. 2.7)
Рис. 2.7. График функции f(x)
Характеристики графика свидетельствуют, что функция непрерывна, и существуют и знакопостоянны на этом отрезке (т. е. функция монотонна и не меняет направление выпуклости) и что корень уравнения лежит в интервале [0,5;1,5], причем единственный. Таким образом, можем применить все вышеперечисленные методы. После этого находим корень с точностью до с помощью встроенной функции системы Mathcad
5. Выписываем точное решение и сравниваем полученные результаты ручного и машинного счета. Определяем погрешность:
6. Определяем с помощью компьютера значение корня методом половинного деления с точностью . По формуле, чтобы удовлетворить погрешности , для начального отрезка единичной длины необходимо провести шагов.
Выписываем автоматически вычисленное по этой формуле в соответствующем разделе количество шагов и таблицу 2.4, содержащую первые и последние три шага получившейся матрицы приближений корня методом половинного деления.
Таблица 2.4
Отыскание корня методом половинного деления
N | … | |||||||
0,5 | 1,0 | 1,0 | 1,125 | … | 1,131836 | 1,131866 | 1,131882 | |
1,5 | 1,5 | 1,25 | 1,25 | … | 1,131897 | 1,131897 | 1,131897 |
Получим корень , абсолютная погрешность которого .
7. Получим на компьютере значение корня методом Ньютона с точностью . Для этого вводим в начале соответсвующего раздела и с помощью запрограммированной формулы (2.7) определяем погрешность . Следовательно, увеличивая N на единицу, вводим в программу и получаем . То есть требуемая точность достигнута (если это не так продолжаем увеличивать N на единицу).
. Выписываем получившуюся таблицу 2.5 для .
Таблица 2.5
Отыскание корня методом Ньютона
N | |||||
1,5 | 1,221958726 | 1,139300286 | 1,131948438 | 1,131892063 |
Получим приближенный корень , абсолютная погрешность которого .
8. Вычисляем на компьютере значение корня методом хорд с точностью . Для этого вводим в начале соответсвующего раздела и с помощью запрограммированной формулы (2.10) получим и оценим погрешность . Следовательно, увеличиваем N на единицу, вводим , для которого . Увеличиваем N еще на единицу, вводим , для которого . Т. е. требуемая точность достигнута (если это не так, продолжаем увеличивать N на единицу).
Выписываем первые и последние два шага из получившейся таблицы для .
Таблица 2.6
Отыскание корня методом хорд
N | … | |||||
0,5 | 0,855440054 | 1,035664929 | … | 1,131891898 | 1,131892012 |
Получим корень , абсолютная погрешность которого .
9. Вычисляем на компьютере значение корня комбинированным методом с точностью . Для этого вводим в начале соответсвующего раздела и получаем по формуле (2.13) погрешность . Следовательно, увеличиваем N на единицу и вводим , для которого . То есть требуемая точность достигнута (если это не так, продолжаем увеличивать N на единицу).
Выписываем получившуюся таблицу 2.7 для .
Таблица 2.7
Отыскание корня комбинированным методом
N | |||||
0,5 | 0,8554400542 | 1,1020813008 | 1,1316586589 | 1,1318920464 | |
1,5 | 1,2219587264 | 1,1393002857 | 1,1319483820 | 1,1318920634 |
Получим корень , абсолютная погрешность которого .
10. Все расчеты оформляются в виде отчета по лабораторной работе.
Вопросы для самоконтроля
1. Уравнение какого типа решается в данной работе?
2. Что называется корнем уравнения ?
3. Как графически решить уравнения ?
4. Перечислите достоинства и недостатки графического метода.
5. В чем состоит этап отделения корней уравнения ?
6. Сколько корней должна иметь функция на начальном отрезке ?
7. Как определить аналитически: возрастает или убывает функция на промежутке?
8. Как определить аналитически: выпукла или вогнута функция на промежутке?
9. Какие условия, наложенные на , гарантируют наличие хотя бы одного корня уравнения на начальном отрезке ?
10. Какие условия, наложенные на , гарантируют наличие не более одного корня уравнения на начальном отрезке ?
11. Привести алгоритм решения уравнения методом половинного деления. Какие условия при этом должны быть наложены на функцию ?
12. Какие условия должны быть наложены на , чтобы уравнение можно было решить методом Ньютона?
13. Как выбирается начальная точка в методе Ньютона?
14. Вывести формулу для вычисления последовательных приближений методом Ньютона, записать формулу оценки погрешности.
15. Какие условия должны быть наложены на , чтобы уравнение можно было решить методом хорд?
16. Как выбирается начальная точка в методе хорд?
17. Вывести формулу для вычисления n последовательных приближений методом хорд, записать формулу оценки погрешности.
18. Какие условия должны быть наложены на , чтобы уравнение можно было решить комбинированным методом?
19. Выписать формулы, по которым уточняются концы начального отрезка комбинированным методом. В зависимости от каких условий осуществляется выбор формул?
20. Указать условие, по которому процесс уточнения отрезка комбинированным методом должен быть прерван? Как затем найти корень?
источники:
http://dit.isuct.ru/IVT/sitanov/Literatura/M869/Pages/Glava1.htm
http://megaobuchalka.ru/7/29269.html
ЛАБОРАТОРНАЯ
РАБОТА №8.
вариант
№ 4
«АНАЛИЗ МОДЕЛЕЙ
В ВИДЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙ»
Задание:
Локализовать
наименьший положительный корень
уравнения 4x
5ln x
5
0 и уточнить его значение с точностью
до ε=0,0001 двумя методами.
1.
методом половинного деления
2.
методом комбинированным: хорд и
касательных.
Теоретическое
введение
-
Уточнение
корней методом половинного деления
Пусть корень
уравнения
отделен на отрезке [a,
b],
т.е. f(a)f(b)<0
и f
’(x)
сохраняет знак (рис. 2.6.).
В качестве начального
приближения корня возьмем точку c0
– середину отрезка:
.
Если f(с0)=0,
то c0
– искомый корень уравнения, если
,
то из двух отрезков [a,
c0]
и [c0,
b]
выбираем тот, на концах которого функция
принимает значение разных знаков.
Новый отрезок
опять делим пополам и далее поступаем
аналогично вышеизложенному. Длина
каждого нового отрезка вдвое меньше
длины предыдущего отрезка, т.е. за n
шагов сократится в 2n
раз.
Вычисления
прекращаем, если длина отрезка
станет меньше заданной погрешности
,
т.е.
.
-
Комбинированный
метод хорд и касательных
Методы хорд и
касательных дают приближения корня с
разных сторон. Поэтому их часто применяют
в сочетании друг с другом, тогда уточнение
корня происходит быстрее.
Пусть дано уравнение
f(x)=0,
корень отделен на отрезке [a,
b].
Рассмотрим случай,
когда f
‘(x)
f
’’(x)>0
(рис. 2.13)
В
этом случае метод хорд дает приближенное
значение корня с недостатком (конец b
неподвижен), а метод касательных – с
избытком (за начальное приближение
берем точку b).
Тогда
вычисления следует проводить по формулам:
;
.
Теперь
корень ξ
заключен в интервале [a1,
b1].
Применяя
к этому отрезку комбинированный метод,
получим:
;
и т.д.
;
(2.6)
Если же f
‘(x)
f
’’(x)<0
(рис. 2.14), то рассуждая аналогично, получим
следующие формулы для уточнения корня
уравнения:
;
.
Вычислительный
процесс прекращается, как только
.
Выполнение
работы
-
Локализация корня
Определим графически
корень уравнения, для этого построим
график функции y=4x
5ln x
5
Из
графика следует , что корень находится
на отрезке [0,5;1]
2.
Метод половинного деления
1-ая
итерация
А=
0,5
Б=1
>0.0001
cледовотельно
точность не достигнута, и необходимо
продолжить расчеты
так
как f(a)*F(c)<0
то следующий интервал для расчетов
(0,5;0,75)
2-я
итерация
А=
0,5
Б=0,75
>0.0001
cледовотельно
точность не достигнута, и необходимо
продолжить расчеты
так
как f(a)*F(c)<0
то следующий интервал для расчетов
(0,5;0,625)
3-я
итерация
А=
0,5
Б=0,625
>0.0001
=>точность не достигнута, и необходимо
продолжить расчеты
так
как f(a)*F(c)<0
то следующий интервал для расчетов
(0,5;0,5625)
Ответ:
х=0,58960,0001
-
Метод комбинированный:
хорд и касательных.
(метод касательных)
(метод хорд)
Ответ:x=0,58960,0001
Соседние файлы в предмете Информатика
- #
- #
- #
- #
- #
- #
- #
Ответы с готовыми решениями:
Найти корень уравнения
Помогите пожалуйста найти корень уравнения в excel x^3-0.3*x^2-1.2=0 на отрезке .
Найти корень уравнения
Найти корень уравнения:
Y=(23x-x^3)/(cos(2x+6)) Через подбор параметра
найти корень уравнения
Найти приблизительное значение корня уравнения используя Эксэль. Описать действия решения….
Найти второй корень уравнения метода хорд excel
Уравнение такое — 3^(x−1)−2−x=0
Один я нашёл, второй не получается.
Помогите исправить.
(Лист 2).
Лабораторная работа № 1.8. Решение нелинейных уравнений заданным методом
(4 – 7 балла)
1.Цель работы
получить представление об итерационных методах определения корней нелинейного скалярного уравнения;
научиться использовать электронные таблицы и средства Excel для определения интервалов существования корней скалярного уравнения и последующего их вычисления с заданной точностью.
2.Необходимые программные и технические средства
- Персональный компьютер.
- Тип операционной системы – Windows XP и выше.
- MS Office версии 97-2003 и выше.
^
3.Общие сведения
Разнообразные проблемы механики, физики, техники сводятся к вопросу о нахождении корней многочлена, причем, иногда достаточно высоких степеней. Точные решения известны для квадратных уравнений, кубических (формула Кардано) и уравнений 4-й степени (метод Феррари). Для уравнений выше 5-й степени не существует формул для выражения корней многочлена. Однако в технических приложениях обычно достаточно знать лишь приближенные значения корней с некоторой заранее заданной точностью. В общем же случае надежд на простое аналитическое решение нет. Более того, доказано, что даже алгебраическое уравнение выше четвертой степени неразрешимо в элементарных функциях. Поэтому решение уравнения проводят численно в два этапа (здесь разговор идет лишь о вещественных корнях уравнения). На первом этапе производится отделение корней – поиск интервалов, в которых содержится только по одному корню. Второй этап решения связан с уточнением корня в выбранном интервале (определением значения корня с заданной точностью).
В общем виде уравнение n-й степени выглядит следующим образом:
где n − некоторое положительное число,
− произвольные числа, причем старший коэффициент должен быть не равен нулю.
Выражение
называется многочленом (полиномом) n
-й степени от неизвестного x
.
Если при некотором x
= x
0 выполняется равенство
, то x
0 называется корнем многочлена .
4.Задание
Задано уравнение f(x)=0. Требуется найти все его корни тремя способами:
1. найти корень с погрешностью eps=0,0001 методом половинного деления (дихотомии) — локализовать один корень уравнения табличным методом и построить график функции в области этого корня;
2. найти корень с помощью инструмента «Подбор параметра»;
3. найти корень с помощью инструмента «Поиск решения».
Варианты заданий:
- х 6 +2х 5 +10х 3 -9х 2 +15х-17,5=0
- х 5 -2,8х 4 +3х 3 -3х 2 +4,4х-5=0
- х 6 +6,5х 5 -14х 4 +14х 3 -17х 2 +21х-22,5=0
- х 6 +10,5х 5 -24х 4 +28х 3 -29х 2 +39х-45=0
- х 5 -1,8х 4 -1,9х 3 -2,3х 2 +2,8х-3=0
- х 6 +10,5х 5 -18х 4 +22х 3 -17х 2 +31х-37,5=0
- х 5 -3х 4 +3,2х 3 -3,5х 2 +4,6х-5=0
- х 6 +7,5х 5 -18х 4 +20х 3 -11х 2 +19х-22,5=0
- х 5 -2х 4 +2,9х 3 -2,44х 2 +4,2х-5=0
- х 6 +9х 5 -18х 4 +19х 3 -19х 2 +30х-35=0
- х 5 -2,6х 4 +2,82х 3 -3,41х 2 +4,12х-3,23=0
- х 6 +6,5х 5 -20х 4 +21х 3 -21х 2 +31х-32,5=0
- х 5 -4х 4 +4х 3 -4,33х 2 +6х-6,67=0
- х 6 +3,5х 5 -14х 4 +14х 3 -17х 2 +21х-22,5=0
- х 5 -1,6х 4 +2,5х 3 -2,7х 2 +3,6х-4=0
- х 6 +8,5х 5 -16х 4 +19х 3 -15х 2 +27х-32,5=0
- х 6 +4,5х 5 -18х 4 +22х 3 -17х 2 +31х-37,5=0
- х 5 -2х 4 +2,09х 3 -2,52х 2 +3х-3,26=0
- х 6 +9,5х 5 -20х 4 +22х 3 -25х 2 +32х-35=0
- х 5 -2х 4 +2,25х 3 -2,58х 2 +3,25х-3,54=0
- х 4 -3х 3 +20х 2 +44х+54=0
- (cos(x)-3sin(x)) 2 -e x =0
- 2cos(x)+2x 2 =1
- ln(x+1)=x 2 +1+5cos(x) 2
- 3cos(x) 2 +2,3sin(x)=0,5ln(x-0.5)
^
5.Порядок выполнения
Прочитайте и уясните материалы разделов лекционного курса «Информатика», относящихся к теме работы.
Ознакомьтесь с общими сведениями о предмете лабораторной работы (см. выше в описании данной работы) и рекомендуемыми дополнительными материалами.
Уясните цель работы.
Подготовьте необходимые программные и технические средства (см. выше в описании данной работы).
Приступайте к выполнению работы:
Действительными корнями многочлена будут абсциссы точек пересечения его графика с осью Х
и только они.
Число положительных корней многочлена равно числу перемен знаков в системе коэффициентов этого многочлена (коэффициенты, равные нулю, не учитываются) или меньше этого числа на четное число.
Число отрицательных корней многочлена равно числу сохранения знаков в системе коэффициентов этого многочлена или меньше этого числа на четное число.
Если многочлен не имеет отрицательных коэффициентов, то многочлен не имеет положительных корней.
О
трезок
локализации всех корней многочлена определяется по выражению:
Для границы a формула справедлива если
Для отыскания корней многочлена с помощью электронной таблицы MS Excel необходимо выполнить следующие шаги:
Провести табулирование заданного многочлена на интервале .
Выявить интервалы локализации каждого корня многочлена (перемена знака в значении ). При необходимости, следует использовать табуляцию многочлена, неоднократно уменьшая шаг табуляции для более точных оценок.
После локализации корней произвести их уточнение.
При последующем уточнении корня на обнаруженном интервале не надейтесь никогда найти точное
значение и добиться обращения функции в нуль при использовании калькулятора или компьютера, где сами числа представлены ограниченным числом знаков. Здесь критерием может служить приемлемая абсолютная
или относительная погрешность
корня. Если корень близок к нулю, то лишь относительная погрешность даст необходимое число значащих цифр. Если же он весьма велик по абсолютной величине, то критерий абсолютной погрешности часто дает совершенно излишние верные цифры. Для функций, быстро изменяющихся в окрестности корня, может быть привлечен и критерий: абсолютная величина значения функции
не превышает заданной допустимой погрешности.
Пример 1
Найти все действительные корни уравнения:
f(x)
= х 5
+ 2х
4
+ 5х
3
+ 8х
2
– 7х – 3 = 0
, где а 5 = 1, а 4 = 2, а 3 = 5, а 2 = 8, а 1 = −7, а 0 = −3.
Число сохраненных знаков
= 4 (в уравнение отрицательных корней 4 или 2).
^
Число перемены знаков
= 1 (в уравнение один положительный корень).
О
пределяем отрезок , на котором существуют корни уравнения
Выполняем приближенное табулирование функции на отрезке [−9; 9] с шагом 1.
Определяем, что функция меняет знак на отрезке [−3; 1].
Производим табулирование функции на отрезке [−3; 1] с шагом 0,1.
Строим график функции.
Используя, таблицу и график функции определяем положение корней уравнения (на рис. 1. отрезки локализации корней выделены желтым цветом).
Из таблицы и графика видно, что многочлен f(x) содержит 3 корня, находящихся в границах отрезков: 1 корень [-2,1; -2]; 2 корень [-0.4; -0,3]; 3 корень .
^
Уточнение корней методом половинного деления (дихотомии)
Самым простейшим из методов уточнения корней является метод половинного деления
, или метод дихотомии
, предназначенный для нахождения корней уравнений, представленных в виде f(x)=
0.
Пусть непрерывная функция f(x)
на концах отрезка [a,b
] имеет значения разных знаков, т.е. f(a)×f(b)
≤ 0 (рис. 2), тогда на отрезке имеется хотя бы один корень.
Возьмем середину отрезка с=(a+b)
/ 2. Если f(a)×f(с)
≤ 0, то корень явно принадлежит отрезку от a
до (a+b
) / 2 и в противном случае от (a+b
) / 2 до b
.
Поэтому берем подходящий из этих отрезков, вычисляем значение функции в его середине и т.д. до тех пор, пока длина очередного отрезка не окажется меньше заданной предельной абсолютной погрешности (b-a
) ε.
Так как каждое очередное вычисление середины отрезка c
и значения функции f(c)
сужает интервал поиска вдвое, то при исходном отрезке [a,b
] и предельной погрешности ε
количество вычислений n
определяется условием (b-a
)/2n
ε, или n
~ log
2((b-a
)/ε
). Например, при исходном единичном интервале и точности порядка 6 знаков (ε ~ 10 -6) после десятичной точки достаточно провести 20 вычислений (итераций) значений функции.
С точки зрения машинной реализации этот метод наиболее прост и используется во многих стандартных программных средствах, хотя существуют и другие более эффективные по затратам времени методы.
Вычислительная процедура в Excel может быть реализована так
В ячейки вносим следующие формулы:
В ячейку А2 − а (левая граница интервала локализации корня);
В ячейку В2 − b (правая граница интервала локализации корня);
В ячейку С2 − =(А2+В2)/2;
В ячейку D2 − =f
(A2)*f
(C2);
В ячейку F2 − 0,0001 (абсолютная погрешность);
В ячейку А3 − =ЕСЛИ(D2
В ячейку B3 − =ЕСЛИ(D2
В ячейку D3 − =f
(A3)*f
(C3);
В ячейку Е3 − =ЕСЛИ(ABS(B3-A3)>$F$2;”продолжаем”;”конец”);
После этого выделяются ячейки А3:Е3 и автозаполнением
буксируются вниз до появления в столбце Е сообщения “конец”. Вычисленный корень с заданной точностью будет находиться в конце столбца F.
Вернемся к примеру, и с помощью метода половинного деления уточним значения корней в выделенных отрезках.
Первый корень находится внутри отрезка = [-2,1; -2] расположенного по адресу А2:В2. Заполняем рабочий лист формулами (рис. 4) и с заданной точностью 0,0001 определяем его значение (рис. 5). Ответ находится в ячейке С12 и равен X 1 = -2,073.
Границы отрезка второго корня находящегося внутри отрезка = [-0,4; -0,3] подставляем в таблицу по адресу А2:В2. Определяем его значение (рис. 6). Ответ находится в ячейке С12 и равен X 2 = -0,328.
Границы отрезка третьего корня находящегося внутри отрезка = подставляем в таблицу по адресу А2:В2. Определяем его значение (рис. 7). Ответ находится в ячейке С12 и равен X 3 = 0,7893.
Как и предполагалось, имеется три корня, два из которых отрицательные (Х 1 = -2,073; Х 2 = -0,32808; Х 3 = 0,789307).
^
Уточнение корней средством “Подбор параметра”
Обширную группу методов уточнения корня представляют итерационные методы
– методы последовательных приближений. Здесь в отличие от метода дихотомии задается не начальный интервал местонахождения корня, а его начальное приближение.
Когда желаемый результат вычислений по формуле известен (подстановка значения корня в уравнение делает его равным нулю), но неизвестны значения, необходимые для получения этого результата, можно воспользоваться средством Подбор
параметр
а.
Для этого выбирается команда Подбор
параметра
в меню Серви
с
. При подборе параметра MS Excel изменяет значение в одной конкретной ячейке до тех пор, пока вычисления по формуле, ссылающейся на эту ячейку, не дадут нужного результата.
Когда задаются условия для применения средства ^
Подбор параметра
, в одной ячейке обычно вводится формула, а переменная, которая используется в формуле (с некоторым стартовым значением), задана в другой ячейке.
В формуле можно применять больше одной переменной, но средство ^
Подбор параметра
позволяет работать только с одной переменной зараз. Для поиска решения в средстве Подбор параметра
применяется итеративный
алгоритм
. Это означает, что функция сначала проверяет заданное исходное значение параметра и проверяет, дает ли это значение нужный результат. Если исходное значение параметра не дает желаемого результата, средство перебирает другие значения, пока не будет найдено решение.
Поскольку поиск точного решения в некоторых задачах может занять много времени, поэтому MS Excel пытается найти компромисс, устанавливая определенные ограничения по точности решения или максимальному количеству итераций.
Средство ^
Подбор параметра
вызывается командой Сервис | Подбор параметра
(рис.8).
В окне диалога Подбор параметра
в поле Установить в ячейке
введем ссылку на ячейку с формулой, в поле Значение
− ожидаемый результат, в поле Изменяя значение ячейки
− ссылку на ячейку, в которой будет храниться значение подбираемого параметра (содержимое этой ячейки не может быть формулой).
Пример 2
Вычислить корень уравнения f(x) = -5х + 6 = 0
с помощью средства ^
Подбор параметра
В ячейку В2 введем любое число, например, 0.
В ячейку В3 введем формулу =-5*В2+6.
Вызовем диалоговое окно Подбор параметра и заполним соответствующие поля.
После нажатия на кнопку ^
ОК
Excel выведет окно диалога Результат подбора параметра.
Если подобранное значение необходимо сохранить, то нажмите на ОК
, и результат будет сохранен в ячейке, заданной ранее в поле Изменяя значения ячейки
.
Для восстановления значения, которое было в ячейке В2 до использования команды ^
Подбор параметра
, нажмите кнопку Отмена
.
Как видно из примера в ячейке B2 установилось точное значение корня уравнения
Х
= 1,2.
При подборе параметра Excel использует итерационный (циклический) процесс. Количество итераций и точность устанавливаются в меню Сервис | Параметры … |
вкладка Вычисления,
в которой задается Предельное число итераций
(по умолчанию 100) и Относительная погрешность
(по умолчанию 0,001).
Если Excel выполняет сложную задачу подбора параметра, можно нажать кнопку ^
Пауза
в окне диалога Результат подбора параметра
и прервать вычисление, а затем нажать кнопку Шаг
, чтобы выполнить очередную итерацию и посмотреть результат. При решении задачи в пошаговом режиме появляется кнопка Продолжить
− для возврата в обычный режим подбора параметра.
Пример 3
Возьмем в качестве примера все тоже квадратное уравнение
f
(x
)
= Х 5 + 2Х 4 + 5Х 3 + 8Х 2 − 7Х – 3 = 0.
Для нахождения корней уравнения с помощью средства ^
Подбор параметра
выполним следующие действия:
В таблице функции (рис.1) выявляем интервалы локализации корней уравнения (перемена знака в значении функции): первый интервал ячейки Е20:Е21, значение (-1,2698 и 3); второй интервал ячейки Е37:Е38, значение (0,80096 и -0,3012); третий интервал ячейки Е48:Е49, значение (-1,6167 и 0,22688);
В каждом интервале выбираем то значение функции, которое ближе к 0 и составляем пары ячеек «аргумент-значение»: первый корень D20:E20; второй корень D38:E38; третий корень D49:E49.
Уточняем значения корней средством ^
Подбор параметра
(рис. 10, 11, 12).
Рис. 10. Корень уравнения Х 1 = -2,073 |
Рис. 11. Корень уравнения Х 2 = -0,32804 |
Рис. 12. Корень уравнения Х 3 = 0,78934 |
Ответ: Х
1 = -2,073; Х
2 = -0,32804; Х
3 = 0,78934.
Значения корней уравнения, полученные приближением методом половинного деления: Х1 = -2,073; Х2 = -0,32808; Х3 = 0,789307.
Определение значения корней скалярного уравнения с заданной степенью точности с помощью инструмента ^
Поиск решения
В качестве примера возьмем тоже уравнение: f(x)
= Х 5
+ 2Х
4
+ 5Х
3
+ 8Х
2
− 7Х – 3 = 0.
Для более точного определения корня в каждом из выделенных диапазонов следует воспользоваться командой ^
Сервис | Поиск решения
. Для этого в ячейку, например, H8 введем формулу для вычисления f(x), а начальное приближение поместим в ячейку G8. Назовем их соответственно Целевая ячейка и Корень. В ячейку G8 введем первоначально значение, принадлежащее первому выделенному диапазону. Возьмем его на середине интервала равным –3,76 (можно эту ячейку оставить пустой). В ячейку H8 введем формулу =G8^5+2*G8^4+5*G8^3+8*G8^2-7*G8-3.
После выбора команды Сервис
| Поиск решения
появится диалог, в котором в поле Установить целевую ячейку
введем $H$8. Затем выберем кнопку Равной значению 0
.
В поле Изменяя ячейки
введем $G$8. В окно Ограничения
с помощью кнопки Добавить
следует указать диапазон поиска корня следующим образом:
- Для левой границы первого интервала –2,1 (оно находится в ячейке D20) $G$8 >= $D$20.
- Для правой границы первого интервала –2 (оно находится в ячейке D21) $G$8
На рис. 13 показан результат выполненных действий, описанных выше, а на рис. 14 диалог, появляющийся после нажатия кнопки Добавить
. Такой же диалог появляется при выборе кнопки Изменить
.
Выбор кнопки Параметры
приводит к появлению диалога (рис. 15), в котором можно задать параметры поиска.
Поле ^
Предельное число итераций
позволяет назначить число «циклов» поиска решения. Значения 100, принятого по умолчанию, достаточно для большинства задач.
Относительная погрешность обеспечивает назначение величины f зад в признаке достижения решения f к =(f k +1 – f k)/f k
Флажок ^
Линейная модель
используется, если задача является задачей линейного программирования. В нашем случае его устанавливать не надо.
Флажок Показывать результаты итераций
позволяет приостанавливать процесс поиска после каждой итерации для анализа процесса поиска. При этом выводится диалоговое окно Текущее состояние поиска
, выбор в котором кнопки Продолжить
позволяет выполнять следующую итерацию. Результаты, полученные на каждой итерации, выводятся в ячейке G8.
Выбор метода решения зависит от типа нелинейности.
Отметим, что задачи решения нелинейных уравнений и методы безусловной оптимизации тесно связаны. Поэтому после нажатия кнопки Выполнить
по окончании поиска появится сообщение, представленное на рис. 16.
Если в верхней части этого окна будет выведено сообщение ^
Р
ешение не найдено
, следует в ячейке H8 использовать формулу, вычисляющую либо |f(x)|, либо (f (x)) 2 . Затем в окне Поиск решения
(рис.13) выбрать переключатель Равной минимальному значению
.
С помощью диалогового окна ^
Результаты поиска решения
можно просмотреть отчеты трех типов: результаты, устойчивость, пределы. Отчеты каждого типа вызываются по следующему алгоритму:
- Курсор на тип вызываемого отчета.
- ОК. (На экране вызванный отчет на новом листе, на ярлычке которого указано название отчета).
- Курсор на ярлычок с названием отчета. (На экране вызванный отчет).
Поиск решения по двум остальным интервалам проведите самостоятельно по описанной выше схеме.
^
6.Оформление результатов
Лабораторная работа 1.8 требует оформления результатов по всем пунктам задания на листе под именем «18» в своей книге Excel «Л.р. по Excel».
^
7.Формулировка выводов
Достигнута ли цель работы?
Роль и возможности инструментов MS Excel для решения скалярного уравнения с заданной степенью точности.
^
Подбор параметров
.
Назначение и особенности инструмента Поиск решения
.
Особенности выполнения математических расчетов и задания целевой ячейки.
^
8.Порядок защиты
- Ответить на вопросы:
- Какое количество действительных корней имеет уравнение n степени?
- Что такое отрезок локализации корня?
- Что значит локализовать корень?
- В чём заключается идея решения уравнений методом деления отрезка пополам?
- Как можно оценить погрешность вычисления корня методом деления отрезка пополам?
- Как с помощью инструмента «Подбор параметра» найти значение корня?
- Уточнение корней методом половинного деления (дихотомии).
- Метод Подбор параметра
. - Метод Поиск решения
.
Вопрос: Нахождение корней уравнения методом деления отрезка пополам
Добрый день,что не так с 3-ьим корнем,никак не хочет выводится.Сверху — 3 корня через подбор параметра.Снизу — методом половинного деления. Округление 0,001 Уравнение x^3-2*x^2-x+2 Кто-нибудь может подправить или дать полезный совет,что не так?
Ответ:
furymaxim
, скобки пропустили
Вопрос: Дешифрование методом Плэйфера в MS Excel
Пожалуйста Подскажите как сделать дешифратор в EXCEL с помощью формул. Или скажите с помощью какой формулы можно с генерировать алфавит
Ответ:
В ячейку А1
Code | ||
|
И протянуть вниз
Вопрос: Тормозит файл-таблица Excel
Доброго времени суток, уважаемые коллеги!
Очень нужна ваша помощь, уже перепробовал все найденные и известные мне методы по уменьшению объема файла. Вроде бы все лишнее там вычистил.
Не смотря на это при работе с таблицей идут тормоза и подвисает, причем они переменны но стабильны (то тормозит, то не тормозит).
мне кажется, что возможно это из-за выпадающего списка с фотографиями, заметил что по мере увеличения выпадающих списков с фотками — увеличиваются и тормоза. Но странно, таблицы все маленькие, галлерея с фотками тоже не большая.
Ответ:
Проблему решил! Просто установив excel 2016 для Mac — вообще никаких тормозов, пока все работает отлично, но не уверен что дальше не столкнусь с этим опять!
Тем не менее, проблема актуальна, т.к. решение не через установку другой версии excel, возможно, кому-то еще пригодится
p.s. предыдущая версия excel была 2011 для Mac
Вопрос: Office 2007 как установить excel 2010
всем привет.
может название темы не совсем точно передает суть(((, но….
у меня win xp sp3 office 2007 Года и excel 2007года.
в excel то ли 2010 то ли 2013 есть функция диаграмм в виде карт стран или континентов powerview или как там точно. там еще карты бин используются.
есть ли какие надстройки для excel2007 чтобы такие диаграммы можно было. если нет, то в каком excel есть эта функция и есть ли возможность установить 2 excel на 1 комп. к примеру 2007 и 2010 на win xp sp3 если функция диаграмм с картами стран есть в 2010????
спасибо.
Ответ:
так а в 2010 excel это есть?? и если есть то как установить excel 2010 не удаляя мой офис 2007???
Добавлено через 3 часа 10 минут
ща смотрел похожие темы. нашел про libreoffice. программа такая как офис только бесплатная. мб у кого-нибудь есть карта Республики Беларусь для это программ????. там расширение geoOOo.
Вопрос: Получение выборки из Excel
Мне требуется на основе данных из Excel файла создать презентацию в PowerPoint.
Ни с тем, ни с другим раньше не работал. Поэтому проверьте алгоритм (наброски):
Получаю с помощью запросов необходимые выборки,
Связываю результаты выборок с шаблоном (пока не читал, как программно создаётся презентация)
Создаю, собственно, презентацию.
И всё это прописываю в макросе.
1. Последовательность правильная?
2. Как мне работать с данными полученными с помощью запросов? Записать их временно; результат каждого запроса на отдельном листе, а после создания файла-презентации закрыть БЕЗ ИЗМЕНЕНИЙ файл Excel? Или как-то по другому?
3. Как правильно написать подобный запрос?
Мой набросок не работает:
Запись результатов запроса из первого листа на второй.
4. Как этот запрос запустить
Что-то типа этого?
Добавлено через 2 часа 42 минуты
Или такое возможно только через временную БД Access?
Ответ:
В смысле сюда? На форум? — Пожалуйста… Дело то не в данных, а в запросах (способах обработки). В Access я это делаю, в Excel — не могу. Например подсчитать продажи для 3-х производителей с самыми большими продажами (ТОП 3), а остальных — суммировать. Я так понимаю — это не автоматизировать… Руками — Да, можно сделать.
Вопрос: Как добавить имена вложений Outlook в Excel с последующим сохранением их в указанной папке
Добрый день всем гуру Excel-я.
Благодаря этому форуму мне получилось наладить документооборот в Excel (точнее регистрацию входящих-исходящих писем) в более-менее автоматизированном виде.
В приложенном файле следующие основные макросы:
1. «Первое_MailSave» — прописывает письма из папки входящие Outlook
2. «Второе_в_шаблон» — выдает входящий номер и выводит данные в определенный шаблон (одобренный руководством в плане удобочитаемости)
3. «Завершение_Печать» — сохраняет лист шаблона в формате pdf в папке с входящим номером и пускает на печать.
Т.е. счастье есть, теперь полная обработка 10 писем занимает 3-4 минуты, а не 30-40.
Проблема с обработкой вложений:
1. Как не в ручную прописывать кол-во вложений
в письме, а автоматом с выводом в ячейку E4 листа «data» количества + 1 (само письмо)
2. Как в листе «Шаблон» в В5 перечислить все вложения по именам
3. Что добавить в макрос «Завершение_Печать», чтобы вложения сохранялись
в новосозданную папку с самим письмом.
Все данные из письма забираются, а вот с вложением так и не придумал как(см.код)
Код Visual Basic | ||
|
Поиски в интернете все ссылаются на макросы для outlook, но регистрация и создание необходимых директорий у меня происходит в excel, соответственно все переменные в нем же.
С одной стороны, у меня три разных вопроса, но, мне кажется, что оптимальнее будет реализовать все три вопроса в одном макросе.
С уважением, Лев
Ответ:
В итоге получился полный и автоматизированный документооборот.
Для переноса писем с вложениями в excel и соотв. папки
Код Visual Basic | ||
|
Ответ:
Строго в модуль книги ThisWorkbook(ЭтаКнига)
личной книги макросов Personal.xls(xlsb)
Visual Basic | ||
|
Пусть корень уравнения (1)
отделен на отрезке
. Требуется найти значение корня с точностью ε
.
«Процедура уточнения положения корня заключается в построении последовательности вложенных друг в друга отрезков, каждый из которых содержит корень уравнения. Для этого находится середина текущего интервала неопределенности (6)
:
В в качестве следующего интервала неопределенности из двух возможных выбирается тот, на концах которых функция F(x)=0
имеет разные знаки»[8
]. «Точность будет достигнута, если:
Корень уравнения вычисляется по формуле x=(a n +b n)/2 (7)
«[1
].
Пусть дана задача следующего характера:
Уточнить корни уравнения cos(2x)+x-5=0
методом половинного деления с точностью до 0,00001, используя:
1. Mathcad
;
Для того чтобы уточнить корни уравнения cos(2x)+x-5=0 методом половинного деления, используя Excel, необходимо выполнить следующие действия:
1. Заполнить ячейки A1:H1 последовательно следующим образом: a, b, c=(a+b)/2, f(a), f(b), f(c), |b-a|<=2*e, e.
2. Ввести в ячейку A2 число 5, в ячейку B2 — число 6.
3. В ячейку B2 ввести формулу: =(A2+B2)/2.
4. В ячейку D2 ввести формулу: =cos(2*A2)+A2-5, скопировать эту формулу в ячейки E2:F2.
5. Ввести в ячейку G2 формулу: =ЕСЛИ(ABS(B2-A2)<=2*$H$2;C2;»-«).
6. Ввести в ячейку H2 число 0,00001.
7. В ячейку A3 ввести формулу: =ЕСЛИ(D2*F2<0;A2;C2).
8. В ячейку B3 ввести формулу: =ЕСЛИ(D2*F2<0;C2;B2).
9. Диапазон ячеек C2:G2 скопировать в диапазон ячеек C3:G3.
10. Выделить диапазон ячеек A3:G3 и с помощью маркера заполнения заполнить все нижестоящие ячейки до получения результата в одной из ячеек столбца G (это ячейки A3:G53).
В итоге получаем следующее:
Ответ: Корень уравнения cos(2x)+x-5=0 равен 5,32977.
- Метод хорд
Берілген әдісті шешу үшін y=F(x) функциясын құру керек
» Для реализации данного метода, нужно построить исходную функциюy=F(x)
и найти значения функции на концах отрезка F(a)
и F(b)
. Затем провести хорду М 1 M 2
c концами в точкахМ 1 (a, F(a))
и M 2 (b, F(b)). Абсцисса точки пересечения хорды М 1 M 2
с осью OX это и есть приближенный кореньx 1
. Далее найти точкуM 3 (X 1 ,F(x 1))
, построить следующую хорду и найти второй приближенный корень x 2
. И так далее. В зависимости от поведения функции возможны два случая
:
Для первого случая
(Рис. 1) справедлива следующая формула (8)
:
и справедливо неравенство: F(a)*F»»(a)>0, где x 0 =b.
Для второго случая
(Рис. 2) справедлива следующая формула (9)
:
и справедливо неравенство: F(b)*F»»(b)>0
, где x 0 =a
.
Условия сходимости метода секущих аналогичны условиям сходимости метода Ньютона, т. е.»[1
]
Пусть дана задача:
Уточнить корни уравнения cos(2x)+x-5=0
методом хорд с точностью до 0,00001, используя:
1. Mathcad
;
Для того чтобы уточнить корни уравнения cos(2x)+x-5=0 методом хорд, используя Excel, необходимо выполнить следующие действия:
1. Выбрать одну из двух предложенных формул для решения задачи, для этого:
o Найти производную первого порядка от функции f(x)=cos(2x)+x-5. Она будет выглядеть следующим образом: f1(x)=-2sin(2x)+1.
o Найти производную второго порядка от функции f(x)=cos(2x)+x-5. Она будет выглядеть следующим образом: f2(x)=-4cos(2x).
o Заполнить ячейки следующим образом:
В ячейку A1 ввести a.
В ячейку A2 ввести цифру 5.
В ячейку B1 ввести b.
В ячейку B2 ввести цифру 6.
В ячейку C1 ввести f(x)=cos(2x)+x-5.
В ячейку C2 ввести формулу =COS(2*A2)+A2-5.
В ячейку D1 ввести f1(x)=-2sin(2x)+1.
В ячейку E1 ввести f2(x)=-4cos(2x).
В ячейку E2 ввести формулу =-4*COS(2*A2).
В ячейку F1 ввести Выбор формулы.
В ячейку F2 ввести формулу =ЕСЛИ(C2*E2>0;»Воспользоваться формулой 8″;»Воспользоваться формулой 9″).
В ячейку G1 ввести e.
В ячейку G2 ввести цифру 0,00001.
o В итоге получается следующее:
2. Исходя из того, что выбрана формула 9, в Excel необходимо выполнить следующие действия:
o В ячейку A4 ввести xn.
o В ячейку B4 ввести f(xn).
o В ячейку C4 ввести b-xn.
o В ячейку D4 ввести f(xn)*(b-xn).
o В ячейку E4 ввести f(b).
o В ячейку F4 ввести f(b)-f(xn).
o В ячейку G4 ввести xn-f(xn)*(b-xn)/f(b)-f(xn).
o В ячейку H4 ввести |f(xn)|<=e.
o В ячейку A5 ввести цифру 5.
o В ячейку B5 ввести формулу =COS(2*A5)+A5-5.
o В ячейку C5 ввести формулу =$B$2-A5.
o В ячейку D5 ввести формулу =B5*C5.
o В ячейку E5 ввести формулу =COS(2*$B$2)+$B$2-5.
o В ячейку F5 ввести формулу =$E$5-B5.
o В ячейку G5 ввести формулу =A5-(B5*C5/F5).
o В ячейку H5 ввести формулу =ЕСЛИ(ABS(B5)<=$G$2;A5;»-«).
o В ячейку A6 ввести формулу =G5.
o Выделить диапазон ячеек B5:D5 и скопировать его методом протягивания в диапазон ячеек B6:D6.
o Выделить диапазон ячеек F5:H5 и скопировать его методом протягивания в диапазон ячеек F6:H6.
o Выделить диапазон ячеек A6:H6 и скопировать его методом протягивания в диапазон ячеек ниже до получения результата в одной из ячеек столбца H (A6:H9).
В итоге получаем следующее:
Ответ: Корень уравнения cos(2x)+x-5=0 равен 5,32976.
Электронная
таблица
Microsoft
Excel
.
Средства и методы решения уравнений.
Цель
работы:
Освоить численный метод
решения уравнения ивстроенные средства
решения уравнений..
Содержание
1
Численный метод решения нелинейных уравнений
.
1
1.1
Область локализации корней
.
1
1.2
Критерии сходимости при решении уравнений
.
2
1.3
Метод дихотомии (половинного деления)
3
Пример
решения уравнения методом дихотомии
.
4
2
Решение уравнений, используя “Подбор параметра”
.
6
2.1
Пример решения уравнения, используя “Подбор параметра”
.
6
3
Решение уравнений и систем уравнений, используя надстройку “Поиск решения”
.
9
3.1
Пример решения уравнения, используя надстройку “Поиск решения”
.
10
Задание
1. Решение уравнений численным методом
..
12
Задания
2. Решение уравнений встроенными средствами “Подбор параметра” и “Поиск
решения”
12
Контрольные
вопросы
..
13
1 Численный метод решения нелинейных уравнений
1.1 Область локализации корней
В общем виде любое
уравнение одной переменной принято записывать так , при этом корнем (решением) называется такое значение x*,
что оказывается верным тождеством. Уравнение может иметь
один, несколько (включая бесконечное число) или ни одного корня. Как легко видеть,
для действительных корней задача отыскания решения уравнения легко
интерпретируется графически: корень есть такое значение независимой переменной,
при котором происходит пересечение графика функции, стоящей в левой части
уравнения
f
(
x
)
, с осью абсцисс.
Например
,
для уравнения выполним
преобразование и приведем его к виду f(x)=0
т.е. . График этой функции представлен на рисунке 1. Очевидно, что
данное уравнение имеет два действительных корня – один на отрезке [-1, 0] , а
второй – .
Рисунок 1. График функции
Таким образом, можно приблизительно
определять область локализации корней
уравнения. Заметим, что отделить корень можно не единственным образом: если
корень отделён на каком-либо отрезке, то годится и любой меньший отрезок,
содержащий этот корень. Вообще говоря, чем меньше отрезок, тем лучше, но при
этом не следует забывать о том, что на отделение корня на меньших отрезках
также тратятся вычислительные усилия, и, быть может, весьма значительные. Таким
образом, часто для начала довольствуются весьма широким отрезком, на котором
корень отделён.
Некоторые виды
уравнений допускают аналитическое решение. Например, степенные
алгебраические уравнения степени n
приn
≤ 4. Однако, в общем виде, аналитическое решение
, как правило,
отсутствует. В этом случае, применяются численные методы
. Все численные методы решения уравнений
представляют собой
последовательного приближения к корню
уравнения. То есть, выбирается начальное приближение к корню
x
0
и затем с помощью итерационной формулы
генерируется последовательность
x
1 ,
x
2 , …,
x k
сходящаяся
к корню уравнения
.
1.2 Критерии сходимости при решении уравнений
Ø
Абсолютная погрешность — абсолютное
изменение приближения на соседних шагах итерации
Ø
Относительная погрешность — относительное
изменение приближения на соседних шагах итерации
Ø
Близость к нулю вычисленного значения левой
части уравнения (иногда это значение называют невязкой
уравнения,
так как для корня невязка равна нулю)
1.3
Метод половинного деления
(метод дихотомии)
Метод половинного
деления основан на последовательном делении отрезка локализации корня пополам.
Для этого выбирается
начальное приближение к отрезку [
a
,
b
], такое, что
f
(
a
)
×
f
(
b
)<0
, затем
определяется знак функции в точке — середине отрезка [
a
,
b
]. Если он
противоположен знаку функции в точке a
,
то корень локализован на отрезке [
a
,
c
], если же нет – то на отрезке [
c
,
b
].
Схема метода
дихотомии приведен
на рис
у
нке 2.
Рисунок 2. Последовательное деление отрезка
пополам и приближение к корню
Алгоритм
метода дихотомии можно записать так:
1.
представить решаемое уравнение в виде
2.
выбрать a, b и вычислить
3.
если
f(a)
×
f(с
)<0,
то
a=a; b = c
иначе
a = c; b=b
4.
если критерий сходимости не выполнен, то перейти к п. 2
Пример решения уравнения методом
дихотомии
Найти решение заданного уравнения методом
дихотомии с точностью до 10 -5 .
Пример создания расчетной схемы на основе
метода дихотомии на примере уравнения:
на отрезке
Данный метод заключается в проверке на
каждой итерации условия:
если
f
(
a
)
×
f
(с)<0 и выбор
соответствующего отрезка для следующей итерации.
a)
|
b)
|
Рисунок 3. Последовательность итераций метода дихотомии
при поиске корня уравнения на отрезке
a
)
схема
расчета
(зависимые ячейки); b)
режимотображения формул;
Для нашего примера итерационная
последовательность для нахождения решения принимает вид:
Точность до пятой значащей цифры достигается
за 20 итераций.
Скорость сходимости этого метода является
линейной.
При выполнении начального условия он
сходится к решению всегда.
Метод половинного деления удобен при решении
физически реальных уравнений, когда заранее известен отрезок локализации
решения уравнения.
2 Решение уравнений
, используя “Подбор
параметра
”
Используя
возможности Excel можно находить корни нелинейного уравнения вида f(x)=
0 в допустимой области определения
переменной. Последовательность операцийнахождения корней следующая:
1.
Производится табулирование функции в диапазоне
вероятного существования корней;
2.
По таблице фиксируются ближайшие приближения к
значениям корней;
3.
Используя средство Excel Подбор параметра,
вычисляются корни уравнения с заданной точностью.
При
подборе параметра Excel использует итерационный (циклический) процесс.
Количество итераций и точность устанавливаются в меню Сервис/Параметры/вкладка Вычисления
. Если Excel выполняет сложную
задачу подбора параметра, можно нажать кнопку Пауза
в окне диалога Результат
подбора параметра
и прервать вычисление, а затем нажать кнопку Шаг
, чтобы выполнить очередную итерацию
и просмотреть результат. При решении задачи в пошаговом режиме появляется
кнопка П
родолжить
— для возврата в обычный
режим подбора параметра.
2.1 Пример решения уравнения,
используя “Подбор параметра”
Например
, найдем все корни уравнения 2x 3 -15sin(x)+0,5x-5=0
на отрезке [-3 ;
3].
Для
локализации начальных приближений необходимо определить интервалы значений Х,
внутри которых значение функции пересекает ось абсцисс, т.е. функция меняет
знак. С этой целью табулируем функцию на отрезке [–3; 3] с шагом 0,2, получим
табличные значения функции. Из полученной таблицы находим, что значение функции
трижды пересекает ось Х, следовательно, исходное уравнение имеет на заданном
отрезке все три корня.
Рисунок 4. Поиск приближенных значений
корней уравнения
Выполните
команду меню Сервис/Параметры,
во
вкладке Вычисления
установите
относительную погрешность вычисленийE=0,00001, а число итераций N=1000, установите флажок Итерации.
Выполните
команду меню Сервис/Подбор параметра
.
В диалоговом окне (рисунок 9) заполните следующие поля:
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ
ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ
ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«САМАРСКИЙ ГОСУДАРСТВЕННЫЙ
АРХИТЕКТУРНО-СТРОИТЕЛЬНЫЙ УНИВЕРСИТЕТ»
Кафедра прикладной математики и вычислительной техники
Excel
и
Mathcad
МЕТОДИЧЕСКИЕ УКАЗАНИЯ
к выполнению лабораторных работ
по дисциплине «Вычислительная математика»
Решение нелинейных уравнений в
Excel и
Mathcad
:
Метод. указ. / Сост. , — Самара: СГАСУ, 20с.
Методические указания разработаны в соответствии с Государственным образовательным стандартом изучения дисциплины «Вычислительная математика».
Рассмотрена реализация численных методов при решении нелинейных уравнений и систем уравнений в Excel и MathCad. Приведены варианты заданий для индивидуального выполнения и вопросы для самоконтроля и тестирования.
Предназначены для студентов специальности 230201 – «Информационные системы и технологии» всех форм обучения.
Рецензент к. ф-м. н.
Ó , составление, 2012
ã СГАСУ, 2012
1.2 Отделение корней |
|
1.5 Метод хорд |
|
1.6 Метод Ньютона (касательных) |
|
1.7 Комбинированный метод |
|
1.8 Метод итераций |
|
2.2 Решение систем нелинейных уравнений методом Ньютона |
|
3 Задания к лабораторным работам |
|
Лабораторная № 1. Отделение корней и стандартные инструменты решения нелинейного уравнения |
|
Лабораторная № 2. Сравнение методов уточнения корней нелинейного уравнения |
|
Лабораторная № 3. Решение систем нелинейных уравнений |
|
Лабораторная № 4. Программирование методов решения нелинейных уравнений и систем |
|
4 Вопросы и тесты для самоконтроля |
|
1 Решение нелинейного уравнения
1.1 Общие сведения о решении нелинейного уравнения
Как правило, нелинейное уравнения общего вида f(х)=0
невозможно решить аналитически. Для практических задач достаточно найти приближенное значение x
, в определенном смысле близкое к точному решению уравнения хточн
.
В большинстве случаев поиск приближенного решения включает два этапа. На первом этапе
отделяют
корни, т. е. находят такие отрезки, внутри которых находится строго один корень. На втором этапе
уточняют
корень на одном из таких отрезков, т. е. находят его значение с требуемой точностью.
Достигнутая точность может оцениваться либо «по функции» (в найденной точке x
, функция достаточно близка к 0, т. е. выполняется условие |f(x)|≤
e
f
, где e
f
требуемая точность по оси ординат), либо «по аргументу» (найден достаточно маленький отрезок [
a,
b]
, внутри которого находится корень, т. е. |
b–
a|≤
e
x
, где e
x
требуемая точность по оси абсцисс).
1.2 Отделение корней
Отделение корней может производиться сочетанием графического
и аналитического
исследования функции. Такое исследование опирается на теорему Вейерштрасса, в соответствии с которой для непрерывной на отрезке [
a,
b]
функции f(х
) и любого числа y
, отвечающего условию f(a)≤y≤
f(b)
, существует на этом отрезке точка x
, в которой функция равна y
. Следовательно, для непрерывной функции достаточно найти отрезок, на концах которого функция имеет разные знаки, и можно быть уверенным, что на этом отрезке есть корень уравнения f(х)=0
.
Для ряда методов уточнения желательно, чтобы найденный на первом этапе отрезок содержал только один корень уравнения. Это условие выполняется, если функция на отрезке монотонна. Монотонность, можно проверить либо по графику функции, либо по знаку производной.
Пример
Найти с точностью до целых все
корни нелинейного уравнения y(x)=
x3 ‑ 10
x + 7=0
а) построив таблицу и б) построив график. Найти корень уравнения на выделенном отрезке, используя опции «Подбор параметра» и «Поиск решения».
Решение
Создадим в Excel таблицу, содержащую аргументы и значения функции и по ней построим точечную диаграмму
. На рисунке 1 приведен снимок решения.
На графике видно, что уравнение имеет три корня, принадлежащие отрезкам [-4, -3], и . Эти отрезки можно выявить и наблюдая за сменой знаков функции в таблице. По построенному графику можно сделать вывод, что на указанных отрезках функция f
(x
)
монотонна и, следовательно, на каждом из них содержится только по одному корню.
Такой же анализ может быть выполнен и в пакете Mathcad. Для этого достаточно набрать определение функции f
(x
)
, используя оператор присваивания (:=) и естественные общепринятые обозначения математических операций и стандартных функций, задать цикл для изменения аргумента, например, а затем вывести на экран таблицу значений функции (расположенными в одной строке командами x
=
f
(x
)=
) и график. Цикл можно задать, например, командой x
:=-5,-4.5…5
. Шаг цикла формируется путем задания начального и следующего за ним значений переменной, а перед конечным значением переменной ставится точка с запятой, которая будет визуально отображена на экране в виде многоточия.
https://pandia.ru/text/78/157/images/image002_56.jpg» width=»640″ height=»334″>
Рисунок 1 – Таблица и график для отделения корней нелинейного уравнения
1.3 Уточнение корней стандартными средствами Excel и Mathcad
Во всех методах уточнения корней необходимо задать начальное приближение, которое затем и будет уточняться. Если уравнение имеет несколько корней, в зависимости от выбранного начального приближения будет найден один из них. При неудачно выбранном начальном приближении решение может и не быть найдено. Если в результате первого этапа расчетов уже выделен отрезок, содержащий единственный корень уравнения, в качестве начального приближения можно взять любую точку этого отрезка.
В Excel для уточнения значений корней можно использовать опции «Подбор параметра» и «Поиск решения». Пример оформления решения приведен на рисунках 2 и 3.
https://pandia.ru/text/78/157/images/image004_31.jpg» width=»501″ height=»175 src=»>
Рисунок 3 – Результаты использования средств решения уравнения в
Excel
В Mathcad для уточнения корней уравнения можно использовать функцию root
(….)
или блок решения
. Пример использования функции root(…) приведен на рисунке 4, а блока решения на рисунке 5. Следует обратить внимание, что в блоке решения (после заголовка блока Given
) между левой и правой частями уравнения должен стоять жирный знак равенства
(тождества), который можно получить выбором из соответствующей палитры инструментов, либо нажатием одновременно клавиши Ctrl
и =
.
243″ height=»31″>
Рисунок 5 – Решение уравнения с использованием блока решения в
Mathcad
Как видим, каждый стандартный инструмент находит решение уравнения с определенной точностью. Эта точность зависит от метода, используемого в пакете и, в определенной степени, настроек пакета. Управлять точностью результата здесь достаточно сложно, а часто и невозможно.
В то же время, очень просто построить собственную таблицу или написать программу, реализующие один из методов уточнения корней. Здесь можно использовать критерии точности расчета, задаваемые пользователем. При этом достигается и понимание процесса расчетов без опоры на принцип Митрофанушки: «Извозчик есть, довезет».
Далее рассмотрены несколько наиболее распространенных методов. Отметим очевидный момент: при прочих равных условиях тот метод
уточнения корней будет более эффективен, в котором результат с той же погрешностью найден с меньшим
числом вычислений функции f(x)
(при этом достигается и максимальная точность при одинаковом числе вычислений функции).
1.4 Метод деления отрезка пополам
В этом методе на каждом шаге отрезок делится на две равные части. Затем сравнивают знаки функции на концах каждой из двух половинок (например, по знаку произведения значений функций на концах), определяют ту из них, в которой содержится решение (знаки функции на концах должны быть разные), и. сужают отрезок, перенося в найденную точку его границу (а
или b
). Условием окончания служит малость отрезка, где содержится корень («точность по x
»), либо близость к 0 значения функции в средине отрезка («точность по y»). Решением уравнения считают середину отрезка, найденного на последнем шаге.
Пример
. Построить таблицу для уточнения корня уравнения x
3
–10
x
+7=0
на отрезке [-4, -3]
методом деления отрезка пополам. Определить сколько шагов надо сделать методом деления отрезка пополам и какая при этом достигается точность по х,
для достижения точности по y
, равной 0,1; 0,01; 0, 001.
Решение
Для решения можно использовать табличный процессор Excel, позволяющий автоматически продолжать строки. На первом шаге заносим в таблицу значения левого и правого концов выбранного начального отрезка и вычисляем значение середины отрезка с
=(a
+b
)/2, а затем вводим формулу для вычисления функции в точке a
(f
(a
)) и растягиваем (копируем) её для вычисления f
(c
) и f
(b
). В последнем столбца вычисляем выражение (b
—a
)/2, характеризующего степень точности вычислений. Все набранные формулы можно скопировать во вторую строку таблицы.
На втором шаге нужно автоматизировать процесс поиска той половины отрезка, где содержится корень. Для этого испльзуется логическая функция ЕСЛИ (Меню
: ВставкаФункцияЛогические). Для нового левого края отрезка мы проверяем истинность условия f
(a
)*f
(c
)>0, если оно верно, то мы в качестве нового значения левого конца отрезка берем число c
a
,
c
a
. Аналогично, для нового правого края отрезка мы проверяем истинность условия f
(c
)*
f
(b
)>0, если оно верно, то мы в качестве нового значения правого конца отрезка берем число c
(т. к. это условие показывает, что корня на отрезке [c
,
b
] нет), иначе оставляем значение b
.
Вторую строку таблицы можно продолжить (скопировать) на необходимое число последующих строк.
Итерационный процесс завершается, когда очередное значение в последнем столбце становится меньшим, чем заданный показатель точности ex. При этом, значение середины отрезка в последнем приближении, принимается в качестве приближенного значения искомого корня нелинейного уравнения. На рисунке 6 приведен снимок решения. Для построения аналогичного процесса в Mathcad можно использовать бланк, подобный приведенному на рисунке 7. Число шагов N может варьироваться до достижения в таблице результатов требуемой точности. При этом таблица будет автоматически удлиняться или укорачиваться.
Итак, одним из трех корней нелинейного уравнения x
3 – 10x
+ 7=0, найденным с точностью e=0,0001, является x
= — 3,46686. Как мы видим, он действительно принадлежит отрезку [-4; -3].
https://pandia.ru/text/78/157/images/image018_6.jpg» width=»563″ height=»552 src=»>
Рисунок 7 – Уточнение корня методом деления отрезка пополам в
Mathcad
1.5 Метод хорд
В этом методе нелинейная функция f(x)
на отделенном интервале [а, b
] заменяется линейной – уравнением хорды, т. е. прямой соединяющей граничные точки графика на отрезке. Условие применимости метода – монотонность функции на начальном отрезке, обеспечивающая единственность корня на этом отрезке. Расчет по методу хорд аналогичен расчету методом деления отрезка пополам, но теперь на каждом шаге новая точка x
внутри отрезка [a
,
b
] рассчитывается по любой из следующих формул:
(х) > 0
), или правая его граница: x0 = b
(если f(b) f»(х)>0
). Расчет нового приближения на следующем шаге i
+1
производится по формуле:
https://pandia.ru/text/78/157/images/image021_4.jpg» width=»596″ height=»265 src=»>
Рисунок 8 – Уточнение корня методом касательных в E
xcel
Расчеты в Mathcad выполняются аналогично. При этом значительное облегчение доставляет наличие в этом пакете оператора, автоматически вычисляющего производную функции.
Наиболее трудоемким элементом расчетов по методу Ньютона является вычисление производной на каждом шаге.
При определенных условиях может использоваться упрощенный метод Ньютона
, в котором производная вычисляется только один раз – в начальной точке. При этом используется видоизмененная формула
.
Естественно, что упрощенный метод, как правило, требует большего числа шагов.
Если вычисление производной связано с серьезными трудностями (например, если функция задана не аналитическим выражением, а вычисляющей ее значения программой) используется модифицированный метод Ньютона, получивший название – метод секущих
. Здесь производная приближенно вычисляется по значениям функции в двух последовательных точках, то есть используется формула
.
В методе секущих необходимо задаться не одной, а двумя начальными точками – x
0
и x
1
. Точка x1
обычно задается сдвигом x0
к другой границе отрезка на малую величину, например, на 0.01.
1.7 Комбинированный метод
Можно показать, что если на начальном отрезке у функции f(x)
сохраняются неизменными знаки первой и второй производных, то методы хорд и Ньютона приближаются к корню с разных. В комбинированном методе для повышения эффективности на каждом шаге использует оба алгоритма одновременно. При этом интервал, где содержится корень, сокращается с обеих сторон, что обусловливает другое условие окончания поиска. Поиск можно прекратить, как только в середине интервала, полученного на очередном шаге значение функции станет по модулю меньшим, чем предварительно заданной погрешности e
f
.
Если, в соответствии со сформулированным выше правилом, метод Ньютона применяется к правой границе отрезка, для вычислений используются формулы:
https://pandia.ru/text/78/157/images/image025_10.gif» width=»107″ height=»45 src=»>.
Если метод Ньютона применяется к левой границе, – в предыдущих формулах меняются местами обозначения a
и b
.
1.8 Метод итераций
Для применения этого метода исходное уравнение f(x)=0
преобразуют к виду: x
=y
(х)
. Затем выбирают начальное значение х0
и подставляют его в левую часть уравнения, получая, в общем случае, x
1
=
y
(х0)
¹
х0
¹
y
(х1)
, поскольку х0
взято произвольно и не является корнем уравнения. Полученное значение х1
рассматривают как очередное приближение к корню. Его снова подставляют в правую часть уравнения и получают следующее значение х2=
y
(х1)
). Расчет продолжают по формуле хi+1=
y
(хi)
. Получающаяся таким образом последовательность: х0, х1, х2, х3 х4,…
при определенных условиях сходиться к корню хточн
.
Можно показать, что итерационный процесс сходится при условии
|y
’
(x
) | < 1 на [a
, b
].
Существуют различные способы преобразования уравнения f(x)
= 0 к виду y
(х)
= х
, причем в конкретном случае одни из них приведут к сходящемуся, а другие – к расходящемуся процессу вычислений.
Один из способов, заключается в применении формулы
https://pandia.ru/text/78/157/images/image027_10.gif» width=»188″ height=»44 src=»>
где М
= max |y
’
(x
)| на [a
, b
].
2 Решение систем нелинейных уравнений
2.1 Общие сведения о решении систем нелинейных уравнений
Систему n
нелинейных уравнений с n
неизвестными x1, x2
, …, xn
записывают в виде:
где F1, F2
,…, Fn
– функции независимых переменных, среди которых есть нелинейные.
Как и в случае систем линейных уравнений, решением системы является такой вектор X
*, который при подстановке обращает одновременно все уравнения системы в тождества.
https://pandia.ru/text/78/157/images/image030_8.gif» width=»191″ height=»56″>
Начальные значения x
0
и y
0
определяются графически. Для нахождения каждого последующего приближения (xi
+1
,
yi
+1
)
используют вектор значений функций и матрицу значений их первых производных, рассчитанные в предыдущей точке (xi
,
yi
)
.
https://pandia.ru/text/78/157/images/image032_5.gif» width=»276″ height=»63 src=»>
Для расчета новых приближений на шаге i+1
используется матричная формула
https://pandia.ru/text/78/157/images/image034_4.gif» width=»303″ height=»59 src=»>.
Приведенные формулы особенно легко записать в Mathcad, где имеются операторы для вычисления производных и действий с матрицами. Однако при правильном использовании матричных операций эти формулы достаточно просто записываются и в Excel. Правда, здесь придется заранее получить формулы для вычисления производных. Для аналитического вычисления производных также может быть использован Mathcad.
2.3 Решение систем нелинейных уравнений методами итераций
Для реализации этих методов исходную систему уравнений необходимо путем алгебраических преобразований явно выразить каждую переменную через остальные. Для случая двух уравнений с двумя неизвестными новая система будет иметь вид
https://pandia.ru/text/78/157/images/image036_5.gif» width=»114″ height=»57 src=»>.
Если одно из решений системы и начальные значения x
0
и y
0
лежат в области D
, задаваемой неравенствами: a
≤ x
≤ b
, c
≤ y
≤ d
, то расчет по методу простых итераций сходится при выполнении в области D
соотношений:
https://pandia.ru/text/78/157/images/image038_5.gif» width=»75 height=48″ height=»48″>< 1.
В методе итераций Зейделя
для каждого расчета используют уже найденные наиболее точные значения каждой переменной. Для рассматриваемого случая двух переменных такая логика приводит к формулам
0 » style=»border-collapse:collapse;border:none»>
Инструмент (опция)
Начальное приближение
Корень
x
f(x)
3.Отсортировать полученные результаты по точности решения.
Лабораторная работа
Отделение корней
нелинейного уравнения
Пусть имеется
нелинейное уравнение .
Требуется найти
корни этого уравнения. Численный процесс приближенного решения поставленной
задачи разделяют два этапа: отделение корня и уточнение корня.
Для отделения
корня необходимо определить промежуток аргумента , где
содержится один и только один корень уравнения. Одна из точек этого промежутка
принимается за начальное приближение корня. В зависимости от метода, который
предполагается использовать для уточнения корня, требуется определение
некоторых свойств отделенного корня и поведения функции на отрезке отделения.
Например, при использовании метода деления пополам, необходимо и достаточно
установить лишь непрерывность функции на отрезке отделения.
Этап отделения
корня уравнения алгоритмизирован только для некоторых классов уравнений
(наиболее известным из которых является класс алгебраических уравнений),
поэтому отделение корней нелинейных уравнений, обычно, выполняется «вручную» с
использованием всей возможной информации о функции . Часто
применяется графический метод отделения действительных корней, обладающий
большой наглядностью.
Методы
отделения корней
Отделение корней
во многих случая можно произвести графически. Учитывая, что действительные
корни уравнения F(x)=0 – это есть точки пересечения графика
функции y=F(x) с осью абсцисс y=0, нужно построить
график функции y=F(x) и на оси OX отметить отрезки,
содержащие по одному корню. Но часто для упрощения построения графика функции y=F(x)
исходное уравнение заменяют равносильным ему уравнением f1(x)=f2(x).
Далее строятся графики функций y1=f1(x)
и y2=f2(x), а затем по оси OX
отмечаются отрезки, локализующие абсциссы точек пересечения двух графиков.
На практике
данный способ реализуется следующим образом: например, требуется отделить корни
уравнения cos(2x)+x-5=0 графически на отрезке [–10;10],
используя Excel.
1 способ
Построим график функции f(x)=cos(2x)+x-5
в декартовой системе координат. Для этого нужно:
1.
Ввести в
ячейку A1 текст х.
2.
Ввести в
ячейку B1 текст y=cos(2x)+x-5.
3.
Ввести в
ячейку А2 число -10, а в ячейку А3 число -9.
4.
Выделить
ячейки А2 и А3.
5.
Навести
указатель «мыши» на маркер заполнения в правом нижнем углу рамки, охватывающий
выделенный диапазон. Нажать левую кнопку «мыши» и перетащить маркер так, чтобы
рамка охватила диапазон ячеек А2:А22.
6.
Ячейки
автоматически заполняются цифрами :
7.
Ввести в
ячейку В2 формулу =COS(2*A2)+A2-5.
8.
Методом
протягивания заполнить диапазон ячеек В3:В22.
9.
Вызвать
«Мастер диаграмм» и выбрать диаграмму график (первый вид), нажать
«далее».
10.
Указать
диапазон данных, для этого щелкнуть кнопку в поле «Диапазон» и выбрать диапазон
данных В2:В22.
11.
Выбрать
вкладку ряд, указать имя ряда, щелкнув кнопку в поле «ряд» и выбрав В1.
12.
В поле
«подписи по оси Х», щелкнуть кнопку и выбрать диапазон А2:А22, нажать «далее».
13.
Подписать
названия осей x и y соответственно, нажать «далее».
14.
Вывести
диаграмму на том же листе, что и таблица, нажать кнопку «готово».
В итоге получаем следующее
(рисунок 1):
Рисунок 1 – Локализация корня
Анализируя полученное
изображение графика, можно сказать, что уравнение cos(2x)+x-5=0
имеет один корень – это видно из пересечения графика функции y=cos(2x)+x-5
с осью OX. Можно выбрать отрезок, содержащий данный корень:[5;6] – отрезок локализации.
2 способ
Для подтверждения полученных
данных, можно решить эту же задачу вторым способом. Для этого необходимо
уравнение cos(2x)+x-5=0 преобразовать к виду: cos(2x)=5-x.
Затем следует каждую часть уравнения рассмотреть как отдельную функцию. Т. е. y1=cos(2x)
и y2=5-x. Для решения этой задачи в Excel необходимо
выполнить следующие действия:
1.
Вести в
ячейки А1:C1 соответственно текст: «x», «y1=cos(2x)»,
«y2=5-x».
2.
A2:A22
заполнить так же как при решении задачи первым способом.
3.
В В2
ввести формулу =COS(2*A2).
4.
Методом
протягивания заполнить диапазон ячеек В3:В22.
5.
В С2
ввести =5-A2.
6.
Методом
протягивания заполнить диапазон ячеек С3:С22.
7.
С помощью
Мастера диаграмм выбрать график (первый вид).
8.
В данном
случае диапазон данных следует указывать для построения двух графиков. Для
этого нужно нажать кнопку в поле «Диапазон» и выделить ячейки В2:В22, затем
нажать Ctrl (на клавиатуре) и выделить следующий диапазон C2:C22.
9.
Перейти
на вкладку ряд, где выбрать именем ряда 1 ячейку В1, а именем ряда 2 ячейку С2.
10.
Подписать
ось x , выбрав диапазон А2:А22.
11.
Подписать
соответственно оси x и y.
12.
Поместить
диаграмму на имеющемся листе.
Результат представлен на
рисунке 2: Анализируя
полученный результат, можно сказать, что точка пересечения двух графиков
попадает на тот же самый отрезок локализации [5;6], что и при решении задачи
первым способом.
Рисунок 2 – Локализация корня
Аналитический способ отделения
корней
Аналитический
способ отделения корней основан на следующей теореме, известной из курса математического
анализа.
ТЕОРЕМА: Если непрерывная на функция ,
определяющая уравнение , на концах отрезка принимает значения разных знаков, т.е. , то на этом отрезке содержится, по
крайней мере, один корень уравнения. Если же функция непрерывна
и дифференцируема и ее производная сохраняет знак внутри отрезка , то на этом отрезке находится только один
корень уравнения.
В случае, когда
на концах интервала функция имеет одинаковые знаки, на этом интервале корни
либо отсутствуют, либо их четное число.
Для отделения корней
аналитическим способом выбирается отрезок , на
котором находятся все интересующие вычислителя корни уравнения. Причем на
отрезке функция F(x) определена, непрерывна и F(a)*F(b)<0.
Требуется указать все частичные отрезки ,
содержащие по одному корню.
Будем вычислять
значение функции F(x),
начиная с точки x=a, двигаясь вправо с некоторым шагом h. Если F(x)*F(x+h)<0,
то на отрезке [x;x+h] существует корень (рисунок 3).
Рисунок 3 –
Аналитический способ локализации корней
Если F(xk)=0,
xk-точный корень.
Доказательство существования
и единственности корня на отрезке.
В качестве примера рассмотрим
функцию f(x)=cos(2x)+x-5.
1. Ввести в ячейки А1, В1 и С1
соответственно «x», «y=cos(2x)+x-5» и «ответ».
2. В А2 и А3 ввести граничные
значения отрезка изоляции.
3. В В2 ввести формулу
=COS(2*A2)+A2-5 и методом протягивания заполнить В3.
4. В С2 ввести формулу
=ЕСЛИ(B2*B3<0;»корень существует»;»корень не
существует»).
Таким образом, на отрезке
изоляции корень существует:
Рисунок 4 – Проверка существования корня на отрезке
Для доказательства
единственности корня на отрезке изоляции необходимо выполнить следующие
действия:
1.
Продолжить
работу в том же документе MS Excel.
2.
Заполнить
D1 и E1 соответственно: «y’=-sin(2x)*2+1» и «ответ» (причем выражение
y’=-sin(2x)*2+1 – это производная первого порядка от функции y=cos(2x)+x-5).
3.
Ввести в
D2 формулу =-SIN(2*A2)*2+1 и методом протягивания заполнить D3.
4.
Ввести в
E2 =ЕСЛИ(D2*D3>0;»корень на данном отрезке единственный»;»Корень
не единственный»).
В результате получаем
(рисунок 5):
Рисунок 5 –
Доказательство единственности корня на отрезке
Таким образом доказано
существование и единственность корня на отрезке изоляции.
Рассмотрим
решение задачи отделения корней уравнения
cos(2x)+x-5=0 аналитическим
способом с шагом 1 на отрезке [-10;10].
Чтобы отделить корни
уравнения аналитическим способом с помощью Excel, необходимо выполнить
следующее:
1.
Заполнить
ячейки A1:D1 соответственно: «x», «y=cos(2x)+x-5»,
«h», «ответ».
2.
В С2
ввести значение 1.
3.
Ввести в
А2 значение -10.
4.
Ввести в
А3 =A2+$C$2 и методом протягивания заполнить ячейки А4:А22.
5.
В В2
ввести =COS(2*A2)+A2-5 и методом протягивания заполнить диапазон В3:В22.
6.
В С3 ввести формулу
=ЕСЛИ(B2*B3<0;»Корень на отрезке существует»;ЕСЛИ(B3=0;»точный
корень»;»-«)) и методом протягивания заполнить диапазон ячеек
С4:С22.
В результате получаем
следующее (рисунок 6):
Рисунок 6 –
Отделение корня
Следующий пример (рисунок
7) демонстрирует отделение нескольких корней. Пусть исследуется функция cos(x)=0,1x на интервале [–10;10] с шагом
1.
Табулирование
функции и построение графика осуществляется как в предыдущих примерах. Видно,
что на заданном отрезке имеем 7 корней, находящихся внутри отрезков: [-10;-9];
[-9;-8]; [-5;-4]; [-2;-1]; [1;2]; [5;6]; [7;8].
Рисунок 7 – Отделение корней
Обратим внимание
на то, что надежность рассмотренного алгоритма отделения корней уравнения
зависит как от характера функции F(x), так и от выбранной величины шага h.
Для повышения надежности следует выбирать при отделении корней достаточно малые
значения h.
Задание
1. Выполнить отделение корней
следующих функций:
№ п/п |
Уравнение |
A |
B |
1 |
tg(x) = 1/x |
0 |
n/2 |
2 |
e -x = x |
0 |
1 |
3 |
ln(x) = 1/x |
1 |
2 |
4 |
2 +ln(x) = 1/x |
0 |
1 |
5 |
x — x3 + 1 =0 |
1 |
2 |
2. Выполнить индивидуальные задания
Вариант 1
1 |
ctg(x) = —x2 |
1,6 |
4,5 |
2 |
2ln(x)+sin(x) =e x |
0 |
2 |
3 |
lg(x) = 2 x-x3 |
0 |
10 |
4 |
cos(x) = 1/x |
0 |
p/2 |
5 |
cos(x) = ln(1+x) |
0 |
p/2 |
Вариант 2
1 |
tg(x) = 1/x—x2 |
1,6 |
4,5 |
2 |
2ln(x) =e x |
0 |
2 |
3 |
lg(x) = sin(x) |
0 |
10 |
4 |
cos(x)+x2 = 1/x |
0 |
p/2 |
5 |
cos2(x) = ln(1+x2) |
0 |
p/2 |
Вариант 3
1 |
cos2(x) = x |
0 |
p/2 |
2 |
1 — 3 x + x3=0 |
0 |
1 |
3 |
1 — 3x + x4=0 |
0 |
1 |
4 |
1 — 3 x + x5=0 |
0 |
1 |
5 |
tg(x) = 1/x2 |
0 |
p/2 |
Вариант 4
1 |
ln(x) = sin(x) |
1 |
3 |
2 |
e — x = sin(x) |
0 |
p/2 |
3 |
e x = 1 /sin(x) |
0 |
p/2 |
4 |
e — x = x2 |
0 |
1 |
5 |
2 + ln(x) = 1/x2 |
0 |
1 |
Вариант 5
1 |
ln(x) = Sin (x) |
0 |
p/2 |
2 |
x — x3 + 2 =0 |
1 |
2 |
3 |
x +5 = x3 |
1 |
2 |
4 |
x 2— 0,5 x-2=0 |
0 |
0,5 |
5 |
ln(x -1)+ x-2=0 |
1 |
3 |
Вариант 6
1 |
ln(x+3) = Sin (x) |
0 |
p/2 |
2 |
2x — x3 + 3 =0 |
1 |
2 |
3 |
x +8 = x3+x2 |
1 |
2 |
4 |
x — 0,5 x2+4=0 |
0 |
0,5 |
5 |
(x -1)2— x=0 |
1 |
3 |
Вариант 7
1 |
x +0,5 = e -x |
0 |
1 |
2 |
2 — x + x3=0 |
-2 |
0 |
3 |
sin(x) = 1/x |
0 |
p/2 |
4 |
sin(x) = x/2 |
п/2 |
п |
5 |
ln(x) = e -x |
0 |
2 |
Вариант 8
1 |
lg(x) = e — x |
0 |
1 |
2 |
cos(x) = x |
0 |
p/2 |
3 |
cos(x) = ln(x) |
0 |
p/2 |
4 |
cos(x) = tg(x) |
0 |
p/2 |
5 |
cos(x) = x3 |
0 |
p/2 |
Вариант 9
1 |
1 — 5 x + x3=0 |
0 |
1 |
2 |
1 — 5 x + x4=0 |
0 |
1 |
3 |
1 — 3 x + x5=0 |
1 |
2 |
4 |
4cos(x) = x+x2 |
0 |
p/2 |
5 |
cos(x) = ln(x)-x2 |
0 |
p/2 |
Вариант 10
1 |
ln(x) = cos (x) |
0 |
p/2 |
2 |
x2 — x3 + 2 =0 |
1 |
2 |
3 |
cosx +5 = x3 |
1 |
2 |
4 |
x 0,5— x=0 |
0 |
0,5 |
5 |
(x -1)2— x=5 |
1 |
3 |
Вариант 11
1 |
ln(x+3) =sin (x)-x |
0 |
p/2 |
2 |
2x — x3 + 3 =ln(x) |
1 |
2 |
3 |
cos(x) +2 = x3+x2 |
1 |
2 |
4 |
3x — 0,5 x2+cosx=0 |
0 |
0,5 |
5 |
(x -1)2— x+tg(x-1)=0 |
1 |
3 |
Вариант 12
1 |
x2 +0,5 = e -x |
0 |
1 |
2 |
2 – sin(x) + x3=0 |
-2 |
0 |
3 |
sin(x) = 1/x-x2 |
0 |
p/2 |
4 |
sin(x) = x/2+cos(x) |
п/2 |
п |
5 |
ln(x)-x = e -x |
0 |
2 |
Вариант 13
1 |
lg(x)+2sin2(x) = e — x |
0 |
1 |
2 |
cos(x+p/2) = x+x2 |
0 |
p/2 |
3 |
cos(x) = ln(x)-x2 |
0 |
p/2 |
4 |
cos(x) = |
0 |
p/2 |
5 |
cos(x) = x3-log(x2) |
0 |
p/2 |
Вариант 14
1 |
1 — 5 cos(x) + x3=0 |
0 |
1 |
2 |
1 — 5 tg(x) + x4=0 |
0 |
1 |
3 |
1 – 3( x-2)2 + x5=0 |
1 |
2 |
4 |
4cos(x) = cos(x)+x2 |
0 |
p/2 |
5 |
cos(x) = ln(x)-x2+x |
0 |
p/2 |