MATLAB第4讲-MATLAB 数值计算.ppt

数值线性代数问题 4.3 函数和数值积分 一、函数的绘图及分析 1、绘制函数曲线 格式： fplot(‘函数名’,lims，’s’) 功能： 绘制指定函数的曲线， lims为x,y 轴的最小最大值, s可指定线形 ■函数和数值积分库(funfun) ■特殊函数库（specfun） 例题 4.3 函数和数值积分 2、函数的简易绘图 格式：ezplot(‘函数名’,lims) 例如：ezplot(‘sin’) ezplot(x^2 - x - 1) ; ezplot(sin(t),cos(t)); (画一椭圆) 除 ezplot外，还有： ezpolar, ezplot3， ezmesh , ezsurf等。 例题 4.3 函数和数值积分 3、求函数极小值 格式：fmin(‘函数名’,x1,x2) 求函数在x=[x1 x2]之间的极小值 4、求函数零点 格式：fzero(‘函数名’,x0)。 求函数在x0附近的过零点 。 例如：fzero(‘humps’,1)=1.2995 4.3 函数和数值积分 二、特殊函数 特殊的数学函数， 如：贝塞尔函数、误差函数等。 4.3 函数和数值积分 三、函数的数值积分 1、定积分（一维数值积分） 格式：quad(‘函数名’,x1,x2) （或 quad8 ，高阶方法） 对函数在区间[x1,x2]内的定积分 例如： quad(‘humps’,0,1)=29.8583 利用定积分可以求不定积分的数值解 数值积分 [y, n]=quad( F, a, b, tol); [y, n]=quad8(F, a, b, tol); [y, n]=quadl(F, a, b, tol); 其中F为被积函数， a, b为积分上下限， tol 为误差限 y 为积分结果，n 为被积分函数的调用次数。 tol 的缺省值对 quad 函数为 1e-3 对 quad8 和 quadl 函数为 1e-6 quad8 和 quadl 函数使用的算法比 quad 高得多 例子：求无穷定积分 f=inline(1/sqrt(2*pi)*exp(-x.^2/2),x); [y,kk]=quad(f,-8,8) y = 1.00000197533430 kk = 81 [y1,kk1]=quad8(f,-8,8) y1 = 1.00000000000023 kk1 = 161 [y,kk]=quad(f,-15,15) y = 0.99999920879563 kk = 89 [y1,kk1]=quad8(f,-15,15) y1 = 0.99999999999999 kk1 = 769 该无穷定积分的理论值为 1 4.3 函数和数值积分 三、函数的数值积分 2、二重数值积分 格式： dblquad(‘函数名’,x1,x2,y1,y2) 例子：双重积分 用函数 dblquad(F, xmin, xmax, ymin, ymax, tol) 求双重积分 f=inline(exp(-x.^2/2).*sin(x.^2+y),x,y); y=dblquad(f,-2,2,-1,1) y = 1.57450259485169 4.3 函数和数值积分 四、常微分方程的求解 (ODE ,Ordinary Differential Equation ) 格式：[x,y]=ode23(‘函数名’,x0,xn,y0) (2,3阶) 或 ode45 (4,5阶) 函数名为微分方程名（含x,y） x0,xn为x区间，y0为y初值。 4.3 函数和数值积分 四、常微分方程的求解 例如：求微分方程： dy/dx=-3y+2x, y(1)=2 区间为[1,3]的解 首先建立函数：f=myde(x,y) (myde.m) f= -3*y+2*x [x,y]=