基于matlab开发的多元散射校正和变量标准化Matlab处理程序，可以对建模前的原始数据进行校正、处理.rar

function [C,E,B,M]=MSC(RR) %定义多元线性校正函数，输入RR是一个包含多个样品光谱值的矩阵，输出是四个参数，RR矩阵是横坐标代表波长，坐标代表 C=mean(RR); %C是对所有原始光谱矩阵RR求平均，得到一个所有样品光谱的平均值，把该平均光谱作为基准光谱 RRT=RR'; %对 [m,n]=size(RRT); %m是波长数，n是样品数 for i=1:n X=RRT(:,i);%特定空间位置的波值 b=mlr(C',X); B(i)=b(1); M(i)=b(2); end; for j=1:n %大循环是对n个样品 for i=1:m %小循环是对n个样品的每个光谱中的任意一个波长下的反射值进行校正 A(i,j)=(RRT(i,j)-B(j))/M(j); end; end; E=A'; function B=mlr(XX,D) %定义mlr函数，多元线性回归函数，输入所有样品平均光谱XX矩阵和一个样品光谱，求平均光谱与一个样品光谱的回归系数。 [k,l]=size(XX); %输入是XX矩阵，求输入矩阵维数 z=ones(k,1); %生成常数项 XX=[z XX]; %输入矩阵中加入常数项 % 求回归系数regression B=(inv(XX'*XX))*XX'*D;