[理学]matlab教程.ppt

注1：fzero只能求零点附近变号的根，试以1.1为初 值，用fzero和fsolve求解(x-1)2=0, 看看发生了什 么？ 注2：fzero和fsolve只能求实根，试用它们解 X2+x+1=0,看看发生了什么？ -----用solve求符号代数方程的根 syms a b c x; F=a*x^2+b*x+c; solve(F) ans = [ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))] [x,y] = solve(x^2 + x*y + y = 3,x^2 - 4*x + 3 = 0) x = [ 1 3] y = [ 1 -3/2] 二 验证洛必达法则 syms a b x; f=a^x-b^x;g=x;L=limit(f/g,x,0) 例11 df=diff(f,x);dg=diff(g,x);L1=limit(df/dg,x,0) L = log(a)-log(b) L1= log(a)-log(b) L=L1 三. 函数的单调性 例12 解： 求可导函数的单调区间与极值，就是求导函数的正负区间与正负区间的分界点。 先求导函数的零点，再画出函数图像，根据图像可直观地看出函数的单调区间与极值。 syms x; f=x^3-6*x^2+9*x+3; df=diff(f,x); s=solve(df) ezplot(f,[0,4]) s = [ 1] [ 3] 3 ) 3 ( 7 ) 1 ( , ) ( ) 3 , 1 ( ) ( ) , 3 ( ) 1 , ( min max = = ˉ - +￥ -￥ f f x f x f ，在区间 与 在区间 可建立一个名为 dandiao.m 文件，用来求函数的单调区间。 disp(输入函数(自变量为x)); syms x; f=input(f(x)=); df=diff(f); s=solve(df); s=eval(s) a=[]; for i=1:length(s); a(i)=s(i); end ezplot(f,[min(a)-1,max(a)+1]) 四.Taylor展开 函数f(x)在a点的Taylor展开为： Matlab命令： taylor(f,n,a,x)表示自变量为x的函数f在a点展开为n-1阶的幂级数 例13 f=xcosx的展开式 syms x f=x*cos(x) T1=taylor(f,1) T2=taylor(f,2) T3=taylor(f,3) T4=taylor(f,4) T5=taylor(f,5) T6=taylor(f,6) T7=taylor(f,7) 图形比较： x=-4:0.1:4; f=x.*cos(x); T2=x; T4=x-1/2*x.^3; T6=x-1/2*x.^3+1/24*x.^5; plot(x,f,x,T2,:,x,T4,*,x,T6,o) Taylor级数逼近分析界面 在命令窗口运行taylortool,将引出逼近分析界面。 五 常微分方程求解 1 常微分方程（组）符号求解 dsolve(‘eqn’,var) eqn是常微分方程,var是变量,默认是t. dsolve(‘eqn1’,‘eqn2’,…,‘eqnm’,var) 是m个方程,var是变量,默认是t. 例14 求下列常微分方程的通解 输入命令： Y1=dsolve(‘x*Dy*log(x)+y=a*x*(log(x)+1)’,‘x’) Y2=dsolve(‘D2y+2*Dy+5*y=sin(2*x)’,‘x’) Y3=dsolve(‘D3y+D2y-2*Dy=x*(exp(x)+4)’,‘x’) 例15 求下列常微分方程的通解 输入命令： s=dsolve(Dx=y, Dy=-x) y=s.y,x=s.x dx y dt dy x dt = = - 2 常微分方程特解 dsolve(eqn, condition1‘, …, condition n,var) eqn是常微分方程，var是变量，默认是t, condittion为条件。 例16 求下列常微分方程的特解 输入命令： y=dsolve(D2y=cos(2*x)-y,y(0)=1,Dy(0)=0,x) simplify(y) 练习： 2.用taylor命令展开以下函数的前8项，并做图比较 1 . 确定下列函数的单调区间 1). 7 18 6 2 2 3 - - - = x x x y ； 2 ) . ) 1 ln( 2 x