使用方法如下:
I=trapz(x,y)
其中x和y分別是自變數和對應的值,例如有函式y=x^3-2x-3,為了計算在[0,1]上的積分,可以這麼做:
>> format compact
>> x=0:0.05:1;
>> y=x.^3-2.*x-3;
>> I=trapz(x,y)
I =
-3.7494
這個函式是可以直接使用經典積分理論計算的,精確值為 -15/4=-3.75,誤差為0.016%。
擴充套件資料:
注意事項
MATLAB中的trapz()函式是基於復化梯形公式設計編寫的,其一般呼叫格式為:
I=trpaz(x,y,dim)
其中x,y是觀測資料,x可以為行向量或列向量,y可以為向量或矩陣,y的行數應等於x向量的元素個數;dim表示按維進行求積,若dim=1(預設值),則按行求積,若dim=2,則按列求積。
如:計算函式y=x^3-2x-3,為了計算在[0,1]上的積分
x=0:0.05:1;
y=x.^3-2.*x-3;
trapz(x,y)
ans =
使用方法如下:
I=trapz(x,y)
其中x和y分別是自變數和對應的值,例如有函式y=x^3-2x-3,為了計算在[0,1]上的積分,可以這麼做:
>> format compact
>> x=0:0.05:1;
>> y=x.^3-2.*x-3;
>> I=trapz(x,y)
I =
-3.7494
這個函式是可以直接使用經典積分理論計算的,精確值為 -15/4=-3.75,誤差為0.016%。
擴充套件資料:
注意事項
MATLAB中的trapz()函式是基於復化梯形公式設計編寫的,其一般呼叫格式為:
I=trpaz(x,y,dim)
其中x,y是觀測資料,x可以為行向量或列向量,y可以為向量或矩陣,y的行數應等於x向量的元素個數;dim表示按維進行求積,若dim=1(預設值),則按行求積,若dim=2,則按列求積。
如:計算函式y=x^3-2x-3,為了計算在[0,1]上的積分
x=0:0.05:1;
y=x.^3-2.*x-3;
trapz(x,y)
ans =
-3.7494