MATLAB语句求帮忙修改
l = 0.1; ; %装蔗糖溶液的玻璃管长度
yx = 5.225; %待测蔗糖溶液的旋光度
x = [0 20 40 60 100]; %已知蔗糖溶液的浓度
y = [0 1.250 2.825 4.100 6.875]; %测得蔗糖溶液的旋光度
coef = polyfit(x,y,l); %对x,y进行一元线性回归
k = coef(1), b = coed(2); %k为正比例系数,b为截距
r = corrcoef(x,y); %计算x、y的线性相关系数
a = k/1; %蔗糖溶液的旋光率
stda = sqrt((1/r 2-1)/3)*a %旋光率的标准误差
xx = yx/k; %待测蔗糖溶液的浓度
rexx = abs(xx-80)/80; %待测蔗糖溶液浓度的相对误差
xbest = [0 20 40 60 100 120];
ybest = k*xbest+b; %xbest、ybest对已知浓度和旋光度进行拟合的数据对
axis([0,120,0,8]); %纵横坐标轴的取值范围
plot(x, y, ‘*’, xbest, ybest); %坐标系中描点和根据描点进行直线拟合 展开
function xuanguang %编制函数名称为xuanguang
l = 2; %装蔗糖溶液的玻璃管长度
yx = 5.225; %待测蔗糖溶液的旋光度
x = [0 20 40 60 100]; %已知蔗糖溶液的浓度
y = [0 1.250 2.825 4.100 6.875]; %测得蔗糖溶液的旋光度
coef = polyfit(x,y,l); %对x,y进行一元线性回归
k = coef(1), b = coef(2); %k为正比例系数,b为截距
r = corrcoef(x,y); %计算x、y的线性相关系数
r = r(2);
a = k/1; %蔗糖溶液的旋光率
stda = sqrt((1/r^2-1)/3)*a %旋光率的标准误差
xx = yx/k; %待测蔗糖溶液的浓度
rexx = abs(xx-80)/80; %待测蔗糖溶液浓度的相对误差
xbest = [0 20 40 60 100 120];
ybest = k*xbest+b; %xbest、ybest对已知浓度和旋光度进行拟合的数据对
axis([0,120,0,8]); %纵横坐标轴的取值范围
plot(x, y, '*', xbest, ybest);
改了两个地方,第二行的I=2,指的是拟合方程的阶数,以及r=r(2),是获取相关系数。然后就是一些中英文标点符号之类的地方,不知道是复制网页的原因还是您原来程序中的问题,你注意一下就好。
结果如下: