微分學(xué)基本定理

1、實(shí)驗(yàn)二 微分學(xué)基本定理實(shí)驗(yàn)的目的1、掌握用matlab驗(yàn)證微分基本定理的正確性；2、掌握利用matlab對(duì)函數(shù)的Taylor展開；3、利用函數(shù)的Taylor展開式做近似計(jì)算。實(shí)驗(yàn)的基本理論與方法1、羅爾(Rolle)定理：如果函數(shù)在閉區(qū)間上連續(xù)，在開區(qū)間可導(dǎo)，且在區(qū)間端點(diǎn)的函數(shù)值相等，即，那么在內(nèi)至少有一點(diǎn)，使得函數(shù)在該點(diǎn)的導(dǎo)數(shù)等于0：。2、拉格朗日(Lagrange)中值定理：如果函數(shù)在閉區(qū)間上連續(xù)，在開區(qū)間可導(dǎo)，那么在內(nèi)至少有一點(diǎn)，使等式成立。3、柯西(Cauchy)中值定理：如果函數(shù)及在閉區(qū)間上連續(xù)，在開區(qū)間內(nèi)可導(dǎo)，且在內(nèi)每一點(diǎn)處均不為零，那么在內(nèi)至少有一點(diǎn)，使等式成立。4、泰勒(Tay

2、lor)中值定理：若函數(shù)在含的某個(gè)區(qū)間內(nèi)具有直到階導(dǎo)數(shù)，則對(duì)任意，有，其中，(Lagrange余項(xiàng))，或（Peano余項(xiàng)）。5、常用函數(shù)的麥克勞林（Maclaurin）展開式：實(shí)驗(yàn)使用的函數(shù)與命令1、求導(dǎo)指令diff；2、繪圖指令plot；3、用符號(hào)運(yùn)算求解一般方程指令： solve(f)：解方程f；solve(f，x )：對(duì)變量x解方程f。4、泰勒(Taylor) 展式指令：taylor(f,n)：返回函數(shù)f的n-1階麥克勞林（Maclaurin）展開式，其中f為自變量v的符號(hào)函數(shù)；taylor(f,n,v,a) ：返回函數(shù)f在a點(diǎn)的n-1階泰勒(Taylor)展開式，其中f為自變量v的符號(hào)

3、函數(shù)，n, v,與 a 可以省略即taylor(f)指令， 此時(shí)a=0、n默認(rèn)為6、v為f中所確定的符號(hào)變量。 5、符號(hào)替換指令：subs(S,old,new) 用新的符號(hào)變量或數(shù)值變量替換表達(dá)式S中舊的符號(hào)變量或表示一個(gè)變量名的字符串。詳細(xì)的應(yīng)用可見Matlab help。實(shí)驗(yàn)指導(dǎo)例1 驗(yàn)證拉格朗日(Lagrange)中值定理對(duì)函數(shù)在區(qū)間上的正確性。解：為解題方便，首先求出函數(shù)的導(dǎo)數(shù)：syms x;y=4*x3-5*x2+x-2;dy=diff(y)輸出結(jié)果：dy=12*x2-10*x+1下面利用dy驗(yàn)證拉格朗日(Lagrange)中值定理，現(xiàn)取a=0，b=2進(jìn)行驗(yàn)證。a=0;ya=4*a3

4、-5*a2+a-2; %求出a=0時(shí)函數(shù)值b=2;yb=4*b3-5*b2+b-2; %求出b=2時(shí)函數(shù)值eq=12*x2-10*x+1-(yb-ya)/(b-a)=0;eq=dy-(yb-ya)/(b-a)=0;mm=solve(eq,x);%求符號(hào)解mm=subs(mm,a,b,yb,ya,-1,1,yb,ya)%代值,求數(shù)值解輸出結(jié)果：mm =-0.40411.2374從實(shí)驗(yàn)結(jié)果中可以觀察到mm =1.2374，從而驗(yàn)證拉格朗日(Lagrange)中值定理。例2 設(shè)函數(shù)，說明有幾個(gè)實(shí)根，指出他們所在的區(qū)間，并求出所有根。解：由羅爾(Rolle)定理易知有三個(gè)根，且分別在區(qū)間內(nèi)，函數(shù)與的圖

