数值分析(第二版)最小二乘法拟合及练习三4.(即课本57页第36题)

  • 格式:doc
  • 大小:28.00 KB
  • 文档页数:2

练习三 4.

x=[-1.00 -0.75 -0.50 -0.25 0 0.25 0.50 0.75 1.00];

>> y=[-0.2209 0.3295 0.8826 1.4392 2.0003 2.5645 3.1334 3.7061 4.2836];

>> a=polyfit(x,y,1)

a =

2.2516 2.0131

>> x1=[-1.00:0.1:1.00];

>> y1=a(2)+a(1)*x1;

>> plot(x,y,'*')

>> hold on

>> plot(x1,y1,'-r')

-1-0.8-0.6-0.4-0.200.20.40.60.81-0.500.511.522.533.544.5

结论:满足条件的一次拟合多项式为:2.0131+2.2516x

x=[-1.00 -0.75 -0.50 -0.25 0 0.25 0.50 0.75 1.00]; >> y=[-0.2209 0.3295 0.8826 1.4392 2.0003 2.5645 3.1334 3.7061 4.2836];

>> a=polyfit(x,y,2)

a =

0.0313 2.2516 2.0001

>> x1=[-1.00:0.1:1.00];

>> y1=a(3)+a(2)*x1+a(1)*x1.^2;

>> plot(x,y,'*')

>> hold on

>> plot(x1,y1,'-r')

-1-0.8-0.6-0.4-0.200.20.40.60.81-0.500.511.522.533.544.5

结论:满足条件的二次拟合多项式为:2.0001+2.2516x+0.0313x^2