Учебная работа № 1413. Расчет двойного интеграла при помощи метода Симпсона

1 Звезда2 Звезды3 Звезды4 Звезды5 Звезд (6 оценок, среднее: 4,67 из 5)
Загрузка...
Контрольные рефераты

Учебная работа № 1413. Расчет двойного интеграла при помощи метода Симпсона

/*****************************************************************************

* .FILE : numeric.c

* .TITLE : Расчет двойного интеграла при помощи метода Симпсона (парабол)

* .DESCR : Курсовой проект по численным методам (1994 год)

* :

* : Выполнил: Петренко В.С. (гр. ПС301)

* : Проверил: Панюков А.В.

* :

* .NOTE : NOT FOR RENTAL OR SALE.

* : FEDERAL LAW PROVIDES SEVERE CIVIL & CRIMINAL PENALTIES FOR

* : UNAUTHORIZED DUPLICATION OR DISTRIBUTION.

* :

* : (C) ’94 by P$P

*****************************************************************************/

#include <math.h>

#include <stdio.h>

/*****************************************************************************

* .NAME : m_Simpson

* .TITLE : Расчет интеграла методом Симпсона (парабол)

* .DESCR :

* :

* .PARAMS : double m_Simpson (double (*func) (double, double),

* : double t_fix, double t_limit, int N);

* : double (*func) (double, double) подынтегральная фя

* : double t_fix фиксированный первый аргумент

* : double t_limit верхний предел интегрирования,

* : нижний равен t_limit

* : int N число точек разбиения

* .RETURN : Значение вычисленного интеграла

*****************************************************************************/

double m_Simpson (double (*func) (double, double),

double t_fix, double t_limit, int N)

{

double sum1 = 0; /* ¬ */

double sum2 = 0; /* ¦ временные переменные */

double sum3 = 0; /* */

double sum; /* конечный результат */

double h = (2 * t_limit) / N; /* шаг сетки */

int i; /* временная */

sum1 = (*func) (t_fix, t_limit) + (*func) (t_fix, +t_limit);

for (i = 1; i <= N 1; i++) sum2 += (*func) (t_fix, t_limit + (i * h));

sum2 *= 2;

for (i = 1; i <= N; i++) sum3 += (*func) (t_fix, t_limit + ((i0.5) * h));

sum3 *= 4;

sum = sum1 + sum2 + sum3;

sum = (h / 6) * sum;

return sum;

}

/* */

/* Глобальные переменные */

/* ~~~~~~~~~~~~~~~~~~~~~ */

#define PI 3.1415926536 /* число П */

double k; /* параметр функции задается пользователем */

int N_MAX; /* число узлов сетки разбиения */

double (*currFunc) (double, double); /* выбранная пользователем функция */

double f1 (double x, double y);

double f2 (double x, double y);

double f3 (double x, double y);

/*****************************************************************************

* .NAME : double F (double dummy, double t)

* .TITLE : Вычисляет внутренний интеграл (G (t)).

* : См. текст курсового проекта.

* .DESCR : П/2

* :

* : первая вычисляемая функция ¦ G (t + П/2) * sin (t + П/2) dt

* : t

* : П/2

* : где G (t) = ¦ currFunc (t, tau) dtau

* :

* : t

* .PARAMS : double F (double dummy, double t);

* : double dummy фиктивный первый аргумент, при вызове этой

* : функции он не используется, т.к. она

* : функция одного аргумента

* : double действительный второй аргумент

* .RETURN : Значение функции: G (t) * sin (k * t);

*****************************************************************************/

double F (double dummy, double t)

{

double G;

t = t + PI / 2; /* сдвижка начала координат, чтобы пределы */

/* были симметричны (в нашем случае на П/2) */

G = m_Simpson (currFunc, t, t, N_MAX);

return G * sin (k * t);

}

/*****************************************************************************

* .NAME : main

* .TITLE : Основная диалоговая функция.

* .DESCR : Запрашивается интересующая пользователя функция,

* : параметр k и число узлов сетки N_MAX.

* : Выводит на экран вычисленное значение интеграла и

* : два справочных значения П и П/2.

* :

* .PARAMS : void main (void);

* .RETURN :

*****************************************************************************/

void main (void)

{

double integral; /* значение вычисленного интеграла */

int selection; /* номер выбранной функции */

/* массив доступных функций */

double (*functions []) (double, double) = { f1, f2, f3 };

printf («\n Вычисление интеграла методом Симпсона (парабол) «);

printf («\n ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ «);

printf («\n «);

printf («\n I = ¦¦ sin k(x + y) f (x, y) dx dy «);

printf («\n «);

printf («\n D «);

printf («\n где D = { (x, y): x, y >= 0; x + y <= П }, f Е C (D)»);

printf («\n»);

printf («\nДля какой функции рассчитывать: «);

printf («\n~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ «);

printf («\n 1) f (x, y) = 0.5 * cos (y) «);

printf («\n 0; k != 1 «);

printf («\n ¦ sin x * sin (kx) dx => ¦ «);

printf («\n L П/2; k = 1 «);

printf («\n «);

printf («\n 2) f (x, y) = 0.5 sin (y) «);

printf («\n «);

printf («\n ¦ x * sin (kx) dx =====> П; k = 1 «);

printf («\n «);

printf («\n «);

printf («\n 3) f (x, y) = sqrt (x * x + y * y)»);

printf («\n»);

do

{

printf («Ваш выбор: «);

scanf («%d», &selection);

} while (!(1 <= selection && selection <= 3));

printf («Параметр k: «);

scanf («%lg», &k);

do

{

printf («Число узлов сетки N: «);

scanf («%d», &N_MAX);

} while (!(N_MAX > 0));

printf («\n»);

printf («\n Расчет интеграла …»);

currFunc = functions [selection 1]; /* текущая функция */

integral = m_Simpson (F, 0, PI / 2, N_MAX); /* вычисляем интеграл */

printf («\n Значение интеграла равно: %.12lg», integral); /* вывод */

printf («\n Величины: П = %.12lg; П/2 = %.12lg», PI, PI / 2);

}

/*****************************************************************************

* .FILE : func.c

* .TITLE : Содержит функции пользователя, которые можно

* : изменять без перекомпиляции основной программы

* .DESCR : После изменения этого модуля его необходимо перекомпилировать

* : и слинковать с numeric.obj

* :

* :

* .NOTE : NOT FOR RENTAL OR SALE.

* : FEDERAL LAW PROVIDES SEVERE CIVIL & CRIMINAL PENALTIES FOR

* : UNAUTHORIZED DUPLICATION OR DISTRIBUTION.

* :

* : (C) ’94 by P$P

*****************************************************************************/

#include <math.h>

/* выбираемая пользователем функция No.1 */

double f1 (double x, double y)

{ return 0.5 * cos (y); }

/* выбираемая пользователем функция No.2 */

double f2 (double x, double y)

{ return 0.5 sin (y); }

/* выбираемая пользователем функция No.3 */

double f3 (double x, double y)

{ return sqrt (x * x + y * y); }

Учебная работа № 1413. Расчет двойного интеграла при помощи метода Симпсона