綜合實驗一定積分的近似計算.ppt

1、數(shù)學(xué)實驗,綜合實驗一 定積分的近似計算,問題背景和實驗?zāi)康?實驗二、定積分的近似計算,定積分計算的基本公式是牛頓萊布尼茲公式。但當被積函數(shù)的原函數(shù)不知道時，如何計算？這時就需要利用近似計算。特別是在許多實際應(yīng)用中，被積函數(shù)甚至沒有解析表達式，而是一條實驗記錄曲線，或一組離散的采樣值，此時只能用近似方法計算定積分。,本實驗主要研究定積分的三種近似計算算法：矩形法、梯形法和拋物線法。同時介紹 Matlab 計算定積分的相關(guān)函數(shù)。,矩形法,定積分的定義：,實驗二、定積分的近似計算,矩形法,n 充分大，x 充分小,定積分的近似：,左點法,右點法,中點法,點 可以任意選取，常見的取法有： 左端點 ，右端

2、點 和中點 。,步長,節(jié)點,右點法：,中點法：,左點法：,左點法、右點法和中點法,解：,矩形法舉例,a=0, b=1, n=100,例：用不同的矩形法計算下面的定積分 ( 取 n=100 )， 并比較這三種方法的相對誤差。,左點法：,右點法：,中點法：,(i = 0,1,2,.,100),理論值：,左點法相對誤差：,誤差分析,矩形法舉例,右點法相對誤差：,中點法相對誤差：,不同的方法有不同的計算精度,有沒有更好的近似計算定積分的方法 ?,定積分幾何意義,曲邊小梯形的面積可以由直邊小梯形的面積來近似,整個曲邊梯形的面積：,梯形法,如果我們 n 等分區(qū)間 a,b，即令：,則,=,梯形公式,梯形法,

3、梯形公式與中點公式有什么區(qū)別 ?,解：,=,例：用梯形法計算下面定積分 ( 取 n=100 )， 并計算相對誤差,梯形法舉例,a=0, b=1, n=100, f (x) = 1/( 1+x2 ),相對誤差：,2n 等分區(qū)間 a,b ，得,該直線用拋物線代替，計算精度是否會更好？,計算每個節(jié)點上的函數(shù)值：,拋物線法,在區(qū)間 x0, x2 上，用過以下三點,的拋物線來近似原函數(shù) f (x) 。,設(shè)過以上三點的拋物線方程為：,則在區(qū)間 x0, x2 上，有,y = x2 + x + = p1(x),拋物線法,同理可得：,相加即得：,拋物線法,整理后可得：,或,辛普森 (Simpson) 公式,拋物

4、線法公式,拋物線法,=,例：用拋物線法計算下面定積分 ( 取 n=100 )， 并計算相對誤差,解：,a=0, b=1, n=100, yi = f (xi) = 1/( 1+xi2 ),相對誤差：,拋物線法,梯形法：trapz,trapz(x,y)x 為分割點（節(jié)點）組成的向量，y 為被積函數(shù)在節(jié)點上的函數(shù)值組成的向量。,Matlab 近似計算定積分的相關(guān)函數(shù),Matlab 計算定積分函數(shù)介紹,前面的做法,例：用梯形法計算下面定積分 ( 取 n=100),解：,a=0, b=1, n=100, yi = f (xi) = 1/( 1+xi2 ), x=0:1/100:1; y=1./(1+x

5、.2); trapz(x, y),trapz函數(shù),trapz(x,1./(1+x.2),trapz 舉例,quad(f,a,b,tol) f = f(x) 為被積函數(shù)，a,b 為積分區(qū)間，tol 為計算精度,將自變量看成是向量,拋物線法：quad,拋物線法,解：, quad(1./(1+x.2),0,1), quad(1./(1+x.2),0,1,10e-10), quad(1./(1+x.2),0,1,10e-16),函數(shù)表達式一定要用 單引號 括起來！ 涉及的運算一定要用 數(shù)組運算！,例：用 quad 計算定積分：,quad 舉例,拋物線法計算二重積分： dblquad,dblquad(f

6、,a,b,c,d,tol),tol 為計算精度，若不指定，則缺省精度為 10-6,f(x,y) 可以由 inline 定義，或通過一個函數(shù)句柄傳遞,a,b 是第一積分變量的積分區(qū)間，c,d 是第二積分變量 的積分區(qū)間,按字母順序，大寫字母排在小寫字母的前面,二重積分的計算, f=inline(4*x*y+3*y2); I=dblquad(f, -1,1,0,2),f(x,y) 中關(guān)于第一自變量的運算是數(shù)組運算， 即把 x 看成是向量，y 看成是標量。 也可以全部采用數(shù)組運算,例2：計算二重積分, dblquad(inline(4*x*y+3*x2),-1,1,0,2), dblquad(inl

7、ine(4*x*y+3*x.2),-1,1,0,2),X,例1：計算二重積分,dblquad 舉例,例：計算二重積分, dblquad(x,y)4*x*y+3*x.2 , -1,1, 0, 2),指定 x、y 分別是第一和第二積分變量, dblquad(inline(4*x*y+3*x.2), -1,1, 0, 2),被積函數(shù) f (x,y) 的另一種定義方法：匿名函數(shù),dblquad 舉例,int(f,a,b) 計算 f 關(guān)于默認自變量 的定積分，積分區(qū)間為a,b。,int(f) 計算 f 關(guān)于默認自變量 的不定積分。,int(f,v,a,b) 計算函數(shù) f 關(guān)于自變量 v 的定積分，積分區(qū)

8、間為 a, b,int(f,v) 計算函數(shù) f 關(guān)于自變量 v 的不定積分,findsym(f,1),符號積分： int,int 符號積分, syms x y; f=y*sin(x); int(f,x) int(f,y) int(f) int(a+b),ans=-y*cos(x),ans=1/2*y2*sin(x),ans=-y*cos(x),ans=a*b+1/2*b2,例：指出下面各條語句的輸出結(jié)果,int 舉例,例：用 int 函數(shù)計算定積分：,解：, syms x; f=1/(1+x2); int(f,x,0,1), f=sym(1/(1+x2); int(f,x,0,1), int(1/(1+x2),x,0,1),或, int(1/(1+x2),0,1),或,或,int 舉例,double(a) 將 a 轉(zhuǎn)化為雙精度型，若 a 是字符，則取對應(yīng)的 ASCII 碼, a=3; double(a) double(a),例：,ans = 3,ans = 97,其它相關(guān)函數(shù), x=1:0.001:2; y=exp(x.(-2); trapz(x,y),梯形法：,拋物線法：, quad(exp(x.(-2),1,2,10e-10),符號積分法：, syms x int(exp(x(-2),x,1,2),例 1：用 Matlab 函數(shù)近