用矩形积分近似求解定积分的动画示意图

解答过程中的发现：

• imwrite的帮助文档说明了如何保存gif，很详细
• 开始定义的syms函数，将被积函数传入作图的函数中，但有些运算不支持syms类型，于是就再转成double型，发现这个过程还挺耗时的！！
• 后来就改成了函数句柄的形式，计算得到的值直接就是double型

还未解决的问题：

• 关于text放置的位置，可能会出现矩形和text重合的情况，但还未找到方法如何自动寻找合适的位置

主程序

clear;clc
fun = @(x)x.*x; % 被积函数
a = 0; % 积分下限
b = 5; % 积分上限
n = 11:5:51; % 区间个数 = n(i) - 1
filename = 'test.gif'; % 保存为 test.gif
delayTime = 1; % 每个动画图像显示之间指定了 delayTime 秒的时滞

for i = 1:length(n)
    myPlot(fun, a, b, n(i));
    
    drawnow;
    F=getframe(gcf);
    I=frame2im(F);
    [I,map]=rgb2ind(I,256);
    if i == 1
        imwrite(I,map,filename,'gif', 'DelayTime',delayTime);
        % 'LoopCount',Inf 使动画连续循环
    else
        imwrite(I,map,filename,'gif','WriteMode','append','DelayTime',delayTime);
    end
end

作图函数

function myPlot(f, a, b, n)
%%% 实现右矩形积分的示意图

% f - 被积函数 function
% a - 积分下限
% b - 积分上限
% n - 划分的区间个数

val = integral(f, a, b); % 计算真实的积分

x = linspace(a, b, n); % 生成a到b等间隔的n个数
y = f(x);

srcX = a:0.001:b;
srcY = f(srcX);

delta = x(2) - x(1); % 间隔
z = sum(delta.*y(2:end));  % 右矩形积分的值

error = z - val; % 近似值 - 准确值

minX = min(srcX);
maxX = max(srcX);
minY = min(srcY);
maxY = max(srcY);

clf
axis([minX, maxX, minY, maxY + 3]);
for i = 1:length(x) - 1
    rectangle('Position',[x(i), 0, delta, y(i + 1)],...
        'facecolor',[255, 203, 153]./255); % 255,203,153 矩形的填充颜色
    hold on
end
hold on
%%% 画图
% 画积分曲线
plot(srcX, srcY, 'red', 'linewidth', 2);

tmp = char(f);
strTitle = strcat('y = ', 32, tmp(5:end));
title(strTitle);

% 创建 textbox
str1 = strcat('误差：', num2str(error));
str2 = strcat('近似值：', num2str(z));
str3 = strcat('真实值：', num2str(val));
str4 = strcat('区间个数：', num2str(length(x) - 1));

annotation('textbox',...
    [0.4 0.67 0.28 0.19],...
    'String',{str1, str2, str3, str4},...
    'LineStyle','none',...
    'FontWeight','bold',...
    'FontSize',14,...
    'FontName','宋体',...
    'FitBoxToText','off');