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

UptoLike

ВУЗ: 

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

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

Рубрика: 

505
1300 print “x(“;i;”),y(“;i;”) “;
1400 input x(i),y(i)
1700 next i
2000 s1=N
2100 s2=0;s3=0;s4=0;s5=0
2200 for i=1 to N
2300 s2=s2+x(i)
2400 s3=s3+y(i)
2500 s4=s4+x(i)⋅x(i)
2600 s5=s5+y(i)⋅x(i)
2610 s6=s6+y(i)⋅x(i)
2620 s7=s2/N
2630 s8=s3/N
2700 next i
3000 d1=s1⋅s4-s2⋅s2
3100 d2=s3⋅s4-s5⋅s2
3200 d3=s1⋅s5-s2⋅s3
4000 a=d2/d1
4100 b=d3/d1
5000 print “Уравнение регрессии”
5100 print “ y = “;: print using “####.####”;a;: print”+ “;: print using
“####.####”; b;:print “ x ”
6000 s=0
6100 for i=1 to N
6200 s=s+(y(i)-(a⋅b⋅x(i)))^2
6300 next i
6500 print “Сумма квадратов отклонений =“;: print using “####.####”;s
6550 z1=s5/s1-s7⋅s8
6560 z2=sqr(s4/s1-s7⋅s7)
6570 z3=sqr(s6/s1-s8^2)
6580 r=z1/(z2⋅z3)
6581 if r>1 then r=2-r
6600 print “Коэффициент корреляции =“;: print using “#.####”;r
6650 print “Хотите сравнить исходные значения y с рассчитанными
значениями ? (Y/N)”
6660 a$=input$(1)
6670 if a$=“Y” or a$=“y” goto 6700
6680 if a$=“n” or a$=“N” then 9999 else 6650
6700 print
6750 print “ x y(измер) Y(теор) “
6800 for i=1 to N
506
6850 print using “####.####”; x(i),y(i),a+b⋅x(i)
6900 next i
7040 goto 9999
9999 end
Приложение III
Программа расчета параболической функции
y = a
0
⋅
x
n
+ a
1
⋅
x
n-1
+ ... + a
n
1 rem Расчет параболы
100 cls
1100 input “Число пар значений “; N
1110 if N<4 goto 7500
1120 dim x(N)
1130 dim y(N)
1200 for i=1 to N
1300 print “x(“;i;”),y(“;i;”) “;
1400 input x(i), y(i)
1700 next i
2000 s1=N
2100 s2=0;s3=0;s4=0;s5=0;s6=0;s7=0;s8=0;s9=0
2200 for i=1 to N
2250 s2=s2+y(i)
2300 s3=s3+x(i)
2400 s4=s4+x(i)⋅x(i)
2500 s5=s5+(x(i)⋅x(i)⋅x(i))
2600 s6=s6+(x(i)⋅x(i)⋅x(i)⋅x(i))
2700 s7=s7+x(i)⋅y(i)
2800 s8=s8+(x(i)⋅x(i))⋅y(i)
2820 o(i)=(y(i)-s2/s1)^2
2830 s9=s9+o(i)
2900 next i
3100 L1=(s4⋅s4-s3⋅s5)
3200 L2=(s3⋅s4-s1⋅s5)
3300 L3=(s4⋅s7-s2⋅s5)
3500 L4=s4⋅s5-s3⋅s6
3600 L5=s4^2-s1⋅s6
3700 L6=s4⋅s8-s2⋅s6
3800 c=(L3⋅L4-L1⋅L6)/(L2⋅L4-L1⋅L5)
3810 b=L3/L1-(L2/L1)⋅c
3811 a=s2/s4-(s3/s4)⋅b-(s1/s4)⋅c 
1300 print “x(“;i;”),y(“;i;”) “;                                           6850 print using “####.####”; x(i),y(i),a+b⋅x(i)
1400 input x(i),y(i)                                                       6900 next i
1700 next i                                                                7040 goto 9999
2000 s1=N                                                                  9999 end
2100 s2=0;s3=0;s4=0;s5=0
2200 for i=1 to N                                                                                                                   Приложение III
2300 s2=s2+x(i)                                                                             Программа расчета параболической функции
2400 s3=s3+y(i)                                                                                      y = a0⋅xn + a1⋅xn-1 + ... + an
2500 s4=s4+x(i)⋅x(i)                                                          1 rem Расчет параболы
2600 s5=s5+y(i)⋅x(i)                                                        100 cls
2610 s6=s6+y(i)⋅x(i)                                                       1100 input “Число пар значений “; N
2620 s7=s2/N                                                               1110 if N<4 goto 7500
2630 s8=s3/N                                                               1120 dim x(N)
2700 next i                                                                1130 dim y(N)
3000 d1=s1⋅s4-s2⋅s2                                                        1200 for i=1 to N
3100 d2=s3⋅s4-s5⋅s2                                                        1300 print “x(“;i;”),y(“;i;”) “;
3200 d3=s1⋅s5-s2⋅s3                                                        1400 input x(i), y(i)
4000 a=d2/d1                                                               1700 next i
4100 b=d3/d1                                                               2000 s1=N
5000 print “Уравнение регрессии”                                           2100 s2=0;s3=0;s4=0;s5=0;s6=0;s7=0;s8=0;s9=0
5100 print “ y = “;: print using “####.####”;a;: print”+ “;: print using   2200 for i=1 to N
“####.####”; b;:print “ x ”                                                2250 s2=s2+y(i)
6000 s=0                                                                   2300 s3=s3+x(i)
6100 for i=1 to N                                                          2400 s4=s4+x(i)⋅x(i)
6200 s=s+(y(i)-(a⋅b⋅x(i)))^2                                               2500 s5=s5+(x(i)⋅x(i)⋅x(i))
6300 next i                                                                2600 s6=s6+(x(i)⋅x(i)⋅x(i)⋅x(i))
6500 print “Сумма квадратов отклонений =“;: print using “####.####”;s      2700 s7=s7+x(i)⋅y(i)
6550 z1=s5/s1-s7⋅s8                                                        2800 s8=s8+(x(i)⋅x(i))⋅y(i)
6560 z2=sqr(s4/s1-s7⋅s7)                                                   2820 o(i)=(y(i)-s2/s1)^2
6570 z3=sqr(s6/s1-s8^2)                                                    2830 s9=s9+o(i)
6580 r=z1/(z2⋅z3)                                                          2900 next i
6581 if r>1 then r=2-r                                                     3100 L1=(s4⋅s4-s3⋅s5)
6600 print “Коэффициент корреляции =“;: print using “#.####”;r             3200 L2=(s3⋅s4-s1⋅s5)
6650 print “Хотите сравнить исходные значения y с рассчитанными            3300 L3=(s4⋅s7-s2⋅s5)
значениями ? (Y/N)”                                                        3500 L4=s4⋅s5-s3⋅s6
6660 a$=input$(1)                                                          3600 L5=s4^2-s1⋅s6
6670 if a$=“Y” or a$=“y” goto 6700                                         3700 L6=s4⋅s8-s2⋅s6
6680 if a$=“n” or a$=“N” then 9999 else 6650                               3800 c=(L3⋅L4-L1⋅L6)/(L2⋅L4-L1⋅L5)
6700 print                                                                 3810 b=L3/L1-(L2/L1)⋅c
6750 print “ x y(измер) Y(теор) “                                          3811 a=s2/s4-(s3/s4)⋅b-(s1/s4)⋅c
6800 for i=1 to N

                                  505                                                                          506

Страницы