rh.m

来自「用matlab仿真对光谱进行定标和很好的一个随机数发生器」· M 代码 · 共 16 行

function [a,xi,yi,RHX,ERR]=RH(x,y,xi,m)
a=polyfit(x,y,m);
x2=50:0.05:70;
y2=a(1)*x2.^2+a(2)*x2+a(3);
plot(x,y,'*',x2,y2,'-r')
grid on
grid minor
yi=spline(x,y,xi);
RH=1.096776*1.0e+007;
RHX=zeros(3,1);
ERR=zeros(3,1);
for k=1:3
    RHX(k)=1.0e+009/yi(k)/(0.25-1/(6-k)^2);
    ERR(k)=abs(RHX(k)-RH)/RH;
end