2.3.2.1. Phương phỏp bỡnh phương tối thiểu nghịch đảo (ILS)
Thuật toỏn ILS giải trong Matlab như sau:
- Nhập ma trận nồng độ C (40x4) của 40 dung dịch chuẩn chứa 4 dạng As phõn tớch vào phần mềm
- Nhập ma trận độ hấp thụ quang A (40x5) - Tớnh ma trận hệ số hồi qui theo cụng thức:
P = inv(A'*A)*A*C
- Nhập ma trận giỏ trị Ax (k*5) của k mẫu cần định phõn - Tớnh nồng độ cỏc dạng theo cụng thức
Nguyễn Thị Thu Hằng CH ƯƠNG 2. THựC NGHIệM
2.3.2.2. Phương phỏp hồi qui cấu tử chớnh (PCR)
Thuật toỏn PCR giải trong Matlab như sau:
- Nhập ma trận nồng độ C (40x4) và ma trận tớn hiệu đo A của 40 dung dịch chuẩn chứa 4 dạng As cần phõn tớch
- Bỡnh phương tập số liệu chứa biến phụ thuộc:
D = A’*A
- Sử dụng một trong 3 hàm tớnh PC để xỏc định cỏc PC. Cõu lệnh sau sử dụng hàm SVD:
[V S] = svd(D)
- Tớnh ma trận phần trăm phương sai của cỏc PC
d = diag(S)/sum(diag(S))*100
- Từ giỏ trị phần trăm phương sai của cỏc PC, căn cứ vào yờu cầu cụ thể của bài toỏn để quyết định chọn số PC làm cơ sở cho khụng gian mới của tập số liệu (n):
f = V(:,1:n)
- Chuyển đổi tập số liệu ban đầu và tớnh ma trận hệ số hồi qui:
Aj = A*f
F = inv(Aj'*Aj)*Aj'*C Fj=f*F
- Nhập ma trận biến phụ thuộc của k mẫu cần định phõn Ax(k*5) và tớnh nồng độ cỏc dạng As trong mẫu theo cụng thức:
Cx=Ax*Fj
Cỏc thao tỏc tớnh sai số và hiệu suất thu hồi sử dụng cỏc cõu lệnh tớnh toỏn thụng thường trờn ma trận trong phần mềm.
Nguyễn Thị Thu Hằng CHƯƠNG 3. KếT QUả Và THảO LUậN
CHƯƠNG 3. KẾT QUẢ VÀ THẢO LUẬN
3.1. NGHIấN CỨU CÁC ĐIỀU KIỆN TỐI ƯU XÁC ĐỊNH HÀM LƯỢNG As(III) BẰNG PHƯƠNG PHÁP HVG – AAS