test_fitting.m

来自「数值分析中的数值逼近方法以及曲线拟合算法」· M 代码 · 共 16 行

function Test_fitting()
clc;
clear;
x=[0 0.25 0.5 0.75 1 1.25 1.5]';
y=[1 1.284 1.6487 2.117 2.7183 3.4903 4.4817]';
t=0:0.01:2;
[Fy1 a1]=Fitting(x,y,3);
[Fy2 a2]=Fitting(x,y,4);
Poly1=polyval(fliplr(a1'),t);
Poly2=polyval(fliplr(a2'),t);
plot(x,y,'r*');
hold on;
plot(t,Poly1,'b');
hold on;
plot(t,Poly2,'g');
legend('原始参数','3次曲线拟合','4次曲线拟合');