MATLAB基础教程-第6章-曲线拟合.ppt

科学实验中，常常需要从一组测量数据中找到实验规律的数学表达式，例如经验公式。第六章曲线拟合温度t(0C)20.532.751.073.095.7电阻R(?)7658268739421032已知热敏电阻数据：求60℃时的电阻R。设R=at+ba,b为待定系数第六章曲线拟合6.1线性最小二乘法线性最小二乘法是曲线拟合最常用的解法。其基本思路为：第六章曲线拟合第一步:先选定一组函数r1(x),r2(x),…rm(x),mn,令f(x)=a1r1(x)+a2r2(x)+…+amrm(x)（1）其中a1,a2,…am为待定系数。6.1线性最小二乘法第六章曲线拟合第二步：确定a1,a2,…am的准则（最小二乘准则）：使n个点（xi,yi)与曲线y=f(x)的距离?i的平方和最小。线性最小二乘法是曲线拟合最常用的解法。其基本思路为：6.1线性最小二乘法第六章曲线拟合记问题归结为，求a1,a2,…am使J(a1,a2,…am)最小。线性最小二乘法是曲线拟合最常用的解法。其基本思路为：第六章曲线拟合6.2线性最小二乘法的求解超定方程组：方程个数大于未知量个数的方程组超定方程一般是不存在解的矛盾方程组。如果有向量a使得达到最小，则称a为上述超定方程的最小二乘解。第六章曲线拟合6.2线性最小二乘法的求解在MATLAB中，执行a=R\y命令，a即为超定方程的最小二乘解所以，曲线拟合的最小二乘法要解决的问题，实际上就是求以下超定方程组的最小二乘解的问题。其中Ra=y（3）第六章曲线拟合6.3用MATLAB进行曲线拟合方法一：多项式f(x)=a1xm+…+amx+am+1拟合可用polyfit函数：方法二：建立超定方程组也可得到最小二乘意义下的解（a为向量），用a=polyfit(x,y,m)输出拟合多项式系数a=[a1,…am,am+1](数组)输入同维的数组X，Y拟合多项式次数第六章曲线拟合6.3用MATLAB进行曲线拟合即要求出二次多项式:中的使得:例6_1对下面一组数据作二次多项式拟合第六章曲线拟合6.3用MATLAB进行曲线拟合1）输入以下命令：x=0:0.1:1;y=[-0.4471.9783.286.167.087.347.669.569.489.3011.2];R=[(x.^2);x;ones(size(x))]’；A=R\y解法1．用解超定方程的方法2）计算结果：A=-9.810820.1293-0.0317第六章曲线拟合解法2．用polyfit函数6.3用MATLAB进行曲线拟合1）输入以下命令：x=0:0.1:1;y=[-0.4471.9783.286.167.087.347.669.569.489.3011.2];A=polyfit(x,y,2)z=polyval(A,x);plot(x,y,k+,x,z,r)%作出数据点和拟合曲线的图形2）计算结果：Ａ=-9.810820.1293-0.0317第六章曲线拟合例6_2用函数g(x)=c1+c2x+c3sin(πx)+c4sin(2πx)，拟合下面的数据，并作出图形6.4曲线拟合实例i123456789xi0.10.20.30.40.50.60.70.80.9yi0.00002.12203.02443.25683.13992.85792.51402.16391.8358第六章曲线拟合练习1、在某个低温过程中，函数y依赖于温度θ（℃）的实验数据如下：θ1234y0.81.51.82.0已知经验公式的形式为y=aθ+bθ2，试用最小二乘法求出a，b2、用形如aex+bsin(x)+cln(x)+dcos