Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Учебники 80153

.pdf
Скачиваний:
1
Добавлен:
01.05.2022
Размер:
632.5 Кб
Скачать

пеции имеет вид I

i

 

h

y

4y

y

. Применяя эту форму-

 

 

3

i 1

i

i 1

 

лу на каждом элементарном отрезке, получим составную квад-

ратурную формулу Симпсона:

 

 

 

I S

 

 

 

h

y 4y y

 

 

y

 

4y y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

3

 

0

1

 

2

 

 

 

 

 

 

2

 

 

 

3

 

4

 

 

 

 

 

 

 

 

 

 

y

n 2

4y

n 1

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

h

y 4y 2y

2

 

4y 2y 2y

n 2

4y

n 1

y

n

 

 

 

3

0

1

 

 

 

 

 

 

 

3

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

n 2

 

 

 

 

 

 

 

 

 

n 2 1

 

 

 

 

 

 

 

 

 

 

 

 

y y 4

 

y

2i 1

2

 

y

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

0

 

 

n

 

 

 

 

 

 

 

 

2i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

 

 

 

 

р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

а

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

y4

 

 

 

 

 

 

 

 

 

ар

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yп

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

x1

 

 

x2

 

 

 

 

 

 

 

x3

 

 

 

 

x4

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если проводить параболы через точки с абсциссами xi 1, xi 1 h2, xi и соответствующими ординатами yi 1, yi 12 , yi , то формула Симпсона примет вид:

 

 

 

h

 

 

n

 

n 1

 

 

I S

C

 

 

y y 4

y

2

 

y

.

 

 

 

6

 

0

n

i 1 2

 

i

 

 

 

 

 

i 1

 

i 1

 

 

Отметим, что формула Симпсона имеет четвертый порядок точности.

В курсе «Теория вероятностей и математическая статистика» дифференцирование широко используется для нахождения плотности распределения вероятностей при известной функции распределения вероятностей случайной величины

pX x

dFX x

 

и

случайного вектора

непрерывного

типа

dx

 

 

 

 

 

 

 

 

 

2F

(x, y)

 

 

pXY (x, y)

 

XY

 

 

. Интегрирование

используется

для

 

x y

 

 

 

 

 

 

 

 

нахождения функции распределения вероятностей случайной

x

 

 

 

 

величины FX x

pX t dt

и

случайного вектора непре-

 

 

 

 

 

x

y

 

 

рывного типа FXY (x,

y)

 

pXY (s, t)dsdt , для проверки ус-

 

 

 

 

 

 

 

 

 

ловия нормировки в одномерном

pX x dx FX 1 и

 

 

 

 

 

 

 

 

 

 

двумерном случае

dx pXY (x,

y)dy 1, для нахождения ве-

 

 

 

 

 

роятности попадания случайной величины в определенный ин-

x2

тервал P x1 X x2 F x2 F x1 pX x dx или в оп-

x1

ределенную область P (X, Y) D pXY (x, y)dxdy , для оп-

D

ределения плотности распределения вероятностей отдельных компонент случайного вектора:

pX (x) pXY (x, y)dy и

pY (y) pXY (x, y)dx .

 

 

39

40

Задание к выполнению лабораторной работы.

1)Найти аналитические выражения для первой и второй производной функции f (x) своего варианта (табл. 1).

2)Вычислить точные и приближенные значения первой и второй производной функции f (x) в середине отрезка локали-

зации корня (см. лабораторную работу № 1).

 

3)

Найти аналитические

выражения для производных

функции

g(x, y) :

g

 

,

g

,

2g

,

 

2g

,

2g

. Задать функцию

x

y

x2

y2

x y

 

 

 

 

 

 

 

в виде g(x, y) arctg

x

, где n – номер варианта.

 

 

 

 

 

ny

 

 

 

 

 

 

 

 

4) Найти аналитическое выражение f (x)dx, точное зна-

1 n

 

 

 

 

 

 

 

 

 

 

 

 

1 n

чение f (x)dx, приближенное значение f (x)dx, точные

0

0

 

0

значения

f (x)dx

и

f (x)dx.

Задать

функцию в

виде

 

0

 

 

 

 

 

 

 

f (x) nxe nx2

, n – номер варианта.

 

 

5) Найти аналитическое выражение

f (x, y)dxdy , точ-

 

 

 

 

 

1 n

x 25

 

 

