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

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
*****************************************
Result
>> [a,xi,yi,RHX,ERR]=RH(x,y,xi,2)

a =

  1.0e+003 *

    0.0007   -0.0712    2.1526


xi =

   50.5000   56.3000   65.7000


yi =

  431.7617  487.2209  655.4012


RHX =

  1.0e+007 *

    1.1029
    1.0946
    1.0986


ERR =

    0.0056
    0.0019
    0.0016