Lqp trình thù tue giai hen phuang trìn hn dn bang phuang phăp khù Gauss

Một phần của tài liệu Nghiên cứu phương pháp hấp thụ quang phân tử xác định đồng thời các chất có phổ hấp thụ xen phủ nhau dựa trên thuật toán lọc Kalman (Trang 61)

O buóc song X2, chùng ta lai băt dău su dung C, ^j,^ lăm giă tri uóc

3.1.3.Lqp trình thù tue giai hen phuang trìn hn dn bang phuang phăp khù Gauss

3.1.3. Lqp trình thù tue giai he n phuang trình n dn bang phuang phăp khù Gauss khù Gauss

Sau khi thu dugc he n phuong trình n ăn, ta phăi lăp thù tue giăi he phuong trình tuyĩn tinh n ăn ò trĩn theo phuang phăp khù Gauss sĩ thu dugc căc giă tri CJ, Q , . . ., C^ lă nóng dò cùa căc cău tu co phó hăp thu quang phăn tu xen phù nhau trong hòn hgp.

He phuang trình tuyĩn tfnh n phuong trình n ăn so (2.11) dugc giăi theo phuong phăp khù Gauss. Phuong phăp khù Gauss giăi he phuang trình tuyĩn tinh dugc lăp trình tuong dói de, tuy nhiĩn viĩc tfnh toăn khă phùc tap nĩn khóng thfch hgp cho viĩc giăi bang taỵ Nguyĩn tăc cùa phuong phăp lă căc phuang trình ban dău phăi dugc biĩn dói sao cho ma tran he so chùa giă tri bang 1 tai duòng chĩo chinh vă bang 0 tai căc vi tri phfa duói vă ò ben trai duòng chĩo chfnh.

Giă su co he n phuong trình tuyĩn tfnh n ăn so

aiix,+ai2X2 + . . . + ai„x„=â„,ij (3.1.1) ^21^1 + a22X2 + . . . + a2nX„=a2(„^ij (3.1.2)

^niXi + a,2X2 + . . . + âx„=a„(„^ij (3.1.n) Căch biĩh dói nhu sau: Thoat tiĩn ta khù x, bang căch chia că hai ve

cùa phuong trình (3.1.1) cho â i :

x, + b,2X2 + . . . + bi„x,-bi(„,ij (3.2.1)

trong do ^\k/^u = ^ik vói k = 2 ^ (n+1).

Thay giă tri cùa Xi ò phuong trình (3..2.1) văo phuong trình (3.1..2) ta co:

(^22 - a2iai2)x2 + . . . + (a2n - a2ia„)x, = a2(n^i) - ^ii^n^i) (3.2.2)

Duói dang tóng quăt khi thay văo phuong trình (3.12.n) ta co:

(an2 - a„iai2)x2 + . . . + (a„, - âjâJx, = a„(,^ij - ^^â^^^^^ (3.2.n)

Day chfnh lă he (n-1) phuang trình mói khòng chùa ăn x^ {tu phuang trình (3.2.2) dĩn phuang trình (3.2.n)} Tiĩp tue thuc hiĩn viĩc khù X2 trong he (n-1) phuang trình mói năy bang căch chia că hai ve eùa căc phuong trình cho he so eùa X2 lă (a„2 - aniai2) ta thu dugc he (n-1) phuong trình mói chùa ăn lă X2, X3, . . . x^ vói he so eùa X2 lă 1. Tinh toăn tuong tu ta sĩ dugc he (n-2) phuang trình khóng chùa ăn Xj vă X2 v.v...

Sau khi thuc hiĩn phĩp khù Gauss, he phuong trình (3.1.1) -^ (3.1.n) sĩ chuyĩn thănh dang

Xi + b,2X2 + . . . + b„X„ = b,(,,,j (3.3.1) X2 + . • . + b2,x„ = b2(,,ij (3.3.2)

Xn = bn(n+i) (3.3.n) Thay thỉ giă tri x„ ò phuong trình (3.3.n) văo phuong trình ò ngay trĩn

nò sĩ tìm dugc ăn so tiĩp theo x^^.^ vă tiĩp tue nhu văy sĩ tìm dugc căc giă tri Xj con lai [5, 28, 29].

Sau khi lap dugc chuang trình xăc dinh dóng thòi căc cău tu su dung sai so tuong dói cùa phĩp do quang (phu lue 1), chùng tói dă nghiĩn cùu kiĩm chùng phuong phăp vă chuong trình dă lap bang căch ăp dung tinh toăn xăc dinh dóng thòi căc cău tu co phó hăp thu quang phăn tu xen phù nhaụ Kĩt qua

dugc chi ra a phăn 3.2.

3.1.4. Quy trình chung xăc dinh dóng thòi căc cău tă co pho hăp thu quang phăn tu xen phù nhau theo phuang phăp su dung sai so tuang dói

+ Pha căc dung dich chuan cùa căc cău tu riĩng rĩ cùng nhu căc dung dich hón hgp, lua chgn thuóc thù thfch hgp de tao phùc măụ Khăo săt diĩu

kiĩn tói uu cho su tao phùc mău (khoăng buóc song, pH, lugng du thuóc thù, khoăng nóng dò tuyĩn tinh,. . . )

+ Do phó hăp thu quang phăn tu eùa căc dung dich phùc mău chùa căc cău tu riĩng rỉ cùng nhu cùa hón hgp ò khoăng buóc song thfch hgp, so lieu dugc ghi thănh căc file co phăn mò ròng lă ".txt"; ".dat".

+ Su dung file so lieu góm: so cău tu trong hòn hgp; so buóc song do do hăp thu quang; giă tri nóng dò cùa căc cău tu trong căc mău riĩng rĩ vă trong hòn hgp; dò hăp thu quang cùa căc cău tu riĩng rĩ vă cùa hòn hgp. Tiĩn hănh tinh toăn theo chuong trình dă lăp. Kĩt qua nóng dò eùa căc cău tu vă sai so tuong dói cùa phĩp xăc dinh dugc hiĩn thi lĩn man hình vă ghi ra file KETQUẠTXT.

Một phần của tài liệu Nghiên cứu phương pháp hấp thụ quang phân tử xác định đồng thời các chất có phổ hấp thụ xen phủ nhau dựa trên thuật toán lọc Kalman (Trang 61)