Системы семи дифференциальных уравнений
Поможем в ✍️ написании учебной работы
Поможем с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой

 

Sub DU()

x1=1 'начальные условия при t=0

x2=0 'начальные условия при t=0

x3=0 'начальные условия при t=0

x4=0 'начальные условия при t=0

x5=0 ' начальные условия при t=0

x6=0 ' начальные условия при t=0

x7=0 ' начальные условия при t=0

Sheets("1").Cells(k+2;2).Value=x1

Sheets("1").Cells(k+2;3).Value=x2

Sheets("1").Cells(k+2;4).Value=x3

Sheets("1").Cells(k+2;5).Value=x4

Sheets("1").Cells(k+2;6).Value=x5

Sheets("1").Cells(k+2;7).Value=x6

Sheets("1").Cells(k+2;8).Value=x7

dt=30/50

a12=Sheets("1").Cells(5;9).Value ' инт. потока

a13=Sheets("1").Cells(5;10).Value ' инт. потока

a21=Sheets("1").Cells(5;11).Value ' инт. потока

a23=Sheets("1").Cells(5;12).Value ' инт. потока

a34=Sheets("1").Cells(5;13).Value ' инт. потока

a45=Sheets("1").Cells(5;14).Value ' инт. потока

a52=Sheets("1").Cells(5;15).Value ' инт. потока

a26=Sheets("1").Cells(5;16).Value ' инт. потока

a62=Sheets("1").Cells(5;17).Value ' инт. потока

a67=Sheets("1").Cells(5;18).Value ' инт. потока

a72=Sheets("1").Cells(5;19).Value ' инт. потока

For k = 0 To 50

k1=One(x1;x2;x3;x4;x5;x6;x7;a12;a13;a21)*dt

m1=Two(x1;x2;x3;x4;x5;x6;x7;a12;a26;a21;a23;a52;a62;a72)*dt

n1=Three(x1;x2;x3;x4;x5;x6;x7;a13;a23;a34)*dt

o1=Four(x1;x2;x3;x4;x5;x6;x7;a34;a45)*dt

p1=Five(x1;x2;x3;x4;x5;x6;x7;a45;a52)*dt

r1=Six(x1;x2;x3;x4;x5;x6;x7;a26;a67;a62)*dt

s1=Seven(x1;x2;x3;x4;x5;x6;x7;a67;a72)*dt

k2=One(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1; x6+0,5*r1;x7+0,5*s1;a12;a13;a21)*dt

m2=Two(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1; x6+0,5*r1;x7+0,5*s1;a12;a26;a21;a23;a52;a62;a72)*dt

n2=Three(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7+0,5*s1;a13;a23;a34)*dt

o1=Four(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7+0,5*s1;a34;a45)*dt

p1=Five(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7+0,5*s1;a45;a52)*dt

r1=Six(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7+0,5*s1;a26;a67;a62)*dt

s1=Seven(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7+0,5*s1;a67;a72)*dt

k3=One(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s2;a12;a13;a21)*dt

m3=Two(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s2;a12;a26;a21;a23;a52;a62;a72)*dt

n3=Three(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s2;a13;a23;a34)*dt

o3=Four(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s2;a34;a45)*dt

p3=Five(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s2;a45;a52)*dt

r3=Six(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2; x6+0,5*r2;x7+0,5*s2;a26;a67;a62)*dt

s3=Seven(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s2;a67;a72)*dt

k4=One(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a12;a13;a21)*dt

m4=Two(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a12;a26;a21;a23;a52;a62;a72)*dt

n4=Three(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a13;a23;a34)*dt

o4=Four(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a34;a45)*dt

p4=Five(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a45;a52)*dt

r4=Six(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a26;a67;a62)*dt

s4=Seven(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a67;a72)*dt

x1=x1+(k1+2*k2+2*k3+k4)/6

x2=x2+(m1+2*m2+2*m3+m4)/6

