loadmeasurelight.m

自己写的计算色坐标的小程序

function  measurelight=loadmeasurelight
%measurelight=rand(77,1);
light=[4.074306e-004  4.388927e-004  5.166482e-004  6.152719e-004  7.143472e-004  8.069272e-004  9.071190e-004   1.027211e-003   1.150631e-003   1.269596e-003   1.396638e-003   1.525937e-003   1.660326e-003   1.808515e-003   1.971032e-003   2.131645e-003   2.303256e-003   2.496001e-003   2.664299e-003   2.887398e-003   3.091158e-003   3.311144e-003   3.539080e-003   3.778035e-003   3.999974e-003   4.247623e-003   4.483671e-003   4.763809e-003   5.030733e-003   5.307072e-003  5.591081e-003    5.828209e-003   6.081896e-003   6.330627e-003   6.631210e-003   6.922470e-003   7.164997e-003   7.490010e-003   7.748103e-003   8.065600e-003   8.408342e-003   8.709728e-003   9.077325e-003   9.411747e-003   9.747465e-003  1.012576e-002    1.047375e-002   1.081643e-002   1.118476e-002   1.153827e-002   1.192276e-002   1.234816e-002   1.277273e-002   1.314492e-002   1.357903e-002   1.396501e-002   1.436345e-002   1.474644e-002   1.507827e-002   1.545587e-002   1.578515e-002   1.616291e-002   1.660261e-002   1.701633e-002   1.746768e-002   1.789071e-002   1.841099e-002   1.889329e-002   1.931690e-002   1.970791e-002   2.013917e-002   2.060998e-002   2.096251e-002   2.137294e-002   2.173144e-002   2.211609e-002   2.246917e-002   2.270622e-002   2.317161e-002   2.352277e-002   2.393077e-002   2.427152e-002   2.453062e-002   2.482466e-002   2.509615e-002   2.547622e-002   2.560750e-002   2.596728e-002   2.608598e-002   2.624237e-002   2.653779e-002   2.664719e-002   2.693053e-002   2.707939e-002   2.739903e-002   2.754029e-002]
i=375
j=376

measurelight=zeros(1,77);
[l,k]=size(light);
m=0;
for n=1:77
    i=i+5;
    j=j+4;
    m=m+1
    if 2*i==2*j+4
         measurelight(1,n)=(light(1,m)+light(1,m+1))./2;
        
    else
        measurelight(1,n)=light(1,m);
    end
    if i==j+4
        j=j+4;
        i=j;
        measurelight(1,n-1)=light(1,m);
        m=m+1;
        measurelight(1,n)=light(1,m);
    end
end
x=380e-9:4e-9:760e-9;
y=380e-9:5e-9:760e-9;
hold on
plot(x,light,'b')

plot(y,measurelight,'r');
     hold off

%L=length(mea);
%x=1:L;
%p=polyfit(x',mea,10);