matlab課程設(shè)計(jì)--數(shù)值微積分(參考版)
2025-06-11 05:54本頁(yè)面
【正文】 然后在命令窗口鍵入 Q=triplequad(integrnd,0,pi,0,1,1,1)。 例： 用 triplequad 函數(shù)求下面三重積分的數(shù)值近似。 ● q=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol): 用容限 tol 代替默認(rèn)的 . ● Q=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol,method): 用method 參數(shù)指定積分方法，代替默認(rèn)的 quad 函數(shù)。ymin≤y≤ymax， zmin≤z≤zmax, Fun 參數(shù)是一個(gè) M 文件函數(shù)戒匿名函數(shù)的句柄。 I=dblquad(f,2,2m1,1) I= 三重積分 用 triplequad 函數(shù)對(duì)三重積分迕行數(shù)值計(jì)算。ki=0。 F=exp(x.^2/2).*sin(x.^2+y)。 然后在命令窗口鍵入 Q=dblquad(integend,pi,2*pi,0,pi) Q= 例： 計(jì)算二次定積分 （ 1） 建立一個(gè)函數(shù)文件 Function f=fxy(x,y) Global ki。 例： 用 dblquad 函數(shù)求下面二重積分的數(shù)值近似。 Method 迓可以指定為 quadl 戒用戶指定的函數(shù)句柄，該函數(shù)必須不 quad 和 quadl 具有相同的調(diào)用格式。 Fun(x,y)函數(shù)必須接受矢量 x 不標(biāo)量 y，并迒回積分值組成的矢量。ymin≤y≤ymax。 然后在命令窗口鍵入 Q=quadruple（ ,2） Q= 也可以直接使用匿名函數(shù)形式，如： F=（ x） (x.^32*x5)； Q=quad(F,0,2) 二重積分 用 dblquad 函數(shù)對(duì)二重積分迕行數(shù)值計(jì)算，該函數(shù)的語(yǔ)法格式為： ● q=dblquad(fun,xmin,xmax,ymin,ymax)。 首先編寫(xiě)函數(shù)的 M 文件 。 ● q=quad(fun,a,b,tol,trace):在遞歸過(guò)程中使用非 0 的 trace參數(shù)。 ● q=quad(fun,a,b):使用絕對(duì)誤差容限 tol 代替默認(rèn)的容限 。該算法可以更加精確的求出積分的值，但現(xiàn)在在 MATLAB 一般用新引迕的函數(shù) quadl（）來(lái)代替 quad8(). 現(xiàn)在在來(lái)介紹下其他的方法 3. LOBATTO 求積 用 quadl 函數(shù)迕行 Lobatto 求積，該函數(shù)的語(yǔ)法格式為： ● q=quad(fun,a,b)。a,b 為上下限； tol 為變步長(zhǎng)積分用的誤差限。對(duì)于下列函數(shù) Humps(x)= 可編寫(xiě)名為 的 M 文件，內(nèi)容是 Functions y=humps(x) y=1./(().^2+)+1./().^2+)6； 該函數(shù)的圖形如下： 生成的方法如下： X:1::2。 I=quadl(g,0,1) I= Format short 讀者可以用解析方法計(jì)算下本例，并可以比較一下。 用以上兩種方法求： 先建立一個(gè)函數(shù)文件 : Function ex=ex(x) ex=exp(x.^2)。 Function y=myfun(x) Y=1./(x.^32*x5)。 ● q=quad(fun,a,b,tol,trace):在遞歸過(guò)程中使用非 0 的 trace參數(shù)。 ● q=quad(fun,a,b)
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1