5、形如圖2.3所示。M文件程序設(shè)計(jì)如下x=0.5:0.1:4.5;y1=(x-1).*(x-2).*(x-3).*(x-4); hold on圖2.3 函數(shù)曲線plot(x,y1);%繪出原函數(shù)圖形syms x;y=(x-1)*(x-2)*(x-3)*(x-4);dy=diff(y);x=0.5:0.1:4.5;y2=subs(dy,x,x);plot(x,y2,-r); %繪出導(dǎo)函數(shù)圖形從圖中可以直觀地觀察出函數(shù)曲線與y軸相交情形，即可得到方程根的個(gè)數(shù)。下面求解方程的三個(gè)根。由上面的程序可知的導(dǎo)數(shù)為 dy dy = (x-2)*(x-3)*(x-4)+(x-1)*(x-3)*(x-4)+(x-

6、1)*(x-2)*(x-4)+(x-1)*(x-2)*(x-3)利用已求導(dǎo)數(shù)，則求的根M文件程序?yàn)閑q=(x-2)*(x-3)*(x-4)+(x-1)*(x-3)*(x-4)+(x-1)*(x-2)*(x-4)+(x-1)*(x-2)*(x-3)=0;solve(eq,x)%求符號(hào)解輸出結(jié)果ans = 5/2 5/2+1/2*5(1/2) 5/2-1/2*5(1/2)例3 求下列函數(shù)的麥克勞林（Maclaurin）展開式。1) ,；2) ,；3）,。解：利用指令taylor(f,n,v)可以迅速求解：syms xf=Log(1+x);T1=taylor(f,7)f=tan(x);T2=tayl

7、or(f,8)f=1/(5+4*cos(x);T3=taylor(f,8)輸出結(jié)果：T1 = x-1/2*x2+1/3*x3-1/4*x4+1/5*x5-1/6*x6 T2 = x+1/3*x3+2/15*x5+17/315*x7 T3 = 1/9+2/81*x2+5/1458*x4+49/131220*x6例4 求下列函數(shù)在指定點(diǎn)的泰勒(Taylor)展開式。1),；2) ,；解：syms xf=x4-5*x3+x2-3*x+4;T1=taylor(f,5,1.5)f=x2*Log(x);T2=taylor(f,4,1)輸出結(jié)果：T1 = 325/16-81/4*x-8*(x-3/2)2+(

8、x-3/2)3+(x-3/2)4 T2 = x-1+3/2*(x-1)2+1/3*(x-1)3例5 繪出正弦函數(shù)的階麥克勞林（Maclaurin）展開式的圖形，并觀察展開式圖形逼近正弦函數(shù)的變化。解：為更好觀察，我們把階麥克勞林（Maclaurin）展開式的圖形與正弦函數(shù)的圖形繪到同一窗口。如圖2.4所示，M文件程序如下圖2.4hold onx=0:0.1:3/2*pi;y=sin(x);plot(x,y,-r); %sin(x)圖形for n=1:10 syms x f=sin(x); T=taylor(f,n); x=0:0.1:(2*n)/(n+5)*pi; %系數(shù)(2*n)/(n+5)

9、為了控制圖形的顯示區(qū)域 y=subs(T,x,x); %繪出n=1.10的展開式函數(shù)圖形 plot(x,y);end從圖2.4中可以看到，當(dāng)?shù)闹抵饾u增大，函數(shù)的展開式圖形越逼近原函數(shù)圖形。實(shí)際上，matlab提供了一個(gè)圖形用戶界面，用戶可以方便處理Taylor展開式圖形。圖2.5例6 Taylor級(jí)數(shù)計(jì)算器的使用。函數(shù) taylortool格式 taylortool %該命令生成圖形用戶界面，顯示缺省函數(shù)f=x*cos(x)在區(qū)間-2*pi，2*pi內(nèi)的圖形，同時(shí)顯示函數(shù)f的前N=7項(xiàng)的Taylor多項(xiàng)式級(jí)數(shù)和(在a=0附近的)圖形通過更改f(x)項(xiàng)可得不同的函數(shù)圖形。taylortool(f) %對(duì)指定的函數(shù)f，用圖形用戶界面顯示出Taylor展開式。（圖2.5）例 taylortool( cos(x)再通過改變相關(guān)的參量，可得如圖2.6。圖2.6 函數(shù)cos(x)的taylortool界面實(shí)驗(yàn)內(nèi)容與練習(xí)1、驗(yàn)證柯西(Cauc