函数逼近与最小二乘法.docxVIP

函数逼近与最小二乘法 实验目的： 熟悉Mablab编程； 2?学习最小二乘法及程序设计算法 实验题目： 1.对于给函数1 1. 对于给函数 1 在区间【T, 1］上取 Xi二-1+0. 2i (i二0, 1-10),试求 3 1 + 25% 次曲线拟合，试画出你和曲线并打印出方程。 有实验给出数据表 x 0.0 0.1 0.2 0.3 0.5 0.8 1.0 y 1.0 0.41 0.50 0.61 0.91 2.02 2.46 试求3次、4次多项式的曲线拟合，再根据数据曲线形状，求一个另外函数的拟合曲线，用 3.给泄数据点 3.给泄数据点 Xi 0 0.5 0.6 0.7 0.8 0.9 1.0 yj 1 1.75 1.96 2.19 2.44 2.71 3.00 3必)如表所示。 用最小二乘法求拟合数据的二次多项式，并求平方误差。 (3)实验内容 l.i = 0:10; x = -l+0.2*i; y= l./(l+25*x.A2); p=polyfit(x,y,3)； s=vpa(poly2sym(p)) f = polyval(p,x); plot( x, f, x, y, fo *) 0.000000000000000018136579684810800950071230066733*xT - 0.57518273581808898597955703735352*xA2 ? 0.00000000000000021309598658452197357889078781603*x + 0.4841249248489066792089374757力93 Q Figure 1 II 口 I 回 2. :(i)三次拟合，程序如下 sym x; xO二[0.0 0.1 0.2 0.3 0.5 0.8 1.0]; y0=[1.0 0.41 0.50 0.61 0.91 2.02 2.46]; cc=polyfit(xO,yO,3); S3=cc(l)+cc(2)*x+cc(3)*xA2+cc(4)*xA3 %三次拟合多项式 xx=xO ⑴:O.l:xO(length(xO)); yy=polyval(cc,xx); plotfx^yy；-); hold on; plot(xO,yO/x*); xlabel(x); ylabel(Y); 运行结果 S3 = -7455778416425083/1125899906842624+1803512222945437/140737488355328*x-655705 280524945/140737488355328*xA2+4172976892199509/4503599627370496*xA3 图像如下 ⑵4次多项式拟合 sym x; x0=[0.0 0.1 0.2 0.3 0.5 0.8 1.0]; y0=[1.0 0.41 0.50 0.61 0.91 2.02 2.46]; cc=polyfit(x0,y0,4); S3=cc(l)+cc(2)*x+cc(3)*xA2+cc(4)*xA3+cc(5)*xA4 xx=x0(l):0.1:x0(length(x0)); yy=polyval(cc,xx); plot(xxzyy；r); hold on; plot(xOzyOzx); xlabel(x); ylabel(y); 运行结果 S3 = 3248542900396215/1125899906842624-3471944732519151/281474976710656*x+4580931 990070637/281474976710656*xA2-5965836931688425/1125899906842624*xA364650307/9007199254740992*xA4 图像如下 (3)另一个拟合曲线， 新建一个M-file 程序如下： function [C,L]=lagran(x,y) w=len gth(x); n=w-l; L=zeros(w,w); for k=l:n+l V=l; for j=l:n+l if k~=j V=conv(V,poly(x(j)))/(x(k)-x(j)); end end L(k,：)=V; end C=y*L 在命令窗口屮输入以下的命令： x=[0.0 0.1 0.2 0.3 0.5 0.8 1.0]; y=[1.0 0.41 0.50 0.61 0.91 2.02 2.46]; cc=polyfit(x,y,4); xx=x(l):0.1:x(length(x)); yy=polyval(cc,xx); plot(xx,yy,T); hold on; plot(x,y；x); xlabel(x);