数学实验4-2.pptVIP

4.3 函数的极限与连续性 求极限的格式 limit(f,x,a) %表示x→a时f的极限 limit(f,a) %同上 limit(f) %表示x→0时f的极限 limit(f,x,a,’right’) %x→a+0时f的右极限 limit(f,x,a,’left’) %x→a-0时f的左极限 例 计算 syms x h a Y1=limit(sin(x)/x) Y2=limit(sin(x)/x,inf) Y3=limit((x-2)/(x^2-4),2) Y4=limit(1/x,x,0,’left’) Y5=limit(1/x,x,0,’right’) Y6=limit((sin(x+h)-sin(x))/h,h,0) Y7=limit((1+a/x)sin(x),x,a) 3 连续性 例34 判断函数f(x)在x0点的连续性 通用程序 lxx.m Syms x; F=input(‘f=’) X0=input(‘x0=‘) F1=limit(f,x,x0) X=x0 F2=eval(f); If f1==f2 A=‘函数在x0连续’ Else A=‘函数在x0不连续’ end 4.4 导数与微分 3．求微分 算法：dy=ydx： 程序如下： syms x dx f=input(f=) fl=diff(f,x) df=fl*dx 4 求平面曲线在x0点的切线及法线方程 已知曲线参数方程x=x(t)，y=y(t)对应x0的参数值为t0,则切线与法线： 通用程序如下：编辑qxfx.m文件： syms x y t x0=input(x=); y0=input(y=); yd=diff(y0,t) t=input(t0=) yl=eval(y0)+eval(yd)*(x-eval(x0)) y2=eval(y0)-1/eval(yd)*(x-eval(x0)) 例36 求曲线x=cost,y=sint在t0=pi/4的切线及法线方程． 程序： qxfx x=cos(t) y=sin(t) t0=pi/4 4.5 泰勒展开 1．命令格式两种 (1)taylor(f) ％默认在x=0点展开6项 (2) taylor(f,n,x0)．％在x=x0点展开n项 例37 将e^x在x=0点展开5项；再在x=3点展开5项 程序： syms x f=exp(x) yl=taylor(f,x,5) y2=taylor(f,x,5,3) 例38 在x0=0点展开sin(x)成7项． 程序： syms x y t f1=taylor(sin(x),x,7) 例39 在x0=1点展开sin(x)成6项． 程序： syms x f2=taylor(sin(x),x,6,1) 例40 在x0=1点展开log(X)成5项． 程序如下： syms x f3=taylor(log(x),x,5,1) 例41 在x0=0点展开，f=exp(x)+2*cos(x)成4项． 程序如下： syms x f=exp(x)+2*cos(x) f4=taylor(f,x,4) 4.6 求一元函数最小值 1．命令格式 x=fminbnd(f，a，b) x=flminbnd(f，a，b,options) x=fminbnd(f，a，b，options，p1，p2) [x，fval]=fminbnd(…) [x，fval，exitflag]=fminbnd(…) [x，fval，exitflag，output]=frninbnd(…) optnew = optimset(options,TolX,1e-44); [x,fv,ef,out] = fzero(@sin,1,optnew) [x,fv,ef,out] = fminbnd(@sin,1,optnew) 例42 求f=x^3-x^2-x+1在(一2，2)内的极小值与极大值． 程序如下： syms x; f=x^3-x^2-x+1; [xl,minf]=fminbnd(f,-2,2) [x2,maxf]=fminbnd(-x^3+x^2+x-1,-2,2) maxf=-maxf 例43 求y=2*e^(-x)*sinx的极值 程序： syms x; f=2*exp(-x)*sin(x); ezplot(f,[0,8]) [x1,minf]=fminbnd(2*exp(-x)*sin(