Исследование точности численного интегрирования
Информация - Компьютеры, программирование
Другие материалы по предмету Компьютеры, программирование
Министерство общего и профессионального образования РФ.
Уральский государственный технический университет УПИ
Кафедра “Технология и средства связи”
"Исследование точности численного интегрирования"
"Research of Accuracy of Numerical Integration"
Отчет
по лабораторной работе
дисциплины
"Информатика",
третий семестр
Преподаватель: Болтаев А.В.
Студенты: Степанов А.Г
Черепанов К.А.
Группа: Р-207
Екатеринбург
2000
Содержание
- Задание исследования3
- Подробное описание задачи и способы ее решения3
- Результаты исследований4
- Сравнение результатов12
- Список библиографических источников13
- Текст программы13
Задание исследования
Провести исследование внутренней сходимости численного интегрирования методом Симпсона и трапеций различных функций, задаваемых с помощью языка С.
Подробное описание задачи и способы ее решения
Необходимо провести исследования так называемой внутренней сходимости численного интегрирования методами Симсона и трапеций различных функций, задаваемых с помощью функций языка С. Предполагается, что отрезок интегрирования [a,b] разбит на n равных частей системой точек (сеткой).
Контроль внутренней сходимости заключается в циклическом вычислении приближенных значений интеграла для удваимого по сравнению со значением на предыдущем прохождении цикла числа n. Отношения абсолютной величины разности этих значений к абсолютной величине предыдущего приближенного значения принимается в качестве критерия достижения точности интеграла.
Построить зависимости количеств итераций от различных величин критерия точности.
Построить обратные зависимости критерия точноти от количества итераций.
Повторить все вышеуказанные исследования для случая, когда при вычислении критерия точности разность значений интеграла относится не к предыдущему значению, а к точному значению аналитически вычисленного интеграла.
Исследовать влияние увеличения верхнего предела интегрирования на точность (при прочих неизменных условиях)
Метод трапеций
, где
Метод Симпсона
, где
Результаты исследований
Таблица и график зависимости количества итераций от различных значений критерия точности
Для
Критерий точностиКоличество итераций-0,167663114-0,151891616-0,004693112-0,002653111-0,000263910-0,00017092-0,00012979-0,00005573-0,0000258-0,00001984-0,00000965-0,000003860150,000005270,07108913
Критерий точностиКоличество итераций-0,112727116-0,075028815-0,054067714-0,002141512-0,000571111-0,00004589-0,00003812-0,00001913-0,0000084-0,0000045-0,00000197-0,000000260,00000580,0002983100,016437713
Критерий точностиКоличество итераций-0,006670913-0,004236714-0,000356110-0,00000165-0,00000140,000000530,000000660,000000920,000000970,000022380,00005690,0002782110,0003474120,005293160,005326715
Критерий точностиКритерий точности-61,446979512-5,7140473-1,021575513-0,72414332-0,51211174-0,322264311-0,21636147-0,15366299-0,0930261140,0353183160,057059150,169737150,2025534100,250472860,62025928
Критерий точностиКоличество итераций-0,011930816-0,000783413-0,00000793-0,00000414-0,00000377-0,00000275-0,00000276-0,0000028-0,000001620,0000003100,000006290,0000385110,0000802120,0005452150,001668914
Критерий точностиКоличество итераций-0,002628616-0,001241614-0,00001183-0,00001074-0,00000465-0,00000469-0,00000286-0,00000217-0,000000520,0000011100,000001880,0000023110,000058120,0001049130,002792815
Таблица и график зависимости значений критерия точности от количества итераций
Для функции
По отношению к предыдущему значениюПо отношению к аналитическому значениюКритерий точностиКоличество итерацийКритерий точностиКоличество итераций-0,00017092-0,00019322-0,00005573-0,00006293-0,00001984-0,00002244-0,00000965-0,00001085-0,00000386-0,000004360,000005270,00000587-0,0000258-0,00002838-0,00012979-0,00014669-0,000263910-0,000298310-0,002653111-0,00299811-0,004693112-0,0052891120,071089130,079740313-0,167663114-0,201436514015015-0,151891616-0,151891616
Для функции
По отношению к предыдущему значениюПо отношению к аналитическому значениюКритерий точностиКоличество итерацийКритерий точностиКоличество итераций-0,00003812-0,00006662-0,00001913-0,00003353-0,0000084-0,00001414-0,0000045-0,00000695-0,00000026-0,00000046-0,00000197-0,000003370,00000580,00000888-0,00004589-0,000080290,0002983100,00052210-0,000571111-0,000999711-0,002141512-0,0037465120,0164377130,028695513-0,054067714-0,095937814-0,075028815-0,125933115-0,112727116-0,175012416
Сравнение результатов
Таблица сравнительных результатов
Метод трапеции n=1000000Метод Симпсона
n =1000000Аналитический результатФункцияПределы4,50514754,52401834,49980967f(x)=1/x0,1…..91,74914621,75007611,791756469f(x)=1/x*x0,3…..51,99918851,99995052f(x)=sin(x)0…….?-0,00005120,0000030f(x)=sin(2*x)0…….?0,28571570,28569350,285714285f(x)=sin(7*x)0…....?0,22220530,22221330,222222222f(x)=sin(9*x)0…....?Таблица влияния увеличения верхнего предела на точность интегрирования
Аналитическое значениеПрактическое значениеВерхний пределПогрешность4,499809674,52179969-0,021989934,6051701864,62496910-0,0197988144,7874917434,803941212-0,0164494574,9416424234,