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

UptoLike

ВУЗ: 

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

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

Рубрика: 

265
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
3812 m1=o:m3=0
3820 for i=1 to N
3850 s10(i)=(a⋅x(i)^2+b⋅x(i)+c)
3860 m1=m1+s10(i)
3870 m2= m1/s1
3875 o1(i)=(s10(i)-m2)^2
3880 m3=m3+o1(i)
3881 next i
3890 R1=m3/s9
3891 R=sqr(R1)
3892 if R<1 goto 3900
3893 R=2-R
3900 cls:print “Уравнение параболы ”
4000 print “Y=“ using “####.#####”;a;:print”⋅x^2+”
using “####.#####”;b;:print”x+”using “####.#####”;c
4100 print “Коэффициент корреляции R =“
using “#.####”;R
4150 print “Хотите сравнить исходные значения y
с рассчитанными значениями Y (Y/N)”;
4200 a$=input$(1)
4250 if a$=“N” goto 7500
4251 if a$=“n” goto 7500
4253 if a$=“Y” or a$=“y” then 4400 else 4150
266
4400 print “ x y(измер) Y(теор) “
4450 for i=1 to N
4500 print using “####.#####”;x(i),y(i),(a⋅x(i)^2+
b⋅x(i)+c)
4550 next i
7500 end
Приложение IV
Программа для решения системы
из n уравнений методом Гаусса
0 rem “Программа Гаусса-Жордана”
1rem **********************************************
2rem ***Решение системы из n линейных уравнений***
3rem*** c n неизвестными методом Гаусса-Жордана***
4rem **********************************************
40000 cls:screen 0
50000 dim a(20,18)
50100 input “Число неизвестных”;N
50101 if N<2 goto 63999
50200 print “Ввод расширенной матрицы
коэффициентов”
50300 for i=1 to N
50400 for j=1 to N
50500 print “a(“;i;”,”;j;”)= “;:input a(i,j)
50600 next j
50700 print “b(“;i;”)= “;:input a(i,N+1)
50800 next i: print
50900 for s=1 to N
51000 for t=s to N
51100 if a(t,s)><0 then 51300 
                            265                                                     266


2830 s9=s9+o(i)                                       4400 print “ x y(измер) Y(теор) “
2900 next i                                           4450 for i=1 to N
3100 L1=(s4⋅s4-s3⋅s5)                                 4500 print using “####.#####”;x(i),y(i),(a⋅x(i)^2+
3200 L2=(s3⋅s4-s1⋅s5)                                      b⋅x(i)+c)
3300 L3=(s4⋅s7-s2⋅s5)                                 4550 next i
3500 L4=s4⋅s5-s3⋅s6                                   7500 end
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
                                                                                           Приложение IV
3811 a=s2/s4-(s3/s4)⋅b-(s1/s4)⋅c
                                                                     Программа для решения системы
3812 m1=o:m3=0
                                                                      из n уравнений методом Гаусса
3820 for i=1 to N
3850 s10(i)=(a⋅x(i)^2+b⋅x(i)+c)
                                                       0 rem “Программа Гаусса-Жордана”
3860 m1=m1+s10(i)
                                                      1rem **********************************************
3870 m2= m1/s1
                                                      2rem ***Решение системы из n линейных уравнений***
3875 o1(i)=(s10(i)-m2)^2
                                                      3rem*** c n неизвестными методом Гаусса-Жордана***
3880 m3=m3+o1(i)
                                                      4rem **********************************************
3881 next i
                                                      40000 cls:screen 0
3890 R1=m3/s9
                                                      50000 dim a(20,18)
3891 R=sqr(R1)
                                                      50100 input “Число неизвестных”;N
3892 if R<1 goto 3900
                                                      50101 if N<2 goto 63999
3893 R=2-R
                                                      50200 print “Ввод расширенной матрицы
3900 cls:print “Уравнение параболы ”
                                                            коэффициентов”
4000 print “Y=“ using “####.#####”;a;:print”⋅x^2+”    50300 for i=1 to N
using “####.#####”;b;:print”x+”using “####.#####”;c   50400 for j=1 to N
4100 print “Коэффициент корреляции R =“               50500 print “a(“;i;”,”;j;”)= “;:input a(i,j)
      using “#.####”;R                                50600 next j
4150 print “Хотите сравнить исходные значения y       50700 print “b(“;i;”)= “;:input a(i,N+1)
     с рассчитанными значениями Y (Y/N)”;             50800 next i: print
4200 a$=input$(1)                                     50900 for s=1 to N
4250 if a$=“N” goto 7500                              51000 for t=s to N
4251 if a$=“n” goto 7500
                                                      51100 if a(t,s)><0 then 51300
4253 if a$=“Y” or a$=“y” then 4400 else 4150

Страницы