matlab積分

1、一.數值積分的實現方法1變步長辛普生法基于變步長辛普生法，MATLAB給出了quad函數來求定積分。該函數的調用格式為：[In]=quad(fnameabtoltrace)其中fname是被積函數名。a和b分別是定積分的下限和上限。tol用來控制積分精度，缺省時取tol=0.001。trace控制是否展現積分過程，若取非0則展現積分過程，取0則不展現，缺省時取trace=0。返回參數I即定積分值，n為被積函數的調用次數。例81求定積分。

2、(1)建立被積函數文件fesin.m。functionf=fesin(x)f=exp(0.5x).sin(xpi6)(2)調用數值積分函數quad求定積分。[Sn]=quad(fesin03pi)S=0.9008n=772牛頓－柯特斯法基于牛頓－柯特斯法，MATLAB給出了quad8函數來求定積分。該函數的調用格式為：[In]=quad8(fnameabtoltrace)其中參數的含義和quad函數相似，只是tol的缺省值取106。該函

3、數可以更精確地求出定積分的值，且一般情況下函數調用的步數明顯小于quad函數，從而保證能以更高的效率求出所需的定積分值。(1)被積函數文件fx.m。functionf=fx(x)f=x.sin(x).(1cos(x).cos(x))(2)調用函數quad8求定積分。I=quad8(fx0pi)I=2.4674分別用quad函數和quad8函數求定積分的近似值，并在相同的積分精度下，比較函數的調用次數。調用函數quad求定積分：fmatl

4、ongfx=inline(exp(x))[In]=quad(fx12.51e10)I=0.28579444254766n=65調用函數quad8求定積分：fmatlongfx=inline(exp(x))[In]=quad8(fx12.51e10)DV=diff(V)%計算V的一階差分例87用不同的方法求函數f(x)的數值導數，并在同一個坐標系中做出f(x)的圖像。程序如下：f=inline(sqrt(x.^32x.^2x12)(x5)

5、.^(16)5x2)g=inline((3x.^24x1).sqrt(x.^32x.^2x12)216.(x5).^(56)5)x=3:0.01:3p=polyfit(xf(x)5)%用5次多項式p擬合f(x)dp=polyder(p)%對擬合多項式p求導數dpdpx=polyval(dpx)%求dp在假設點的函數值dx=diff(f([x3.01]))0.01%直接對f(x)求數值導數gx=g(x)%求函數f的導函數g在假設點的導數p

6、lot(xdpxxdx.xgx)%作圖Matlab數值積分的一些做法Filedunder:未分類—franz@9:31pm今天是元宵節(jié)，突然來講Matlab的確很奇怪，連我自己也這樣感覺由于Matlab對于計算過程的細節(jié)要求定義詳細，往往讓人覺得使用不方便與同樣很流行的Mathematic相比，Matlab更象是程序，而不是象Mathematic那么直觀的寫作業(yè)Matlab同樣提供符號積分（不定積分）和數值積分（定積分）兩大功能符號積分