3853

Обчислення означених інтегралів

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

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

Обчислення означених інтегралів Мета роботи: вивчити методи наближених обчислень і запрограмувати алгоритми обчислення означених інтегралів. Короткі теоретичні відомості Формули прямокутників. Нехай на відрізку задана неперервна функція . Потрібно о...

Украинкский

2012-11-09

79 KB

8 чел.

Обчислення означених інтегралів

Мета роботи: вивчити методи наближених обчислень і запрограмувати алгоритми обчислення означених інтегралів.

Короткі теоретичні відомості

Формули прямокутників.

Нехай на відрізку задана неперервна функція . Потрібно обчислити інтеграл

Розіб’ємо відрізок  на n  рівних частин точками , i=0,1,…n-1, довжина кожної з яких дорівнює  . Через  позначимо значення функції  в точках  і складемо суми

 або  

            Кожна з цих сум є інтегральною сумою для  на відрізку і тому наближено виражають означений інтеграл:

                                                       (1)

                                                                  (1/)

            Ці формули називаються формулами прямокутників.

2. Формула трапецій.

                                          (3)

            3. Формула парабол (Сімпсона).

4. Формула трьох восьмих:

  Якщо в формулі Ньютона-Котеса взяти n = 3, тобто функцію f(x) замінити інтерполяційним багаточленом третього степеня, побудованим за значення функції f(x) у точках x0=a, x1=a+h, x2=a+2h, x3=b, h=(b-a )/3. то одержимо таку квадратурну формулу:

де

Ця квадратурна формула називається малою квадратурною формулою трьох восьмих. Використовуючи цю формулу, легко записати велику квадратурну формулу трьох восьмих.

Варіант 17

Обчислити інтеграл методом прямокутників, трапецій, парабол, трьох восьмих, Монте-Карло.

Заданий інтеграл обчислити наближено та точно.

1.

2.

3.

Текст програми:

//---------------------------------------------------------------------------

#include <vcl.h>

#include <math.h>

#pragma hdrstop

#include "IntegralUn.h"

//---------------------------------------------------------------------------

#pragma package(smart_init)

#pragma resource "*.dfm"

TForm1 *Form1;

//---------------------------------------------------------------------------

__fastcall TForm1::TForm1(TComponent* Owner)

 : TForm(Owner)

{  }

//---------------------------------------------------------------------------

float __fastcall TForm1::func(float arg1, float arg2=0)

{

 switch(RadioGroup1->ItemIndex)

 {

   case 0: return (sin(arg1/17))/(cos(arg1/17));

   case 1: return arg2*(sin(arg1/17))/(cos(arg1/17));

   case 2: return atan(17*sin(arg1));

 }

}

void __fastcall TForm1::Edit1KeyPress(TObject *Sender, char &Key)

{

 if ((Key < '0' || Key > '9') && Key != 8 && Key != ',' && Key != '-') Key= 0;

}

//---------------------------------------------------------------------------

void __fastcall TForm1::Button1Click(TObject *Sender)

{

 int i, part = StrToFloat(Edit5->Text);

 float xa = StrToFloat(Edit1->Text),

       xb = StrToFloat(Edit2->Text),

       ya = StrToFloat(Edit4->Text),

       yb = StrToFloat(Edit3->Text),

       S=0;

 int k1=0,k2=0;

 float xh;

if(RadioGroup1->ItemIndex!=1)

{

 //ліві

 xh=(xb-xa)/part;

 for (i=0;i<part;i++)

   S+=func(xa+i*xh)*xh;

 LabeledEdit1->Text=FloatToStrF(S,ffFixed,14,6);

 //праві

 S=0;

 for (i=1;i<=part;i++)

   S+=func(xa+i*xh)*xh;

 LabeledEdit2->Text=FloatToStrF(S,ffFixed,14,6);

 //середні

 S=0;

 for (i=0;i<part;i++)

   S+=func(((xa+(i+1)*xh)+(xa+i*xh))/2)*xh;

 LabeledEdit3->Text=FloatToStrF(S,ffFixed,14,6);

 //трапеції

 S=0;

 for (i=1;i<part;i++)

   S+=func(xa+i*xh);

 S+=(func(xa)+func(xb))/2;

 S*=(xb-xa)/part;

 LabeledEdit4->Text=FloatToStrF(S,ffFixed,14,6);

 //парабола

 S=0;

 xh=(xb-xa)/(2*part);

 for (i=1;i<2*part;i++)

 {

   if (i%2) S+=4*func(xa+i*xh);

   else S+=2*func(xa+i*xh);

 }

 S+=func(xa)+func(xb);

 S*=xh/3;

 LabeledEdit5->Text=FloatToStrF(S,ffFixed,14,6);

 //три-восьмих

 S=0;

 xh=(xb-xa)/(3*part);

 for (i=1;i<3*part;i++)

 {

   if (i%3) S+=3*func(xa+i*xh);

   else S+=2*func(xa+i*xh);

 }

 S+=func(xa)+func(xb);

 S*=3/(float)8*xh;

 LabeledEdit6->Text=FloatToStrF(S,ffFixed,14,6);

 //Монте-Карло

   srand(time(NULL));

   rand();

 S=0;

 part=1000000;

 xh=(xb-xa)/(3*part);

 for (i=0;i<part;i++)

 {

   float x = (float)rand()/(float)RAND_MAX*xb+xa,

         y = (float)rand()/(float)RAND_MAX*100-50;

   if ((func(x)>=0) && (y>=0) && (y<=func(x))) k1++;

   if ((func(x)<0)  && (y<0) && (y>func(x))) k2++;

 }

 S=(xb-xa)*100*k1/part-(xb-xa)*100*k2/part;

 LabeledEdit7->Text=FloatToStrF(S,ffFixed,14,6);

}

else

{

 //Монте-Карло

   srand(time(NULL));

   rand();

 S=0;

 part=1000000;

 k1=0,k2=0;

 xh=(xb-xa)/(3*part);

 for (i=0;i<part;i++)

 {

   float x = (float)rand()/(float)RAND_MAX*xb+xa,

         y = (float)rand()/(float)RAND_MAX*yb+ya,

         z = (float)rand()/(float)RAND_MAX*100-50;

   if ((func(x,y)>=0) && (z>=0) && (z<=func(x,y))) k1++;

   if ((func(x,y)<0)  && (z<0) && (z>func(x,y))) k2++;

 }

 S=(xb-xa)*(yb-ya)*100*k1/part-(xb-xa)*(yb-ya)*100*k2/part;

 LabeledEdit7->Text=FloatToStrF(S,ffFixed,14,6);

}

}

