matlab作业重积分

发布 2022-09-20 19:59:28 阅读 6254

《matlab语言》课程**。

matlab在重积分计算中的应用。

姓名:学号:

专业:班级:

指导老师:学院:

完成日期:matlab在重积分中的应用。

摘要]高等数学课程是理工科专业中非常重要的基础课程,重积分是《高等数学》相对难学的部分,并且计算复杂,而matlab 软件在求解重积分的数值解方面有较大优势,既能进行数值求解,又能绘制有关曲线,非常方便实用。而且可以极大地提高同学们的学习兴趣,培养同学们利用matlab解决实际问题的能力。

关键词]matlab 高等数学重积分计算。

一、 问题的提出。

matlab是矩阵实验室(matrix laboratory)的简称,是 1984 年美国 mathworks 公司推出的一套高性能的数值计算和可视化数学软件,被誉为“巨人肩上的工具”。 现已成为国际公认的最优秀的科技应用软件之一。由于使用matlab编程运算与人进行科学计算的思路和表达方式完全一致,所以不像学习其他高级语言——如basic、fortran和c语言等那样难以掌握,用matlab编写程序犹如在演算纸上排列出公式与求解问题,所以又被称为演算纸式科学算法语言。

在这个环境下,对所要求解的问题,只需简单的列出数学表达式,其结果便以数值或图形方式显示出来。在英美等发达国家的理工类大学里,matlab 软件是大学生必须掌握的一种基本工具;在研究设计单位和工业部门,它更是研究和解决计算问题的标准软件,是工程技术人员必备的软件。

众所周知,高等数学已经被大部分高校列为重要的公共基础课,高等数学中许多重要方法,如求极限、求导数、求不定积分、求定积分、解常微分方程、向量运算、求偏导数、计算重积分、级数展开等,只靠笔算是难以完成的。为提高同学们用高等数学解决实际问题的能力,使同学们能在理解、掌握数学理论知识的同时,迅捷地计算出繁杂的数**算结果,而不必去考虑用什么算法以及如何实现等问题,用matlab求高等数学中问题,能快捷、准确地得出解,显示出matlab在数学计算上的优越性。

本**讲述matlab在重积分中的应用,重积分是高等数学教学的重点和难点,特别是重积分的计算比较复杂,其中涉及积分区域的确定、交换积分次序等; matlab 软件在求解重积分的数值解方面有较大优势,其函数库中包含了求解简单重积分的函数; 对复杂的积分,可以先利用 matlab 绘制出积分区域,确定积分区域类型,再利用相应函数进行求解; 针对二重积分和三重积分的几种不同情况利用 matlab 进行了实例说明,并给出了相应的程序; matlab 操作和演示较为方便,在辅助重积分教学方面优势明显。为提高读者用高等数学解决实际问题的能力,使学生能在理解、掌握数学理论知识的同时,迅捷地计算出繁杂的数**算结果,而不必去考虑用什么算法以及如何实现等问题,可以把课程中的计算用数学软件来解决,让学生亲自上机实践,可以提高学习兴趣,提高解题效率和学习效果。本文阐述 matlab 在高等数学微积分中的应用,以加强学生对学习数学的兴趣和能力。

二、重积分的有关理论。

定积分:积分是微分的无限和,函数在区间上的积分定义为。

其中从几何意义上说,对于上非负函数,记分值是曲线与直线及轴所围的曲边梯形的面积。有界连续(或几何处处连续)函数的积分总是存在的。

二重积分的定义为。

当非负时,积分值表示曲顶柱体的体积。二重积分的计算主要是转换为两次单积分来解决,无论是解析方法还是数值方法,如何实现这种转换,是解决问题的关键。

重积分:重积分的数值计算可通过若干次单积分的组合实现,如对于二重积分。

先化为二次计分利用梯形法,先将区间等分,利用梯形积分公式可得。

再将区间等分,利用梯形积分公式可得。

微积分基本定理(newton-leibniz公式):在上连续,且,则有。这个公式表明导数与积分是一对互逆运算,它也提供了求积分的解析方法:

为了求的定积分,需要找到一个函数,使的导数正好是,我们称是的原函数或不定积分。不定积分的求法有学多数学技巧,常用的有换元积分和分部积分法。从理论上讲,可积函数的原函数总是存在的,但很多被积函数的原函数不能用初等函数表示,也就是说这些积分不能用解析方法求解,需用数值积分法解决。

在应用问题中,常常是利用微分进行分析,而问题最终归结为微分的和(即积分)。一些更复杂的问题是含微分的方程,不能直接积分求解。多元函数的积分称为多重积分。

三、matlab在重积分计算中的应用。

1、matlab在一重积分计算中的应用。

不定积分的计算:在matlab中,int函数用于求符号函数的不定积分,有两种调用格式:(1) int(f):

没有指定积分变量和积分阶数时,系统按findsym函数指示的默认变量对被积函数或符号表达式f求不定积。(2) int(f,v):以v为自变量,对被积函数或符号表达式f求不定积分。

求符号函数的积分也可使用int(f,v,a,b),其中a、b分别表示积分的下限和上限。该函数求被积函数f在区间[a、b]上的定积分。a和b可以是两个具体的数,也可以是一个符号表达式还可以是无穷(inf)。

若为一元函数求积分,则变量v省略;若为多元函数求积分。则需指明变量为v。若是求不定积分,则上下限a,b省略,当函数f关于变量x在闭区间[a,b]上可积时,函数返回一个定积分结果。

