前言:

在第本系列第六篇中,我们讲解了一下相对来讲最普通的不定积分和定积分的问题,相信大家对积分已经有了大概的了解了,可是仍有一些小问题,比如如何计算多个变量的积分,比如二重积分、甚至三重三重积分呢?本讲中,我将做详细讲解。

本期内容

首先复习一下一元的定积分、然后讲解二重积分,如果可能,讲解三重积分。

复习定积分(单变量)

在之前,先复习一下之前学习的内容。使用matlab求解定积分的步骤大概如下:

  1. 定义符号变量(syms关键字)
  2. 定义内联函数(inline函数)
  3. 使用matlab内置函数进行计算定积分或者不定积分(使用int函数,具体求解定积分还是不定积分,根据参数的数量决定)

我们都知道:定积分是求解一个图形与坐标轴围成的面积。而今天我们需要深入理解一下:定积分是求解两条曲线之间围成的面积,如果只有一个公式,比如x,我们可以理解成 x-0,y=x 这个函数减去了 y=0这个常数函数,或者说,上面函数与x轴围成的面积减去了下面函数与x轴围成的面积(并且这个面积在x轴上方是正的,下方是负的),这样我们就可以推广,假设我们现在有两条曲线,一个是sinx,一个是 ,我们求在0-pi的范围内,两个图形围成的面积,容易知道,两个图形大概是这样围起来的:

PS: 上图的代码和讲解如下

%% 目标图像
clc;clear all;  % 清除原先的东西x = 0:0.01:pi;  % x的范围是0-pi,以0.01为一个小分割点细分出一个范围
y1 = sin(x);    % sinx
y2 = exp(x);    % e^xfigure;         % 弄一个figure的框体
hold on;        % 保持这个窗体,保证后面的内容都画在了上面,hold off解除subplot(221);   % 开辟子窗体,221:共2行2列,里面的第一个
fill(x,y1,'r'); % 画图并填充颜色,变量是x,函数是y1,颜色是r,红色
title('y=sinx');% 设置标题,必需先fill,才能添加标题,否则会失败subplot(222);   % 开辟子窗体,222,共2行2列,里面第二个
fill(x,y2,'b'); % b 蓝色
title('y=e^x');subplot(223);   % 开辟子窗体,223,共2行2列,里面第三个
fill([x, fliplr(x)],[y1,fliplr(y2)],'g');% 两个函数直接填充颜色,每一个方括号代表一个函数的内容
title('两个函数中间的部分');
hold off;

OK,那么我们就可以很容易的使用int()函数算出这样一个定积分,代码如下:

%% 对上面的函数进行计算定积分
syms x;
f = 'exp(x) - sin(x)';
res_int = int(f, x, 0, pi);
disp(res_int);

上面左边是我手算结果,右边是matlab计算结果,可见我算的还是很准的哈哈哈哈哈哈哈。

二重积分

上面复习了定积分的计算和一般概念,下面我们看二重积分,先了解下二重积分是什么吧:

二重积分的一个概念是求体积(还有其他的,比如薄片质量等,在这里先不谈,先说体积的问题),那么明确概念其实就不难了,定积分求面积、二重积分求体积,那么这个体积是哪里的体积呢?

如图,粉红色的是z=f(x,y)的函数图像,下面是地面的阴影部分,红色的是边界,这个体积就是粉色的顶部与底面阴影之间的一个曲顶柱体的体积,那这个体积怎么求呢?

在定积分中,我们把整个面积分成了无数个小面积,计算后累加,同样的,在二重积分中,我们也采取同样的方法,我们把这个柱体切成一个一个的小柱体,就像图中的那个黄色柱体一样,假设我们把整个底面分成了无数份,每一份的底面积都无限接近于0,假设这个底面积是d ,那么这个地方的高就是f(x,y),所以这个小柱体的体积就是 f(x,y)d ,我们只需要将这些小体积全部加起来,就能得到整个曲顶柱体的体积公式了:

并且我们很容易发现,积分的区间其实就是底面,被积函数是曲面的公式。公式中各个符号在下图中体现:

现在,计算一个例子给大家看,假设底面如图所示(灰色阴影部分):

先积分x(Y型)时区域可表示为:Dx = {(x,y)|0<=y<=1, y<=x<=1}。

先积分y(X型)时区域可表示为:Dy = {(x,y)|0<=y<=x, 0<=x<=1}。

下面我们以X型为例,使用matlab计算曲面为 f(x,y) = xcos(y)的曲顶柱体的体积:

详细的代码以及注释见下方:

%% 计算二重积分
clear all;
f = @(x,y)x.*cos(y); % 定义一个函数句柄,两个变量分别是x,y
% 这样定义函数句柄时,需要在所有的运算前加一个.(表示对应元素各自计算,不按照矩阵的规则)
ymax = @(x) x;  % 根据0<=y<=x处,右方的y=x反解出y,
res = integral2(f,0,1,0,ymax);
% 计算二重积分的函数
% integral2(fun, xmin, xmax, ymin, ymax)
disp(res);

最后,让我们欣赏一下这个曲面的形状吧:

代码如下:

%% 看一看这个曲面吧
[x,y] = meshgrid(0:0.001:1, 0:0.001:1);
% 设置x和y的区间,以及细分程度
z = (x.*cos(y));  % 使用一个z来代表这个函数
mesh(x,y,z)  % 画出二维曲面图像
title('x*cos(y)')

那么,今天就先到这里吧,三重积分留到下次在讲解,谢谢观看!

今日小结

