加权最小二乘法matlab

加权最小二乘法是一种优化方法，它可以用来拟合一组数据点到一个函数模型。如果每个数据点的权值不同，那么我们应该如何求解加权最小二乘法呢？下面是在MATLAB中实现加权最小二乘法的一个例子：

假设我们有一组数据点$(x_1,y_1), (x_2,y_2), \dots ,(x_n,y_n)$，每个数据点都有一个权值$w_i$。我们想要用一个函数模型$f(x)=a_1x+a_2$拟合这些数据点。我们的目标是找到$a_1$和$a_2$，使得加权残差平方和最小：

$$\sum_{i=1}^{n}w_i(y_i-f(x_i))^2$$

为了求解这个问题，我们可以将其转化为一个线性最小二乘问题。我们定义一个向量$\mathbf{y}$，其中每个元素为$y_i$，一个矩阵$\mathbf{X}$，其中每一行为$[x_i,1]$，还有一个对角矩阵$\mathbf{W}$，其中对角线元素为$w_i$。然后，我们可以将加权残差平方和表示为：

$$(\mathbf{y}-\mathbf{X}\mathbf{a})^T\mathbf{W}(\mathbf{y}-\mathbf{X}\mathbf{a})$$

其中$\mathbf{a}$是一个包含$a_1$和$a_2$的向量。我们的目标是最小化这个式子，即：

$$\mathbf{a}=(\mathbf{X}^T\mathbf{W}\mathbf{X})^{-1}\mathbf{X}^T\mathbf{W}\mathbf{y}$$

下面是在MATLAB中实现加权最小二乘法的代码：

% 输入数据点和权值
x = [1,2,3,4,5];
y = [1.2,3.5,4.6,5.9,7.1];
w = [1,3,2,4,1];

% 构造矩阵X和向量y
X = [x',ones(size(x'))];
y = y';

% 构造对角矩阵W
W = diag(w);

% 计算最小二乘解
a = inv(X'*W*X)*X'*W*y;

% 绘制结果
xx = 0:0.1:6;
yy = a(1)*xx+a(2);
plot(x,y,'o',xx,yy)

在这个例子中，我们输入了数据点$x$和$y$以及权值$w$。然后，我们构造了矩阵$\mathbf{X}$和向量$\mathbf{y}$以及对角矩阵$\mathbf{W}$。最后，我们使用公式$\mathbf{a}=(\mathbf{X}^T\mathbf{W}\mathbf{X})^{-1}\mathbf{X}^T\mathbf{W}\mathbf{y}$计算最小二乘解$a_1$和$a_2$。最后，我们使用这些参数绘制了拟合曲线。