nihe.m

来自「一个处理化学数据的MATLAB程序」· M 代码 · 共 23 行

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:400);      
YJD=YY1(1:400)';
a=polyfit(XJD,YJD,2);

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

P=polyval(a,XJD')
plot(P,'g')