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

House

MATLAB源程序代码分享：MATLAB梯形法求离散数据点的数值积分
% ~) i" C. d: k/ v! }%% MATLAB trapz 函数, 使用梯形法, 求解离散数据点的数值积分
6 ^0 x3 ~5 w6 Aclear;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
Q=trapz(x,y)               % 计算积分结果：离散数据点所连成的折线与 x 轴所围面积
2 d. X$ u4 s4 X- G
%% 动画演示, 梯形法求解数值积分的计算过程
; N, c: G8 C, o1 Z3 {* Qplot(x,y,'*','MarkerSize',10,'Color',[0.2 0.7 0.2])  % 绘制离散的数据点
axis([-4 8 0 6])
xlabel('横坐标 x')
) o! z7 i. ?+ P* k  n+ jylabel('纵坐标 y')
/ R5 w3 m3 Z0 H) e0 H+ mtitle('梯形法求解数值积分的动画演示')
3 t; B& F4 V$ r. qhold on
( |! T3 V; t2 F. h8 l
' E( @4 i' E' n) Bfor k=1:length(x)-1
    pause(1)
    plot(x(k:k+1),y(k:k+1),'--','LineWidth',2,'Color',[0.3 0.7 1])  % 将离散的数据点, 使用直线依次相连
end
, e8 d  Y1 }2 t
6 E9 w/ f0 [8 A& I: y6 O! wfor k=1:length(x)
2 s% P$ C+ v# X    pause(1)
5 ^$ a3 F3 L% D( C    plot([x(k) x(k)],[0 y(k)],'--','LineWidth',2,'Color',[0.5 0.2 0.6])  % 绘制各离散数据点与横轴的垂线, 得到各个小梯形
end

for k=1:length(x)-1
    pause(1)
# G" I1 Z; H: {5 o8 _% v) s) }    patch([x(k) x(k) x(k+1) x(k+1)],[0 y(k) y(k+1) 0],[1 0.8 0.1])  % 依次给每一个小梯形面积着色, 表示各个小梯形面积的求和过程
  {  p8 Y6 r! h2 g9 Dend

fanichicl

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