Компьютерная обработка информации (в пакете MathCAD). Мартьянова А.Е. - 107 стр.

UptoLike

ВУЗ: 

АГТУ | Астрахань

Составители: 

Рубрика: 

Мартьянова А.Е. Компьютерная обработка информации
107
5.6.2. Решение жесткой системы ОДУ I с помощью
функций stiffb и stiffr
acc 0.000000001:= kmax 5:= s10:=
ZZ stiffb y t1, t2, acc, D, J, kmax, s,():=
ZZ1 stiffr y t1, t2, acc, D, J, kmax, s,():=
ZZ
0
12.718
36.521
92.637
100
1
0.295
0.033
2.01 10
4−
×
1.029 10
4−
×
0
2.757 10
5−
×
3.019 10
6−
×
1.827 10
8−
×
9.355 10
9−
×
0
0.705
0.967
1
1
⎛
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎟
⎠
=
ZZ1
0
10.107
20.295
30.327
100
1
0.377
0.146
0.058
1.029 10
4−
×
0
3.549 10
5−
×
1.347 10
5−
×
5.326 10
6−
×
9.355 10
9−
×
0
0.623
0.854
0.942
1
⎛
⎜
⎜
⎜
⎜
⎜
⎜
⎜
⎝
⎞
⎟
⎟
⎟
⎟
⎟
⎠
=
Z
0
〈〉
()
100
100=
Z
1
〈〉
()
100
1.029 10
4−
×=
Z1
0
〈〉
()
100
100=
Z1
1
〈〉
()
100
1.029 10
4−
×=
Z
2
〈〉
()
100
9.355 10
9−
×=
Z
3
〈〉
()
100
1=
Z1
2
〈〉
()
100
9.356 10
9−
×=
Z1
3
〈〉
()
100
1=
Мартьянова А.Е. Компьютерная обработка информации                                             107



    5.6.2. Решение жесткой системы ОДУ I с помощью
    функций stiffb и stiffr
   acc := 0.000000001            kmax := 5         s := 10
   ZZ := stiffb ( y , t1 , t2 , acc , D , J , kmax , s)
   ZZ1 := stiffr ( y , t1 , t2 , acc , D , J , kmax , s)

        ⎛ 0               1              0          0 ⎞
        ⎜
                                             −5
        ⎜ 12.718       0.295       2.757 × 10     0.705 ⎟
        ⎜                                    −6
                                                        ⎟
   ZZ = ⎜ 36.521       0.033       3.019 × 10     0.967 ⎟
        ⎜                     −4             −8         ⎟
        ⎜   92.637  2.01 ×  10     1.827 × 10       1   ⎟
        ⎜                      −4            −9
        ⎝ 100 1.029 × 10           9.355 × 10       1 ⎠
          ⎛ 0              1              0          0 ⎞
          ⎜
                                               −5
          ⎜ 10.107       0.377      3.549 × 10     0.623 ⎟
          ⎜                                    −5
                                                          ⎟
   ZZ1 = ⎜ 20.295        0.146      1.347 × 10     0.854 ⎟
          ⎜                                    −6         ⎟
          ⎜  30.327      0.058      5.326 × 10     0.942  ⎟
          ⎜                     −4             −9
          ⎝ 100 1.029 × 10          9.355 × 10       1 ⎠

    ( Z 〈0〉 ) 100 = 100                                      ( Z 〈1〉 ) 100 = 1.029 × 10− 4
    ( Z1 〈0〉 ) 100 = 100                                     ( Z1 〈1〉 ) 100 = 1.029 × 10− 4
    ( Z 〈2〉 ) 100 = 9.355 × 10− 9                            ( Z 〈3〉 ) 100 = 1
    ( Z1 〈2〉 ) 100 = 9.356 × 10− 9                           ( Z1 〈3〉 ) 100 = 1

Страницы