当a、b中有一个是inf时,函数返回一个广义积分。当a,b中有一个符号表达式时,函数返回一个符号函数。

整理的matlab求函数的积分基本语句:

1)int(函数) %计算不定积分;

2)int(函数,变量名) %计算不定积分;

3)int(函数) %计算定积分;

4)int(函数变量名) %计算定积分.

用matlab命令计算定积分的相关函数及简介。

1) sum(a):求数组a的和。

2) format long:长格式,即屏幕显示15位有效数字。

3) double():若输入的是字符则转化为相应的asxii码;若输入的是整型数值则转化为相应的的实型数值。

4) trapz():梯形法求数值积分。

格式:trapz(x,y):其中x为带有步长的积分区间;y为数值形式的运算。

5) quad():抛物线法求数值积分。

格式:quad(fun,a,b):注意此处的fun是函数,并且为数值形式的,所以使用*、\等运算时要在其前加上小数点,即。*、等。

6) fblquad(‘fun’,a,b,c,d):矩形区域二重数值积分,fun表示被积函数的m函数名,a,b分别为x的上、下限,c,d分别为y的上、下限。

例1、计算.

解输入命令。

> syms x %定义符号变量x

int(x^2*log(x计算原积分。

ans =x^3*(log(x) -1/3))/3

注意设置符号变量。

例2、计算下列不定积分:;

解首先建立函数向量。

> syms x %定义符号变量x

y=[sqrt(a^2-x^2),(x-1)/(3*x-1)^(1/3),x^2*asin(x)];定义y的内容。

int(y,x) %对y积分可得对y的每个分量积分的结果。

ans = (a^2*asin(x/(a^2)^(1/2)))2 + x*(a^2 - x^2)^(1/2))/2, (3*x - 1)^(2/3)*(3*x - 6))/15, (x^3*asin(x))/3 + 1 - x^2)^(1/2)*(x^2 + 2))/9]

例3、计算不定积分。

输入程序如下:

> x=sym('x定义符号变量x

f=(3-x^2)^3将所给表达式用f表示。

int(f求表达式f的不定积分。

ans = x^7/7 + 9*x^5)/5 - 9*x^3 + 27*x

例4、求函数f(x)=1/(2x+2x+3),g(x)=1/(2x+2x-3)在负无穷到正无穷的积分。

输入程序如下:

> syms x定义一个符号变量x

f=1/(x^2+2*x+3将f(x)表达式用f表示。

g=1/(x^2+2*x-3将g(x)表达式用g表示

intf=int(f,-inf,inf) %求从负无穷到正无穷f的积分 intg=int(g,-inf,inf) %求从负无穷到正无穷g的积分。

intf =

pi*2^(1/2))/2

结果说明f(x)在整个数轴上可积分,而g(x)在整个数轴上不可积分。

例5、 计算。

解输入命令:

> int(exp(x),0,1) %计算原方程从0到1的积分。

ans =exp(1) -1

例6、 求定积分。

int(exp(x)*(1+sin(x)),x,0,1) %计算f(x)从0到1的积分。

ans =exp(1)-1/2*exp(1)*cos(1)+1/2*exp(1)*sin(1)-1/2

或者。y=exp(x)*(1+sin(x));将f(x)表达式用f表示。

int(y,0,1)

ans= exp(1)-1/2*exp(1)*cos(1)+1/2*exp(1)*sin(1)-1/2

说明:当实际问题中需要求出具体数值时,可以应用vpa命令求出任意精度得结果)

vpa(exp(1)-1/2*exp(1)*cos(1)+1/2*exp(1)*sin(1)-1/2,5)

ans=2.6276

例7、 计算。

解输入命令

> int(abs(x-1),0,2)ans =

例8、判别广义积分、与的敛散性,收敛时计算积分值。

解对第一个积分输入命令syms p real;int(1/x^p,x,1,inf)

得结果ans =limit(-1/(p-1)*x^(-p+1)+1/(p-1),x = inf).由结果看出当时,x^(-p+1)为无穷,当时,ans=1/(p-1),这与课本例题是一致的。

对第二个积分输入命令int(1/(2*pi)^(1/2)*exp(-x^2/2),-inf,inf)

得结果:ans=7186705221432913/18014398509481984*2^(1/2)*pi^(1/2)

matlab函数微积分课程设计

matlab在数学分析 函数微积分 中应用的研究。摘要。matlab是现如今最流行和广泛的数学软件,它将数值分析 矩阵计算 科学数据可视化以及非线性动态系统的建模和 等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究 工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,因此...

matlab作业

2011029170002王柳。a 一个问题的病态性如何,与求解它的算法有关系。错 b 无论问题是否病态,好的算法都会得到它好的近似解。错 c 计算中使用更高的精度,可以改善问题的病态性。错 d 用一个稳定的算法计算一个良态问题,一定会得到它好的近似解。对 e 浮点数在整个数轴上是均匀分布。错 f ...

matlab作业

matlab语言 第3次作业 字符串,单元数组和结构体 专业 海洋技术 海洋测绘方向 姓名 张体强学号 1026222 1 如何将一个char 数据类型的向量转化为相应的double 型数据类型的数据向量。从式1 到8,判断这些语句是否正确。如果它们正确,那么将产生什么结果?这题不要在电脑中做。1....