Студопедия

Главная страница Случайная страница

КАТЕГОРИИ:

АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатикаИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханикаОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторикаСоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансыХимияЧерчениеЭкологияЭкономикаЭлектроника






Способ 2.






С помощью команды inttrap. Интегрирование экспериментальных данных с помощью трапецеидальной интерполяции (метод трапеций).

При вычислении интеграла между соседними узлами функция интерполируется линейно. Этот метод вычислений называется методом трапеций.

Вычислим интеграл от той же самой функции

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

--> y=exp((x-3).^2/8)

--> v=inttrap(x, y)

Получаем:

 

v =

4.8017553

Или

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= inttrap (x, y)

Получаем:

 

v =

4.80172

В первых двух способах вычисляется площадь заштрихованных фигур.

Способ 3. С помощью команды integrate. Это интегрирование по квадратуре. Может задаваться требуемая точность вычислений.

Синтаксис [x]=integrate(expr, v, x0, x1 [, ea [, er]])

Параметры

expr: подынтегральная функция;

v: переменная интегрирования;

x0, x1: пределы интегрирования;

ea, er: действительные числа. Первое из них – еа - предельная абсолютная ошибка. По умолчанию принимает значение 0. Второе - er -, - предельная относительная ошибка. По умолчанию принимает значение 1.d-8.

Пример: вычислить .

--> integrate('exp((x-3)^2/8)', 'x', 1, 5)

ans =

4.7798306

Можно это сделать и так.

Набираем и сохраняем в окне редактора под именем a.sci файл

function g=a(x)

g=exp((x-3).^2/8);

endfunction

Загружаем этот файл в среду Scilab (Load into Scilab).

Далее в строке ввода набираем:

--> integrate('a', 'x', 1, 5)

Получаем:

ans =

4.7798306


Поделиться с друзьями:

mylektsii.su - Мои Лекции - 2015-2024 год. (0.005 сек.)Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав Пожаловаться на материал