MATLAB中通过最小二乘法逼近函数

我们使用最小二乘法逼近函数

1.基函数

基函数我们选用基本的 1，$x$，$x^2$ , $x^3$ 函数

2.被逼近函数

我们设为
$$y=2x^3+x+2$$
并加上一个扰动

3.结论

我们主要使用的最小二乘的逼近公式为：
$$A^T\ast W\ast A=C\ast B$$
A为基函数节点处的值所构成的列向量
W为权重矩阵，此处为单位阵
C为我们要求的系数
B为函数在节点处基函数与被逼近函数值的乘积矩阵

4.实现代码

%A
a0 = ones(21,1);
a1 = linspace(0,1,21)';
a2 = a1.^2;
a3 = a1.^3;
A = [a0,a1,a2,a3]
%W权重矩阵
w = a0;
W = diag(w);
%正负一交替的数列
signs = (-1).^(1:21)
%函数值
b = (2*a1'.^3 + a1' + 2 + signs*0.0001)'
%利用Normal_matrix*coefficient = B 求解系数
B = [b'*a0;b'*a1;b'*a2;b'*a3]
Normal_matrix = A'*W*A;
%系数
coefficient = inv(Normal_matrix)*B
%函数值
coefficient_ = [coefficient(4);coefficient(3);coefficient(2);coefficient(1)]
y_pred = A*coefficient_
%%绘制图形
figure(1)
plot(a1,b,'r');
hold on 
plot(a1,y_pred,'b')
axis([0 1 1 6])
%其中红线为原始函数，蓝色的线为逼近的三次函数。

5.图像结果