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

Практикум

.pdf
Скачиваний:
68
Добавлен:
25.04.2014
Размер:
726 Кб
Скачать

) !(

( 0 3 3 ( "

3 0 4 ( ( 0 3

3 % 3 (

( ( ( 4 ( f " 0 f ( " 0 ( ( 0

( ( ( 0 >0 ˆ f

0 ( / 0 0 3/ 3 0 3 3 0

3 % ( 0 ( 0 Z

f ( 3 0 3/ 3 0 0 0

0 % 0 3/ % %

f 5 5Z 0 0 3/

0 ( 00 ( 5Z / 0 %

g

% % % 0 / g

3 0 0 ( 0 4 ( (

( ( g

3 ( ( % 0 5 0 ( g

0 0 3 ( 0 ( " 0 3 3 3 (

00 ( / / ( ( 0 0

4 3 / 3 0 0 3/ f 5Z 3 / 3

," -5 .#$A,A$"% 0 x y" Z (

R\

/ 3" x1" 5Z g 3 % 3

f ( y1 0 3 3 0 0 3/ 3

(

' !O"! , @= #$% $" !O! "O"" ," -5 .#$A,A$"% -5=+#$A,A; ;A$"A,"A; ;%

0 0 3/ 3 0 % 3 / " ( 5Z %

% 0 T % 4 0 5 3 3 ( % /

. 4 4 0 3/ 00 ( / /

% % 0 3 % / T " 0

" 5 3 0 % / c %

%d 0 / 00 ( / 0 3/ V 3 0 % 00 ( / " 0 ( 0 ( "

( " 0 3 00 ( /

.. ! (

! Š W ( 0 5 / 5 -=5,: +" 3

0 ( 0 00 ( 5Z / ( (

4 (( / 3 Z f /

0 ( b 0 ( 0 0 ( Z f / 3

p( 0 3 /

!0 ( QQ ( JJ ( U 3

0 cLD8:iBd x1" x2" xn ( / f / 0 ( c>?8EBd aij " i = 1, . . . n" j = 0, . . . , N c N ˆ 0 0 ( d

/ 5

 

 

a10(x − x1)N

+ a11(x − x1)N −1 + . . .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ a1 N

1

(x

x1) + a1N , x < x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

+ a (x x )

N

1

+ . . .

 

 

 

 

 

 

a (x x )

 

 

 

 

 

 

 

 

 

 

 

20

 

2

 

 

21

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ a2 N 1(x x2) + a2N , x2 < x < x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (x) =

 

 

 

 

N

 

 

 

· · ·

 

N

 

 

1

 

 

 

 

 

 

· · ·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xk )

 

+ . . .

 

 

 

 

 

 

ak0(x xk ) + ak1(x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1(x

 

xk ) + akN , xk < x < xk+1

 

 

 

 

 

 

 

 

 

+ ak N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

· · ·

 

 

 

 

 

 

 

 

 

N

1

 

· · ·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an 1 0(x xn 1) + an 1 1(x xn 1) + . . .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ an 1 N 1(x xn 1) + an 1 N , xn 1 < x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ 3( " ( QQ ( ( 3

V 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

--b 5#lA$%

3 /

l

 

( $ c0

-=5,b 5 3 0 ( d

: >.*#l% ( 3 0 5 / l c 3

JJ ( d

: +#l% 3 0 5 / l" Z 5Z 5 3

x0

: @17 0 3 ( 0 / / 3( "

0 (

: @17#lA @%" @ ˆ 3 " ( / 5 l ( @ > 8 % 3% 0 f / 3

: @17#lA •% 3 (( % / l •

!(

V 3 0 3 0 3/ 0 S31 % ( 0

3/ 5 / {xj , yj } % 3 ! Z

/ " 3 0 . ( ( 4 ˆ

@ -.

/ 3 @ -. 0 3/ 0 S31 0

0 3/ / R= {xj }" [= {yj } (

4 f ( (3 ( (

S

j @ -.#RA[%

! Z ( / @ -. 3 3 3 0 ( S31" 0

5Z ( / {xj , yj } 35Z ( ( ( 3("

0 3 0 ( 5 0 j 3 QQ (

0 ( / 5 sin 0 ( [0, π] 0 03

3Z ( (

b ( R" Z

R 5 - @.#!A- A(%

! ( 3 / sin f % % 0 ( (

[

[ #R%

) 0 0 ( 0

j @ -.#RA[%

( 0 j / sin 0 ( [−π/2, 3π/2]

$ E - O! !"O - F -5=+#$A--b 5#jA$%A$A #$%ARA[A; ;%

!

b 0 % 0 3( /

T 3 0 / 0 ( 0 3/

! 0 ( 0 3/ 0 3 3 0

( /

! +#" 0 : ! :

S1

3 !

! +#%

!

3 W 0 3 0 3 0 5 c 0 3 5d / 5

0 3 5

4 0 5 " 0 j 3

( % %

!/ @ -. 0 ( 5 % 5 3

3 0 0 0 % x0 xN (

3 0 0 0 3/ % 0 ( 0 "

0 % 0 R ( ( ( 0 ( (

! +#& ! P3

x0 x1 x2 x3

dP3(x0) = dS(x0) . dx dx

- 0 0 " 35Z ( ( 3(

( 0 3 / @ -. 0 ( 5" / 3 @ -.

3 W R ( ( ! ( (

; =+ H =+;" ;-.* => @;" ;@=1-5.+.;" ; .@= >;" ;b * + = 5;

! +#+ 5 Z.5- @ -.

: "5

) ' '

% ( 0 0 3/ 5 4 /

( 3 0 0 ( 0

( 0 Fj ( 0 0 5

/ δij " " 0 F0 3 3( 0 3/

F0(x0) = 1, F0(x1) = 0, . . . , F0(xN ) = 0.

W

( 0 Fk 3 3( 0 3/

Fk (x0) = 0, . . . , Fk (xk−1) = 0, Fk (xk ) = 1, Fk (xk+1) = 0,

. . . , F0(xN ) = 0.

- ( 0 0 " 0 3/ 0

/ f (x) 3 5 ((

 

N

 

 

j

 

S(x) =

f (xj )Fj (x).

cR Sd

 

=0

 

T 0 ( % 0 / / 3( " ( QQ

( ( 0 S

! +#. 3 : ! " "5 "

7-

! +#/ 3 : ! " "5 "

! *7-+ & : @17

! GG

*7-+ !

7-

! +#2 ; ! " !

f (x) = cos(5x) exp(0.3x2)

 

f (x) = sin(5x) exp(0.3x2)

H6 CI !

! " "5

" ) !

J

!( Z / 00 ( /

/ / 0 ( 00 ( / 3 3 (

R

( / " 0 f 0 3 3 0 3( JD?Q8D@G kHA@?D V 3 f 0 / 6??=B

( 5 3 (

t:BA> IA@@A;M 00 ( / c dg v:@: K@:@AB@A>B 0 ( %

( t:BA> IA@@A;M " 5Z 0 3 00 ( / 0 " f ( 0 ( 0

3( S c 3 00 ( / 3d S ! ( 0 5

c 3 00 ( / 3d W c 3 00 ( / 3d

( v:@: K@:@AB@A>B ( 0 4

%

8 ( 0 0 5Z 3

t" 0 o

 

 

R

[

Y

S

U " (!

[X

R

SXS

[ \

S X

S SX

 

 

 

 

 

 

 

( 0

E!A !A&!A3!AG!A"!!F

a E ! 34A ! (&A ! "4"A ! 3!9A " !(4A " ("4F -5=+# Aa ;=;%

' 0 ( (" 0 33 ( % (

c 3 % ( 0 3( %d

0 ( 6??=B { t:BA> IA@@A;M" ( 0 0

3/ ! f ( ( ( % 0

3/ ( 0 3/ " 5Z 3

3 0 5Z / 00 ( / 0 3 3

(

2 0 0 ( K7?F 8‹<:@A?;B c 3d" (0 ( 0 0 3/ 0 0 ( /

' /" 0 ( 6??=B { v:@: K@:@AB@A>B" ( 0

3 ( % 0 ( %" 0 % ( | m T ( 0 0 ( f ( " ( 5

5Z 0 3 " 0 ( ( ( (" ( ( ( % ( ( ( ("

( ( ( ( ( %

! +#6 25 0?

! /

f 2K!

-6

46

C6

L6

-C6

1C6

v M

--71

--C4

--L-

--NL

-1-O

-11N

; ! ! )

! +#; 5 /

t 0?

6

16

76

P6

N6

-66

p

P-6N*1+

1 44L*4+

L 4L7*4+

-OO-*7+

7 L4C*7+

-6-4*C+

; ! ! )

[

/ 0

b % 3 4 ( (

( 0 ( V 0 ( r f % f (x) = 0 % 0 x = F (x) !

( ( ( ( 4 3 %

( ( ' % ( ( " 3 0 ( 3

( ( 0 3 0 % 3

( " f %

,

1 !

 

 

 

 

 

 

' ( 3

 

 

 

 

 

 

 

 

 

 

2

' ( " ( ( S

 

0 ( 3 0 x = 1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 0

0 3

) x2

4 " f x 4 ( !

0 3 ( ( ( 4 x = 1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 " 03 0

0 0 (

) x2

( 4 " f x 4 ( (

 

3 0 0 ( ( 0 ( 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

;

1

1

;

 

 

1

3

;

1

5

;

1

13

;

1

27

;

. . .

 

 

 

 

 

 

 

32

 

 

2

 

4

 

 

 

8

 

16

 

 

 

64

 

 

V ( ( 0 ( Z 5 5Z

0 ((

" 7 H !

TZ 5. #7 V .- % X #HU"!!% $ # 7%

X

: $P V 7$

.5 .$

. >

H H "

. >

$ 3 0 0 ( % 3 0 ( % 3%

0 4 ) ( ( 0 (( 0 0 3

3 Z

0 0 3 0 (( ' 0 4 3 0 ( " 3

3 3 f (x) = 0 Z 0 [a, b]"

3 (

! .#" ; " :

C1 ! "5

*

:=*1 + Z.$+

::3 !9.334: 7@@ 7 ::3 !9.334: 7@>

( )

V 3 % 3 3 f (x) = 0 0 (( %

( ( (

:\ @+ = 7 .@+#:A !A7!A.-A1$% u ' x O

u ;:@ ; Bx y^ i{ d y^ i{ y d ;:@ 1; u !A7! x i ^x i x ' xB Bi i e u .- x M Bxe

u 1$ iB de ^Bx M Bd x { u ' M O

u H M Bd d ' x {

Y