命令窗口輸入被積函數

2

1

2

0

t

t

e

dt

?

。

輸入應為：。

Step2

：執行函數。輸入形式為

mymulNewtonCotes(ft,a,b,m,n)

；

其中

ft

—被積函數，此體重，已經在第一步賦值；

a

—積分下限，本題中為

0

；

b

—積分上限，本題中為

1

；

m

—將區間

[a,b]

等分的子區間數量，本題可選為

10

；

n

—采用的

Newton-Cotes

公式的階數，必須滿足

n<8

，否則積分沒法保證穩

定性。

當

n=1

時，即為復化梯形公式；

n=2

時，即為復化復化辛普森公式。

所以，分別輸入

mymulNewtonCotes(ft,0,1,10,1)

和

mymulNewtonCotes(ft,0,1,10,2)

就可以

得到兩種方法的積分計算結果。

2

、計算結果

而根據積分運算，可得：

說明復化梯形和復化辛普森公式計算出的結果基本一致，與實際結果相符。

3

、程序代碼

functionyy=mymulNewtonCotes(ft,a,b,m,n)

%

復化

Newton-Cotes

數值積分公式，即在每個子區間上使用

Newton-Cotes

公式，然后求和

,

%

參考的輸入形式為

mymulNewtonCotes(ft,0,1,10,2)

%

參數說明

:

%ft

——

被積函數，此題中

ft=@(t)t.*exp(t^2/2)

%a

——

積分下限

%b

——

積分上限

%m

——

將區間

[a,b]

等分的子區間數量

%n

——

采用的

Newton-Cotes

公式的階數，必須滿足

n<8

，否則積分沒法保證穩定性

%(1)n=1

時為復化梯形公式

總結

以上是生活随笔為你收集整理的复化科特斯公式matlab_matlab实现复化NewtonCotes公式求积分的程序应用和代码的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。