插值与拟合精要.pptVIP

* 零件加工问题 已知某机翼断面的下轮廓线上的部分数据，见下表，现在需要得到x每改变0.1时y的坐标，试完成加工所需数据，画出曲线，并求13≤x≤15范围内y的最小值。 x 0 3 5 7 9 11 12 13 14 15 y 0 1.2 1.7 2.0 2.1 2.0 1.8 1.2 1.0 1.6 * 5 曲线拟合 最小二乘法的基本思想：构造函数 使得 达到最小，其中 为待定系数， 称为基函数。 令 ，可求得 。 * 1. 通过机理分析建立数学模型来确定 f(x)（机理模型）； + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + f=a1+a2x 2. 将数据 (xi,yi) i=1, …,n 作图，通过直观判断确定 f(x)（经验模型）： f=a1+a2x+a3x2 f=a1+a2x+a3x2 f=a1+a2/x f=a1exp(a2x) f=a1exp(a2x) * 数据拟合的MATLAB实现 多项式拟合： a =polyfit(x,y,n） 多项式的最高阶数 所拟合的数据向量 拟合多项式的系数向量（降幂排列） * 例：已知一组温度t和电阻R的数据见下表，拟合电阻与温度之间的关系R=at+b，并预测60℃时的电阻有多大。 t 20.5 32.7 51.0 73.0 95.7 R 765 826 873 942 1032 （1）根据数据进行拟合 t=[20.5,32.7,51.0,73.0,95.7];r=[765,826,873,942,1032]; p=polyfit(t,r,1);a=p(1),b=p(2), 运行得 a = 3.3987，b = 702.0968 （2）作出拟合曲线，可检验拟合效果 y=polyval(p,t); plot(t,r,+,t,y,r),xlabel(t),ylabel(R) （3）求出60℃的电阻 y=polyval(p,60) 运行得 y =906.0212 * 例：已知实验数据，见下表，求x与y的拟合曲线。 x 7.5 6.8 5.1 4.5 3.6 3.4 y 359.26 165.60 59.17 41.66 25.92 22.37 x 2.6 2.5 2.1 1.5 2.7 3.6 y 13.47 12.87 11.87 6.69 14.87 24.22 （1）画出数据点 x=[7.5,6.8,5.1,4.5,3.6,3.4,2.6,2.5,2.1,1.5,2.7,3.6]; y=[359.26,165.60,59.17,41.66,25.92,22.37,13.47,12.87, 11.87,6.69,14.87,24.22]; plot(x,y,r*) * （2）将非线性拟合转化为线性拟合 Y=log(y);p=polyfit(x,Y,1);B=p(1);A=p(2);b=B,a=exp(A) 运行得 b = 0.6241，a = 2.7039 （3）检验拟合效果 X=8:-0.01:1;Y=a*exp(b.*X); plot(x,y,r*,X,Y,b) legend(数据点,拟合曲线) * 非线性最小二乘拟合函数： x =lsqcurvefit(‘fun’,x0,x1,y1） 拟合函数中待定参数的初值，可任意取值 所拟合的函数 拟合的数据向量 * x =lsqnonlin(‘fun’,x0,x1,y1） 拟合函数中待定参数的初值，可任意取值 所拟合的函数 拟合的数据向量 * 例：用下表的数据拟合c(t)=a+be﹣0.02kt中的参数a,b,k。 t 100 200 300 400 500 600 700 800 900 1000 c 4.54 4.99 5.35 5.65 5.90 6.10 6.26 6.39 6.50 6.59 （1）编写M文件 fun1.m function f=fun1(x,t) f=x(1)+x(2)*exp(-0.02*x(3)*t)%其中 x(1)=a; x(2)=b；x(3)=k; （2）输入命令 t=100:100:1000 c=1e-03*[4.54,4.99,5.35,5.65,5.90,6.10, 6.26,6.39,6.50,6.59]; x0=[0.2,0.05,0.05]; x=lsqcurvefit (fun1,x0,t,c) f= fun1(x,t) * （3）运算结果为： f =0.0