% 最小二乘法 % 注释的地方需要修改 clear xi=[0;1;2;3]; % 数据点xi yi=[1;0.5;0.5;0.25]; % 数据点yi A=[xi.^0, xi.^2]; % x的线性无关函数族 b=yi.^(-1); % y的函数 % 这里的x函数族和y函数判断方法见同目录下的“最小二乘法.png” a=(A'*A)\(A'*b); vpa(a,10) % 结果有效数字位数