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

House

MATLAB源程序代码分享：MATLAB梯形法求离散数据点的数值积分
& @+ R( D; Y7 v2 }4 {. n! u1 P+ @' D%% MATLAB trapz 函数, 使用梯形法, 求解离散数据点的数值积分
5 O4 v0 d7 H7 k9 hclear;clc;close all
/ {2 i+ r' v# u6 L& W7 sx=[-3 -1 0.5 2 3 6 7.5];   % 离散数据点的横坐标 x
3 F6 [( F; F$ x! Q$ ^& P. Ry=[1 2 5 5.5 1.5 3 1];     % 离散数据点的纵坐标 y
! @% D! w  |! D+ TQ=trapz(x,y)               % 计算积分结果：离散数据点所连成的折线与 x 轴所围面积

1 z0 y3 d- C) }* i: X%% 动画演示, 梯形法求解数值积分的计算过程
plot(x,y,'*','MarkerSize',10,'Color',[0.2 0.7 0.2])  % 绘制离散的数据点
axis([-4 8 0 6])
  B; o: G5 A* @* ?! O& j4 Y" W9 sxlabel('横坐标 x')
ylabel('纵坐标 y')
title('梯形法求解数值积分的动画演示')
, }6 j( [! Q2 j9 ^9 ], e  q) e& q- _hold on

! E0 I$ f3 R1 p7 v/ W) F9 Tfor k=1:length(x)-1
    pause(1)
7 V5 F/ b0 W& G  v) P    plot(x(k:k+1),y(k:k+1),'--','LineWidth',2,'Color',[0.3 0.7 1])  % 将离散的数据点, 使用直线依次相连
& Y, B/ P- z' j. _end

9 I+ v) N, G1 _" W: N) xfor k=1:length(x)
    pause(1)
    plot([x(k) x(k)],[0 y(k)],'--','LineWidth',2,'Color',[0.5 0.2 0.6])  % 绘制各离散数据点与横轴的垂线, 得到各个小梯形
+ l1 H' G. H1 o. u2 {$ yend

for k=1:length(x)-1
    pause(1)
( ^) B' R$ d) f% g# g8 v& U    patch([x(k) x(k) x(k+1) x(k+1)],[0 y(k) y(k+1) 0],[1 0.8 0.1])  % 依次给每一个小梯形面积着色, 表示各个小梯形面积的求和过程
. Q- L3 o, ], e) W' R- Cend
# T% m; ~  {" N/ \& l9 u

fanichicl

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