66521

Вычисление определенных интегралов

Лабораторная работа

Информатика, кибернетика и программирование

То, насколько точно методом Монте-Карло будет вычислен интеграл, зависит от количества поставленных точек и количества точек попавших в область интегрирования, поэтому при вычислении каждый раз значение интеграла будет отличаться от предыдущего.

Русский

2014-08-22

172 KB

4 чел.

Выполнил: Марудо А.В., 2 курс, 3 группа

Проверил: Шапочкина Ирина Викторовна

Лабораторная работа #3(вариант #26)

Тема: Вычисление определенных интегралов.

Цель: Вычислить значения интегралов на основании ряда значений подынтегральной функции методами Монте-Карло и Ньютона-Котеса.

Задание 1

Условие: Методом Монте-Карло найти объем тела, ограниченного снизу плоскостью , с боков плоскостями , , , , а сверху — эллиптическим параболоидом . Сравнить полученное значение с результатом, вычисленным средствами MathCAD.

Ход работы: 

Интеграл принимается, как площадь под графиком.

Пусть имеется n точек, которые «бросают» на область интегрирования. Координаты точек, принадлежащих области, выбираются случайным образом. Тогда

Код программы:

 begin

   E:=0;

   sum:=0;

   n:=StrToInt(Edit1.Text);

   randomize;

   //вычисляем значение интеграла

   for i:=1 to n do

    begin

     x:=random*3;

     y:=random*7;

     if y<=7 then

      begin

       inc(E);

       sum:=sum+(sqr(x)+sqr(y));

      end;

    end;

   II:=sum*3*7/E;

   Label1.Caption:=floattostr(II);

 end;

Полученные результаты:

Результаты MathCAD:

Вывод:

То, насколько точно методом Монте-Карло будет вычислен интеграл, зависит от количества поставленных точек и количества точек попавших в область интегрирования, поэтому при вычислении каждый раз значение интеграла будет отличаться от предыдущего. Но при большом количестве поставленных точек, получается результат близкий к истинному.

Задание 2

Условие: Вычислить значение определенного интеграла по квадратурным формулам Ньютона-Котеса с точностью , определяя шаг интегрирования по оценки остаточного члена:

Ход работы:

Вычислять будем по формулам закрытого типа, т.к. a и b являются узлами интерполяции. Воспользуемся формулой Симпсона для четного кол-ва узлов.

Код программы:

//функция под интегралом

   function f(x0: real): real;

   begin

     Result:=ln(3/2* exp(x0+0.3) )/ln(3* exp(-5*sqr(x0)+0.8) );

   end;

//интегрирование по формулам Ньютона-Котеса

  function integral(a1, b1: real; n1: integer): real;    

  var h1, sigma1, sigma2: real;

      i: integer;

  begin

   sigma1:=0;

   sigma2:=0;

   h1:=(b1-a1)/n1;

   for i:=1 to (n1-1) do

    begin

     if (i div 2)=(i/2) then

       begin

        sigma2:=sigma2+f(a1+i*h1);

       end

      else

       begin

        sigma1:=sigma1+f(a1+i*h1);

       end;

    end;

//формула Симпсона для четного количества узлов

   Result:=(h1/3)*(f(a1)+f(b1)+4*sigma1+2*sigma2);

  end;

begin

xa:=-0.21;

xb:=0.35;

e:=0.0005;

 q:=false;

//первоначально количество точек (кол-во в процессе интегрирования должно быть четн.)

 n:=2;  

while q=false do

  begin

//первая формула Рунге

  R:=(integral(xa,xb,n)-integral(xa,xb,n*2))/(2*2-1);  

  if abs(R)>e then

    n:=n*2

   else

    begin

     q:=true;

     S:=integral(xa,xb,n);

//вторая форулма Рунге

     S:=S+R;               

    end;

 end;

Полученные результаты:

Блок-схема метода Симпсона

Вывод:

Полученное значение интеграла совпадает со значением вычисленным в MathCAD с точностью до погрешности.

БГУ

Физический факультет

2011/2012 учебный год

Минск

PAGE   \* MERGEFORMAT 2


 

А также другие работы, которые могут Вас заинтересовать

