matlab trapz二重积分函数_小小知识点（十二）利用MATLAB计算定积分

一重定积分

1. Z = trapz(X,Y,dim)

梯形数值积分，通过已知参数x,y按dim维使用梯形公式进行积分

%举例说明1

clc

clear all

% int(sin(x),0,pi)

x=0:pi/100:pi; %积分区间

y=sin(x); %被积函数

z = trapz(x,y) %计算方式一

z = pi/100*trapz(y) %计算方式二

运行结果

被积函数曲线

2、[q,fcnt]= quad(fun,a,b,tol,trace,p1,p2...)

自适应simpson公式数值积分，适用于精度要求低，积分限[a,b]必须是有限的,被积函数平滑性较差的数值积分.

[q,fcnt] = quadl(fun,a,b,tol,trace,p1,p2...)

自适应龙贝格数值积分，适用于精度要求高，积分限[a,b]必须是有限的，被积函数曲线比较平滑的数值积分

%举例说明2

% 被积函数1/(x^3-2*x-p)，其中参数p=5，积分区间为[0,2]

clc

clear all

F = @(x,n)1./(x.^3-2*x-n); %被积函数

Q1 = quad(@(x)F(x,5),0,2) %计算方式一

Q1 = quad(F,0,2,[],[],5) %计算方式二