matlab課程設(shè)計--數(shù)值微積分-閱讀頁
2025-01-28 16:35本頁面
【正文】 8函數(shù)進行自適應(yīng)遞歸Simpson求積?！袂蠛瘮?shù)fun從a到b的積分近似，函數(shù)x=fun(x)應(yīng)該接受矢量變量x并返回結(jié)果矢量y。 q=quad(fun,a,b):?！瘛? Function y=myfun(x) Y=1./(x.^32*x5)。 用以上兩種方法求：然后在MATLAB命令窗口，輸入命令： Format long I=quad(‘ex’,0,1) I= I=quadl(‘ex’,0,1) I= 也可以不建立關(guān)于被積函數(shù)文件，而使用語句函數(shù)求解，命令如下： q=inline(‘exp(x.^2)’)。同樣。Polt(x,humps(x))現(xiàn)在需要求humps從0到1積分，可使用下面的命令： q=quad(‘humps’,0,1) q= 注意quad的一般調(diào)用格式為 [y,n]=quad(F,a,b,tol)其中F為描述被積分的字符變量，該函數(shù)的一般格式為y=F(x)。 MATLAB中常見的一元函數(shù)數(shù)值積分指令如下 表2 其中quad()與quad8()這兩者的調(diào)用格式完全一致，只不過在函數(shù)quad8（）中tol的默認值較低。 q=quad(fun,a,b)。●Tol的值更大時，需要的計算次數(shù)少，計算快。 q=quad(fun,a,b,tol,trace):在遞歸過程中使用非0的trace參數(shù)。 []=quad(…): 返回函數(shù)計算次數(shù) 例： 用quadl函數(shù)求下面積分的數(shù)值積分。 Function y=myfun(x) Y=1./(x.^32*x5)。 q=dblquad(fun,xmin,xmax,ymin,ymax)。ymin≤y≤ymax。Fun(x,y)函數(shù)必須接受矢量x與標量y，并返回積分值組成的矢量。 q=dblquad(fun,xmin,xmax,ymin,ymax,tol):.●Method還可以指定為quadl或用戶指定的函數(shù)句柄，該函數(shù)必須與quad和quadl具有相同的調(diào)用格式。使用MATLAB提供的dblquad函數(shù)就直接可以求出上述積分的數(shù)值解，該函數(shù)的調(diào)用格式為：I=dblquad(f,a,b,c,d,)該函數(shù)求f(x,y)在[a,b][c,d]區(qū)域上的二重積分。,代碼為： Function z=integrnd(x,y) z=y*sin(x)+x*cos(y)。 例： 計算二次定積分 （1）Ki=ki+I。(2)ki=0。 I=dblquad(f,2,2m1,1) I= 三重積分 用triplequad函數(shù)對三重積分進行數(shù)值計算。 q=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax): 調(diào)用quad函數(shù)對fun(x,y)函數(shù)進行三重積分運算，積分區(qū)間為xmian≤x≤xmax。Fun(x,y)函數(shù)必須接受矢量x與標量yz并返回積分值組成的矢量。 q=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol): .●Method還可以指定為quadl或用戶指定的函數(shù)句柄，該函數(shù)必須與quad和quadl具有相同的調(diào)用格式。 然后在命令窗口鍵入 Q=triplequad(integrnd,0,pi,0,1,1,1)
高考資料相關(guān)推薦
鄂ICP備17016276號-1