CHƯƠNG II : THỰC NGHIỆM
2.4. Chương trình máy tính của các phương pháp hồi quy đa biến
2.4.4. 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 nồng độ X0 (mxk) của m mẫu chuẩn chứa k cấu tử (m hàng, k cột)
Chun ngành hóa phân tích 39 Trường ĐHKHTN
+ 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 [4, 8] %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):
f = V(:,1:n);
% Chuyen doi tap so lieu ban dau va tinh ma tran he so hoi qui voi 3 cau tu chinh:
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;
Chun ngành hóa phân tích 40 Trường ĐHKHTN
Saiso=(X0ktra-Xktra)*100./X0ktra ;
% TINH NONG DO CUA CHAT TRONG MAU BAT KI. %Nhap ma tran do hap thu quang cua mau thuc: Y
X=Y*Fj;
% TINH LOD, LOQ CỦA 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;
- 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
- Kích chuột vào giá trị Saiso, X trong WORKSPACE thu được các dữ liệu mong muốn.
Chun ngành hóa phân tích 41 Trường ĐHKHTN