* 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 (25x2) của 25 dung dịch chuẩn chứa 2 dạng Sb phân tích vào phần mềm.
Vũ Thị Thảo Thực nghiệm
- Nhập ma trận độ hấp thụ quang A (25x2) - 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
Cx = Ax*P
* 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 25 dung dịch chuẩn chứa 2 dạng Sb 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:
Vũ Thị Thảo Thực nghiệm
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.