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

MATHCAD_ИНОСТР_НОВАЯ2

.pdf
Скачиваний:
57
Добавлен:
02.02.2015
Размер:
1 Мб
Скачать
☆

ТЕМА 8 Нахождение корней уравнения f (x) 0

8.1. Краткие теоретические сведения

Для того, чтобы найти решение уравнения с одним неизвестным графическим способом, нужно, перенеся все члены его в левую часть, представить это уравнение в виде f (x) 0. После этого необходимо построить график функции y f (x). Абсциссы точек пересечения этого графика с осью X равны корням исходного уравнения. Если таких точек нет, то уравнение не имеет решений.

8.2 Задачи для выполнения лабораторных работ

Листинг 8.1. Найти корни уравнения sin(x) cos(x) 0, в интервале

x

Найти корни уравнения

y(x) : sin(x) cos(x)

x : 3.14, 2.9 .. 3.14

Исходное уравнение

Задание диапазона изменения переменной

Построение графика для определения количества корней уравнения

2

1

y(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

2

0

2

4

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

x

Задание начального приближения и нахождение корней уравнения

x :

1

root(y(x),x)

 

0.785

x :

2.2

root(y(x),x)

 

2.356

51

Пояснение

Функция root используется для решения одного уравнения, например, f(x)=0

root (выражение, переменная)

 

 

левая часть

имя переменной, относительно

yравнения f(x)=0

которой решается уравнение

Варианты задач представлены в табл.8.1.

Таблица 8.1. Варианты задач

№ вар.

Уравнение

Диапазон изменения переменной x

1

x3 4x2 3x 7 0

5 x 5

2

x3 2x2 6x 3 0

5 x 5

3

x sin(x) 1 0

 

 

x

 

5

 

5

 

 

4

8cos(x) x 6

 

 

x

 

5

 

5

 

 

5

2 ln(x 7) 5 sin(x) 0

 

 

x

 

5

 

5

 

 

6

4 cos(x) 0,3x 0

 

 

x

 

5

 

6

 

 

7

2x2 5 2x

5 x 5

8

5x 45x2 170 0

5

x 5

9

5x 45x 70 0

5 x 5

10

3x 1 4 x 0

5 x 5

11

2x 3x 2 0

5 x 5

12

4x 1 9x 0

2 x 5

52

Нахождение корней полинома

 

Для нахождения

корней полинома

(выражения, имеющего вид

vnxn ...v2x2 v1x v0 )

предпочтительнее

использовать

функцию

polyroots(v), которая в отличие от функции root не требует начального приближения. Эта функция возвращает сразу все корни, как вещественные, так и комплексные.

Аргумент v представляет собой вектор с перечнем коэффициентов полинома, начиная с младшего. При записи нельзя игнорировать нулевые коэффициенты (т.е. если полином третьей степени и отсутствует член со второй степенью x, то необходимо указать 0).

Листинг 8.2. Найти корни полинома 2x3 11x 5 0

Листинг 8.3. Найти корни полинома f (x) x4 6x3 12x2 10x 3

В данном примере имеются комплексные корни.

53

ТЕМА 9 Нахождение экстремумов функции

9.1. Краткие теоретические сведения

Исследование функции с помощью производной.

Теорема 1 (Необходимое условие экстремума функции). В точке экстремума производная обращается в нуль.

Теорема 2 (Достаточное условие экстремума функции).

Если в некоторой точке производная обращается в нуль и, кроме того, производная, проходя через неё, меняет свой знак, то в этой точке функция достигает экстремума.

Точки экстремума надо искать прежде всего среди корней производной!

Листинг 9.1. Найти экстремумы функции (min, max)

T(x)

3sin(x)3

2,3

на отрезке 3 x 7

h 0,1

x 2cos(x)3

 

 

 

 

54

 

Найти экстремумы функции

 

3 sin(x)3

2.3

 

 

 

 

 

T(x) :

 

 

Исходное уравнение

 

x 2 cos(x)3

 

 

 

 

 

Определение производной

 

 

 

 

T1(x) : d

T(x)

 

 

 

 

 

 

dx

 

 

 

 

 

 

9 cos(x) sin(x)2

6 sin(x) cos(x)2

1 3 sin(x)3 2.3

T1(x)

x 2 cos(x)3

 

 

 

x 2 cos(x)3 2

x : 3,3.1.. 7

Задание диапазона изменения переменной

Построение графика для определения

 

размещения и количества экстремумов

 

 

1

 

 

 

 

 

 

 

0.5

 

 

 

 

 

 

T(x)

 

 

 

 

 

 

 

T1(x)

3

4

5

6

7

8

 

0.5

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

x

 

 

 

 

Задание начального приближения и вычисление значения экстремума

x : 4.6

x1 : root(T1(x),x)

 

x1 4.696

T(x1) 0.149

x : 6.2

x1 : root(T1(x),x)

 

x1 6.157

T(x1) 0.546

x : 6.7

x1 : root(T1(x),x)

 

x1 6.727

T(x1) 0.483

 

 

 

 

 

55

 

 

9.2. Задачи для выполнения лабораторных работ

Найти экстремумы функции (min, max)

Таблица 9.1. Варианты задач

№ вар.

Функция

 

 

 

Диапазон изменения

 

 

 

 

 

 

 

 

 

переменной x

1

Y(x) lg(x)3 sin(x)2 lg(x)

1 x 10

h 0,1

2

B(x) lg(x)cos(x)2

sin(x)

1 x 10

h 0,1

3

Y(x) ln(x)5 sin(x)2 x 20

1 x 10

h 0,1

4

Z(x) sin(x)3 0,3

xcos(x) 5

1 x 8

h 0,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

T(x) sin(x)

3

0,6

x

2

1 x 10

h 0,1

 

 

 

 

 

6

F(x) 5cos(x) 5sin(x)2

3

0 x 10

h 0,1

7

C(x) (sin(x)3 0,7)

x

1 x 10

h 0,1

8

S(x) (cos(x)3 0,2)

x

1 x 10

h 0,1

9

Z(x) (cos(x)3 0,5x)

(xsin(x) 7)

1 x 7

h 0,1

10

T(x)

3sin(x)

3

2,3

 

3 x 9

h 0,1

 

 

 

 

 

 

x 2cos(x)3

 

 

 

 

 

 

 

 

 

 

11

F(x) (sin(x)3 0,5)

(xsin(x) 7)

1 x 10

h 0,1

12

G(x) (cos(x)3 0,5)

(x cos(x))

1 x 10

h 0,1

2-ой способ нахождения локального экстремума

Для поиска локальных экстремумов имеются две встроенные функции:

Minimize( f ,x1,...xM )– вектор значений аргументов, при которых функция f достигает минимума;

Maximize( f ,x1,...xM )– вектор значений аргументов, при которых функция f достигает максимума;

f – функция;

x1,...xM – аргументы, по которым производится

минимизация (максимизация).

56

В качестве примера рассмотрим задачу численного поиска экстремумов полинома четвертой степени

f (x) x4 3x3 2x2 x 0.5 в интервале 0.7 x 1.8

57

ТЕМА 10 Линейная и полиномиальная аппроксимация

10.1. Краткие теоретические сведения

ЗАДАЧА. По данным експеримента (или табличным данным) определить функциональную зависимость. Выбрать вид кривой и определить коэффициент функциональной зависимости, которая описывает экспериментальные данные с ошибкой ±5%.

Линейная регрессия (1 способ)

Приближение данных (xi,yi ) осуществляется линейной функцией y(x) b0 b1 x. На координатной плоскости (x,y) линейная функция, как известно, представляется прямой линией. Еще линейную регрессию часто называют методом наименьших квадратов.

В соответствии с методом наименьших квадратов коэффициенты b0 и b1 находят из условия минимизации суммы квадратов отклонений между экспериментальными и вычисленными по расчетному уравнению значениями, т.е.

N

 

2

min

или, что равносильно:

 

F (yi yi)

 

i1

N

F(yi (b0 b1xi ))2 min,

i1

где xi – экспериментальные значения фактора (входного параметра); yi – экспериментальные значения выходного параметра;

yi – вычисляемые по расчетному уравнению yi b0 b1xi значения выходного параметра;

N– размер выборки (количество экспериментальных точек)

N N N N N

N (yixi ) xi yi

b

i 1

i 1

i 1

1

N

 

N

2

 

N xi2

 

xi

 

i 1

i 1

 

b

yi

b

xi

i 1

i 1

N

N

0

1

ЗАДАЧА. Определить функциональную зависимость вязкости

от

температуры ºС для воздуха

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Температура,

0

100

200

300

400

500

600

700

800

 

(°C)

 

 

 

 

 

 

 

 

 

 

 

вязкость

0,0171

0,0218

0,0259

0,0294

0,0320

0,0357

0,0384

0,0411

0,0437

 

58

Линейная аппроксимация - способ 1

 

ORIGIN : 1

 

 

N : 9

 

 

 

 

0

 

 

 

 

 

0.0171

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

0.0218

 

 

 

 

 

 

 

 

 

 

 

 

200

 

 

 

 

 

0.0259

 

 

 

300

 

 

 

 

 

0.0294

 

X :

400

 

 

 

Y :

0.0320

 

 

 

 

 

 

 

 

 

 

 

 

 

500

 

 

 

 

 

0.0357

 

 

600

 

 

 

 

 

0.0384

 

 

 

 

 

 

 

 

 

 

 

 

 

700

 

 

 

 

 

0.0411

 

 

 

 

 

 

 

 

 

 

 

 

800

 

 

 

 

 

0.0437

 

 

N

 

 

 

 

N

 

N

N

syx :

Yi Xi

 

sx :

Xi

sy : Yi

sx2 : Xi 2

i 1

 

 

 

 

i 1

 

i 1

i 1

b1 :

N syx sx sy

 

b1 3.26 10 5

 

 

 

 

 

 

 

 

 

N sx2 sx2

 

 

 

 

 

 

 

 

b0 :

sy

b1

sx

 

b0 0.019

 

 

 

 

 

 

 

 

N

N

 

 

 

 

 

 

 

 

 

YR(X) : b0 b1 X

corr(Y,YR(X)) 0.995

 

 

 

0.019

 

 

 

 

 

 

 

 

 

 

 

 

0.022

0.0250.028

YR(X) 0.0320.035

0.0380.041

0.045

0.05

 

 

 

 

0.04

 

 

 

 

Y

 

 

 

 

YR(X) 0.03

 

 

 

 

0.02

 

 

 

 

0.01

200

400

600

800

0

 

 

X

 

 

59

Линейная регрессия (2 способ)

c использованием функций intercept(x, y), slope(x, y).

intercept(x, y) – коэффициент b0линейной регрессии;

slope(x, y) – коэффициент b1 линейной регрессии:

x– вектор действительных данных аргумента;

y– вектор действительных данных значений того же размера.

60

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