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

House

MATLAB源程序代码分享：MATLAB梯形法求离散数据点的数值积分
%% MATLAB trapz 函数, 使用梯形法, 求解离散数据点的数值积分
1 M; y+ u$ z4 K# ^clear;clc;close all
" e+ g+ G4 |5 T  i; a" M4 Ux=[-3 -1 0.5 2 3 6 7.5];   % 离散数据点的横坐标 x
y=[1 2 5 5.5 1.5 3 1];     % 离散数据点的纵坐标 y
- q9 e9 q( Y  {4 Q! s* h. }Q=trapz(x,y)               % 计算积分结果：离散数据点所连成的折线与 x 轴所围面积

; N* _7 \$ m8 t" P, ?%% 动画演示, 梯形法求解数值积分的计算过程
6 l, u% v! U7 [8 i% j. f8 n$ w- {plot(x,y,'*','MarkerSize',10,'Color',[0.2 0.7 0.2])  % 绘制离散的数据点
axis([-4 8 0 6])
xlabel('横坐标 x')
. r% ~# U# `# q1 Jylabel('纵坐标 y')
title('梯形法求解数值积分的动画演示')
9 `% ^( g3 D' h! \/ U( f2 d. A8 ghold on
4 N# V& U' r+ F6 X: h7 w0 I
  R- h0 s% O# t8 Vfor k=1:length(x)-1
    pause(1)
    plot(x(k:k+1),y(k:k+1),'--','LineWidth',2,'Color',[0.3 0.7 1])  % 将离散的数据点, 使用直线依次相连
9 B1 U4 Z$ |/ j0 q/ z7 M9 kend
8 f& B" b$ }# o/ M9 I- v
, M" W" @/ k7 I& T" dfor k=1:length(x)
    pause(1)
' l3 [9 c) R+ s+ v+ b    plot([x(k) x(k)],[0 y(k)],'--','LineWidth',2,'Color',[0.5 0.2 0.6])  % 绘制各离散数据点与横轴的垂线, 得到各个小梯形
end
6 x# f* R+ J0 i. O& D1 ?$ C& F
for k=1:length(x)-1
, k6 C1 b3 X$ g/ u! z4 U& C    pause(1)
1 G* [" f* U, M8 w9 Y' ?    patch([x(k) x(k) x(k+1) x(k+1)],[0 y(k) y(k+1) 0],[1 0.8 0.1])  % 依次给每一个小梯形面积着色, 表示各个小梯形面积的求和过程
end
3 N% U$ p% |. {6 u( J8 b" z7 o

fanichicl

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