ное и приближенное значения

 

dx f (x, y)dy. Задать функ-

 

 

 

 

 

1

0

 

 

цию в виде

f (x, y) nx2

xy2

n, n – номер варианта.

 

Порядок выполнения лабораторной работы

 

1) Создаем функцию f (x), определяем функцию

f1(x)

как производную f

 

с помощью встроенных функций diff

(x)

и define, определяем функцию f 2(x) как производную f (x)

спомощью встроенных функций diff и define.

2)Вычисляем точные значения первой и второй производной функции f (x) в середине отрезка локализации корня.

Определяем шаг h как длину отрезка локализации, деленную на 10. Создаем выражение f1p(x) для вычисления приближенного значения первой производной согласно формуле (3.3) и вычисляем приближенное значение первой производной функции f (x) в середине отрезка локализации корня. Создаем выражение f 2p(x) для вычисления приближенного значения первой производной согласно формуле (3.4) и вычисляем при-

ближенное значение второй производной функции

f (x)

в се-

редине отрезка локализации корня.

 

 

 

 

 

 

 

3) Создаем функцию

g(x, y)

и с помощью встроенной

функции

diff

находим аналитические выражения для произ-

 

g1x(x, y)

 

g

 

g1y(x, y)

g

g2x(x, y)

2g

водных:

 

 

,

 

,

 

,

 

 

 

2g

 

 

x

 

 

2g

 

 

y

 

 

x2

g2y(x, y)

 

, g2xy(x, y)

.

 

 

 

 

 

 

y2

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4) Создаем функцию

f (x). Находим аналитическое вы-

ражение

f (x)dx

 

с

 

помощью

встроенной

функции

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

integrate( f (x),x), а

точные значения

f (x)dx с

помощью

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

встроенной функции integrate( f (x),x,a,b) . Приближенное зна-

b

чение f (x)dx находим с помощью встроенной функции

a

romberg( f (x),x,a,b) из пакета romberg .

41

42

5) Создаем функцию f (x, y). Находим аналитическое

выражение f (x, y)dxdy dx f (x, y)dy с помощью встро-

енной функции integrate(integrate( f (x),y),x) . Точное значение

b y2(x)

интеграла dx f (x, y)dy находим с помощью встроенной

a y1(x)

функции integrate(integrate( f (x), y, y1(x), y2(x)),x,a,b), а при-

ближенное значение находим с помощью встроенной функции romberg(romberg( f (x), y, y1(x),y2(x)),x,a,b).

Пример выполнения лабораторной работы

1) (%i) f(x):=log(cos(x^2+1)); (%i) define(f1(x),diff(f(x),x));

(%i) define(f2(x),diff(f(x),x,2));

2)(%i) f1(0.25); -0.897466078632706 (%i) f2(0.25); -4.645309676827188 (%i) h:0.05;

(%i) f1p(x):=(f(x+h)-f(x-h))/(2*h); (%i) f1p(0.25);

-0.903569482862924

(%i) f2p(x):=(f(x-h)-2*f(x)+f(x+h))/(h^2); (%i) f2p(0.25);

-4.66663672974179

3)(%i) g(x,y):=atan(x/(25*y));

(%i) g1x(x,y):=diff(g(x,y),x,1); (%i) g1x(x,y);

(%i) g1y(x,y):=diff(g(x,y),y,1); (%i) g1y(x,y);

(%i) g2x:diff(g(x,y),x,2);

(%i) g2y:diff(g(x,y),y,2);

(%i) g2xy:diff(g1x(x,y),y);

4) (%i) f(x):=25*x*exp(-25*x^2); (%i) integrate(f(x),x);

(%i) integrate(f(x),x,0,1/25);

43

44

(%i) float(%); 0.0196052804238384 (%i) load(romberg);

(%i) romberg(f(x),x,0,1/25); 0.0196053203562632

(%i) integrate(f(x),x,0,inf);

(%i) integrate(f(x),x,minf,0);

1

2

5) (%i) f(x,y):=25*x^2+x*y^2/25; (%i) integrate(integrate(f(x,y),y),x);

(%i) integrate(integrate(f(x,y),y,0,x/25), x,1,1+25);

(%i) float(%); 114245.7777546667

(%i) romberg(romberg(f(x,y),y,0,x/25),x, 1,1+25);

114245.7820949445

Контрольные вопросы

1. В каких случаях применяется численное дифференцирование?

