如果您无法下载资料，请参考说明：
1、部分资料下载需要金币，请确保您的账户上有足够的金币
2、已购买过的文档，再次下载不重复扣费
3、资料包下载后请先用软件解压，在使用对应软件打开
程序：x=-1:0.2:1;f=1./(1+25*x.^2);y=polyval(f,x);%计算原函数每个x值所对应的函数值p=polyfit(x,y,3)%对（x，y）进行三次拟合并输出三次多项式系数z=polyval(p,x);%计算拟合后的多项式对应x点的函数值plot(x,y,'r',x,z,'b')%画图输出：p=1.82221.2090-0.3619-0.1140则三次曲线拟合的方程为：2、3次和4次多项式拟合：程序：x=[0.00.10.20.30.50.81.0];y=[1.00.410.500.610.912.022.46];%3次拟合%p3=polyfit(x,y,3)xi=0:0.1:1.0;yi=polyval(p3,xi);subplot(1,2,1);plot(x,y,'*',xi,yi,'r');xlabel('x');ylabel('y');%4次拟合%p4=polyfit(x,y,4)xi=0:0.1:1.0;yi=polyval(p4,xi);subplot(1,2,2);plot(x,y,'*',xi,yi,'r');xlabel('x');ylabel('y');输出：p3=-6.622112.8147-4.65910.9266p4=2.8853-12.334816.2747-5.29870.9427则3次拟合多项式为：4次拟合多项式为：另一函数拟合：定义函数：function[C,L]=lagran(x,y)%x为n个节点的横坐标所组成的向量，y为纵坐标所组成的向%C为所求的牛顿插值多项式的系数构成w=length(x);n=w-1;L=zeros(w,w);fork=1:n+1V=1;forj=1:n+1ifk~=jV=conv(V,poly(x(j)))/(x(k)-x(j));%conv求积，poly(x)将该多项式的系数赋给向量endendL(k,:)=V;endC=y*L输入命令：x=[0.00.10.20.30.50.81.0];y=[1.00.410.500.610.912.022.46];cc=polyfit(x,y,4);xx=x(1):0.1:x(length(x));yy=polyval(cc,xx);plot(xx,yy,'r');holdon;plot(x,y,'x');xlabel('x');ylabel('y');x=[0.00.10.20.30.50.81.0];y=[0.940.580.470.521.002.002.46];%y中的值是根据上面两种拟合曲线而得到的估计数据，不是真实数据[C,L]=lagran(x,y);xx=0:0.01:1.0;yy=polyval(C,xx);holdonplot(xx,yy,'b',x,y,'.');输出：C=40.6746-110.2183110.3671-57.326423.4994-5.47640.9400
立即下载
