Химические методы анализа - 255 стр.

UptoLike

ВУЗ: 

ВСГУТУ | Улан-Удэ

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

Рубрика: 

509
56300 next t
56400 return
63999 end
Приложение 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
510
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 
56300 next t                                                            255 H2=H2+O1(j)⋅O2(j)
56400 return                                                            260 H3=H3+O1(j)⋅O3(j)
63999 end                                                               265 H4=H4+O1(j)⋅O4(j)
                                                                        270 W1=W1+O1(j)⋅O(j)
                                                   Приложение V         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)
 10 rem Множественная регрессия
                                                                        295 W2=W2+O2(j)⋅O(j)
 15 cls: screen 0,0
                                                                        300 H8=H8+O3(j)^2
 20 M$= “! X1 X2 X3             X4    Y     Y(расч) !”
                                                                        305 H9=H9+O3(j)⋅O4(j)
 30 input “Количество переменных: “;N
                                                                        310 W3=W3+O3(j)⋅O(j)
 31 if N<3 goto 610
                                                                        315 H=H+O4(j)^2
 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)                                                       320 W4=W4+O4(j)⋅O(j)
 40 Z1=):Z2=):Z3=0:Z4=0:Z5=0                                            325 next
 60 for i=1 to N                                                        330 L1=W4⋅H1-H4⋅W1
 65 print: print “Значение номер “; i                                  335 L2=H1⋅H-H4^2
 70 input “Значение Y”; Y(i)                                            340 L3=H7⋅H1-H2⋅H4
 85 Z5=Z5+Y(i)                                                          345 L4=W2⋅H1-H2⋅W1
 90 input “Значение Х1”; X1(i)                                          350 L5=H2⋅H3-H1⋅H6
105 Z1=Z1+X1(i)                                                         355 L6=H9⋅H1-H3⋅H4
110 input “Значение Х2”;Х2(i)                                           360 L7=W3⋅H1-W1⋅H3
125 Z2=Z2+X2(i)                                                         370 L9=H8⋅H1-H3^2
130 input “Значение Х3”;Х3(i)                                           375 L=H1⋅H5-H2^2
145 Z3=Z3+X3(i)                                                         380 M1=L9⋅L-L5^2
150 input “Значение Х4”;Х4(i)                                           390 M2=L6⋅L+L5⋅L3
165 Z4=Z4+X4(i)                                                         395 M3=L7⋅L+L4⋅L5
166 next i
                                                                        400 M4=L1⋅L-L3⋅L4
167 print “Все верно (Y/N) ?”: ASD$=input$(1)
                                                                        405 M5=L2⋅L-L3^2
168 if ASD$=“Y” or ASD$=“y” goto 170
169 if ASD$=“N”or ASD$=“n” then gosub 1000 else goto 167                407 if S4=0 then 415 else D=(M4⋅M1-M2⋅M3)/(M5⋅M1-M2^2)
170 S1=Z1/N: S2=Z2/N: S3=Z3/N: S4=Z4/N                                  410 goto 420
180 S=Z5/N                                                              415 D=0
210 for j=1 to N                                                        420 G1=M3-M2⋅D
220 O(j)=Y(j)-S                                                         421 G=G1/M1
225 O1(j)=X1(j)-S1                                                      425 B=(L4+L5⋅G-L3⋅D)/L
230 O2(j)=X2(j)-S2                                                      430 A=(W1-H2⋅B-H3⋅G-H4⋅D)/H1
235 O3(j)=X3(j)-S3                                                      440 for i=1 to N
240 O4(j)=X4(j)-S4                                                      450 K(i)=S+A⋅O1(i)+B⋅O2(i)+G⋅O3(i)+D⋅O4(i)
250 H1=H1+o1(j)^2                                                       460 V=V+K(i)
                                                                        470 next

                                509                                                                   510

Страницы