2.Запишите простейшие формулы для вычисления производных первого и второго порядка.

3.Дайте геометрическую интерпретацию формул для производной первого порядка.

4.Каков порядок точности этих формул?

5.Почему приближенное дифференцирование является некорректной операцией?

6.Как реализуется вычисление производных с помощью стандартных средств пакета Maxima?

7.В каких случаях применяется численное интегрирова-

ние?

8.Как разбивается площадь исходной криволинейной трапеции на сумму площадей элементарных криволинейных трапеций?

9.Как вычисляется площадь элементарной криволинейной трапеции в методе левых прямоугольников?

10.Как вычисляется площадь элементарной криволинейной трапеции в методе правых прямоугольников?

11.Как вычисляется площадь элементарной криволинейной трапеции в методе центральных прямоугольников?

12.Как вычисляется площадь элементарной криволинейной трапеции в методе трапеций?

13.Как вычисляется площадь элементарной криволинейной трапеции в методе Симпсона?

14.Какой порядок точности имеют методы прямоугольников, трапеций, Симпсона?

15.Как реализуется точное вычисление интегралов с помощью стандартных средств пакета Maxima?

16.Как реализуется точное вычисление интегралов с помощью стандартных средств пакета Maxima?

17.В каких случаях используется дифференцирование в курсе «Теория вероятностей и математическая статистика»?

18.В каких случаях используется интегрирование в курсе «Теория вероятностей и математическая статистика»?

45

46

Форма отчетности: Отчет должен содержать название лабораторной работы, задание для конкретного варианта, все командные строки, результаты расчетов аналитические и численные (в десятичных дробях три цифры после запятой).

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1.Плис А.И. Лабораторный практикум по высшей математике / А.И. Плис, Н.А. Сливина. М.: Высш. шк., 1994.

2.Амосов А.А. Вычислительные методы для инженеров / А.А. Амосов, Ю.А. Дубинский, Н.В. Копченова. М.: Высш.

шк., 1994.

3.Бабенко К.И. Основы численного анализа / К.И. Бабен-

ко. М.: Наука, 1986.

4.Бахвалов Н.С. Численные методы / Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков. М.: Наука, 1987.

5.Боглаев Ю.П. Вычислительная математика и программирование / Ю.П. Боглаев. М.: Высш. шк., 1990.

6.Волков Е.А. Численные методы / Е.А. Волков. М.:

Наука, 1987.

7.Демидович Б.П. Основы вычислительной математики / Б.П. Демидович, И.А. Марон. М.: Наука, 1966.

8.Иванов В.В. Методы вычислений на ЭВМ / В.В. Иванов. Киев: Наукова думка, 1986.

9.Калиткин Н.Н. Численные методы / Н.Н. Калиткин. М.:

Наука, 1978.

10.Марчук Г.И. Методы вычислительной математики / Г.И. Марчук. М.: Наука, 1989.

11.Самарский А.А. Численные методы / А.А. Самарский, А.В. Гулин. М.: Наука, 1989.

12.Турчак Л.И. Основы численных методов / Л.И. Турчак. М.: Наука, 1987.

13.Воробьева Г.Н. Практикум по вычислительной математике / Г.Н. Воробьева, А.Н. Данилова. М.: Высш. шк., 1990.

14.Гмурман В.Е. Теория вероятностей и математическая статистика / В.Е. Гмурман. М.: Высш. шк., 1998.

15.Вентцель Е.С. Теория вероятностей / Е.С. Вентцель.

М.: Наука, 1999.

16.Кремер Н.Ш. Теория вероятностей и математическая статистика / Н.Ш. Кремер. М.: ЮНИТИ-ДАНА, 2004.

17.Кибзун А.И. Теория вероятностей и математическая статистика. Базовый курс с примерами и задачами / А.И. Кибзун, Е.Р. Горяинова, А.В. Наумов. М.: Физматлит, 2002.

18.Андронов А.М. Теория вероятностей и математическая статистика / А.М. Андронов, Е.А. Копытов, Л.Я. Гринглаз. СПб.: Питер, 2004.

19.Письменный Д.Т. Конспект лекций по теории вероятностей и математической статистике / Д.Т. Письменный. М.: Айрис-пресс, 2004.

20.Дубровская А.П. Теория вероятностей и элементы математической статистики / А.П. Дубровская, В.И. Минаков. Воронеж: ВПИ, 1993.

