实验四 插值问题离散数据的拟合

《实验四 插值问题离散数据的拟合》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、计算方法实验四插值问题实验五离散数据的曲线拟合班级：自动控制093班宿舍：八公寓323宿舍姓名：庄倩倩200908801胡瑞200908803刘银梅200908809李艳芳200908812丁春平200908816姚小娇2009088192011年6月14日实验四插值问题一、拉格朗日插值实验程序：functionf=agui_lagrange(x0,y0,x)n=length(x0);m=length(x);formatlongs=0.0;fork=1:np=1.0;forj=1:nifj~=kp=p*(

2、x-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endf=s;end实验结果：二、牛顿插值实验程序：functionf=aguiNewton(x,y,x0)symst;if(length(x)==length(y))n=length(x);c(1:n)=0.0;elsedisp('xºÍyµÄάÊý²»ÏàµÈ£¡');return;endf=y(1);y1=0;l=1;for(i=1:n-1)for(j=i+1:n)y1(j)=(y(j)-y(i))/(x(j)-x(

3、i));endc(i)=y1(i+1);l=l*(t-x(i));f=f+c(i)*l;simplify(f);y=y1;if(i==n-1)if(nargin==3)f=subs(f,'t',x0);elsef=collect(f);f=vpa(f,6);endendend实验结果：实验五离散数据的曲线拟合实验程序：functionf=agui_fit(x,y,m)A=zeros(m+1,m+1);fori=0:mforj=0:mA(i+1,j+1)=sum(x.^(i+j));endb(i+1)=su

4、m(x.^i.*y);endc=Ab';f=c';实验结果：