//---------------------------------------------------------------------------

Результат виконання програми:

Висновок: На цій лабораторній роботі я вивчив методи наближених обчислень і запрограмувати алгоритми обчислення означених інтегралів.


 

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

35073. Теория государства и права. Сборник статей 2.92 MB
  и права является как бы деятельностью второго порядка способом методом регуляции. Матузов вслед за многими другими авторами в своей работе призывает разграничивать понятия правовая коллизия и правовой конфликт указывая на то что конфликт это противоречие между людьми в отличие от коллизий представляющих собой противоречия между нормами права способные служить причиной поводом для конфликта. Она появилась еще во времена зарождения самого права.
35074. МОТИВАЦИЯ И СТИМУЛИРОВАНИЕ ПЕРСОНАЛА 2.48 MB
  В ней рассмотрены основные концепции и психологические аспекты мотивации труда методы стимулирования работников и мотивационные ресурсы управления организацией. Исследования мотивации труда проведенные автором на предприятиях различных форм собственности и сфер деятельности позволили ему предложить новую систему диагностики мотивации труда работников. Они позволят руководителям выбрать из многообразия ресурсов методов и схем управления наиболее близкую для себя систему формирования мотивационного механизма и учтя специфические...
35075. СПРАВОЧНИК ПО ОТДЕЛКЕ МЕБЕЛИ 6.11 MB
  При первом способе необходимое по расчету количество красителя 01 05 к массе лака рабочей вязкости предварительно растворяют в растворителе. Некоторые красители изменяют цвет под действием перекисных соединений содержащихся в полиэфирных лаках [43]. Ее изготовляют из суховальцованных паст СВП спирторастворимых красителей и нитролака. Грунтовку НЦ 0140 можно применять при отделке как полиэфирными так и нитроцеллюлозными лаками.
35076. Экономическое содержание налогов 33.93 KB
  а Государство и налоги б Влияние налогов на экономику и социальную жизнь . Налоги как основной источник образования государственных финансов через изъятие в виде обязательных платежей части общественного продукта известны с незапамятных времен. В разные времена и в разных странах налоги собирались различными способами. По этой причине налоги в истории всегда ассоциировались с несправедливостью.
35077. Змішана економіка та її типи 229 KB
  Змішана економіка та її типи План Сутність і ознаки змішаної економікиумови і принципи виникнення. Особливості становлення ринкової економіки в Україні. Моделі і шляхи переходу до ринкової економіки. Сутність і ознаки змішаної економікиумови і принципи виникнення.
35078. Франчайзинг, проблемы его развития в России 216.5 KB
  Развитие франчайзинга как системы организации бизнеса насчитывает уже полтора века с тех пор когда Зингер впервые применил его для сбыта швейных машин. В России интерес к этому экономическому инструменту возродился около двадцати лет назад когда практически одновременно возникло несколько франчайзинговых систем отечественных и с участием иностранного капитала. Однако надежды на быстрое развитие франчайзинга не оправдались. В 2011 году когда российское правительство делает упор на развитие инновационных идей и переходу от поставщика...
35079. СВОЙСТВА ЛИЧНОСТИ 235.5 KB
  Из истории учений о типах темперамента. Физиологическая основа темперамента. Психологическая характеристика типов темперамента. Несмотря на то что предпринимались неоднократные и постоянные попытки исследовать проблему темперамента до сих пор эта проблема относится к разряду спорных и до конца не решенных проблем современной психологической науки.
35080. Инвестиции. Понятие и классификация инвестиций 103 KB
  Понятие и классификация инвестиций Экономическая деятельность отдельных хозяйствующих субъектов и страны в целом в значительной мере характеризуется объемом и формами осуществляемых инвестиций. Источником прироста капитала и движущим мотивом осуществления инвестиций является получаемая от них прибыль. Чистые инвестиции представляют собой сумму валовых инвестиций уменьшенную на сумму амортизационных отчислений в определенном периоде. Факторы влияющие на объемы инвестиций довольно многочисленны.
35081. Теория государства и права. Конспект лекций 1.58 MB
  Предмет и методология теории государства и права. Место и роль теории государства и права в системе других наук. Основные категории и структура курса теории государства и права.