Исследовать влияние числа разбиений k на вид решения yi (x) .
Таким образом, студенты, решающие 4-й вариант, должны написать программу по решению уравнения
|
|
y |
y |
x |
ln(yx) |
yx3 |
|
yx1/10 |
y , |
на участке |
x0 |
4 |
1 |
0,2 , |
xn |
4 |
1 |
2 |
2,2 |
для начального зна- |
|
|
|
|
|
|
|
25 |
|
25 |
|
|
|
|
|
|
|
|
|
|
чения y0 |
4 4 |
|
при k = 3·4 + 20 = 32. |
|
|
|
25 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В отчете необходимо представить постановку задачи, алгоритм, программу и результаты, содержащие табулированную функцию y = f(x) для трех различных шагов разбиения. Для выбора первого шага (k) воспользоваться вышеприведенной формулой. Два других шага выбрать самостоятельно. Построить графики y = f(x) для участка [x0, xn] для различных шагов разбиения.
5.9.Задачи по разработке технологии численного решения уравнений в частных производных в MS EXCEL
Постановка задачи. Разработать алгоритм и программу для решения уравнения
где T(t, x, y) – температура в момент t для координат области x, y; q – коэффициент температуропроводности.
Вид двухмерной плоскости задается на основе рис. 5.1.
В прямоугольной области среды располагаются две отливки прямоугольной формы из разных металлов. Размеры прямоугольников, области расположения отливок, начальные значения температуры клеток областей (отливок и среды), коэффициент температуропроводности выбираются по данным, представленным в табл. 5.9. В табл. 5.9 значения a и b – это высота и ширина области среды, a1
и b1 – это высота и ширина отливки 1, a2 и b2 – это высота и ширина отливки 2 (размерность значений – число клеток в MS EXCEL), q1, q2, q – коэффициенты температуропроводности соответственно отливки 1, отливки 2, области среды. Начальная температура отливки 1 вычисляется по формуле
T1 300 50 n ;
начальная температура отливки 2
T2 700 n2 ;
начальная температура области среды
|
|
|
|
|
|
|
T |
10 n , |
|
|
где n – номер варианта. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Среда |
b2 |
|
|
|
T |
|
b1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отливка 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отливка 1 |
|
|
|
|
|
T |
a |
|
|
|
|
|
T1 |
|
a2 |
|
T2 |
|
а1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 5.1
Рассчитать температурное поле охлаждения отливок для трех случаев:
–для заданных коэффициентов температуропроводности (q1, q2, q3);
–коэффициенты температуропроводности повышены на 50 %;
–коэффициенты температуропроводности повышены на 75 %.
Таблица 5.9
Номер варианта, n |
a1 |
b1 |
a2 |
b2 |
a |
b |
q1 |
q2 |
q |
1 |
8 |
25 |
16 |
10 |
20 |
50 |
30 |
45 |
70 |
2 |
10 |
30 |
18 |
8 |
21 |
49 |
31 |
44 |
69 |
3 |
10 |
40 |
17 |
6 |
22 |
48 |
32 |
43 |
68 |
4 |
12 |
32 |
20 |
12 |
23 |
47 |
33 |
42 |
67 |
5 |
14 |
33 |
19 |
10 |
24 |
46 |
34 |
41 |
66 |
6 |
10 |
35 |
21 |
14 |
25 |
45 |
35 |
40 |
65 |
7 |
15 |
37 |
23 |
15 |
26 |
44 |
36 |
39 |
64 |
8 |
8 |
38 |
24 |
13 |
27 |
43 |
37 |
38 |
63 |
9 |
12 |
39 |
26 |
12 |
28 |
42 |
38 |
37 |
62 |
10 |
10 |
25 |
25 |
11 |
29 |
40 |
39 |
3 |
61 |
11 |
11 |
26 |
28 |
10 |
30 |
40 |
40 |
34 |
60 |
12 |
7 |
27 |
15 |
8 |
20 |
41 |
41 |
35 |
59 |
13 |
10 |
28 |
16 |
7 |
21 |
42 |
42 |
33 |
58 |
14 |
12 |
29 |
19 |
10 |
22 |
43 |
43 |
32 |
57 |
15 |
13 |
30 |
21 |
9 |
23 |
44 |
44 |
31 |
56 |
16 |
14 |
25 |
20 |
8 |
24 |
45 |
45 |
30 |
55 |
17 |
12 |
26 |
23 |
6 |
25 |
46 |
46 |
50 |
54 |
18 |
8 |
27 |
17 |
6 |
26 |
47 |
47 |
49 |
53 |
19 |
9 |
28 |
25 |
10 |
27 |
48 |
48 |
48 |
52 |
20 |
10 |
22 |
26 |
13 |
28 |
49 |
49 |
47 |
51 |
21 |
14 |
19 |
25 |
10 |
29 |
50 |
50 |
46 |
50 |
22 |
12 |
25 |
27 |
8 |
30 |
49 |
51 |
45 |
70 |
23 |
7 |
15 |
13 |
6 |
20 |
48 |
52 |
44 |
69 |
24 |
6 |
30 |
14 |
5 |
21 |
47 |
53 |
43 |
68 |
25 |
10 |
22 |
16 |
7 |
22 |
46 |
54 |
42 |
67 |
Так, например, для варианта 2 (n = 2) условия задачи будут задаваться следующим образом.
Шаг 1. Определение размеров прямоугольников и коэффициентов температуропроводности из табл. 5.9:
a1 = 10 клеток, b1 = 30 клеток, a2 = 18 клеток, b2 = 8 клеток, a = 21 клетка, b = 49 клеток,
q1 = 31, q2 = 44, q = 69.
Шаг 2. Определение начальных температур отливок и области среды по формулам
T1 |
300 50 2 |
400 |
С , |
T |
700 n2 700 |
22 |
704 С , |
2 |
|
|
|
T 10 2 12 C.
Шаг 3. Построение в MS EXCEL клеточной структуры среды и отливки. Для того чтобы среда отображалась во весь экран без прокрутки, необходимо уменьшить ширину столбцов, для этого нужно выделить столбцы, которые необходимо уменьшить, выбрать Формат – Столбец – Ширина и в появившемся диалоговом окне задать ширину столбцов, например, 2. Для корректного отображения температуры в клетках необходимо уменьшить шрифт. Вид полученной клеточной структуры представлен на рис. 5.2.
Рис. 5.2
5.10. Задачи по разработке технологии решения системы линейных уравнений в MS EXCEL
Постановка задачи. Разработать алгоритм решения системы линейных уравнений с пятью неизвестными в MS EXCEL, если коэффициенты системы вычисляются по закону, заданному в табл. 5.10. Например, для 1-го варианта коэффициенты вычисляются как
аij= sin(i2 2ij j) , ci tg3 (i 10) ,
где i = 1; 5 (1), j = 1; 5 (1), (1 – шаг изменения переменной).
Система линейных уравнений после вычисления коэффициентов должна иметь вид
a112 x1 a12x2 a13x3 a14x4 a15x5 c1, a21x1 a222 x2 a23x3 a24x4 a25x5 c2, a31x1 a32x2 a332 x3 a34x4 a35x5 c3, a41x1 a42x2 a43x3 a442 x4 a45x5 c4,
a51x1 a52 x2 a53x3 a54 x4 a552 x5 c5.
Решить сформированную систему линейных уравнений, использовав встроенные функции MS EXCEL, предварительно определив коэффициенты при неизвестных j по формулам, представленным в табл. 5.10 по вариантам заданий.
265
5.11.Задачи по разработке технологии решения нелинейных уравнений в MS EXCEL
Постановка задачи. Разработать алгоритм решения нелинейных уравнений вида f(x) = 0, используя метод дихотомии на отрезке [а, b]. Определить число корней на указанном отрезке, используя графические средства MS EXCEL. Найти три любых корня уравнения с точностью ε = 0,01. Варианты заданий представлены в табл. 5.11.
Таблица 5.11
Номер варианта |
|
Функция f(x) |
|
|
a |
b |
1 |
|
|
|
2 |
|
|
|
3 |
4 |
1 |
|
sin 4 x |
|
4 |
|
|
-10 |
10 |
|
|
|
|
|
|
|
|
|
|
|
2 |
5cos x |
4 |
sin 7x |
-5 |
15 |
|
|
|
|
|
|
|
|
|
3 |
4 cos |
1 |
|
|
|
cos 4x |
-11 |
21 |
|
|
|
|
|
|
x |
1 |
|
|
|
|
|
|
|
|
4 |
sin cos2 |
x 4 |
|
|
-14 |
6 |
5 |
cos 3x |
|
2sin 5x |
cos 4x |
-5 |
25 |
|
|
|
|
|
|
|
|
6 |
tg 4x |
|
|
ctg 7x |
|
-16 |
8 |
|
|
|
|
|
|
|
|
7 |
4cos 7x |
|
x2 |
2x |
|
-13 |
40 |
8 |
5sin 11x |
2sin 2 12x |
3 |
-8 |
30 |
9 |
log4 x |
12 |
|
|
|
log5 |
x |
13 |
-17 |
8 |
10 |
cos3 3x |
|
|
2cos2 3x |
|
-20 |
14 |
11 |
sin(3x) sin(6x) |
|
-12 |
18 |
|
|
|
|
|
12 |
tg(5x) tg(7x) |
|
-19 |
5 |
|
|
|
|
|
13 |
ctg(x) ctg(2x) |
|
-4 |
24 |
|
|
|
|
|
|
|
|
14 |
cos3 3x |
|
|
2cos2 3x |
|
-5 |
17 |
15 |
sin 4 3x |
|
|
4sin 2 3x |
|
-10 |
10 |
16 |
tg2 3x |
|
|
|
tg 4x |
|
-5 |
15 |