12541. Работа с программными средствами Internet. Утилиты ping и traceroute 22.64 KB
  ЛАБОРАТОРНАЯ РАБОТА Работа с программными средствами Internet. Утилиты ping и traceroute 1. ЦЕЛЬ РАБОТЫ Исследование вероятностновременных характеристик сети с использованием утилиты ping исследование топологии фрагментов Internet с использованием утилиты traceroute. 2. ЗАДАНИЕ Н
12542. Работа с операционной системой Linux. Утилита netstat 23.83 KB
  ЛАБОРАТОРНАЯ РАБОТА Работа с операционной системой Linux. Утилита netstat 1. ЦЕЛЬ РАБОТЫ С помощью утилиты netstat исследовать состояние локальной IPсети. 2. ЗАДАНИЕ НА ЛАБОРАТОРНУЮ РАБОТУ 2.1. С помощью утилиты netstat получить список соединений открытых на сервере pds.sut.ru.
12543. РАБОТА С СЕТЕВЫМИ УТИЛИТАМИ И ПРОТОКОЛАМИ ПРИКЛАДНОГО УРОВНЯ 276.5 KB
  Работа с сетевыми утилитами и протоколами прикладного уровня Цель работы: освоить приемы работы с сервисными сетевыми утилитами прикладного уровня изучить команды основных утилит получить представление о методах работы с ними под управлением различных операционных ...
12544. Изучение адресации сети Интернет 149 KB
  Сети ЭВМ и ТК Лабораторная работа №1 Изучение адресации сети Интернет Цель работы: изучить адресацию сети Интернет Адресация сети Интернет В стеке протокола TCP/IP адресацию обеспечивает протокол IP. Согласно стандарту IP каждому хосту должно быть присвоено ун...
12545. Исследование протоколов IP-сетей 575.5 KB
  Лабораторная работа по курсу Вычислительные комплексы и сети Аппаратнопрограммные средства телекоммуникаций Исследование протоколов IPсетей Цель работы. Развитие практических навыков работы с протоколами стека ТСР/IP и исследование возможностей протоко...
12546. ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ ИЗГОТОВЛЕНИЯ ОДНОСЛОЙНЫХ ПЕЧАТНЫХ ПЛАТ 826.5 KB
  ЛАБОРАТОРНАЯ РАБОТА N1 ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ ИЗГОТОВЛЕНИЯ ОДНОСЛОЙНЫХ ПЕЧАТНЫХ ПЛАТ СОДЕРЖАНИЕ Цель работы 3 Теоретические сведения 4 Опис
12547. ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ ИЗГОТОВЛЕНИЯ ДВУСЛОЙНЫХ ПЕЧАТНЫХ ПЛАТ 1.15 MB
  ЛАБОРАТОРНАЯ РАБОТА N2 ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ ИЗГОТОВЛЕНИЯ ДВУСЛОЙНЫХ ПЕЧАТНЫХ ПЛАТ МЕТОДИЧЕСКОЕ РУКОВОДСТВО N СОДЕРЖАНИЕ Цель работы 3 Теоретические сведения 4 Описание лабораторного макета 36 Лабораторное задание 36 Порядок выпол...
12548. КОНСТРУКЦИЯ И ТЕХНОЛОГИЯ ПОВЕРХНОСТНОГО МОНТАЖА ИНТЕГРАЛЬНЫХ СХЕМ 14.46 MB
  ЛАБОРАТОРНАЯ РАБОТА № 5 КОНСТРУКЦИЯ И ТЕХНОЛОГИЯ ПОВЕРХНОСТНОГО МОНТАЖА ИНТЕГРАЛЬНЫХ СХЕМ МЕТОДИЧЕСКОЕ РУКОВОДСТВО СОДЕРЖАНИЕ Цель работы Теоретические сведения Описание лабораторного макета Лабораторное задание...
12549. ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ ИЗГОТОВЛЕНИЯ МНОГОСЛОЙНЫХ ПЕЧАТНЫХ ПЛАТ 1.11 MB
  ЛАБОРАТОРНАЯ РАБОТА N3 ТЕХНОЛОГИЧЕСКИЕ ПРОЦЕССЫ ИЗГОТОВЛЕНИЯ МНОГОСЛОЙНЫХ ПЕЧАТНЫХ ПЛАТ МЕТОДИЧЕСКОЕ РУКОВОДСТВО СОДЕРЖАНИЕ Цель работы3 Теоретические сведения4 Описание лабораторного макета58 Лабораторное задание 58 Порядок выполнения работы58 Треб...