x3=x3+(n1+2*n2+2*n3+n4)/6

x4=x4+(o1+2*o2+2*o3+o4)/6

x5=x5+(p1+2*p2+2*p3+p4)/6

x6=x6+(r1+2*r2+2*r3+r4)/6

x7=x7+(s1+2*s2+2*s3+s4)/6

Sheets("1").Cells(k+3;2).Value=x1

Sheets("1").Cells(k+3;3).Value=x2

Sheets("1").Cells(k+3;4).Value=x3

Sheets("1").Cells(k+3;5).Value=x4

Sheets("1").Cells(k+3;6).Value=x5

Sheets("1").Cells(k+3;7).Value=x6

Sheets("1").Cells(k+3;8).Value=x7

Next

End Sub

Function One(x1;x2;x3;x4;x5;x6;x7;a12;a13;a21)'Вер.P1

One=-(a12+a13)*x1+a21*x2

End Function

FunctionTwo(x1;x2;x3;x4;x5;x6;x7;a12;a26;a21;a23;a52;a62;a72)'Вер.P4

Two=a12*x1-(a26+a21+a23)*x2+a52*x5+a62*x6+a72*x7

End Function

Function Three(x1;x2;x3;x4;x5;x6;x7;a13;a23;a34)'Вер.P3

Three=a13*x1+a23*x2-a34*x3

End Function

Function Four(x1;x2;x3;x4;x5;x6;x7;a34;a45)'Вер.Р4

Four=a34*x3-a45*x4

End Function

Function Five(x1;x2;x3;x4;x5;x6;x7;a45;a52)'Вер.Р5

Five=a45*x4-a52*x5

End Function

Function Six(x1;x2;x3;x4;x5;x6;x7;a26;a67;a62)'Вер.Р6

Six=a26*x2-(a67+a62)*x6

End Function

Function Seven(x1;x2;x3;x4;x5;x6;x7;a67;a72)'Вер.Р7

Seven=a67*x6-a72*x7

End Function

 

Приложение 3

Оформление рабочего листа MS EXCEL в задаче 3

  A B C D E F G H I J K L M N O P Q R S
1   P1 P2 P3 P4 P5 P6 R                      
2 1 1 0 0 0 0 0 0 T1 T2 T3 T4 T5 T6 T7 T8 T9    
3 2 9,33E-1 6,61E-02 7,86E-04 4,49E-04 2,49E-05 3,58E-05 1,87E-06 433 4 0,4 0,5 3 2 624 2,6 3,5    
4 3 9,13E-1 8,39E-02 1,04E-03 1,12E-03 1,23E-04 7,14E-05 9,83E-06 12 13 21 23 34 45 52 26 62 67 72
5 4 9,07E-1 8,87E-02 1,11E-03 1,79E-03 2,85E-04 1,00E-04 2,17E-05 0,25 0,002 2,5 0,002 2 0,33 0,5 0,002 0,39 0,39 0,29
6 5 9,04E-1 8,99E-02 1,14E-03 2,41E-03 4,94E-04 1,22E-04 3,60E-05                      
7 6 9,03E-1 9,02E-02 1,14E-03 2,96E-03 7,34E-04 1,39E-04 5,14E-05                      
8 7 9,02E-1 9,03E-02 1,15E-03 3,44E-03 9,92E-04 1,52E-04 6,71E-05                      
9 8 9,01E-1 9,02E-02 1,14E-03 3,86E-03 1,26E-03 1,61E-04 8,26E-05                      
10 9 9,00E-1 9,02E-02 1,14E-03 4,23E-03 1,52E-03 1,68E-04 9,76E-05                      
11 10 8,99E-1 9,02E-02 1,14E-03 4,56E-03 1,78E-03 1,73E-04 1,12E-04                      
12 11 8,98E-1 9,01E-02 1,14E-03 4,84E-03 2,02E-03 1,77E-04 1,25E-04                      
                     
