一次多项式耦合:
y=[4,6.4,8,8.8,9.22,9.5,9.7,9.86,10,10.20,10.32,10.42,10.5,10.55,10.58,10.6]
x=[1:16]
t=ones(1,16)
X=[x;t]
X2=X.'
P=inv(X*X2)
Q=X*y.'
w=P*Q
scatter(x,y)
hold on
plot(x,w(1,1)*x+w(2,1))
二次多项式耦合:
y=[4,6.4,8,8.8,9.22,9.5,9.7,9.86,10,10.20,10.32,10.42,10.5,10.55,10.58,10.6]
x=[1:16]
t=ones(1,16)
X=[x.^2;x;t]
X2=X.'
P=inv(X*X2)
Q=X*y.'
w=P*Q
scatter(x,y)
hold on
plot(x,w(2,1)*x+w(1,1)*x.^2+w(3,1))
三次多项式耦合:
y=[4,6.4,8,8.8,9.22,9.5,9.7,9.86,10,10.20,10.32,10.42,10.5,10.55,10.58,10.6]
x=[1:16]
t=ones(1,16)
X=[x.^3;x.^2;x;t]
X2=X.'
P=inv(X*X2)
Q=X*y.'
w=P*Q
scatter(x,y)
hold on
plot(x,w(2,1)*x.^2+w(1,1)*x.^3+w(3,1)*x+w(4,1))