給定n 個(gè)整數(shù)a ,a , ,an 1 2 組成的序列， a n i | |£ ，1 £ i £ n。如果對(duì)于i £ j ，有
0 = å
=
j
k i
k a ，則稱序列區(qū)間i i j a , a , , a +1 為一個(gè)零和區(qū)間，相應(yīng)的區(qū)間長度為j-i+1。
給定n 個(gè)整數(shù)a ,a , ,an 1 2 組成的序列， a n i | |£ ，1 £ i £ n。如果對(duì)于i £ j ，有
0 = å
=
j
k i
k a ，則稱序列區(qū)間i i j a , a , , a +1 為一個(gè)零和區(qū)間，相應(yīng)的區(qū)間長度為j-i+1。
function [R,k,b] = msc(A)
% 多元散射校正
% 輸入待處理矩陣，通過多元散射校正，求得校正后的矩陣
%% 獲得矩陣行列數(shù)
[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
%% 求得結(jié)果
for i = 1:n
Ai = (A(:,i)-b(i))/k(i);
if i == 1
R = Ai;
else
R = [R,Ai];
end
end
