Page 48 - 数学建模算法与应用

Mathematical Modeling Algorithms and Applications
数学建模算法与应用

其对应的加权最小方差为表达式

例 根据 WLS 数据拟合方法，自行编写使用 WLS 方法拟合数据的 M 函数，

然后使用 WLS 方法进行数据拟合。
解 在 M 文件编辑器中输入下面的程序代码
function [th，err，yi] = polyfits(x，y，N，xi，r)
% x，y： 数据点序列

% N： 多项式拟合的系统
% xi： 加权系数逆矩阵
% r： 权重系数的逆矩阵

M = length(x)；
x = x(：)；
y = y(：)；

% 判断调用函数的格式
else if nargin == 3
xi = x；

r = 1；
end

% 求解系数矩阵
A = zeros(M， N+1)；

A(：， N+1) = ones(M， 1)；
for n = N：-1：1
A(：， n) = A(：， n+1).*x；
end

if length(r) == M
for m = 1：M
A(m，：) = A(m，：)/r(m)；

38