spirmain.m

一个处理化学数据的MATLAB程序

load y1.txt
y1=flipud(y1);
YY1=-y1(1:1000,2);
T=1.0;
for k=1:1000
   X(k)=k*T;
end  
for i=1:1000
   jx(i)=0.0;
end
plot(jx,'r')                       
hold on                           
%plot(YY1,'k')  

XJD=X(1:454);      
YJD=YY1(1:454)';  

[a,P]=spir0(XJD,YJD,454,3);

fid = fopen('a1.txt','w');
fprintf(fid,'%f',a);
fclose(fid);

plot(X,YY1,'k',XJD,P,'g')











valuee=1.0e-5;
%由三点公式计算原图像在X=454处的一阶导数，取步长为h=1
f=-3.0000*YY1(454)+4.0000*YY1(455)-1.0000*YY1(456);
fd=f/2.0000;

%由三点公式计算拟合图像在X=454处的一阶导数，取步长为h=1
p=1.0000*P(452)-4.0000*P(453)+3.0000*P(454);
pd=p/2.0000;

k=1.0000*fd-1.0000*pd;
k=abs(k);
while k>valuee
  YJD=1.000*YJD-1.000*P; 
  YJD=abs(YJD);
  [a,P]=spir0(XJD,YJD,454,3);
  
  fid = fopen('a2.txt','w');
  fprintf(fid,'%f',a);
  fclose(fid);
  
  p=1.0000*P(452)-4.0000*P(453)+3.0000*P(454);
  pd=p/2.0000;
  k=1.0000*fd-1.0000*pd;
end
hold off
plot(X,YY1,'k',XJD,P,'g')