基于MATLAB的三重积分与曲线积分

这篇具有很好参考价值的文章主要介绍了基于MATLAB的三重积分与曲线积分。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

目录

一. 三重定积分的数值求解

例题1

例题2

例题3

二. 第一类曲线积分:对弧长的曲线积分

例题4

例题5

三. 第二类积分:对坐标的曲线积分

例题6

例题7


一. 三重定积分的数值求解

三重定积分的数学表达形式为:

在MATLAB中调用的格式为如下:

在调用时,最后一个@quadl为具体求解一元积分的数值函数,当然也可以修改为@quad或者自己编写的函数,只要调用格式同上式子就行。与求解的精度相关。

例题1

计算三重积分

解:

MATLAB代码:

clc;clear;
triplequad(inline('4*x.*z.*exp(-x.*x.*y-z.*z)',...
    'x','y','z'), 0, 1, 0, pi, 0, pi,1e-7,@quadl)

运行结果:

ans =1.732762284542535

integral3也可以对三重积分进行数值计算。正常的调用格式为:

q=integral3(fun,xmin,xmax,ymin,ymax,zmin,zmax)

如果增加相关的参数限制,调用格式为:

q=integral3(fun,xmin,xmax,ymin,ymax,zmin,zmax,Name,Value)

例题2

在笛卡尔坐标系中对单位球面计算积分,要求在MATLAB代码中使用到匿名函数

基于MATLAB的三重积分与曲线积分

解:

MATLAB代码:

clc;clear;
fun=@(x,y,z)x.*cos(y)+x.^2.*cos(z);
xmin=-1; xmax=1;
ymin=@(x)-sqrt(1-x.^2);   
ymax=@(x)sqrt(1-x.^2);
zmin=@(x,y)-sqrt(1-x.^2-y.^2);
zmax=@(x,y)sqrt(1-x.^2-y.^2);
q=integral3(fun,xmin,xmax,ymin,ymax,zmin,zmax,'Method','tiled')    
%使用 'tiled' 方法计算定积分

运行结果:

q =0.779555454656150

例题3

求函数基于MATLAB的三重积分与曲线积分在区域的积分。其中参数a=2。

解:

MATLAB代码:


clc;clear;
a=2;
f=@(x,y,z)10./(x.^2+y.^2+z.^2+a);
q1=integral3(f,-Inf,0,-100,0,-100,0)
q2=integral3(f,-Inf,0,-100,0,-100,0,'AbsTol',0,'RelTol',1e-9)

运行结果:

q1 =2.734244598320928e+03


q2 =2.734244599944285e+03

MATLAB工具箱中已经拥有了很多计算数值积分的函数,包含integral,integral2,integral3,这些函数都可以计算积分限为无穷大的积分。由于精度限制的原因quad,quadl,quadv,dblquad,triplequad这几个函数可能后续会被MATLAB删除。

二. 第一类曲线积分:对弧长的曲线积分

第一类曲线积分起源于对不均匀分布的空间函数曲线总质量的求取。设空间曲线L的密度函数为f(x,y,z),则其总质量可以利用积分求解如下:

积分中s为曲线上某点的弧长,所以第一类积分又称之为对弧长的曲线积分。

在数学中,弧长积分可以表示为基于MATLAB的三重积分与曲线积分,可以简化记作基于MATLAB的三重积分与曲线积分,由此原来的曲线积分可以转化为与t相关:

基于MATLAB的三重积分与曲线积分

例题4

求解基于MATLAB的三重积分与曲线积分。其中l为为螺线。

解:

MATLAB代码:

clc;clear;
syms t; 
syms a positive; 
x=a*cos(t); y=a*sin(t); z=a*t;
I=int(z^2/(x^2+y^2)*sqrt(diff(x,t)^2+diff(y,t)^2+ diff(z,t)^2),t,0,2*pi);
pretty(I)

运行结果: 

基于MATLAB的三重积分与曲线积分

数学形式结果:。

例题5

基于MATLAB的三重积分与曲线积分,其中l曲线为和的交线。

解:

MATLAB代码:

clc;clear;
X=0:.001:1.2; 
Y1=X; 
Y2=X.^2;
%绘图
plot(X,Y1,X,Y2) %观察图形决定积分的正负

syms x; 
y1=x; 
y2=x^2; 
I1=int((x^2+y2^2)*sqrt(1+diff(y2,x)^2),x,0,1);
I2=int((x^2+y1^2)*sqrt(1+diff(y1,x)^2),x,1,0); 
%注意I1和I2的积分自变量顺序
I=I2+I1

运行结果:

I =(349*5^(1/2))/768 - (2*2^(1/2))/3 - (7*log(5^(1/2) + 2))/512

三. 第二类积分:对坐标的曲线积分

第二类积分起源于变力f(x,y,z)沿曲线l移动时做功的研究,如下:

向量

向量也为向量,利用参数方程可以表示此曲线:

到此,利用此两个向量的点乘可以求解此类积分。

例题6

求曲线积分基于MATLAB的三重积分与曲线积分,其中l正向为圆周基于MATLAB的三重积分与曲线积分

解:

MATLAB积分:

clc;clear;
syms t; 
syms a positive; 
x=a*cos(t); y=a*sin(t);
F=[(x+y)/(x^2+y^2),-(x-y)/(x^2+y^2)]; 
ds=[diff(x,t);diff(y,t)];
I=int(F*ds,t,2*pi,0)  % 正向圆周 

运行结果:I =2*pi 

例题7

