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.1. Sa dò thuqt toăn biĩn dai hem phuang trình năn thănh hen phuang trình n dn vói sai so tuang dói cùa phĩp do quang nhò nhăt
k : = l
y[k]:=0
[k, i]:= x[k,i]+ s|j, i]* sU, k]* a|j]
j : = l
y[k]:=y[k]+ s[j, k]* a|j] * a[j]
Hình 3.1. Sa dò thuăt toăn biĩn dói he m phuang trình n ăn thănh he n phuang
Trong so dò thuăt toăn ò hình 3.1, Căc so liĩu ban dău bao góm: so cău tu (n); so buóc song (m) ; he so hăp thu mol phăn tu cùa căc cău tu ùng vói căc buóc song tuong ùng s|j, i] vă giă tri dò hăp thu quang eùa dung dich chùa hòn hgp n cău tu a|j] dugc nhăp tu mot file so lieụ
Thuc ra giă tri s[j, i] khóng dugc nhap truc tiĩp ma dugc tfnh tu giă tri dò hăp thu quang a[j,i] do dugc cùa tùng cău tu riĩng rĩ ò m buóc song.
Viĩc tinh căc he so x,^; theo phuang trình (2.10) dugc thuc hiĩn bang 3 vòng lăp vă ehi dùng lai khi k = n; i = n vă j = m; căc he so yi, dugc thuc hiĩn bang 2 vòng lap vă chi dùng lai khi k = n vă j = m.
3.1.2. Lăp trình thù tue biĩn dòi he m phuang trình n ăn thănh he n phuang trình n dn văi sai so tuang dói cùa phĩp do quang nhò nhăt (biĩn doi ma tran khòng vuòng thănh ma tran vuòng)
Tu so dò thuăt toăn biĩn dói he m phuong trình n ăn thănh he n phuang trình n ăn (hình 3.1). Tiĩn hănh lăp trình theo ngón ngù PASCAL de tfnh căc
he so x^i, x^2^ . . . Xki, . . . Xj^ vă căc giă tri y^ tuong ùng. Day lă phăn quan
trgng nhăt cùa phuong phăp tinh toăn trong viĩc ăp dung phuong phăp bình phuong tói thiĩu cho he da biĩn. Doan chuang trình dugc lap nhu sau:
for k:=l to n do begin
y[k]:=0;
for j:=l to m do y [k] :=y [k]-he[j ,k] *a[j]*a[j]; end;
for k:=l to n do for i:=l to n do
begin x[k,l]:=0;
for j:=l to m do x[k,1]:=x[k,i]+e[j,1]*e[j,k]* a[j]; end;
Vói n, m, a[j], e[j, k]; e[j, i] da dugc nhap truóc do (hoac tu băn phim hoac tu file so lieu do phó).
Sau khi biĩn dói he m phuong trình n ăn thănh bĩ n phuang trình n ăn. Viĩc giăi he phuang trình năy co nhiĩu căch, chùng tói su dung phuang phăp