实验四 用Mathematica解多元函数微积分问题

基于Matlab软件求解多元函数积分一、多元函数积分的概念及背景多元函数积分是对多元函数在一定区域内求和得到的结果，它类似于一元函数积分，但是需要考虑到多个自变量的情况。

在实际应用中，多元函数积分可以用来计算体积、质心、质量、惯性矩、功与位的转换等问题，因此具有广泛的应用价值。

在Matlab中，多元函数积分可以通过syms工具箱中的int函数来求解。

int函数能够处理一元和多元的定积分，通过指定积分变量和积分区间的方式，可以求解出多元函数在给定区域内的积分结果。

1. 定义多元函数在使用Matlab求解多元函数积分之前，首先需要定义待积的多元函数。

Matlab中可以使用syms函数定义符号变量，再通过这些符号变量来定义多元函数。

我们定义一个二元函数 f(x, y) = x^2 + y^2，可以使用如下代码来定义：syms x yf = x^2 + y^2;2. 求解多元函数积分定义好多元函数后，就可以使用int函数来求解多元函数积分。

int函数的语法格式为：int(F, x_min, x_max, y_min, y_max)其中F为待积的多元函数，x_min和x_max分别为x变量的积分下限和上限，y_min和y_max分别为y变量的积分下限和上限。

我们求解函数f在区域R={(x, y)|0≤x≤1,0≤y≤1}内的积分，可以使用如下代码来求解：result = int(f, 0, 1, 0, 1);3. 显示积分结果可以使用disp函数将求解出的积分结果进行显示。

我们使用如下代码来显示上述求解结果：disp(result)通过上述三个基本步骤，就可以使用Matlab求解多元函数积分了。

三、实例演示下面通过一个实例来演示如何使用Matlab对多元函数积分进行求解。

假设我们要求解函数f(x, y) = x^2 + y^2在区域R={(x, y)|0≤x≤1, 0≤y≤1}内的积分。

我们使用syms函数定义符号变量x和y，并定义函数f：然后，我们使用int函数对函数f在R内进行积分求解：我们通过disp函数来显示求解结果：在Matlab命令窗口中执行以上代码，将得到函数f在区域R内的积分结果为2/3。

472 实习四 定积分以及相关应用问题实习目的1.掌握用Mathematica 求定积分2.用定积分求面积、平面曲线的弧长和旋转体的体积。

实习准备1.定积分的运算在不定积分中加入积分的上下限便成为定积分(definite integral)。

Mathematica 的定积分命令和不定积分的命令相同，但必须指定积分变量的上下限。

