matlab计算三角格网面积,MATLAB计算积分

函数的积分和椭圆的周长

1．正弦函数的积分

[问题]求正弦函数从0到π的积分

y = sin x

当x = 0时，积分为0，画出积分的函数曲线。

[数学模型]

定积分的结果为

ππ00

sin d cos 2S x x x ==-=? 不定积分的结果为

sin d cos I x x x C ==-+?

其中C 是积分常量，由初始条件决定。当x = 0时，积分为I = 0，必有C = 1。结果为

I = -cos x + 1

[算法]根据积分的基本概念，将积分区域分为多份，用矩形法求曲线下的近似面积表示积分的近似值

1()n

i i S f x x ==?∑

矩形法的函数是sum(f)。

用梯形法求曲线下的近似面积表示积分的近似值

1

101[()()]2

n i i i S f x f x x -+==+?∑

梯形法的函数是trapz(f)。

用数值积分的函数是quad 和quadl ，常用使用格式是

S = quad(f,a,b)

其中，f 表示被积函数，a 表示积分的下限，b 表示积分的下限。

用符号的函数是int ，常用使用格式是

S = int(f,a,b)

[程序]zqy4_1.m 如下。

%正弦函数的积分

clear %清除变量

x=linspace(0,pi); %自变量向量

dx=x(2); %间隔

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

s1=sum(y)*dx %矩形法积分

s2=trapz(y)*dx %梯形法积分

f=inline('sin(x)'); %被积的内线函数

s3=quad(f,0,pi) %数值定积分