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

House

MATLAB源程序代码分享：MATLAB梯形法求离散数据点的数值积分
%% MATLAB trapz 函数, 使用梯形法, 求解离散数据点的数值积分
clear;clc;close all
x=[-3 -1 0.5 2 3 6 7.5];   % 离散数据点的横坐标 x
y=[1 2 5 5.5 1.5 3 1];     % 离散数据点的纵坐标 y
- F0 b8 F. M8 U  `4 q+ \# kQ=trapz(x,y)               % 计算积分结果：离散数据点所连成的折线与 x 轴所围面积

%% 动画演示, 梯形法求解数值积分的计算过程
plot(x,y,'*','MarkerSize',10,'Color',[0.2 0.7 0.2])  % 绘制离散的数据点
" ^: j7 o# N" N( ~* Vaxis([-4 8 0 6])
/ z6 z2 m4 `# W3 F* `xlabel('横坐标 x')
ylabel('纵坐标 y')
title('梯形法求解数值积分的动画演示')
hold on
, `$ p* m# S( D( D% Q# v
for k=1:length(x)-1
( _) x% L; g9 A* G" s) S2 k. }0 F% b- I! ^    pause(1)
    plot(x(k:k+1),y(k:k+1),'--','LineWidth',2,'Color',[0.3 0.7 1])  % 将离散的数据点, 使用直线依次相连
end
4 v$ _' ?, N: O; N% C  j2 ], h) g
) O9 t- v8 i% L/ G  gfor k=1:length(x)
    pause(1)
6 j1 x! u4 t( \    plot([x(k) x(k)],[0 y(k)],'--','LineWidth',2,'Color',[0.5 0.2 0.6])  % 绘制各离散数据点与横轴的垂线, 得到各个小梯形
end
8 }7 O7 f' t  R
, y1 J; f. g" b9 y7 ~( zfor k=1:length(x)-1
! f4 w) n! o% ~4 y6 u% w    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])  % 依次给每一个小梯形面积着色, 表示各个小梯形面积的求和过程
end

  G  y- y8 E4 v) O5 x& z

fanichicl

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