CHƯƠNG 2 THỰC NGHIỆM
2.4 Câu lệnh thực hiện hồi quy đa biến trong Matlab
Trong các phương pháp CLS, ILS, PLS, PCR chỉ có PCR phù hợp với quy trình phân tích:
- CLS gặp sai số do mơ hình này tính tốn trên cơ sở tập số liệu thơ ban đầu, kết quả tính tốn cuối cùng là kết quả tính tốn trung bình nên tồn phổ phụ thuộc vào giá trị tuyệt đối của độ hấp thụ quang A mà giá trị này lại thay đổi theo từng phép đo và bị ảnh hưởng rất lớn của nền mẫu và của KBr.
- Phương pháp ILS , sai số lớn do chỉ sử dụng một số giá trị hấp thụ tại một số bước sóng đặc trưng do đó nếu số bước sóng lựa chọn không phù hợp sẽ gây ra nguồn sai số.
- Mơ hình PLS khơng phụ thuộc vào tập số liệu thơ đưa vào mà phụ thuộc vào các mơ hình trung gian lựa chọn. Tuy nhiên, khi áp dụng mơ hình này lại có sai số rất lớn, nguyên nhân do mơ hình này giống với mơ hình ILS chỉ sử dụng một số giá
trị số sóng đặc trưng nên mơ hình trung gian tìm được chưa phù hợp với hệ 3 cấu tử. Do vậy, chúng tôi chỉ sử dụng thuật tốn PCR để phân tích tập số liệu.
Phương pháp hồi qui cấu tử chính (PCR)
* Các bước tính tốn PCR trong phần mềm Matlab: - Khởi động phần mềm MATLAB
- Nhập các ma trận dữ liệu trong cửa sổ WORKSPACE
+ Nhập ma trận hàm lượng chất chuẩn X0 (mxk) của m dung dịch chuẩn chứa k cấu tử (m hàng, k cột)
+ Nhập ma trận tín hiệu phân tích Y0 (mxn) (n là số tín hiệu đo) + Nhập tín hiệu phân tích Y của mẫu cần định phân
- Lưu các dữ liệu vừa nhập vào thành 1 file trong Matlab : PCR.mat [6] %Phuong phap PCR:
load PCR.mat;
%Binh phuong tap so lieu chua bien phu thuoc Y0
D = Y0'*Y0;
% Su dung mot trong 3 ham tinh PC de xac dinh cac PC theo cau lenh % su dung ham SVD
[V S] = svd(D);
% Tinh ma tran phan tram phuong sai cua cac PC
d = diag(S)/sum(diag(S))*100;
% Tu gia tri phan tram phuong sai cua cac PC, can cu vao yeu cau cu the cua % bai toan de quyet dinh so PC lam co so cho khong gian moi cua tap so lieu % (n):
% Chuyen doi tap so lieu ban dau va tinh ma tran he so hoi qui:
Yj = Y0*f;
F = inv(Yj'*Yj)*Yj'*X0; Fj=f*F
% Nhap ma tran bien phu thuoc cua k mau can dinh phan va tinh nong do mau % theo cong thuc:X=Y*Fj
% NEU MUON KIEM TRA DO CHINH XAC CUA PHUONG PHAP % Nhap ma tran do hap thu quang cua mau kiem tra:Yktra
%Tinh nong do mau kiem tra theo PCR: Xktra=Yktra*Fj;
%Tinh sai so giua nong do chuan voi nong do xac dinh duoc tu PCR:
Saiso=(X0ktra-Xktra)*100./X0ktra ;
% TINH LOD, LOQ CUA PHUONG PHAP
%Nhap ma tran do lech chuan cua tin hieu do lap lai mau trang: Z %Tinh toan cac gia tri LOD, LOQ theo PCR:
LOD=(3*Z)*Fj; LOQ=(10*Z)*Fj;
% TINH NONG DO CUA CHAT TRONG MAU BAT KI. %Nhap ma tran do hap thu quang cua mau thuc: Y
X=Y*Fj;
- Lưu lại M-file vừa thực hiện được: PCR.m
- Gọi hàm M-file vừa viết được trong cửa sổ COMMAND WINDOW : >> PCR