(1) Integrate[f,{x,下限,上限}](2) ⎰dx x f b a )(例1 计算定积分⎰-dx xx 151。

解 dx xx In 1:]1[51-=⎰ Out[1]=4-2ArcTan[2]和不定积分一样,除了我们指定的积分变量之外,其它所有符号都被作常数处理.例2 计算定积分⎰+dx e x a x 3220。

解 dx a xExp x In ]3[:]2[220+=⎰ 2726272]2[6aa e e Out ++-= 1 数值积分如果Mathematica 无法解出积分的符号表达式或者定积分的结果过于冗长而失去意义时，我们就可以用数值积分求解。

数值积分只能进行定积分的运算，即必须指定上、下限。

用Mathematica 求解数值积分有两种形式：(1) NIntegrate[f,{x,a,b}] x 从a 到b ，做)(x f 的数值积分。

(2) N[⎰dx x f b a )(] 求定积分表达式的数值例3 求定积分⎰dx x )sin(sin 30π。

解 用Integrate 命令无法求)sin(sin x 的定积分，用NIntegrate 命令即可求得473其数值积分。

In[1]:=NIntegrate[Sin[Sin[x]],{x,0,Pi/3}]Out[1]=0.466185求定积分表达式的数值，也能得到与上式相同的结果。

]dx ]]Sin[Sin[N[:]2[In 3/0x Pi ⎰=Out[2]=0.466185例4 求定积分dx e x 210-⎰的近似值。

mathematica多元函数赋值Mathematica是一种强大的符号计算软件，用于数学、科学和工程领域。

它具有许多功能，其中之一是计算多元函数赋值。

本文将介绍Mathematica中如何进行多元函数赋值，以及它在数学和科学领域中的实际应用。

在Mathematica中定义多元函数非常简单。

以下是一个多元函数的示例：f[x_, y_] := x^2 + y^2;在这个例子中，f是一个两个变量的函数，x和y。

我们可以评估这个函数在某个特定的x和y值上，如下所示：f[2, 3]这会得到13，因为f在x=2和y=3时等于13。

在Mathematica中，我们可以使用赋值运算符来将一个值分配给多个变量。

例如，以下代码片段为x、y和z分配了值：x = 2; y = 3; z = 4;现在，我们可以通过调用f[x, y, z]来计算函数f在这些值上的值。

这将返回29。

除了简单的函数之外，我们还可以在Mathematica中处理更加复杂的多元函数。

我们可以使用内置的函数，如Sin、Cos和Exp，以及一系列高级函数来定义多元函数。

例如，以下代码片段定义了一个三变量函数：g[x_, y_, z_] := Sin[x] + Cos[y] + Exp[z];我们还可以使用Mathematica中的其他功能，如求导和积分，以更深入地了解多元函数。

例如，我们可以使用D 函数来找到函数的偏导数。

以下代码片段将计算函数g在x=1、y=2和z=3的偏导数：D[g[x, y, z], x] /. {x -> 1, y -> 2, z -> 3}这将返回cos（1），因为g在x=1，y=2和z=3处的x 偏导数为cos（1）。

我们还可以使用Mathematica中的Plot3D函数来可视化三元函数。

例如，以下代码片段将创建并显示函数g在x 和y的范围内的3D图表：Plot3D[g[x, y, z], {x, -2π, 2π}, {y, -2π,2π}]这将显示一个美丽的3D函数，显示了函数的值在x、y和z坐标上的变化。

实验九 用Mathematica 软件求函数偏导数与多元函数的极值实验目的:掌握用Mathematica 软件求函数偏导数与全微分、多元函数的极值的语句和方法。

实验过程与要求:教师利用多媒体组织教学,边讲边操作示范。

实验的内容:一、求偏导数在Mathematica 系统中与求一元函数导数类似用D 函数求函数f 的偏导数，基本格式为：D[f ,{变量,n }] 给出对变量的n 阶偏导数.D[f ,变量1,变量2，…] 给出高阶混合偏导数.实验 求y x x z cos sin +=的两个一阶偏导数和四个二阶偏导数.解 In[1]:=Clear[x ,y ]In[2]:=f [x _,y _]:=Sin[x ]+x *Cos[y ]In[3]:=D[f [x,y ],x ]In[4]:=D[f [x,y ],y ]In[5]:=D[f [x,y ],{x ,2}]In[6]:=D[f [x,y ],{y ,2}]In[7]:=D[f [x,y ],x,y ]In[8]:=D[f [x,y ],y,x ]Out[3]=Out[4]=Out[5]=Out[6]=Out[7]=Out[8]=二、求全微分在Mathematica 系统中与求一元函数微分类似用Dt 函数求函数f 的全微分，基本格式为：Dt[f ]实验 求函数206933+-+-+=y x xy y x z 的全微分.解 In[9]:=Dt[x ^3+y ^3-x *y +9x -6y +20]Out[9]=三、求多元函数的极值在Mathematica 系统中与求一元函数极小值类似用FindMinimum 函数求多变量函数f 的极小值，基本格式为：FindMinimum [f ,{x,x 0}，{y, y 0},…]其中{ x 0，y 0，…}为初始值，表示求出的是f 在（x 0，y 0，…）附近的极小值.因此，一般需借助于Plot3D 函数先作出函数的图象，由图象确定初始值，再利用FindMinimum 求出f 在（x 0，y 0，…）附近的极小值.仍用FindMinimum 函数求函数的极大值，基本格式为：FindMinimum [-f,{x,x 0}，{y, y 0},…]其中{ x 0，y 0，…}为初始值，表示求出的是-f 在（x 0，y 0，…）附近的极小值，设为W ，实际上间接地求出了f 在（x 0，y 0，…）附近的极大值，为-W.实验 求函数206922+-+-+=y x xy y x z 的极值.解 In[10]:=Clear[f,x,y ]In[11]:=FindMinimum[x ^2+y ^2+9*x -x *y-6y +20,{x ,-4},{y,-4}]In[12]:=Plot3D[x ^2+y ^2+9*x -x*y-6y +20,{x ,-4,5},{y ,-4,5}]Out[11]=表示z 在x =-4,y =1处取得极小值-1该函数无极大值.图形如图实验)ln()()4()3()2()1(.122222y x y x z yx yx z y x z yx e z xy +-=-+=-=+=求下列函数的偏导数：)ln()()4()arcsin()3()2()1(.2222y x y x z xy z y x z e z yx ++====-求下列函数的全微分：x y x y x z y x y x z 933)2()4)1(.3223322-++-=---=（求二元函数的极值：。