химические методы анализа
.pdf1300 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 |