自用程序：各种光谱数据预处理代码matlab

%#										
%#  function [cdatamectest]=center（dataoptnewdata）;			
%#										
%#  AIM: 	Centering along columns rows or double centering		
%#										
%#  PRINCIPLE:  Removal of the column- row- or overall mean from 		
%#              each column row or both respectively 		 		
%# 				 If a test data set is available it can ONLY be 
%#              column centered using the mean from the calibration
%#              data set.
%#
%#
%#  INPUT:	data: （m x n） matrix with m rows and n variables		
%#				opt: optional							
%#		     1 = column centering					
%#		     2 = row centering						
%#		     3 = double centering					
%#          newdata: （mt x n） test matrix with mt rows and n variables
%#					
%#			 							
%#  OUTPUT:	cdata: （m x n） matrix containing centered data			
%#				me: mean vector overall mean （scalar）
%#              newdata: （mt*n） test matrix centered with the mean of data
%#	
%#										
%#  AUTHOR: 	Andrea Candolfi				 			
%#	    			Copyright（c） 1997 for ChemoAc					
%#          	FABI Vrije Universiteit Brussel            			
%#          	Laarbeeklaan 103 1090 Jette					
%#   										
%# VERSION: 1.2 （25/02/2002）							
%#										
%#  TEST:   	I. Stanimirova	& S. Gourv閚ec & M. Zhang
%#										
	
function [cdatamecnewdata]=center（dataoptnewdata）;

[mn]=size（data）;

if nargin==1;
  opt=[4];
  while opt>3 | opt<=0 
    opt=input（‘column centering（1） row centering（2） double centering（3）:‘）;
  end
end


if opt==1			% column centering 
   me=mean（data）;
   cdata=data-ones（m1）*me;
end

if opt==2			% row centering
   me=mean（data‘）‘;
   cdata=data-me*ones（1n）;
end

if opt==3 	% double centering
   me=mean（mean（data））;
   mej=mean（data‘）;
   mei=mean（data）;
   cdata=data-（ones（m1）*mei）-（ones（n1）*mej）‘+（ones（mn）*me）;
end

if exist（‘newdata‘）==1			% center new data
    [mtn]=size（newdata）;
    
    if opt==1				% column centering 
        me=mean（data）;
        cnewdata=newdata-ones（mt1）*me;
    else
        error（‘Row centering and double centering are impossible to perform on a test set‘）;
    end
    
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2158  2015-08-14 15:04  预处理center.m

     文件     521996  2015-08-14 15:04  预处理d1.fig

     文件         52  2015-08-14 15:04  预处理dc.m

     文件       1934  2015-08-14 15:04  预处理DERIV.M

     文件       1346  2015-08-14 15:04  预处理DETREND.M

     文件       2442  2015-08-14 15:04  预处理EXPSMOOT.M

     文件       3588  2015-08-14 15:04  预处理FF.M

     文件        955  2015-08-14 15:04  预处理GENFACT.M

     文件       1146  2015-08-14 15:04  预处理GRAMPOLY.M

     文件       2030  2015-08-14 15:04  预处理HFOUR.M

     文件      15170  2015-08-14 15:04  预处理hs_err_pid2508.log

     文件      15105  2015-08-14 15:04  预处理hs_err_pid3008.log

     文件      14740  2015-08-14 15:04  预处理hs_err_pid3864.log

     文件       2666  2015-08-14 15:04  预处理KENSTON.M

     文件      10000  2015-08-14 15:04  预处理LOWP.M

     文件       2333  2015-08-14 15:04  预处理MSC.M

     文件       1573  2015-08-14 15:04  预处理
ircor.m

     文件       1167  2015-08-14 15:04  预处理
irmaf.m

     文件       1191  2015-08-14 15:04  预处理
irpca2.m

     文件    3051247  2015-08-14 15:04  预处理pre.mat

     文件       2509  2015-08-14 15:04  预处理RANDSEL.M

     文件        118  2015-08-14 15:04  预处理sg_smooth.m

     文件       2532  2015-08-14 15:04  预处理SMOOTH.M

     文件        940  2015-08-14 15:04  预处理SNV.M

     文件        369  2015-08-14 15:04  预处理statmoments.m

     文件        625  2015-08-14 15:04  预处理statxture.m

     文件     894518  2015-08-14 15:04  预处理untitled.bmp

     文件       3128  2015-08-14 15:04  预处理WAVE.M

     文件       1218  2015-08-14 15:04  预处理WEIGHT.M

     文件         89  2015-08-14 15:04  预处理wenlitezheng.m

............此处省略6个文件信息