二维数据模拟程序.txt

窗口移动多项式平滑演示程序。。。二维数据模拟程序

function  datasimu 

%色谱模拟 
c1=gs(10,3,40);c2=gs(18,4,40);c3=gs(28,3.6,40); 
C=[c1 c2 c3]; 
%光谱模拟 
wav=30; 
s1=0.2*gs(8,3.5,wav)+0.07*gs(20,3,wav); 
s2=0.2*gs(13,4.3,wav)+0.05*gs(24,2.5,wav); 
s3=0.20*gs(10,2.6,wav)+0.15*gs(18,2.5,wav)+0.01; 
S=[s1 s2 s3]; 
%产生二维数据 
X=C*S'; 
%绘图 
 figure;plot(C);title('chromatography profiles'); 
 figure;plot(S);title('spectra'); 
 figure;surf(X);title('Two-way data'); 
%高斯函数 
function  vector=gs(a,b,n) 
x=1:n;vector=exp(-(x-a).^2/(2*b^2)); 
vector=vector'; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%% 拷贝以上所有内容存为m文件直接运行即可 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%