例题4-1.m

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

%例题4-1
n=4;        %数据点数
p=9;        %波长数
x=[20 40 60 80];                %X矩阵（浓度）
y=[0.022	0.051	0.088	0.123
0.029	0.071	0.113	0.162
0.046	0.102	0.168	0.224
0.069	0.148	0.231	0.316
0.096	0.202	0.313	0.423
0.116	0.240	0.374	0.507
0.120	0.248	0.390	0.520
0.115	0.236	0.367	0.482
0.098	0.193	0.315	0.410];  %Y矩阵（吸光度）
for i=1:p
    [a,b,r]=linear(x,y(i,:),n)               %调用一元线性回归函数（按照推导过程设计）
    yc=a+b*x                            %按模型计算每个浓度下的吸光度
end
plot(y,'.-')