2-2 拟合.pptVIP

拟 合;拟 合 问 题 引 例 1;拟 合 问 题 引 例 2;曲 线 拟 合 问 题 的 提 法;拟合与插值的关系;最临近插值、线性插值、样条插值与曲线拟合结果：;曲线拟合问题最常用的解法——线性最小二乘法的基本思路;线性最小二乘法的求解：预备知识;线性最小二乘法的求解;线性最小二乘拟合 f(x)=a1r1(x)+ …+amrm(x)中函数{r1(x), …rm(x)}的选取 ;用MATLAB解拟合问题;用MATLAB作线性最小二乘拟合;即要求 出二次多项式:;1）输入以下命令： x=0:0.1:1; y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; R=[(x.^2) x ones(11,1)]; A=R\y;1）输入以下命令： x=0:0.1:1; y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; A=polyfit(x,y,2) z=polyval(A,x); plot(x,y,k+,x,z,r) %作出数据点和拟合曲线的图形;1. lsqcurvefit 已知数据点： xdata=（xdata1，xdata2，…，xdatan）， ydata=（ydata1，ydata2，…，ydatan） ; 输入格式为: x = lsqcurvefit (‘fun’,x0,xdata,ydata); x =lsqcurvefit (‘fun’,x0,xdata,ydata,options); x = lsqcurvefit (‘fun’,x0,xdata,ydata,options,’grad’); [x, options] = lsqcurvefit (‘fun’,x0,xdata,ydata,…); [x, options,funval] = lsqcurvefit (‘fun’,x0,xdata,ydata,…); [x, options,funval, Jacob] = lsqcurvefit (‘fun’,x0,xdata,ydata,…);; lsqnonlin用以求含参量x（向量）的向量值函数 f(x)=(f1(x),f2(x),…,fn(x))T 中的参量x，使得 最小。 其中 fi（x）=f（x，xdatai，ydatai） =F(x,xdatai)-ydatai ;输入格式为： 1） x=lsqnonlin（‘fun’，x0）； 2） x= lsqnonlin （‘fun’，x0，options）； 3） x= lsqnonlin （‘fun’，x0，options，‘grad’）； 4） [x，options]= lsqnonlin （‘fun’，x0，…）； 5） [x，options，funval]= lsqnonlin （‘fun’，x0，…）；;; 1）编写M-文件 curvefun1.m function f=curvefun1(x,tdata) f=x(1)+x(2)*exp(-0.02*x(3)*tdata); %其中 x(1)=a; x(2)=b；x(3)=k;;f =0.0043 0.0051 0.0056 0.0059 0.0061 0.0062 0.0062 0.0063 0.0063 0.0063 x = 0.0063 -0.0034 0.2542; 解法 2 用命令lsqnonlin f(x)=F(x,tdata,ctada)= x=（a，b，k）;3）运算结果为 f =1.0e-003 *(0.2322 -0.1243 -0.2495 -0.2413 -0.1668 -0.0724 0.0241 0.1159 0.2030 0.2792 x =0.0063 -0.0034 0.2542;练习1： 分别用lsqcurvefit,lsqnonlin求解下列曲线拟合问题;练习2. 用给定的多项式，如y=x3-6x2+5x-3，产生一组数据(xi,yi，i=1,2,…,n),再在yi上添加