博文
X,Y坐标轴都存在误差的数据处理方法
||
pro zhang
READCOL,'E:test.dat',F='F,F,F,F',x, y, C, D
FITEXY, x, y, A, B, X_SIG= C, Y_SIG= D, sigma_A_B, chi_sq, q
print,'Intercept=',A
print,'Slope=',B
print,'Intercept standard deviation=',sigma_A_B[0]
print,'Slope standard deviation=',sigma_A_B[1]
print,'Chi_sq=',chi_sq
print,'q - chi-sq probability=',q
end
READCOL,'E:test.dat',F='F,F,F,F',x, y, C, D
FITEXY, x, y, A, B, X_SIG= C, Y_SIG= D, sigma_A_B, chi_sq, q
print,'Intercept=',A
print,'Slope=',B
print,'Intercept standard deviation=',sigma_A_B[0]
print,'Slope standard deviation=',sigma_A_B[1]
print,'Chi_sq=',chi_sq
print,'q - chi-sq probability=',q
end
天文软件包里面有个程序FITEXY,用于处理X，Y都存在误差时候的线性拟合，并且给出了一系列的参数
仅仅Y轴存在误差建议采用下面的程序
result=linfit(x,y,MEASURE_ERRORS=Yerror,/DOUBLE,SIGMA=variable)
result[0]存放的是截距，result[1]存放的是斜率
linfit是idl自带的程序
带权重的线性拟合建议采用下面的程序
SIXLIN, x, y, a, siga, b, sigb, WEIGHT = Weights
https://blog.sciencenet.cn/blog-783377-613435.html
上一篇：两个有用的天文数据库
下一篇：IDL linux版本安装方法
