第5章数值计算.ppt

最小二乘法曲线拟合 利用lsqcurvefit函数进行任意表达式曲线拟合 原理： min sum {(fun(x,xdata)-ydata).^2} 用法： x=lsqcurvefit(fun,x0,xdata,ydata) x=lsqcurvefit(fun,x0,xdata,ydata,lb,ub) x=lsqcurvefit(fun,x0,xdata,ydata,lb,ub,options) [x,resnorm]=lsqcurvefit(fun,x0,xdata,ydata,...) [x,resnorm,residual,exitflag]=lsqcurvefit(fun,x0,xdata,ydata,...) 实 例 xdata = 0:0.1:10; ydata = 3*sin(xdata)+6; fun = inline(x(1)*sin(xdata)+x(2),x,xdata); x = lsqcurvefit(fun,[0 0], xdata, ydata) yfit=feval(fun,x,xdata); plot(xdata,ydata,ro,xdata,yfit) legend(ydata, yfit) 实 例 Sw Kro Krw 0.290 1 0.000 0.356 0.63 0.017 0.449 0.35 0.069 0.496 0.27 0.113 0.543 0.2 0.165 0.590 0.133 0.227 0.684 0.05 0.372 0.730 0.01 0.456 0.780 0 0.546 实 例 load xiangshen.dat Swg=xiangshen(:,1); Krog=xiangshen(:,2); Krwg=xiangshen(:,3); plot(Swg,Krog,r-o,Swg,Krwg,b-o) legend(Kro,Krw) miuo=9; rouo=0.9; miuw=0.45; rouw=1; load water Sw=WAT; Kro=interp1(Swg,Krog,Sw); Krw=interp1(Swg,Krwg,Sw); fw=Krw*rouw/miuw./(Krw*rouw/miuw+Kro*rouo/miuo); figure ,surf(Sw),colorbar,view(2),axis equal,axis tight figure,surf(fw),colorbar,view(2),axis equal,axis tight 实 例 X=magic(4) X = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 [Y,I]=sort(X) Y = 4 2 3 1 5 7 6 8 9 11 10 12 16 14 15 13 I = 4 1 1 4 2 3 3 2 3 2 2 3 1 4 4 1 结果解释： 对矩阵X的每列按升序排列得到新矩阵Y，并返回序数矩阵。 bvp4c 函数的用法如下： sol = bvp4c(odefun,bcfun,solinit)，odefun代表常微分方程组的函数，bcfun是描述边界条件的函数，solinit是对方程解的猜测解； sol = bvp4c(odefun,bcfun,solinit,options)，使用options结构体来设定解法器的参数。 运行结果如下图所示。 习 题 optimget函数得到目前优化器的参数，其具体使用方法如下： val = optimget(options,param)，返回优化器参数param的值； val= optimget(options,param,default)，返回优化器参数param的值。 4．数值积分 MATLAB提供一些的数值积分函数，如 下表所示。 数值积分函数 （1）一