34 33 8,92E-1 9,00E-02 1,13E-03 6,70E-03 4,34E-03 1,87E-04 2,41E-04                      
35 34 8,92E-1 9,00E-02 1,13E-03 6,71E-03 4,37E-03 1,87E-04 2,42E-04                      
36 35 8,92E-1 9,00E-02 1,13E-03 6,72E-03 4,39E-03 1,87E-04 2,43E-04                      
37 36 8,92E-1 9,00E-02 1,13E-03 6,73E-03 4,40E-03 1,87E-04 2,44E-04                      
38 37 8,92E-1 9,00E-02 1,13E-03 6,74E-03 4,42E-03 1,87E-04 2,45E-04                      
39 38 8,92E-1 9,00E-02 1,13E-03 6,75E-03 4,43E-03 1,87E-04 2,46E-04                      
40 39 8,92E-1 9,00E-02 1,13E-03 6,76E-03 4,45E-03 1,87E-04 2,47E-04                      
41 40 8,91E-1 9,00E-02 1,13E-03 6,76E-03 4,46E-03 1,87E-04 2,47E-04                      
42 41 8,91E-1 9,00E-02 1,13E-03 6,77E-03 4,47E-03 1,87E-04 2,48E-04                      
43 42 8,91E-1 9,00E-02 1,13E-03 6,77E-03 4,47E-03 1,87E-04 2,48E-04                      
44 43 8,91E-1 9,00E-02 1,13E-03 6,77E-03 4,48E-03 1,87E-04 2,49E-04                      
45 44 8,91E-1 9,00E-02 1,13E-03 6,78E-03 4,49E-03 1,87E-04 2,49E-04                      
46 45 8,91E-1 9,00E-02 1,13E-03 6,78E-03 4,49E-03 1,87E-04 2,49E-04                      
47 46 8,91E-01 9,00E-02 1,13E-03 6,78E-03 4,50E-03 1,87E-04 2,50E-04                      
48 47 8,91E-01 9,00E-02 1,13E-03 6,78E-03 4,50E-03 1,87E-04 2,50E-04                      
49 48 8,91E-01 9,00E-02 1,13E-03 6,79E-03 4,51E-03 1,87E-04 2,50E-04                      
50 49 8,91E-01 9,00E-02 1,13E-03 6,79E-03 4,51E-03 1,87E-04 2,50E-04                      
51 50 8,91E-01 9,00E-02 1,13E-03 6,79E-03 4,51E-03 1,87E-04 2,51E-04                      
52 51 8,91E-01 9,00E-02 1,13E-03 6,79E-03 4,51E-03 1,87E-04 2,51E-04                      
53 52 8,91E-01 9,00E-02 1,13E-03 6,79E-03 4,52E-03 1,87E-04 2,51E-04                      

Приложение 4

Оформление рабочего листа MS EXCEL в задаче 4

 


  P6   P7
ОГЛАВЛЕНИЕ

Задача 1. Многофакторный регрессионный и корреляционный анализ... 3

Методические указания к решению задачи 1..................................... 6

Пример решения задачи 1.................................................................. 10

Литература к задаче 1......................................................................... 16

Задача 2. Динамическое программирование ............................................ 17

Методические указания к решению задачи 2................................... 18

Литература к задаче 2 ........................................................................ 20

Задача 3. Марковские случайные процессы ............................................. 20

Методические указания к решению задачи 3 .................................. 24

Литература к задаче 3 ........................................................................ 24

Задача 4. Метод Монте-Карло ................................................................... 25

Методические указания к решению задачи 4 .................................. 28

Последовательность решения задачи 4 .......................................... 30

Литература к задаче 4 ........................................................................ 31

Приложение 1 ............................................................................................. 32

Приложение 2 ............................................................................................. 32

Приложение 3 ............................................................................................. 35

Приложение 4 ............................................................................................. 37

 

 


Дата: 2019-05-28, просмотров: 203.