clear all;
close all;
n=40;
a=5;
b=10;
c=6;
x1=10*rand(1,n)-6;
x2=x1.^2;
f=@(k1,k2)a*k2+b*k1+c;
y=f(x1,x2)+30*(rand(1,n)-0.5);
sample=[x2;x1;y]';
x_mean=mean(sample);
sample=sample-x_mean;
covmat=1/n*(sample'*sample);
[V,D]=eig(covmat);
lambda=diag(D);
[m,i]=min(lambda);
v=V(:,i);
v=v*(-1/v(end))
b=x_mean(end)-v(1:2)'*x_mean(1:2)'
f_new=@(x)v(1)*x.^2+v(2)*x+b;
fplot(f_new,[-6,4]);
hold on
scatter(x1,y);
LeastSquare Method最小二乘法(非线性转线性)用于多项式拟合
最新推荐文章于 2022-08-26 23:52:12 发布