第六讲们 数值计算.pptVIP

第六讲 数值计算;MATLAB提供大量具有强大数值计算功能的函数。 本章着重介绍关于数值计算的函数。;目录;6.1 多项式运算; MATLAB提供了关于多项式的函数： 多项式的值； 多项式的根和微分； 多项式拟合曲线； ;多项式函数;1．多项式表示法;;2．多项式求值;;3．多项式乘法和除法;4．多项式的微积分; [q,d] = polyder(b,a)，返回多项式b/a微分的系数向量。 ;;（2）多项式的积分;5．多项式的根和由根创建多项式;（2）由根创建多项式 函数poly()实现由根创建多项式，其具体使用方法如下： p = poly(r)，输入r是多项式所有根，返回值为多项式的系数向量； p = poly(A)，输入A是方阵，返回值为A的特征多项式的系数向量。;6．多项式曲线拟合;;运行结果如下图所示。 ;6.2 插 值 运 算; 插值是根据已知输入/输出数据集和当前输入估计输出值。MATLAB提供大量的插值函数，如下表所示。 ;6.2.1 一维插值; 函数interp1()实现一维插值，其具体使用方法如下： yi=interp1(x,y,xi)，x，y是已知数据集且具有相同长度的向量； yi = interp1(y,xi)，默认x为1:n，其中n为向量y的长度； yi=interp1(x,y,xi,method)。method用于指定插值的方法。 ;;运行结果如下图所示。 ;6.2.2 二维插值; 函数interp2()实现二维插值,其具体使用方法如下： zi = interp2(x,y,z,xi,yi)，x,y,z为原始数据，返回值zi是插值结果； zi = interp2(z,xi,yi)，若z=n×m，则x=1:n，y=1:m； zi = interp2(x,y,z,xi,yi,method)，method用于指定插值的方法 。 ;;运行结果如下图所示。 ;6.3 有限差分;;运行结果如下图所示。;6.4 功 能 函 数; 函数可以通过以下方式来表示： M文件； 匿名函数； 函数inline()。 ;;2．函数画图; 以函数fplot()为例介绍画图函数的用法，其具体使用方法如下： fplot(function,limits)，function为待画图的函数，limits是横坐标数值范围或横纵坐标数值范围；; fplot(function,limits,LineSpec)，LineSpec指定画图的线条属性； fplot(function,limits,tol)，tol指定画图相对精度； fplot(function,limits,tol,LineSpec)，指定画图的线条属性和画图相对精度。;;运行结果如下图所示。;3．函数最小值和零点;（1）求一元函数最小值; x = fminbnd(fun,x1,x2,options)，使用options选项来指定的优化器的参数； [x,fval] = fminbnd(...)，附加返回函数最小值。 ;（2）求多元函数的最小值; [x,fval] = fminsearch(...)，附加返回函数最小值。;（3）求一元函数的零点; x = fzero(fun,x0,options)，用options指定寻找零点的优化器参数； [x,fval] = fzero(...)，附加自变量为x时的函数值。 ;4．数值积分;（1）一元函数的数值积分; q = quad(fun,a,b,tol)，以绝对误差容限tol计算函数fun在[a b]区间内的定积分； q = quad(fun,a,b,tol,trace)，当trace为非零值时，显示迭代过程的中间值。;（2）矢量数值积分;（3）二重积分; q = dblquad(fun,xmin,xmax,ymin,ymax,tol,method)，用method指定计算一维积分时采用的函数。;6.5 微分方程组数值解; 在MATLAB中，可以计算微分方程数值解，如： 常微分方程组的初值问题； 延迟微分方程的问题； 常微分方程组的边界问题。;6.5.1 常微分方程组的初值问题 ; 在MATLAB中可以计算以下初值问题的数值解。 显式常微分方程组； 线性隐式常微分方程组； 完全隐式常微分方程组。 ;1．显式常微分方程组;常微分方程组解法对比 ;2．设置解法器参数; options= odeset(oldopts,‘name1’, value1,…)，修改原来的解法器options结构体oldopts；