计算基于MATLAB的三重积分与曲线积分,其中l为抛物线。

解:

MATLAB代码:

clc;clear;
syms x; 
y=x^2; 
F=[x^2-2*x*y,y^2-2*x*y]; 
ds=[1; diff(y,x)];
I=int(F*ds,x,-1,1)

运行结果:I =-14/15 文章来源地址https://www.toymoban.com/news/detail-420264.html

到了这里,关于基于MATLAB的三重积分与曲线积分的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 基于MATLAB的双重积分的数值求解

    目录 一. Gauss求积公式 1.1 数学理论 1.2 例题与MATLAB代码 例题1 二. 双重积分问题的数值解 2.1 数学理论 2.2 MATLAB代码与例题 例题2  例题3  例题4 例题5 例题6 三. 小结 我们知道数值的方法求解积分往往精度会差一点。这个时候 Gauss求积公式 可以让数值求解积分有较高的代数精度

    2023年04月26日
    浏览(38)
  • 第一型曲线积分与第一型曲面积分、第二型曲线积分与格林公式

    提示:本文的适用对象为已修过《微积分A1》的非数学系学生,文中题型方法为个人总结,为个人复习使用。部分理解虽然不太严谨,但对于解题的实用性较强。若有疏漏or错误,欢迎批评指正。 对于已经熟知第一型曲线积分和第一型曲面积分定义的朋友们来说,我在这里主

    2024年02月04日
    浏览(48)
  • 宋浩高等数学笔记(十一)曲线积分与曲面积分

            个人认为同济高数乃至数学一中最烧脑的一章。。。重点在于计算方式的掌握,如果理解不了可以暂时不强求,背熟积分公式即可。此外本贴暂时忽略两类曲面积分之间的联系,以及高斯公式的相关内容,日后会尽快更新,争取高效率学习。         在数学中

    2024年02月13日
    浏览(42)
  • 【高等数学笔记】两类曲线积分、曲面积分的转化

    整体思想:局部均匀化,用很小的长度/面积元上一点某个量的数值来代替整个元的数值。 设曲线 Γ Gamma Γ 的参数方程为 x = x ( t ) , y = y ( t ) , z = z ( t ) x=x(t),y=y(t),z=z(t) x = x ( t ) , y = y ( t ) , z = z ( t ) 。令 r = ( x , y , z ) bm r=(x,y,z) r = ( x , y , z ) ,则方程为 r = r ( t ) bm r=bm r(t

    2024年02月04日
    浏览(40)
  • 曲线与曲面积分公式整理

    当f(x,y)=1时,表示曲线L的长度 表示线密度为f(x,y)的曲线质量 沿L运动的变力F=f(x,y)做的功 其中cosα与cosβ是L在(x,y)处的切向量相对于x轴和y轴的方向余弦 其中L是单连通区域D的正向边界 其中L是复连通区域D外部正向边界,l(小写L)是复连通区域D内部正向边界(假设D内只有一个“洞

    2024年02月04日
    浏览(46)
  • 高等数学啃书汇总重难点(十一)曲线积分与曲面积分

    依旧是公式极其复杂恶心的一章,建议是: 掌握两种线面积分的计算套路即可 ,和第8章一样属于同济版教材中最不重要的章节,不会对底层理解做过多考察~ 1.弧长曲线积分的几何意义 2.弧长曲线积分的定义和性质 3.弧长曲线积分的计算方式 4.坐标曲线积分的几何意义 5.坐标

    2024年02月06日
    浏览(37)
  • 基于改进ISODATA算法的负荷场景曲线聚类MATLAB程序

    🍞正在为您运送作品详情 https://mbd.pub/o/bread/ZJqclJ9v 参考文献: 基于机器学习的短期电力负荷预测和负荷曲线聚类研究_张辰睿(硕士论文) 参考其第三章 主要内容: 主要包含四种聚类算法,K-means聚类、ISODATA聚类、L-ISODATA聚类及K-L-ISODATA聚类,并且包含了对聚类场景以及聚

    2024年02月13日
    浏览(39)
  • 使用matlab进行回声处理(三重回声)-数字信号处理课设

    1.录制一段声音信号,作为原音频信号,生成频域与时域图形,观察与分析其时域 与频域图形。 2.对该音频信号进行时域处理,并且实现对该声音信号添加第一层回声。 3.对该音频信号再继续一次进行时域处理,实现对该声音信号添加第二层回声。 4.将三段音频信号进行合成

    2024年02月03日
    浏览(46)
  • 曲线生成 | 基于多项式插值的轨迹规划(附ROS C++/Python/Matlab仿真)

    🔥附C++/Python/Matlab全套代码🔥课程设计、毕业设计、创新竞赛必备!详细介绍全局规划(图搜索、采样法、智能算法等);局部规划(DWA、APF等);曲线优化(贝塞尔曲线、B样条曲线等)。 🚀详情:图解自动驾驶中的运动规划(Motion Planning),附几十种规划算法 多项式插值(polynomial

    2024年02月03日
    浏览(36)
  • 【MATLAB源码-第116期】基于matlab的SM空间调制系统仿真并且和QSM进行对比,输出误码率曲线。

    空间调制(SM): 1. 基本概念:SM是一种MIMO(多输入多输出)通信技术,它利用天线阵列的空间维度来传输信息。 2. 工作原理:    - 在每个时间时隙,SM从多个天线中选择一个天线来发送信号。    - 数据流被分为两部分:一部分用于选择激活的天线,另一部分则用于调

    2024年01月18日
    浏览(44)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包