MATHCAD_ИНОСТР_НОВАЯ2
.pdf
ТЕМА 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
