Химические методы анализа. Танганов Б.Б. - 132 стр.

UptoLike

ВУЗ: 

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

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

Рубрика: 

263
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
6850 print using “####.####”; x(i),y(i),a+b⋅x(i)
6900 next i
7040 goto 9999
9999 end
264
Приложение 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 
                              263                                                        264


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                                                                                              Приложение III
6500 print “Сумма квадратов отклонений =“;:                         Программа расчета параболической функции
     print using “####.####”;s                                                y = a0⋅xn + a1⋅xn-1 + ... + an
6550 z1=s5/s1-s7⋅s8                                          1 rem Расчет параболы
6560 z2=sqr(s4/s1-s7⋅s7)                                   100 cls
6570 z3=sqr(s6/s1-s8^2)                                   1100 input “Число пар значений “; N
6580 r=z1/(z2⋅z3)                                         1110 if N<4 goto 7500
6581 if r>1 then r=2-r                                    1120 dim x(N)
6600 print “Коэффициент корреляции =“;: print using       1130 dim y(N)
    “#.####”;r                                            1200 for i=1 to N
6650 print “Хотите сравнить исходные значения y           1300 print “x(“;i;”),y(“;i;”) “;
    с рассчитанными значениями ? (Y/N)”                   1400 input x(i), y(i)
6660 a$=input$(1)                                         1700 next i
6670 if a$=“Y” or a$=“y” goto 6700                        2000 s1=N
6680 if a$=“n” or a$=“N” then 9999 else 6650              2100 s2=0;s3=0;s4=0;s5=0;s6=0;s7=0;s8=0;s9=0
6700 print                                                2200 for i=1 to N
6750 print “ x y(измер) Y(теор) “                         2250 s2=s2+y(i)
6800 for i=1 to N                                         2300 s3=s3+x(i)
6850 print using “####.####”; x(i),y(i),a+b⋅x(i)          2400 s4=s4+x(i)⋅x(i)
6900 next i                                               2500 s5=s5+(x(i)⋅x(i)⋅x(i))
7040 goto 9999                                            2600 s6=s6+(x(i)⋅x(i)⋅x(i)⋅x(i))
9999 end                                                  2700 s7=s7+x(i)⋅y(i)
                                                          2800 s8=s8+(x(i)⋅x(i))⋅y(i)
                                                          2820 o(i)=(y(i)-s2/s1)^2

Страницы