数学实验课件4..docVIP

（1）diff(f) (2) diff(f,v) (3 diff(f,n) ％关于符号变量求n阶导数 (4)diff(f,v,n) ％关于变量v对f求n阶导数 2．举 例 例35 已知函数 程序如下： syms x a b c f=a*x^3+x^2-b*x-c yl=diff(f) y2=diff(f,a) y3=diff(f,2) y4=diff(f,a,2) 3．求微分 算法：dy=ydx： 程序如下： syms x dx f=input(f=) fl=diff(f,x) df=fl*dx 4 求平面曲线在x0点的切线及法线方程 算法： 已知曲线参数方程，对应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=2的切线及法线方程． 程序： x=cos(t) y=sin(t) t0=pi/4 4.5 泰勒展开 1．命令格式两种 (1)taylor(f) ％默认在x=0点展开6项 (2) taylor(f,n，x0)． ％在x=xO点展开n项 例37 将 在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，p l，p2) [x，fval]=fminbnd(…) [x，fval，exitflag]=fminbnd(…) [x，fval，exitflag，output]--frninbnd(…) 2．举 例 例42求f=在(一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 求的极值 程序： syms x; f=2*exp(-x)*sin(x); ezplot(f,[0,8]) [x1,minf]=fminbnd(2*exp(-x)*sin(x),2,5) [x2,maxf]=fminbnd(-2*exp(-x)*sin(x),0,2) maxf=-maxf 例44 对边长为4 m的正方形铁皮，在四个角剪去相等的小正方形以制成无盖水槽，问如何 剪裁可使水槽容积最大? 设x为减去的小正方形边长，v为水槽容积． 程序如下： syms x f=-(4-2*x)^2*x [x,maxf]=fminbnd(f,0,2) maxf=-maxf 4.7 一元函数积分 1 积分格式 int(f) ％对f关于符号变量求不定积分：f(x)dx int(f,v) ％对f关于变量v，求不定积分： int(f,a,b) ％对f关于符号变量从a到b求定积分：f(x)dx int(f,v,a,b) ％对f关于变量v从a到b求定积分：f(v)dv 2．举例 例45 求f(x)=axsinx 的不定积分 程序如下： syms x a f=a*x