21.Глушко Е.Г. Элементы теории вероятностей и математической статистики / Е.Г. Глушко, А.П. Дубровская. Воронеж: ВГТУ, 2004.

22.Колемаев В.А. Теория вероятностей и математическая статистика / В.А. Колемаев, О.В. Староверов, В.Б. Турундаевский. М.: Высш. шк., 1991.

23.Чистяков В.П. Курс теории вероятностей / В.П. Чистяков. М.: Наука, 1987.

24.Гурский Е.И. Теория вероятностей с элементами математической статистики / Е.И. Гурский. М.: Наука, 1971.

25.Захаров В.К. Теория вероятностей / В.К. Захаров, Б.А. Севостьянов, В.П. Чистяков. М.: Наука, 1983.

26.Гнеденко Б.В. Курс теории вероятностей / Б.В. Гнеденко. М.: Наука, 1988.

47

48

27.Прохоров Ю.В. Теория вероятностей / Ю.В. Прохоров, Ю.А. Розанов. М.: Наука, 1973.

28.Боровков А.А. Курс теории вероятностей / А.А. Боровков. М.: Наука, 1976.

29.Пугачев В.С. Теория вероятностей и математическая статистика / В.С. Пугачев. М.: Наука, 1979.

30.Ивашев-Мусатов О.С. Теория вероятностей и математическая статистика / О.С. Ивашев-Мусатов. М.: Наука, 1979.

31.Дежин В.В. Использование математического моделирования при решении экономических вероятностных задач: методические указания для организации самостоятельной работы по курсу «Теория вероятностей и математическая статистика» студентов направления 080100 «Экономика» профиля «Экономика предприятий и организаций» очной формы обучения. Ч. 1 / В.В. Дежин. Воронеж: ВГТУ, 2013. (№ 150-2013)

32.Дежин В.В. Использование математического моделирования при решении экономических вероятностных задач: методические указания для организации самостоятельной работы по курсу «Теория вероятностей и математическая статистика» студентов направления 080100 «Экономика» профиля «Экономика предприятий и организаций» очной формы обучения. Ч. 2 / В.В. Дежин. Воронеж: ВГТУ, 2013. (№ 151-2013)

33.Дежин В.В. Использование математического моделирования при решении экономических вероятностных задач: методические указания для организации самостоятельной работы по курсу «Теория вероятностей и математическая статистика» студентов направления 080100 «Экономика» профиля «Экономика предприятий и организаций» очной формы обучения. Ч. 3 / В.В. Дежин. Воронеж: ВГТУ, 2013. (№ 152-2013)

34.Дежин В.В. Использование математического моделирования при решении экономических вероятностных задач: методические указания для организации самостоятельной работы по курсу «Теория вероятностей и математическая статистика» студентов направления 080100 «Экономика» профиля

«Экономика предприятий и организаций» очной формы обучения. Ч. 4 / В.В. Дежин. Воронеж: ВГТУ, 2013. (№ 153-2013)

35. Дежин В.В. Решение задач по теории вероятностей и математической статистики с экономическим содержанием: методические указания для организации самостоятельной работы по курсу «Математика» для студентов специальности 080502 «Экономика и управление на предприятиях машиностроения» очной формы обучения. Ч. 1 / В.В. Дежин, Т.А. Надеина. Воронеж: ВГТУ, 2009. (№ 588-2009)

СОДЕРЖАНИЕ

Введение в пакет Maxima……………………………...…..1

Лабораторная работа № 1. Табулирование функции на интервале. Построение графика функции y f x .

Решение уравнения f x 0………………….….……………..6

Лабораторная работа № 2. Решение систем линейных уравнений………………………………...……..……14

Лабораторная работа № 3. Дифференцирование и интегрирование…………………….……………...…………..31

Библиографический список……..……………..………....47

49

50

МЕТОДИЧЕСКИЕ УКАЗАНИЯ к выполнению лабораторных работ по курсу

«Теория вероятностей и математическая статистика»

сиспользованием пакета Maxima для студентов направления 080100 «Экономика» профиля «Экономика предприятий

и организаций» очной формы обучения Часть 1

Составитель Дежин Виктор Владимирович

В авторской редакции

Подписано в печать Уч.-изд. л.

ФГБОУ ВПО «Воронежский государственный технический университет»

394026 Воронеж, Московский просп., 14

51

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]