Научная электронная библиотека
Монографии, изданные в издательстве Российской Академии Естествознания

Приложение V. Программа расчета функции по аргументам методом многоуровневого моделирования (ММУМ) параметров химических систем

Программа расчета функции по аргументам методом многоуровневого моделирования (ММУМ) параметров химических систем

10 rem Многоуровневое моделирование

15 cls: screen 0,0

20 M$= “!  X1  X2  X3  X4  Y  Y(расч)  !”

30 input “Количество переменных:  “;N

31 if N<3 goto 610

35 dim (N),X1(N),X2(N),X3(N),X4(N),O(N),O1(N),O2(N),O3(N), O4(N), K(N),VO(N)

40 Z1=):Z2=):Z3=0:Z4=0:Z5=0

60 for i=1 to N

65 print: print “Значение номер  “; i

70 input “Значение Y”; Y(i)

85 Z5=Z5+Y(i)

90 input “Значение Х1”; X1(i)

105 Z1=Z1+X1(i)

110 input “Значение Х2”;Х2(i)

125 Z2=Z2+X2(i)

130 input “Значение Х3”;Х3(i)

145 Z3=Z3+X3(i)

150 input “Значение Х4”;Х4(i)

165 Z4=Z4+X4(i)

166 next i

167 print “Все верно (Y/N) ?”: ASD$=input$(1)

168 if ASD$=“Y” or ASD$=“y” goto 170

169 if ASD$=“N”or ASD$=“n” then gosub 1000 else goto 167

170 S1=Z1/N: S2=Z2/N: S3=Z3/N: S4=Z4/N

180 S=Z5/N

210 for j=1 to N

220 O(j)=Y(j)-S

225 O1(j)=X1(j)-S1

230 O2(j)=X2(j)-S2

235 O3(j)=X3(j)-S3

240 O4(j)=X4(j)-S4

250 H1=H1+o1(j)^2

255 H2=H2+O1(j)×O2(j)

260 H3=H3+O1(j)×O3(j)

265 H4=H4+O1(j)×O4(j)

270 W1=W1+O1(j)×O(j)

275 W=W+O(j)^2

280 H5=H5+O2(j)^2

285 H6=H6+O2(j)×O3(j)

290 H7=H7+O2(j)×O4(j)

295 W2=W2+O2(j)×O(j)

300 H8=H8+O3(j)^2

305 H9=H9+O3(j)×O4(j)

310 W3=W3+O3(j)×O(j)

315 H=H+O4(j)^2

320 W4=W4+O4(j)×O(j)

325 next

330 L1=W4×H1-H4×W1

335 L2=H1×H-H4^2

340 L3=H7×H1-H2×H4

345 L4=W2×H1-H2×W1

350 L5=H2×H3-H1×H6

355 L6=H9×H1-H3×H4

360 L7=W3×H1-W1×H3

370 L9=H8×H1-H3^2

375 L=H1×H5-H2^2

380 M1=L9×L-L5^2

390 M2=L6×L+L5×L3

395 M3=L7×L+L4×L5

400 M4=L1×L-L3×L4

405 M5=L2×L-L3^2

407 if S4=0 then 415 else D=(M4×M1-M2×M3)/(M5×M1-M2^2)

410 goto 420

415 D=0

420 G1=M3-M2×D

421 G=G1/M1

425 B=(L4+L5×G-L3×D)/L

430 A=(W1-H2×B-H3×G-H4×D)/H1

440 for i=1 to N

450 K(i)=S+A×O1(i)+B×O2(i)+G×O3(i)+D×O4(i)

460 V=V+K(i)

470 next

480 Vs=V/N

490 for j=1 to N

500 VO(j)=K(j)-VS

510 R2=R2+VO(j)^2

520 next

529 cls

530 R1=R2/W

540 R=sqr(R1)

545 K1=S-A×S1-B×S2-G×S3-D×S4

581 print: print

582 F$=string$(70,”-”): print F$

583 print M$

584 F$=string$(70,”-”): print F$

585 for i=1 to N

586 print using “#####.##”; X1(i),X2(i),X3(i),X4(i),Y(i),K(i)

587 next

588 print F$

589 print: print

590 print “Y=“ using “####.####”;A;:print “×X1+(“using “####.####”; B;:print”)×X2+(“using “####.####”;G; :print”)×X3+(“using “####.####”;  D;: print”)×X4+(“using “####.####”; K1;:print”)”

600 print “R= “ using “#.####”;R

610 end

1000 rem *****************************************

1050 rem*Подпрограмма обработки ошибок ввода***

rem* данных *********************************

1080 rem******************************************

1100 input “Какое значение не верно (введите его номер)”;i

1200 Z5=Z5-Y(i)

1210 Z1=Z1-X1(i); Z2=Z2-X2(i); Z3=Z3-X3(i); Z4=Z4-X4(i)

1250 print “Значение номер  “,i

1260 input “Значение Y”; Y(i)

1270 Z5=Z5+Y(i)

1280 input “Значение Х1”; X1(i)

1290 Z1=Z1+X1(i)

1300 input “Значение Х2”; X2(i)

1310 Z2=Z2+X2(i)

1320 input “Значение Х3”; X3(i)

1330 Z3=Z3+X3(i)

1340 input “Значение Х4”; X4(i)

1350 Z4=Z4+X4(i)

1360 print “Будут ли меняться другие значения?”

1370 XC$=input$(1)

1380 if XC$=“Y” or XC$=“y” goto 1100

1390 if XC$=“N” or XC$=“n” then return else goto 1360


Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074