MATLAB源程序代码分享：MATLAB梯形法求离散数据点的数值积分

House

MATLAB源程序代码分享：MATLAB梯形法求离散数据点的数值积分
7 W# a' c( a" h7 H%% MATLAB trapz 函数, 使用梯形法, 求解离散数据点的数值积分
3 E  J* D1 g; J3 @" g- @4 uclear;clc;close all
x=[-3 -1 0.5 2 3 6 7.5];   % 离散数据点的横坐标 x
3 j1 ~- @5 R5 e" w* {% Ey=[1 2 5 5.5 1.5 3 1];     % 离散数据点的纵坐标 y
* D) B/ a# ^1 |9 A4 V+ G0 nQ=trapz(x,y)               % 计算积分结果：离散数据点所连成的折线与 x 轴所围面积

%% 动画演示, 梯形法求解数值积分的计算过程
plot(x,y,'*','MarkerSize',10,'Color',[0.2 0.7 0.2])  % 绘制离散的数据点
axis([-4 8 0 6])
$ l1 R, G. `0 _' _; yxlabel('横坐标 x')
/ L5 X* A- V; Q9 c4 E* i2 eylabel('纵坐标 y')
. Z/ W# `6 `: L6 o$ Btitle('梯形法求解数值积分的动画演示')
/ y8 T( s. Q+ O+ p/ C! h. d+ Xhold on
/ _- t8 F  R& `% Z9 s* @  [5 l+ ^
) C! x' ~( b8 r" {* R& cfor k=1:length(x)-1
    pause(1)
    plot(x(k:k+1),y(k:k+1),'--','LineWidth',2,'Color',[0.3 0.7 1])  % 将离散的数据点, 使用直线依次相连
% C9 T, A: o8 o; }- Z8 Aend
( P& `$ n$ q( x: v* {1 Z5 ~5 C
for k=1:length(x)
    pause(1)
    plot([x(k) x(k)],[0 y(k)],'--','LineWidth',2,'Color',[0.5 0.2 0.6])  % 绘制各离散数据点与横轴的垂线, 得到各个小梯形
" P* `+ ~5 o  ]6 |% o8 X% xend
* R1 r3 Z3 U- Q1 C
0 Z8 ^$ m5 s' v, X, Efor k=1:length(x)-1
2 S. s( V8 r- @2 ~, E3 u: `% j    pause(1)
    patch([x(k) x(k) x(k+1) x(k+1)],[0 y(k) y(k+1) 0],[1 0.8 0.1])  % 依次给每一个小梯形面积着色, 表示各个小梯形面积的求和过程
( e4 h. E6 e& C- u) R9 t$ @) Rend

fanichicl

MATLAB源程序代码分享：MATLAB梯形法求离散数据点的数值积分。