实验5 方程(组)求根.pptVIP

[x,f,h]=fsolve(fun,-0.6) x = -0.6180 f = -3.3152e-012 h = 1 例10. 用fsolve求函数 y=xsin(x2 –x-1) 的零点 fun=inline(‘x*sin(x^2-x-1)’,’x’) [x,f,h]=fsolve(fun,-1.6) x = -1.5956 f = 1.4909e-009 h = 1 * 例11. 求下列方程组在原点附近的解 解. 先写一个M文件 function f=fun7(x) f(1)=4*x(1)-x(2)+exp(x(1))/10-1; f(2)=-x(1)+4*x(2)+x(1)^2/8; [x,f,h]=fsolve(@fun7,[0,0]) x =0.2326 0.0565 f=1.0e-006*0.0908 0.1798 h=1 * fun=inline(‘[4*x(1)-x(2)+exp(x(1))/10-1, -x(1)+4*x(2)+x(1)^2/8]’) ; [x,f,h]=fsolve(fun,[0,0]) x =0.2326 0.0565 f=1.0e-006*0.0908 0.1798 h=1 * 二 验证洛必达法则 * syms a b x; f=a^x-b^x;g=x;L=limit(f/g,x,0) 例12 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 * 三. 函数的单调性 例13 解： 求可导函数的单调区间与极值，就是求导函数的正负区间与正负区间的分界点。 先求导函数的零点，再画出函数图像，根据图像可直观地看出函数的单调区间与极值。 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] * * 四.Taylor展开 函数f(x)在a点的Taylor展开为： Matlab命令： taylor(f,n,a,x)表示自变量为x的函数f在a点展开为n-1阶的幂级数 * 例14 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) T1 =0 T2 =x T3 =x T4 =x - x^3/2 T5 =x - x^3/2 T6 =x^5/24 - x^3/2 + x T7 =x^5/24 - x^3/2 + x * 图形比较： 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’) * t2 t6 t4 可以看出阶数越高逼近程度越好 * Taylor级数逼近分析界面 在命令窗口运行taylortool,将引出逼近分析界面。 * 五 常微分方程求解 1 常微分方程（组）符号求解 dsolve(‘eqn’,var): eqn是常微分方程，var是变量，默认是t dsolve(‘eqn1’,’eqn2’, …, ‘eqnm’, var):%是m个方程，var是变量，默认是t * 例15 求下列常微分方程的通解 输入命令： 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’) * 例16 求下列常微分方程的通解 输入命令： s=dsolve(‘Dx=y’,Dy=-x’) y=s.y,x=s.x * 2 常微分方程特解 dsolve(‘eqn’,’condition1’,’…,’condition n’,var):%eqn是常微分方程，var是变量，默认是t,condittion 为条件。 例17 求下列常微分方程的特解 输入命令： y=dsolve(‘D2y=cos(2*x)-y’,’y(0)=1’,’Dy(0)=0’,’x’) simplify(y) * 练习： 2 用taylor命令展开以下函数的前8项，并做图比较。 * 3 用牛顿迭代法求方程xlnx=1的近似根 4 作出函数f(x)=x3e-2x-xsin(x2