. , , , , .., . . - , .
2 n ; 2 n 1. , ( ) ,
(2.19)
2 n 1. . (2.19), k
(2.20)
, n(x).. r() .
.
) r() = 1. = - 1, b = 1. . x k, g k,
xk = ½ (a + b) + ½ (b a)x k,
(2.20,a)
ck = ½ (b a)g k, 1 £ k £ n.
) [ - 1, + 1] r() = 1/(1 2)1/2. n (x).
x k = cos[p(k ½)/ n ], g k = p/ n, 1 £ k £ n. (2.20,)
) 0 £ < ¥; r() = a - , - ¥ < < + ¥ . 2.2.
Pn (x) r() [ a, b ], :
(2.21)
. , (a, b). Pn (x) Pn- 1(x).
, .
(2.22)
(2.23)
, ,
|
|
(2.24)
. 2.1.
n, .
. L 0(x) = 1; L 1(x) = 1; L 2(x) = ½ (3 x 2 1); L 3(x) = = ½ (5 x 3 3 x); L 4(x) = (1/8)(35 x 4 30 x 2 + 3); L 5(x) = (1/8)(63 x 5 70 x 3 + 15 x).
n = 1; x1 = 0; g1 = 2.
n = 2; - x1 = x2 = (1/3)1/2; g1 = g2 = 1.
n = 3; - x1 = x3 = (3/5)1/2; x2 = 0; g1 = g3 = 5/9; g2 = 8/9.
n = 4;
n = 5;
(a=0).
n = 1; x1 = 1; g1 = 1.
n = 2; x1 = 2 (2)1/2, x2 = 2 + (2)1/2; g1 = [2 + (2)1/2]/4; g2 = [2 (2)1/2]/4.
n = 3; x1 = 0,415775, x2 = 2,294280, x3 = 6,289945;
g1 = 0,711093, g2 = 0,278518, g3 = 0,010389.
. H 0(x) = 1; H 1(x) = 2 x; H 2(x) = 4 x 2 2; H 3(x) = 8 x 3 12 x; H 4(x) = 16 x 4 48 x 2 + 12; H 5(x) = 32 x 5 160 x 3 +120 x.
n = 1; x1 = 0; g1 = (p)1/2.
n = 2; - x1 = x2 = 1/(2)1/2; g1 = g2 = (p)1/2/2.
n = 3; - x1 = x3 = (3/2)1/2, x2 = 0; g1 = g3 = (p)1/2/6, g2 = 2(p)1/2/3.
2.1.
Ln (x) | Tn (x) | Un (x) | Hn (x) | |||
a, b | - 1, + 1 | - 1, + 1 | - 1, + 1 | - 1, + 1 | 0, + ¥ | - ¥, + ¥ |
r(x) | (1 x)a(1 + x)b, a, b > - 1 | (1 x 2)-1/2 | (1 x 2)1/2 | x a e - x, a > - 1 | ||
s(x) | 1 x 2 | 1 x 2 | 1 x 2 | 1 x 2 | x | |
Nn | 2a+b+1G(a + n + 1)G(b + n + 1)/[ n!(a + b + 2 n + 1)G(a + b + n + 1)] | 2/(2 n + 1) | p/2 n ¹ 0, p n = 0 | n! ´ ´G(a + n + 1) | (p)1/22 nn! | |
An | n (n +a + b + 1) | n (n + 1) | n 2 | n (n + 2) | n | 2 n |
Bn | (- 1) n /(2 n × n!) | (- 1) n /(2 n ´ ´ n!) | (- 1) n | (- 1) n | ||
an | 2(n + 1)(n + a + b + 1)(2 n + a + b) | n + 1 | ||||
bn | (2 n + a + b)(2 n + a + b + 1)(2 n + a + b + 2) | 2 n + 1 | ||||
cn | (b2 - a2)(2 n + a + b + 1) | 2 n + a + 1 | ||||
dn | 2(n + a)(n + b)(2 n + a + b + 2) | n | n (n + a) | 2 n |
n = 5;
. T 0(x) = 1; T 1(x) = x; T 2(x) = 2 x 2 1;
T 3(x) = 4 x 3 3 x; T 4(x) = 8 x 4 8 x 2 + 1; T 5(x) = 16 x 5 20 x 3 + 5 x;
T 6(x) = 32 x 6 48 x 4 + 18 x 2 1; T 7(x) = 64 x 7 11 x 5 + 56 x 3 7 x.
. U 0(x) = 1; U 1(x) = 2 x; U 2(x) = 4 x 2 1;
U 3(x) = 8 x 3 4 x; U 4(x) = 16 x 4 12 x 2 + 1; U 5(x) = 32 x 5 32 x 3 + 6 x.
. Pn (x) xi r i (1 £ x £ N),
,
l nPn +1(x) = (x an) Pn (x) + bnPn -1(x),
l n .
2.2.
, , .
. -; . .
|
|
[ a, b ] , p ; .
. q, q £ p. , q, . .
2.2.
, = 0. ( - 1 £ £ 0 0 £ £ 2 ), .
. , . [ a, b ] , . , ,
,
(2.25)
, .
.
, . , . Fn = F (xn).
. .
1 , . ,
x = a /(1 t) [ a, ¥) [0, 1]. , .
, "" . ,
.
2 ( ). b,
. ,
. b . ,
(. . 2.1), 4 £ £ 5 . . , .
. 2.1.
3 . , . ,
|
|
(, t ' = t x),
e t ' , . 2.1 , a = 0, ..
(2.26)
. , ,
, ~ 5 % =1, . ~ 5 % ³ 0,5, > 1 1 %.
, f (x) - [ a, b ] (.. , ). .
1 . f (x) = j(x) + y(x), j(x) , y(x) . , . , . , , .
2 . f (x) = j(x)r(x), j(x) , r(x) . r(x) . , . ,
(2.27)
(. . 2.1).
3 , . , (2.27) (xi -1, xi) exp(xi -1/2 )/(1 x 2)1/2, .
(2.28)
2.3. .
,
, 10 -4, , ,
x 2 = p x = p1/2 sin(x 2) 0. , x 2 > p . . -, x = p1/2 , .
-, x = p1/2. sin(x 2) = sin(p - x 2) = -sin(x 2 - p) ,
.. x = p1/2 1.748499467£ x £ 1.796089062. ,
4
f (y) ( h) y, [ a, b ].
|
|
, , .
: 1) , 2) , 3) [ f (y)], 4) .
k | a | b | h | |||
1,0 | 0,5 | 1,0 | 0,05 | 10-5 | ||
- " - | 1,2 | 0,5 | 1,0 | 0,05 | 10-5 | |
- " - | 1,4 | 0,5 | 1,0 | 0,05 | 10-5 | |
- " - | 1,6 | 0,5 | 1,0 | 0,05 | 10-5 | |
- " - | 1,0 | 1,0 | 1,5 | 0,05 | 10-5 | |
- " - | 1,2 | 1,0 | 1,5 | 0,05 | 10-5 | |
- " - | 1,4 | 1,0 | 1,5 | 0,05 | 10-5 | |
- " - | 1,6 | 1,0 | 1,5 | 0,05 | 10-5 | |
- " - | 1,0 | 1,5 | 2,0 | 0,05 | 10-5 | |
- " - | 1,2 | 1,5 | 2,0 | 0,05 | 10-5 | |
- " - | 1,4 | 1,5 | 2,0 | 0,05 | 10-5 | |
- " - | 1,6 | 1,5 | 2,0 | 0,05 | 10-5 | |
0,8 | 0.1 | 10-6 | ||||
-"- | 0,9 | 0.1 | 10-6 | |||
-"- | 1,0 | 0.1 | 10-6 | |||
-"- | 1,1 | 0.1 | 10-6 | |||
-"- | 0,8 | 0.1 | 10-6 | |||
-"- | 0,9 | 0.1 | 10-6 | |||
-"- | 1,0 | 0.1 | 10-6 | |||
-"- | 1,1 | 0.1 | 10-6 | |||
-"- | 0,8 | 0.1 | 10-6 | |||
-"- | 0,9 | 0.1 | 10-6 | |||
-"- | 1,0 | 0.1 | 10-6 | |||
-"- | 1,1 | 0.1 | 10-6 |
k | a | b | h | |||
0,8 | 0,1 | 10-6 | ||||
-"- | 0,85 | 0,1 | 10-6 | |||
-"- | 0,9 | 0,1 | 10-6 | |||
-"- | 0,95 | 0,1 | 10-6 | |||
-"- | 0,8 | 0,1 | 10-6 | |||
-"- | 0,85 | 0,1 | 10-6 | |||
-"- | 0,9 | 0,1 | 10-6 | |||
-"- | 0,95 | 0,1 | 10-6 | |||
-"- | 0,8 | 0,1 | 10-6 | |||
-"- | 0,85 | 0,1 | 10-6 | |||
-"- | 0,9 | 0,1 | 10-6 | |||
-"- | 0,95 | 0,1 | 10-6 | |||
0,6 | 1,0 | 1,3 | 0,03 | 10-6 | ||
-"- | 0,7 | 1,0 | 1,3 | 0,03 | 10-6 | |
-"- | 0,8 | 1,0 | 1,3 | 0,03 | 10-6 | |
-"- | 0,8 | 1,0 | 1,3 | 0,03 | 10-6 | |
-"- | 0,6 | 1,3 | 1,6 | 0,03 | 10-6 | |
-"- | 0,7 | 1,3 | 1,6 | 0,03 | 10-6 | |
-"- | 0,8 | 1,3 | 1,6 | 0,03 | 10-6 | |
-"- | 0,8 | 1,3 | 1,6 | 0,03 | 10-6 | |
-"- | 0,6 | 1,6 | 1,9 | 0,03 | 10-6 | |
-"- | 0,7 | 1,6 | 1,9 | 0,03 | 10-6 | |
-"- | 0,8 | 1,6 | 1,9 | 0,03 | 10-6 | |
-"- | 0,8 | 1,6 | 1,9 | 0,03 | 10-6 |
k | a | b | h | |||
0,4 | 2,5 | 0,05 | 10-6 | |||
-"- | 0,45 | 2,5 | 0,05 | 10-6 | ||
-"- | 0,5 | 2,5 | 0,05 | 10-6 | ||
-"- | 0,55 | 2,5 | 0,05 | 10-6 | ||
-"- | 0,4 | 3,0 | 3,5 | 0,05 | 10-6 | |
-"- | 0,45 | 3,0 | 3,5 | 0,05 | 10-6 | |
-"- | 0,5 | 3,0 | 3,5 | 0,05 | 10-6 | |
-"- | 0,55 | 3,0 | 3,5 | 0,05 | 10-6 | |
-"- | 0,4 | 3,5 | 4,0 | 0,05 | 10-6 | |
-"- | 0,45 | 3,5 | 4,0 | 0,05 | 10-6 | |
-"- | 0,5 | 3,5 | 4,0 | 0,05 | 10-6 | |
-"- | 0,55 | 3,5 | 4,0 | 0,05 | 10-6 | |
0,7 | 0,1 | 0,4 | 0,03 | 10-5 | ||
-"- | 0,9 | 0,1 | 0,4 | 0,03 | 10-5 | |
-"- | 1,1 | 0,1 | 0,4 | 0,03 | 10-5 | |
-"- | 1,3 | 0,1 | 0,4 | 0,03 | 10-5 | |
-"- | 0,7 | 0,4 | 0,7 | 0,03 | 10-5 | |
-"- | 0,9 | 0,4 | 0,7 | 0,03 | 10-5 | |
-"- | 1,1 | 0,4 | 0,7 | 0,03 | 10-5 | |
-"- | 1,3 | 0,4 | 0,7 | 0,03 | 10-5 | |
-"- | 0,7 | 0,7 | 1,0 | 0,03 | 10-5 | |
-"- | 0,9 | 0,7 | 1,0 | 0,03 | 10-5 | |
-"- | 1,1 | 0,7 | 1,0 | 0,03 | 10-5 | |
-"- | 1,3 | 0,7 | 1,0 | 0,03 | 10-5 |
k | a | b | h | |||
1,1 | 4,2 | 4,4 | 0,02 | 10-6 | ||
-"- | 1,2 | 4,2 | 4,4 | 0,02 | 10-6 | |
-"- | 1,3 | 4,2 | 4,4 | 0,02 | 10-6 | |
-"- | 1,4 | 4,2 | 4,4 | 0,02 | 10-6 | |
-"- | 1,1 | 4,4 | 4,6 | 0,02 | 10-6 | |
-"- | 1,2 | 4,4 | 4,6 | 0,02 | 10-6 | |
-"- | 1,3 | 4,4 | 4,6 | 0,02 | 10-6 | |
-"- | 1,4 | 4,4 | 4,6 | 0,02 | 10-6 | |
-"- | 1,1 | 4,6 | 4,8 | 0,02 | 10-6 | |
-"- | 1,2 | 4,6 | 4,8 | 0,02 | 10-6 | |
-"- | 1,3 | 4,6 | 4,8 | 0,02 | 10-6 | |
-"- | 1,4 | 4,6 | 4,8 | 0,02 | 10-6 | |
0,2 | 2,6 | 2,8 | 0,02 | 10-6 | ||
-"- | 0,25 | 2,6 | 2,8 | 0,02 | 10-6 | |
-"- | 0,3 | 2,6 | 2,8 | 0,02 | 10-6 | |
-"- | 0,35 | 2,6 | 2,8 | 0,02 | 10-6 | |
-"- | 0,2 | 2,8 | 3,0 | 0,02 | 10-6 | |
-"- | 0,25 | 2,8 | 3,0 | 0,02 | 10-6 | |
-"- | 0,3 | 2,8 | 3,0 | 0,02 | 10-6 | |
-"- | 0,35 | 2,8 | 3,0 | 0,02 | 10-6 | |
-"- | 0,2 | 3,0 | 3,2 | 0,02 | 10-6 | |
-"- | 0,25 | 3,0 | 3,2 | 0,02 | 10-6 | |
-"- | 0,3 | 3,0 | 3,2 | 0,02 | 10-6 | |
-"- | 0,35 | 3,0 | 3,2 | 0,02 | 10-6 |
|
|
2.4.
y Î [2.6, 3] h = 0.04.
(1 x)-1/2
n , , ,
, f (3) n = 10, 12 14 ( , y =3 ). (. .). , [cos(3 xk)] , .
n = 12 n = 14 , f (y) 12 . xky cos(xky). .
n | k | xk | cos(3 xk) | xk2 | xk2 + 0.3 | |
0.98769 | - 0.98410 | 0.97553 | 1.27553 | 0.24630 | ||
0.89100 | - 0.89220 | 0.79388 | 1.09388 | 0.28720 | ||
0.70711 | - 0.52314 | 0.50000 | 0.80000 | 0.39270 | ||
0.45399 | 0.20731 | 0.20611 | 0.50611 | 0.62073 | ||
0.15643 | 0.89189 | 0.02447 | 0.32447 | 0.96822 | ||
f 10(3) = 0.28818 | ||||||
0.99144 | - 0.98604 | 0.98295 | 1.28295 | 0.20406 | ||
0.92388 | - 0.93234 | 0.85355 | 1.15355 | 0.22695 | ||
0.79336 | - 0.72380 | 0.62942 | 0.92942 | 0.28168 | ||
0.60876 | - 0.25271 | 0.37059 | 0.67059 | 0.39040 | ||
0.38268 | 0.41027 | 0.14644 | 0.44644 | 0.58641 | ||
0.13052 | 0.92432 | 0.01704 | 0.31704 | 0.82576 | ||
f 12(3) = 0.28851 | ||||||
0.99371 | - 0.98715 | 0.98746 | 1.28746 | 0.17430 | ||
0.94388 | - 0.95235 | 0.89091 | 1.19091 | 0.18843 | ||
0.84672 | - 0.82453 | 0.71693 | 1.01693 | 0.22066 | ||
0.70711 | - 0.52314 | 0.50000 | 0.80000 | 0.28050 | ||
0.53203 | - 0.02529 | 0.28306 | 0.58306 | 0.38486 | ||
0.33028 | 0.54798 | 0.10908 | 0.40908 | 0.54854 | ||
0.11196 | 0.94412 | 0.01254 | 0.31254 | 0.71798 | ||
f 14(3) = 0.28854 |
y | xky | cos(xky) | y | xky | cos(xky) |
2.60 | 2.57774 | - 0.84520 | 2.64 | 2.61740 | - 0.86573 |
2.40209 | - 0.73880 | 2.43904 | - 0.76320 | ||
2.06274 | - 0.47234 | 2.09447 | - 0.50006 | ||
1.58278 | - 0.01198 | 1.60713 | - 0.03633 | ||
0.99497 | 0.54453 | 1.01028 | 0.53162 | ||
0.33935 | 0.94297 | 0.34457 | 0.94122 | ||
f (2.6) = 0.62012 | f (2.64) = 0.58406 |
.. .
yi | f(yi) | f(yi,yi+1) | f(yi,yi+1,yi+2) |
2.6 | 0.6201 | ||
-0.9 | |||
2.64 | 0.5841 | 0.15625 | |
-0.8875 | |||
2.68 | 0.5486 | 0.1875 | |
-0.8725 | |||
2.72 | 0.5137 | 0.21875 | |
-0.855 | |||
2.76 | 0.4795 | 0.1875 | |
-0.84 | |||
2.8 | 0.4459 | 0.1875 | |
-0.825 | |||
2.84 | 0.4129 | 0.25 | |
-0.805 | |||
2.88 | 0.3807 | 0.21875 | |
-0.7875 | |||
2.92 | 0.3492 | 0.25 | |
-0.7675 | |||
2.96 | 0.3185 | 0.21875 | |
-0.75 | |||
0.2885 |
3.
, , .. , . , , , () ( ).
, , , . , .
, , . , .
:
A x = b, A x, b ;
() -1;
det| A |;
.
, : , , ?
. , , .
-, (. . 3.1), . , , ( n n n - ). det A = 0, , , . .
( ). . .
det A 0. A x = b x = A -1 b. d E = d(AA -1) = A d A -1 + d AA -1=0,
d x = A -1(d b - d A × x).
, det A ¹ 0 . -1 , (d b d ), . . , det A 0 (. 3.1), , , .
c = || A ||×|| A -1||. , , c ³ 1. , : c ~ 103 104 .
, , . , , .
MATLAB :
k = cond(A) k = condest(A) k = cond(A, p) [k, v] = condest(A)
k = cond(A) , : k = smax/smin.
k = cond(A, p) - k º c = norm(A, p)×norm(inv(A), p), p = 1, 2, inf, 'fro'.
k = condest(A) 1-.
[k, v] = condest(A) , , norm(Av, 1) = norm(A, 1)×norm(v, 1)/k. , v - , k .
k_1 = = rcond(A), , 1-. , k_1 ; , k_1 .
3.1. :
() (b):
>> A = [7.9 5.6 5.7 -7.2; 8.5 -4.8 0.8 3.5; 4.3 4.2 -3.2 9.3; 3.2 -1.4 -8.9 3.3]
A =
7.900000 5.600000 5.700000 -7.200000
8.500000 - 4.800000 0.800000 3.500000
4.300000 4.200000 - 3.200000 9.300000
3.200000 - 1.400000 - 8.900000 3.300000
>> B = [5 7 6 5; 7 10 8 7; 6 8 10 9; 5 7 9 10]
B =
5 7 6 5
7 10 8 7
6 8 10 9
5 7 9 10
:
>> k = cond(A) k_1 = rcond(A) k1 = condest(A)
k = k_1 = k1 =
2.547891 0.417639 5.118441
>> k2 = cond(B) k2_1 = rcond(B) k20 = condest(B)
k2 = k2_1 = k20 =
2.9841e3 3.0882e- 4 4.4480e3
, , , .
Mathcad :
cond1(A) , 1-;
cond2(A) , L 2;
conde(A) , ;
condi(A) , .
Insert/Matrix,
:
cond1(A) = 5.118441, cond2(A) = 2.547891, conde(A) = 5.190408, condi(A) = 5.660286
cond1(B) = 4488, cond2(B) = 2984.092702, conde(B) = 3009,578708, condi(B) = 4488
MATLAB Mathcad, () cond(A) cond2(A), condest(A) cond1(A).
Maple V cond(A, namenorm), Mathcad, namenorm : 1, 2, infinity 'frobenius'. :
cond1(A) cond(A, 1); cond2(A) cond(A, 2); condi(A) cond(A, infinity);
conde(A) cond(A, 'frobenius')
Maple
> A:= matrix(4, 4, [[7.9,5.6,5.7,-7.2],[8.5,-4.8,0.8,3.5],[4.3,4.2,-3.2,9.3],[3.2,-1.4,8.9,
3.3]]);
Maple V Release 5
Maple 6 , .
. , . n £ 200 . () n 103 ¸ 105. .
(3.1)
(3.1,)
. (3.1) , , 1. , .. , .
, .. . , , .
. (k 1)- . :
(3.2)
k -
(3.3)
m - . ,
(3.4)
, k - . . .
, . (3.1). cmk; , .
(3.5)
(3.5) . :
(3.6)
. (3.2) (3.4) , (3.2) : , det A = 0. .
, cmk, . , . , .. k - , , . .
: