您现在的位置是:源码地带 > 资源下载

多元散射校正MSC

  • 上传时间: 2020-03-12
  • 上传用户:15275387185
  • 资源积分:2 下载积分
  • 标      签: MSC 多元 散射 校正

资 源 简 介

function [R,k,b] = msc(A)
%   多元散射校正
%   输入待处理矩阵,通过多元散射校正,求得校正后的矩阵

%%  获得矩阵行列数
[m,n] = size(A);

%%  求平均光谱
M = mean(A,2);

%%  利用最小二乘法求每一列的斜率k和截距b
for i = 1:n
    a = polyfit(M,A(:,i),1);
    if i == 1
        k = a(1);
        b = a(2);
    else
        k = [k,a(1)];
        b = [b,a(2)];
    end
end

%%  求得结果
for i = 1:n
    Ai = (A(:,i)-b(i))/k(i);
    if i == 1
        R = Ai;
    else
        R = [R,Ai];
    end
end

相 关 资 源