Способ 1.

Численное интегрирование

При вычислении определенных интегралов от функций, заданных в виде таблицы или в явном виде (), одним из численных методов используют функции intsplin, inttrap, integrate, intg.

Способ 1.

С помощью команды intsplin. Это интегрирование экспериментальных данных с помощью сплайн-интерполяции. Известно значение интегрируемой функции в дискретных точках (узлах).

Пример: вычислить интеграл от таблично заданной функции

x		1.4	1.8	2.2	2.6		3.4	3.8	4.2	4.6
	1.6487	1.3771	1.1972	1.0833	1.0202		1.0202	1.0833	1.1972	1.3771	1.6487

--> x=1:.4: 5;

--> y=exp((x-3).^2/8)//значения у в таблице получены табулированием этой функции

--> v=intsplin(x, y)

Получаем:

v =

4.7799684

Или:

--> x=[1 1.4 1.8 2.2 2.6 3 3.4 3.8 4.2 4.6 5];

--> y=[1.6487 1.3771 1.1972 1.0833 1.0202 1 1.0202 1.0833 1.1972 1.3771 1.6487];

--> v=intsplin(x, y)

v =

4.7799328

clc

x=[1 1.4 1.8 2.2 2.6 3 3.4 3.8 4.2 4.6 5];

y=[1.6487 1.3771 1.1972 1.0833 1.0202 1 1.0202 1.0833 1.1972 1.3771 1.6487];

plot2d(x, y, -4); //График экспериментальных данных

koeff=splin(x, y);

plot2d(x, y, -3); //Нанесение точек на график

//Построение кубического сплайна

t=min(x): 0.01: max(x);

ptd=interp(t, x, y, koeff);

plot2d(t, ptd);