今天我们学到了如下的知识:

  1. 复习了之前学习的定积分、不定积分的内容
  2. 学习了二重积分的概念,以及应用的大概场景:求曲顶柱体体积,平面薄片质量(底面换位薄片,被积函数即薄片面密度的函数)等。
  3. 学会了如何在一个figure里面绘制多个图像
  4. 学会了如何在matlab中绘制图像、对目标区域填充颜色并命名标题
  5. 学会了怎么使用matlab的integral2计算二重积分
  6. 学会了怎么在matlab中声明(定义)一个函数句柄,或者说创建一个多元函数,以及怎么反解一个函数中某个变量(y = @(x) x  ,此处函数是y=x),其实说白了,反解函数就是定义函数句柄。
  7. 学会了怎么在matlab中绘制二维曲面的图像。

大家下期再见!!

MATLAB新手简明使用教程(八)——高级积分运算、二重积分——新手来看,保证看懂相关推荐

  1. UI设计新手常遇到的问题及解决方法,快来看一看!

    UI是研究人机交互的学科,单从目前它的应用领域来看主要应用于软件.互联网.移动设备.游戏和虚拟现实影音方面这些都是新兴的热门方向和活跃领域.鉴于UI设计行业最近几年的兴起,有许多年轻人都被它好的就业前 ...

  2. MATLAB新手简明使用教程(六)——使用matlab求解定积分和不定积分的问题——新手来看,保证看懂。

    前期回顾 上期的教程虽然简短,但还是有点东西的,上一节中,我们学习了下面的知识点. 在matlab中定义矩阵:A = [1 2 3; 2 3 4; 3 4 5],每行元素空格隔开,两行之间使用分号隔开 ...

  3. MATLAB基础教程(9)——高级积分运算、二重积分

    目录 复习定积分(单变量) 二重积分 今日小结 复习定积分(单变量) 在之前,先复习一下之前学习的内容.使用matlab求解定积分的步骤大概如下: 定义符号变量(syms关键字) 定义内联函数(inl ...

  4. MATLAB新手简明使用教程(七)——使用matlab建立多项式以及求导,商求导乘积求导等——新手来看,保证看懂。

    前期回顾 上一期中,我们学了下面的知识: 定积分的基本概念和一些简单的几何意义. 使用 int 函数计算不定积分. 使用 int 函数计算定积分. 本期内容 本期我打算给大家介绍一下使用matlab对 ...

  5. MATLAB新手简明使用教程(四)——斐波那契数列——新手来看,保证看懂。

    前期回顾 在上一篇我们学习了下面 的知识点: 使用roots函数求一元二次方程的根. 使用solve函数求方程的根. 使用inline函数和ezplot函数一起绘制函数图形. 使用 fzero 函数和 ...

  6. matlab 极限积分,实验二MATLAB中的极限和微分积分运算

    <实验二MATLAB中的极限和微分积分运算>由会员分享,可在线阅读,更多相关<实验二MATLAB中的极限和微分积分运算(28页珍藏版)>请在人人文库网上搜索. 1.实验二 MA ...

  7. matlab 极限积分,实验二matlab中的极限和微分积分运算.ppt

    实验二matlab中的极限和微分积分运算实验二matlab中的极限和微分积分运算 实验二 MATLAB中的极限.微分和积分运算 一.实验目的 熟悉MATLAB软件中关于极限.微分运算和不定积分.定积分 ...

  8. matlab项目实例教程,matlab简明实例教程.doc

    matlab简明实例教程.doc 南京航空航天大学王正盛1MATLAB数学工具软件实例简明教程王正盛编写南京航空航天大学南京航空航天大学王正盛2第一章MATLAB简介MALAB译于矩阵实验室MATri ...

  9. python新手入门教程-Python简明入门教程

    本文实例讲述了Python简明入门教程.分享给大家供大家参考.具体如下: 一.基本概念 1.数 在Python中有4种类型的数――整数.长整数.浮点数和复数. (1)2是一个整数的例子. (2)长整数 ...

最新文章

  1. C 语言编程 — const 关键字
  2. RetinaFace笔记
  3. python九九乘法口诀表-python打出九九乘法口诀表
  4. 为啥总让我“先去博客园其他网站逛逛”?
  5. Pandas数据清洗工具箱
  6. 调用WCF返回Josn的两种方式
  7. linux进程管理fork,Linux -- 进程管理之 fork() 函数
  8. vue中在data中写入img地址 引入问题
  9. Atomic包的4种类型详解
  10. Windows下的Linux开发环境Cygwin
  11. php如何写文件缓存,PHP代码篇(八)--php实现文件缓存的读写
  12. jqgrid使用小记
  13. 多标签文本分类研究进展
  14. CentOS安装并设置MariaDB
  15. Mysql基础-常用sql语句
  16. 企业微信第三方服务商和钉钉ISV开发对比
  17. 网站一直被大流量DDOS攻击,服务器被打死了怎么办?
  18. java获取本月最后一天
  19. 联邦学习纵向逻辑回归公式推导
  20. 定时脚本任务列子(crontab)

热门文章

  1. 动态链接库(DLL)初始化例程失败 解决方法
  2. 利用Django生成PDF文档
  3. 三维空间:点到直线垂足坐标公式推导
  4. ECS+SLB负载均衡实践
  5. 组态王图素制作_组态王教程——很实用(杂项)
  6. 佛山2000坐标-国家2000坐标双向重投影工具
  7. OC 基础 UIButton
  8. 电话号码中间四位隐藏
  9. 医院管理系统完整项目代码以及数据库建表语句分享(包括加盐和拦截器)
  10. 2018年携程笔试分享