例题4-5.m

来自「对“曲线的校正”是很多行业和领域都要面临的问题。本人编写并总结了变量校正的源代码」· M 代码 · 共 51 行

clear
%校正样本吸光度矩阵
y=[0.851	0.654	0.817	0.620	0.830	0.731	0.401	0.598	1.045	1.242	0.912	0.813	1.023	0.826	0.989	0.792
0.881	0.674	0.846	0.639	0.849	0.746	0.401	0.608	1.059	1.266	0.921	0.817	1.027	0.82	0.992	0.785
0.942	0.722	0.902	0.682	0.887	0.777	0.407	0.627	1.102	1.322	0.952	0.842	1.046	0.826	1.006	0.786
0.977	0.750	0.933	0.706	0.900	0.786	0.401	0.628	1.112	1.339	0.954	0.840	1.034	0.806	0.990	0.762
1.001	0.767	0.955	0.721	0.901	0.784	0.387	0.621	1.104	1.338	0.941	0.824	1.004	0.770	0.958	0.724
1.005	0.779	0.950	0.724	0.883	0.770	0.376	0.602	1.095	1.320	0.927	0.814	0.973	0.747	0.918	0.692
0.998	0.778	0.939	0.719	0.866	0.756	0.367	0.587	1.078	1.298	0.910	0.800	0.946	0.726	0.887	0.667
0.996	0.779	0.935	0.718	0.859	0.750	0.364	0.581	1.074	1.291	0.905	0.796	0.937	0.72	0.876	0.659
0.975	0.768	0.912	0.705	0.842	0.739	0.365	0.572	1.066	1.273	0.899	0.796	0.933	0.726	0.870	0.663
0.949	0.750	0.887	0.688	0.830	0.731	0.370	0.569	1.060	1.259	0.898	0.799	0.941	0.742	0.879	0.68
0.921	0.731	0.859	0.669	0.812	0.717	0.370	0.560	1.046	1.236	0.889	0.794	0.937	0.747	0.875	0.685
0.894	0.706	0.837	0.649	0.801	0.707	0.369	0.557	1.029	1.217	0.879	0.785	0.937	0.749	0.880	0.692
0.872	0.675	0.830	0.633	0.813	0.714	0.376	0.573	1.022	1.219	0.881	0.782	0.962	0.765	0.920	0.723
0.870	0.649	0.850	0.628	0.851	0.740	0.386	0.608	1.025	1.247	0.893	0.782	1.005	0.783	0.984	0.762
0.860	0.615	0.861	0.616	0.876	0.754	0.388	0.632	1.011	1.256	0.890	0.767	1.027	0.782	1.028	0.783
0.849	0.592	0.864	0.607	0.900	0.771	0.401	0.658	1.018	1.275	0.904	0.776	1.068	0.811	1.084	0.826
0.850	0.567	0.885	0.602	0.920	0.779	0.389	0.672	0.991	1.274	0.885	0.744	1.061	0.778	1.096	0.813
0.847	0.545	0.898	0.596	0.933	0.782	0.380	0.682	0.968	1.270	0.868	0.717	1.054	0.752	1.105	0.803
0.828	0.517	0.891	0.579	0.931	0.775	0.371	0.682	0.939	1.250	0.846	0.690	1.042	0.731	1.105	0.793]';
%校正浓度矩阵
x=[50	50	50	50	150	150	150	150	250	250	250	250	350	350	350	350
50	10	70	30	70	50	10	50	30	70	30	10	50	10	70	30
20	20	15	15	10	10	5	5	20	20	15	15	10	70	5	5]';
p=19;       %波长数
m=3;        %组分数
n=16;       %校正样本数
%未知样本集光度矩阵
yun=[0.640	0.741	0.902	0.980
0.660	0.756	0.911	0.982
0.706	0.785	0.939	0.999
0.732	0.791	0.940	0.985
0.749	0.787	0.925	0.954
0.759	0.765	0.909	0.923
0.757	0.746	0.892	0.897
0.758	0.739	0.886	0.888
0.747	0.724	0.881	0.885
0.730	0.716	0.881	0.894
0.711	0.701	0.872	0.89
0.688	0.696	0.864	0.891
0.659	0.716	0.869	0.918
0.636	0.763	0.885	0.961
0.606	0.798	0.884	0.984
0.585	0.829	0.902	1.025
0.563	0.856	0.885	1.019
0.543	0.874	0.87	1.012
0.516	0.879	0.85	1.001
0.490	0.869	0.814	0.966]';
ns=4;