高数多元函数微分学教案 第七讲 方向导数与梯度

高数多元函数微分学教案  第七讲  方向导数与梯度
高数多元函数微分学教案  第七讲  方向导数与梯度

第七讲 方向导数与梯度

授课题目:

§8. 7 方向导数与梯度

教学目的与要求:

了解方向导数与梯度的概念及其计算方法。

教学重点与难点:

重点与难点:方向导数的概念及方向导数的计算

讲授内容:

一、 方向导数

1、方向导数的概念

回顾偏导数

现在我们来讨论函数z =f (x , y )在一点P 沿某一方向的变化率问题. 设l 是xOy 平面上以P 0(x 0, y 0)为始点的一条射线, e l =(cos α, cos β)是与l 同方向的单位向量. 射线l 的参数方程为

x =x 0+t cos α, y =y 0+t cos β (t ≥0).

设函数z =f (x , y )在点P 0(x 0, y 0)的某一邻域U (P 0)内有定义, P (x 0+t cos α, y 0+t cos β)为l 上另一点, 且P ∈U (P 0). 如果函数增量f (x 0+t cos α, y 0+t cos β)-f (x 0, y 0)与P 到P 0的距离|PP 0|=t 的比值

t

y x f t y t x f ),()cos ,cos (0000-++βα 当P 沿着l 趋于P 0(即t →t 0+)时的极限存在, 则称此极限为函数f (x , y )在点P 0沿方向l 的方向导数, 记作),(00y x l f

??, 即

),(00y x l f

??t

y x f t y t x f t ),()cos ,cos (lim 00000-++=+→βα. 从方向导数的定义可知, 方向导数

)

,(00y x l f

??就是函数f (x , y )在点P 0(x 0, y 0)处沿方向l 的变化率.

注: ),(),(0000y x f y y x x f -?+?+

))()((),(),(220000y x o y y x f x y x f y x ?+?+?+?=

?x =t cos α, ?y =t cos β,

t y x =?+?22)()(.

2、方向导数的计算:

定理 如果函数z =f (x , y )在点P 0(x 0, y 0)可微分, 那么函数在该点沿任一方向l 的方向导数都存在, 且有 ),(00y x l f

??βαcos ),(cos ),(0000y x f y x f y x +=,

其中cos α, cos β是方向l 的方向余弦.

简要证明: 设?x =t cos α, ?y =t cos β, 则

f (x 0+t cos α, y 0+t cos β)-f (x 0, y 0)=f x (x 0, y 0)t cos α+f y (x 0, y 0)t cos β+o (t ). 所以

t

y x f t y t x f t ),()cos ,cos (lim 00000-+++→βα??sin ),(cos ),(0000y x f y x f y x +=.

这就证明了方向导数的存在, 且其值为

),(00y x l f

??βαcos ),(cos ),(0000y x f y x f y x +=.

讨论: 函数z =f (x , y )在点P 沿x 轴正向和负向, 沿y 轴正向和负向的方向导数如何?

提示:沿x 轴正向时, cos α=1, cos β=0, x

f l f ??=??; 沿x 轴负向时, cos α=-1, cos β=0, x

f l f ??-=??. 沿y 轴正向时, cos α=0, cos β=1, y

f l f ??=??;

沿y 轴负向时, cos α=0, cos β=--1,

y

f l f ??-=??. 课堂练习:习题8-7:2,4 例1 求函数z =xe 2y 在点P (1, 0)沿从点P (1, 0)到点Q (2, -1)的方向的方向导数.

解 这里方向l 即向量→

)1 ,1(-=PQ 的方向, 与l 同向的单位向量为)2

1 ,21(-=l e . 因为函数可微分, 且1)0,1(2)0,1(==??y e x z

, 22)0,1(2)0,1(==??y xe y z ,

所以,所求方向导数为

22)2

1(2211)0,1(-=-?+?=??l z . 对于三元函数f (x , y , z )来说, 它在空间一点P 0(x 0, y 0, z 0)沿e l =(cos α , cos β , cos γ)的方向导数定义为

),,(000z y x l f

??t

z y x f t z t y t x f t ),,()cos ,cos ,cos (lim 0000000-+++=+→γβα. 如果函数f (x , y , z )在点(x 0, y 0, z 0)可微分, 则函数在该点沿着方向e l =(cos α , cos β , cos γ)的方向导数为

),,(000z y x l f

??=f x (x 0, y 0, z 0)cos α+f y (x 0, y 0, z 0)cos β+f z (x 0, y 0, z 0)cos γ.

例2 求f (x , y , z )=xy +yz +zx 在点(1, 1, 2)沿方向l 的方向导数, 其中l 的方向角分别为60?, 45?, 60?.

解 与l 同向的单位向量为

e l =(cos60?, cos 45?, cos60?))2

1 ,22

,21(=. 因为函数可微分, 且

f x (1, 1, 2)=(y +z )|(1, 1, 2)=3,

f y (1, 1, 2)=(x +z )|(1, 1, 2)=3,

,f z (1, 1, 2)=(y +x )|(1, 1)=2,

所以

)235(21212223213)

2,1,1(+=?+?+?=??l

f

. 二、梯度

1、梯度的概念 设函数z =f (x , y )在平面区域D 内具有一阶连续偏导数, 则对于每一点P 0(x 0, y 0)∈D , 都可确定一个向量

f x (x 0, y 0)i +f y (x 0, y 0)j ,

这向量称为函数f (x , y )在点P 0(x 0, y 0)的梯度, 记作grad f (x 0, y 0), 即

grad f (x 0, y 0)= f x (x 0, y 0)i +f y (x 0, y 0)j .

2、梯度与方向导数的关系:

如果函数f (x , y )在点P 0(x 0, y 0)可微分,e l =(cos α , cos β )是与方向l 同方向的单位向量, 则

),(00y x l f

??βαcos ),(cos ),(0000y x f y x f y x +=,

= grad f (x 0, y 0)?e l

=| grad f (x 0, y 0)|?cos(grad f (x 0, y 0),^ e l ).

这一关系式表明了函数在一点的梯度与函数在这点的方向导数间的关系. 特别, 当向量e l 与grad f (x 0, y 0)的夹角θ=0, 即沿梯度方向时, 方向导数),(00y x l f

??取得最大值, 这个最大值就是梯度的模|grad f (x 0, y 0)|. 这就是说: 函数在一点的梯度是个向量, 它的方向是函数在这点的方向导数取得最大值的方向, 它的模就等于方向导数的最大值.

讨论: l

f ??的最大值; 结论: 函数在某点的梯度是这样一个向量, 它的方向与取得最大方向导数的方向一致, 而它的模为方向导数的最大值.

3、梯度与等值线的关系:

我们知道, 一般说来二元函数z =f (x , y )在几何上表示一个曲面, 这曲面被平面z =c (c 是常数)所截得的曲线L 的方程为

?

??==c z y x f z ),(. 这条曲线L 在xOy 面上的投影是一条平面曲线L *, 它在xOy 平面上的方程为

f (x , y )=c .

对于曲线L *上的一切点, 已给函数的函数值都是c , 所以我们称平面曲线L *为函数z =f (x , y )的等值线.

若f x , f y 不同时为零, 则等值线f (x , y )=c 上任一点P 0(x 0, y 0)处的一个单位法向量为

)),(),,(()

,(),(10000002002y x f y x f y x f y x f y x y x +=n . 这表明梯度grad f (x 0, y 0)的方向与等值线上这点的一个法线方向相同, 而沿这个方向的方向导数n

f ??就等于|grad f (x 0, y 0)|, 于是 n n

f y x f ??=),(00grad . 这一关系式表明了函数在一点的梯度与过这点的等值线、方向导数间的关系. 这说是说: 函数在一点的梯度方向与等值线在这点的一个法线方向相同, 它的指向为从数值较低的等值线指向数值较高的等值线, 梯度的模就等于函数在这个法线方向的方向导数.

梯度概念可以推广到三元函数的情形. 设函数f (x , y , z )在空间区域G 内具有一阶连续偏导数, 则对于每一点P 0(x 0, y 0, z 0)∈G , 都可定出一个向量

f x (x 0, y 0, z 0)i +f y (x 0, y 0, z 0)j +f z (x 0, y 0, z 0)k ,

这向量称为函数f (x , y , z )在点P 0(x 0, y 0, z 0)的梯度, 记为grad f (x 0, y 0, z 0), 即

grad f (x 0, y 0, z 0)=f x (x 0, y 0, z 0)i +f y (x 0, y 0, z 0)j +f z (x 0, y 0, z 0)k .

结论: 三元函数的梯度也是这样一个向量, 它的方向与取得最大方

向导数的方向一致, 而它的模为方向导数的最大值.

如果引进曲面

f (x , y , z )=c

为函数的等量面的概念, 则可得函数f (x , y , z )在点P 0(x 0, y 0, z 0)的梯度的方向与过点P 0的等量面 f (x , y , z )=c 在这点的法线的一个方向相同, 且从数值较低的等量面指向数值较高的等量面, 而梯度的模等于函数在这个法线方向的方向导数.

例3 求2

21

y x +grad . 解 这里221),(y x y x f +=. 因为 222)(2y x x x f +-=??, 222)

(2y x y y f +-=??, 所以 221

y x +grad j i 222222)(2)(2y x y y x x +-+-=. 例4 设f (x , y , z )=x 2+y 2+z 2, 求grad f (1, -1, 2).

解 grad f =(f x , f y , f z )=(2x , 2y , 2z ),

于是 grad f (1, -1, 2)=(2, -2, 4).

数量场与向量场: 如果对于空间区域G 内的任一点M , 都有一个确定的数量f (M ), 则称在这空间区域G 内确定了一个数量场(例如温度场、密度场等). 一个数量场可用一个数量函数f (M )来确定, 如果与点M 相对应的是一个向量F (M ), 则称在这空间区域G 内确定了一个向量场(例如力场、速度场等). 一个向量场可用一个向量函数F (M )来确定, 而

F (M )=P (M )i +Q (M )j +R (M )k ,

其中P (M ), Q (M ), R (M )是点M 的数量函数.

利用场的概念, 我们可以说向量函数grad f (M )确定了一个向量场——梯度场, 它是由数量场f (M )产生的. 通常称函数f (M )为这个向量场的势, 而这个向量场又称为势场. 必须注意, 任意一个向量场不一定是势场, 因为它不一定是某个数量函数的梯度场.

例5 试求数量场

r

m 所产生的梯度场, 其中常数m >0,

222z y x r ++=为原点O 与点M (x , y , z)间的距离.

解 32)(r

mx x r r m

r m x -=??-=??, 同理 3)(r my r m y -=??, 3

)(r mz

r m z -=??. 从而 )(2k j i r z r

y r x r m r m ++-=grad . 记k j i e r

z r y r x r ++=, 它是与→OM 同方向的单位向量, 则

r r

m r m e 2-=grad . 上式右端在力学上可解释为, 位于原点O 而质量为m 质点对位于点M 而质量为l 的质点的引力. 这引力的大小与两质点的质量的乘积成正比、而与它们的距平方成反比, 这引力的方向由点M 指向原点. 因此数量场r m 的势场即梯度场grad r m 称为引力场, 而函数r

m 称为引力势. 课外作业: P51-2,4

方向导数和梯度

§ 3 方向导数和梯度 附录:数量场,向量场 数量场:设D 是n R 中的一个区域,f 是定义在D 内的一个实值函数,即R D f →:。则称在D 内有一个数量场f ,或称f 是D 内的数量场。 例如:教室中每一点的温度、位置等;点电荷形成的电位切; 磁铁周围磁力的大小. 等值面:设f 是D 内的一个数量场,称})({C x f D x s =∈= (C 是常数)是数量场f 的等值面,即在S 内每一点x 处,f 所对应的数值是相同的,都等于C. 特别当D 是2R 中的区域时,称S 等值线. 例如:天气预报中的等温面,等压面;地势图上的等高线(海拔相同). 向量场 一、 方向导数 1. 方向导数的定义 三元函数f 在点),,(0000z y x P 的三个偏导数,分别是函数f 在点),,(0000z y x P 沿着平行于坐标轴的直线方向(双向)上的变化率. 函数f 在点0P 沿射线l (单向)方向的变化率,即f 在点0P 沿方向l 的方向导数. 定义1(P124)设三元函数f 在点),,(0000z y x P 的某邻域)(0P ?3R 内有定义 , l 为从点0P 出 发的射线 . ),,(z y x P 为l 上且含于)(0P 内的任一点 , 以ρ表示P 与0P 两点间的距离 . 若极限 ρρρρf P f P f l ?=-++→→000lim )()(lim 存在 , 则称此极限为函数f 在点0P 沿方向l 的方向导数 , 记为0P l f ?? 或)(0P f l 、),,(000z y x f l . 定义1' 设D 是3R 中的一个区域,f 是D 内的一个数量场,D P ∈0,l 是3R 中的一个单位向量,即,1=l 如果t P f tl P f t )()(lim 000-++ →,存在,则称此极限是数量场f 在点0P 沿方向l 的方向导数,记为)(0P l f ??,即t P f tl P f P l f t )()(lim )(0000-+=??+→。也称它是函数f 在点0P 沿方向l 的方向导数,它表示数量场f 在点0P 沿方向l 的变化率。 易见,x f ??、y f ?? 和 z f ??是三元函数f 在点0P 分别沿X 轴正向、Y 轴正向和Z 轴正向的方向导数. 对二元函数),(y x f z =在点),(000y x P , 可仿此定义方向导数 . 根据定义计算方向导数 例1 (P125,有补充)),,(z y x f =32z y x ++. 求f 在点0P ) 1 , 1 , 1 (处沿l 方向的方向导数,

vasp做分子动力学

vasp做分子动力学的好处,由于vasp是近些年开发的比较成熟的软件,在做电子scf速度方面有较好的优势。 缺点:可选系综太少。 尽管如此,对于大多数有关分子动力学的任务还是可以胜任的。 主要使用的系综是NVT和NVE。 下面我将对主要参数进行介绍! 一般做分子动力学的时候都需要较多原子,一般都超过100个。 当原子数多的时候,k点实际就需要较少了。有的时候用一个k点就行,不过这都需要严格的测试。通常超过200个原子的时候,用一个k点,即Gamma点就可以了。 INCAR: EDIFF 一般来说,用1E-4或者1E-5都可以,这个参数只是对第一个离子步的自洽影响大一些,对于长时间的分子动力学的模拟,精度小一点也无所谓,但不能太小。 IBRION=0 分子动力学模拟 IALGO=48 一般用48,对于原子数较多,这个优化方式较好。 NSW=1000 多少个时间步长。 POTIM=3 时间步长,单位fs,通常1到3. ISIF=2 计算外界的压力. NBLOCK= 1 多少个时间步长,写一次CONTCAR,CHG和CHGCAR,PCDAT. KBLOCK=50 NBLOCK*KBLOCK个步长写一次XDATCAR. ISMEAR=-1 费米迪拉克分布. SIGMA =0.05 单位:电子伏 NELMIN=8 一般用6到8,最小的电子scf数.太少的话,收敛的不好. LREAL=A APACO=10 径向分布函数距离,单位是埃. NPACO=200 径向分布函数插的点数. LCHARG=F 尽量不写电荷密度,否则CHG文件太大. TEBEG=300 初始温度. TEEND=300 终态温度。不设的话,等于TEBEG. SMASS -3 NVE ensemble;-1 用来做模拟退火;大于0 NVT 系综。 ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// 1)收敛判据的选择 结构弛豫的判据一般有两种选择:能量和力。这两者是相关的,理想情况下,能量收敛到基态,力也应该是收敛到平衡态的。但是数值计算过程上的差异导致以二者为判据的收敛速度差异很大,力收敛速度绝大部分情况下都慢于能量收敛速度。这是因为力的计算是在能量的基础上进行的,能量对坐标的一阶导数得到力。计算量的增大和误差的传递导致力收敛慢。 到底是以能量为收敛判据,还是以力为收敛判据呢?关心能量的人,觉得以能量

14方向导数与梯度

方向导数与梯度第六节 经常需要研究函数在某点沿某一固定方向的变化率问在实际问题中,z?),xyz?f(实际上就是点的偏导数题,例如我们所学习的函数x?x)x,yP(沿轴方向变化时函数的变化率,由此引入方向导数的概念。一、方向导数PQ方向的变不难看出函数沿我们以二元函数为例介绍方向导数。化率可以用如下极限表示)f(Pf(Q)?lim||PQ0|PQ|?l)f(x,yz?)(x,yP)UP(内有定义,设函数在点的某邻域为一000e),bai?bj?(ae?LP相引射线,方向与,自点向量,其单位向量为exOy),yP(xL为方向向量的直线,由解同,由于面上通过点是且以00L析几何知射线的参数方程为tax??x?0??t?0?。 ?tb??yy?0ta?x?x?0)x,y(QL,则上任意取一点。在 ?y?y?tb?0Qetb)?,(?)?,xx?PQ(?yytaP两点间的距,所以到由于001 t?|?|t|?|PQ||t(a,b)离为 y L Q

P x O )y?f(x,z)P(x,ye的变化率我们可以用函数则函数在点处沿方向 00Q)(P(fQ)?ftP增量到点的比值与点的距离)fy(x,,f(x?tay?tb)?)(P)f(Q?f0000? tt?0?tQlP)时的极限来表示,该极限为函数当趋于点沿直线(即e)?zf(x,yP处沿方向的变化率,称为方向导数。在点l)(xy,?zf)x,Py(的某个邻域内有定义,设函数定义是在点00)b,(?ea 为,如果极限一非零向量,其单位向量l2 )y(x,(x?ta,y?tb)?ff0000lim t?0?t l)yz?f(x,),Py(x的方向导在点处沿方向存在,则称此极限为函数00f?|数,记作,即)x(,y 00l?),yy?tb)?f(xtaf(x?,f?0000|lim?。)yx,( t00?l?0t?f?),yz?f(x|就是函数由方向导数的定义可知,方向导数),y(x00l?l)Py(x,处沿方向在点的变化率。00f?),0?(1e?i|则在地特别,如果,存取,且l)x,y(x?00f?f??f),1j?(0e?|?||则存在取,。如果且,l)yx)(x,y,(),xy(x?x?000000l?f??f|?|。注意,反之未然。))y,(x,(xy y?0000l?方向导数的计算本质上仍然是一元函数导数的计算,因为若令 ?)?tb?ta,y?(t)f(x,则00??)xf(,y)y?f(xta,?tb?)?0(t()0000limlim?, tt??00t?t?3

GIS基本概念集锦

GIS基本概念集锦 1、地理信息系统(geographic information system ,即gis )——一门集计算机科学、信息学、地理学等多门科学为一体的新兴学科,它是在计算机软件和硬件支持下,运用系统工程和信息科学的理论,科学管理和综合分析具有空间内涵的地理数据,以提供对规划、管理、决策和研究所需信息的空间信息系统。gis有以下子系统:数据输入子系统,数据存储和检索子系统,数据操作和分析子系统,报告子系统. 信息系统 非空间的空间的 管理信息系统非地理学的gis cad/cam 其他gis lis 社会经济,人口普查基于非地块,基于地块的 2、比较gis与cad、cac间的异同。 cad——计算机辅助设计,规则图形的生成、编辑与显示系统,与外部描述数据无关。 cac——计算机辅助制图,适合地图制图的专用软件,缺乏空间分析能力。 gis——地理信息系统,集规则图形与地图制图于一身,且有较强的空间分析能力。 3、图层:将空间信息按其几何特征及属性划分成的专题。 4、地理数据采集——实地调查、采样;传统的测量方法,如三角测量法、三边测量法;全球定位系统(g ps);现代遥感技术;生物遥测学;数字摄影技术;人口普查。 5、信息范例——传统的制图方法,称为信息范例,即假定地图本身是一个最终产品,通过使用符号、分类限制的选择等方式交换空间信息的模式。这个范例是传统的透视图方法,由于原始而受到很多限制,地图用户不能轻易获得预分类数据。也就是说,用户只限于处理最终产品,而无法将数据重组为更有效的形式以适应环境或需求的变化。 6、分析范例(整体范例)——存储保存原始数据的属性数据,可根据用户的需求进行数据的显示、重组和分类。整体范例是一种真正的用于制图学和地理学的整体方法。 7、栅格——栅格结构是最简单最直接的空间数据结构,是指将地球表面划分为大小均匀紧密相邻的网格阵列,每个网格作为一个象元或象素由行、列定义,并包含一个代码表示该象素的属性类型或量值,或仅仅包括指向其属性记录的指针。因此,栅格结构是以规则的阵列来表示空间地物或现象分布的数据组织,组织中的每个数据表示地物或现象的非几何属性特征。特点:属性明显,定位隐含,即数据直接记录属性本身,而所在的位置则根据行列号转换为相应的坐标,即定位是根据数据在数据集中的位置得到的,在栅格结构中,点用一个栅格单元表示;线状地物用沿线走向的一组相邻栅格单元表示,每个栅格单元最多只有两个相邻单元在线上;面或区域用记有区域属性的相邻栅格单元的集合表示,每个栅格单元可有多于两个的相邻单元同属一个区域。 8、矢量——它假定地理空间是连续,通过记录坐标的方式尽可能精确地表示点、线、多边形等地理实体,坐标空间设为连续,允许任意位置、长度和面积的精确定义。对于点实体,矢量结构中只记录其在特定坐标系下的坐标和属性代码;对于线实体,用一系列坐标对的连线表示;多边形是指边界完全闭合的空间区域,用一系列坐标对的连线表示。 9、“拓扑”(topology)一词来源于希腊文,它的原意是“形状的研究”。拓扑学是几何学的一个分支,它研究在拓扑变换下能够保持不变的几何属性——拓扑属性(拓扑属性:一个点在一个弧段的端点,一个点在一个区域的边界上;非拓扑属性:两点之间的距离,弧段的长度,区域的周长、面积)。这种结构应包括:唯一标识,多边形标识,外包多边形指针,邻接多边形指针,边界链接,范围(最大和最小x、y坐标值)。地理空间研究中三个重要的拓扑概念(1)连接性:弧段在结点处的相互联接关系;(2)多边形区域定义:多个弧段首尾相连构成了多边形的内部区域;(3)邻接性:通过定义弧段的左右边及其方向性

模糊C均值聚类 及距离函数的优缺点

K-均值聚类分析是一种硬划分,它把每一个待识别的对象严格的划分到某一类中,具有非此即彼的性质。而实际上高光谱值目标在形态和类属方面存在着中介性,没有确定的边界来区分。因此需要考虑各个像元属于各个类别的隶属度问题,进行软划分,从而更好的区分。 设要进行聚类分析的图像像元数N ,图像像元集合{}N x x x X ,...,,21=,其中 {} T p k k k k x x x x ,...,,2 1=,p 为波段数。设把图像分为C 个类别,每个类别的聚类中心 ),...,,(21p i i i i v v v v =,聚类中心集合{}c v v v V ,...,,21=。用ik u 表示像元k x 隶属于以i v 为中心的类别i 的隶属度,定义隶属度矩阵U 为[]N C ik u U ?=。 矩阵U 中每一列的元素表示所对应的像元隶属于C 个类别中各个类的隶属度。满足一下约束条件: ???? ?????≤≤===>∑∑==10,...2,1;,...,2,1. (101) 1ik C i ik N k ik u N k C i u u 对隶属度ik u 进行了迷糊化,ik u 可取0和1之间的任意实数,这样一个像元可以同时隶 属于不同的类别,但其隶属度的总和总是等于1,这符合高光谱像元的实际情况。而属于硬聚类的K-均值聚类,其隶属度具有非此即彼的性质,隶属度ik u 只能取0或1。 定义目标函数J 为 ∑∑==?=N k C i ik m ik m d u V U J 11 2)()(),( 22)(i k ik v x d -=为Euclidean 距离;[)∞∈,1m 为模糊加权指数(当m=1时,同K- 均值的目标函数一致)。最优化的类就是使目标函数取最小值的类,如果一类中的所有点都 贴近于它们的类中心,则目标函数很小。 FCM 算法步骤: (1) 确定聚类数C ,加权指数m ,终止误差ε,最大迭代次数LOOP 。 (2) 初始化隶属度矩阵) 0(U (3) 开始循环,当迭代次数为IT (IT=0,1,2…,C )时,根据) (IT U 计算C-均值向量, 即C i u x u U N k N k m ik k m ik IT ,...,2,1],))(/()([1 1 ) (==∑∑== (4) 对k=1,2,…,N ,按以下公式更新) (IT U 为) 1(+IT U : 若i k v x ≠对所有的i v (i=1,2,…,C)满足,则对此k x 计算 C i v x d d u i k C j m jk ik ik ,...,2,1,,])([1112 =≠=-=-∑ 若对某一个i v ,有k x 满足i k v x =,则对应此k x ,令)(0;1i j u u jk ik ≠==。这样把聚 类中心与样本一致的情形去掉,把隶属度模糊化为0和1之间的实数。 (5) 若ε<+)1() (-IT IT U U 或IT>LOOP ,停止;否则置IT=IT+1,并返回第三步。 FCM 算法允许自由选取聚类个数,每一向量按其指定的隶属度]1,0[∈ik u 聚类到每一聚类中心。FCM 算法是通过最小化目标函数来实现数据聚类的。

数学建模习题答案

数学建模部分课后习题解答 中国地质大学 能源学院 华文静 1.在稳定的椅子问题中,如设椅子的四脚连线呈长方形,结论如何? 解: 模型假设 (1) 椅子四条腿一样长,椅脚与地面接触处视为一点,四脚的连线呈长方形 (2) 地面高度是连续变化的,沿任何方向都不会出现间断(没有像台阶那样的情况), 即从数学角度来看,地面是连续曲面。这个假设相当于给出了椅子能放稳的必要条件 (3) 椅子在任何位置至少有三只脚同时着地。为了保证这一点,要求对于椅脚的间 距和椅腿的长度而言,地面是相对平坦的。因为在地面上椅脚间距和椅腿长度的尺寸大小相当的范围内,如果出现深沟或凸峰(即使是连续变化的),此时三只脚是无法同时着地的。 模型建立 在上述假设下,解决问题的关键在于选择合适的变量,把椅子四只脚同时着地表示出来。首先,引入合适的变量来表示椅子位置的挪动。生活经验告诉我们,要把椅子通过挪动放稳,通常有拖动或转动椅子两种办法,也就是数学上所说的平移与旋转变换。然而,平移椅子后问题的条件没有发生本质变化,所以用平移的办法是不能解决问题的。于是可尝试将椅子就地旋转,并试图在旋转过程中找到一种椅子能放稳的情形。 注意到椅脚连线呈长方形,长方形是中心对称图形,绕它的对称中心旋转180度后,椅子仍在原地。把长方形绕它的对称中心旋转,这可以表示椅子位置的改变。于是,旋转角度θ这一变量就表示了椅子的位置。为此,在平面上建立直角坐标系来解决问题。 设椅脚连线为长方形ABCD,以对角线AC 所在的直线为x 轴,对称中心O 为原点,建立平面直角坐标系。椅子绕O 点沿逆时针方向旋转角度θ后,长方形ABCD 转至A1B1C1D1的位置,这样就可以用旋转角)0(πθθ≤≤表示出椅子绕点O 旋转θ后的位置。 其次,把椅脚是否着地用数学形式表示出来。当椅脚与地面的竖直距离为零时,椅脚就着地了,而当这个距离大于零时,椅脚不着地。由于椅子在不同的位置是θ的函数,因此,椅脚与地面的竖直距离也是θ的函数。 由于椅子有四只脚,因而椅脚与地面的竖直距离有四个,它们都是θ的函数,而由假设(3)可知,椅子在任何位置至少有三只脚同时着地,即这四个函数对于任意的θ,其函数值至少有三个同时为0。因此,只需引入两个距离函数即可。考虑到长方形ABCD 是对称中心图形,绕其对称中心O 沿逆时针方向旋转180度后,长方形位置不变,但A,C 和B,D 对换了。因此,记A ,B 两脚与地面竖直距离之和为)(θf ,C,D 两脚之和为)(θg ,其中[] πθ,0∈,使得)()(00θθg f =成立。 模型求解 如果0)0()0(== g f ,那么结论成立。

(整理)第七节方向导数与梯度

第七节 方向导数与梯度 要求:了解方向导数与梯度的概念,会计算方向导数与梯度方法。 重点:方向导数与梯度的计算。 难点:梯度的几何意义,方向导数与梯度的联系。 作业:习题8-7(60P )2,4,6,8,10 一.方向导数 问题提出:在许多实际问题中,常常需要知道函数),(y x f z =在点(,)P x y 沿任意方向或某个方向的变化率.例如预报某地的风向和风力就必须知道气压在该处沿着哪个方向的变化率,在数学上就是多元函数在一点沿给定方向的方向导数问题. 1.方向导数定义 设函数),(y x f z =在点(,)P x y 的某一邻域内有定义,自P 点引有向直线L ,x 轴正向与直线L 夹角为?,在L 上任取一点'(,)P x x y y +?+?,若'P 沿着L 趋近于P 时,即当0)()(22→?+?= y x ρ时,极限 ρ ρ) ,(),(lim y x f y y x x f -?+?+→ 存在 则称此极限值为函数在点P 沿着L 方向的方向导数.记作 ρ ρ),(),(lim 0y x f y y x x f L f -?+?+=??→. 说明 (1)规定逆时针方向旋转生成的角是正角0>?,顺时针方向旋转生成的角是负角 0

(整理)第七节方向导数与梯度讲课稿

(整理)第七节方向导 数与梯度

第七节 方向导数与梯度 要求:了解方向导数与梯度的概念,会计算方向导数与梯度方法。 重点:方向导数与梯度的计算。 难点:梯度的几何意义,方向导数与梯度的联系。 作业:习题8-7(60P )2,4,6,8,10 一.方向导数 问题提出:在许多实际问题中,常常需要知道函数),(y x f z =在点(,)P x y 沿任意方向或某个方向的变化率.例如预报某地的风向和风力就必须知道气压在该处沿着哪个方向的变化率,在数学上就是多元函数在一点沿给定方向的方向导数问题. 1.方向导数定义 设函数),(y x f z =在点(,)P x y 的某一邻域内有定义,自P 点引有向直线L , x 轴正向与直线L 夹角为?,在L 上任取一点'(,)P x x y y +?+?,若'P 沿着L 趋 近于P 时,即当0)()(22→?+?=y x ρ时,极限 ρ ρ) ,(),(lim y x f y y x x f -?+?+→ 存在 则称此极限值为函数在点P 沿着L 方向的方向导数.记作 ρ ρ),(),(lim 0y x f y y x x f L f -?+?+=??→. 说明 (1)规定逆时针方向旋转生成的角是正角0>?,顺时针方向旋转生成的角是负角0

2.方向导数的计算 定理 若函数),(y x f z =在点(,)P x y 可微分,那么函数),(y x f z =在点 (,)P x y 沿任一方向L 的方向导数都存在,且有计算公式 ??sin cos y f x f L f ??+??=??{},cos ,sin ,f f f f e x y x y ??????????=?=????????????? r . 其中?为x 轴到方向L 的转角,e r 是与L 同方向的单位向量. 证明:因为函数),(y x f z =在点(,)P x y 可微分,所以有 ()f f f x y o x y ρ???= ?+?+??, 上式两边同除以ρ,得 ()() cos sin f f x f y o f f o x y x y ρρ??ρ ρρρρ ???????= ++=++????,则 0lim cos sin f f f f L x y ρ??ρ→????==+??? 例1.求函数y xe z 2=在点(1,0)P 处沿从点(1,0)P 到点)1,2(-Q 的方向的方向导数. 解 这里方向L 即向量{}1,1PQ =-u u u r 的方向,因此x 轴到L 方向的转角4π?=, 又因为y e x z 2=??,y xe y z 22=??,所以在点)0,1(处,1=??x z ,2=??y z , 于是方向导数为 2 2)4sin(2)4cos(1-=-+-?=??ππL z . 另一方法.

空间数据分析教学内容

间数据分析 1. 空间分析:(spatial analysis,SA)是基于地理对性的位置和形态特征的空间数据分析技术,其目的在于提取和传输空间信息,是地理信息系统的主要特征,同时也是评价一个地理信息系统功能的主要指标之一,是各类综合性地学分析模型的基础,为人们建立复杂的空间应用模型提供了基本方法. 2. 空间分析研究对象:空间目标。空间目标基本特征:空间位置、分布、形态、空间关系(度量、方位、拓扑)等。 3. 空间分析根本目标:建立有效地空间数据模型来表达地理实体的时空特性,发展面向应用的时空分析模拟方法,以数字化方式动态的、全局的描述的地理实体和地理现象的空间分布关系,从而反映地理实体的内在规律和变化趋势。GIS空间分析实际是一种对GIS海量地球空间数据的增值操作。 4. ArcGIS9中主要的三种数据组织方式:shapefile,coverage和geodatabase。Shapefile由存储空间数据的dBase表和存储属性数据和存储空间数据与属性数据关系的.shx文件组成。Coverage的空间数据存储在INFO表中,目标合并了二进制文件和INFO表,成为Coverage 要素类。 5. Geodatabase是面向对象的数据模型,能够表示要素的自然行为和要素之间的关系。 6. GIS空间分析的基本原理与方法:根据空间对象的不同特征可以运用不同的空间分析方法,其核心是根据描述空间对象的空间数据分析其位置、属性、运动变化规律以及周围其他对象的相关制约,相互影响关系。方法主要有矢量数据的空间分析,栅格数据的空间分析,空间数据的量算与空间内插,三维空间分析,空间统计分析。 7. 栅格数据在数据处理与分析中通常使用线性代数的二维数字矩阵分析法作为数据分析的数学基础。栅格数据的处理方法有:栅格数据的聚类、聚合分析,复合分析,追踪分析,窗口分析。 8. 栅格数据的聚类与聚合分析区别:聚类是根据设定的聚类条件对原有的数据系统进行有选择的信息提取儿建立的新的栅格数据系统的方法;聚合分析是根据空间分辨率和分类表进行数据类型的合并或转换以实现空间地域的兼并。e.g.:从遥感图像信息中提取某一地物的方法是栅格数据的聚类,而由数字高程模型转换为数字高程分级模型便是空间数据的聚合。 9. ArcGIS9的空间分析功能主要包括:空间分析模块、3D分析模块、地统计分析模块、网络分析模块、跟踪分析模块等。 10.GIS的四种模型:要素模型(矢量类型),场模型(栅格类型),时态模型,网络模型。 11.矢量数据的空间分析方法:空间关系查询;叠置分析;缓冲区分析;泰森多边形分析;网络分析。 12.空间关系查询中要素间的关系有:相邻关系(proximity);包容关系(containment);叠加关系(overlap)。空间关系查询涉及到的:目标层—从其中查询满足条件的要素。选择层—比较此图层中要素与目标层中的关系。选择层是从目标层中得来的。 13.空间关系连接(spatial join)有:根据图层间的关系连接属性表;根据空间位置连接图层属性表。 14.叠置分析:是将代表不同主题的各个数据层面进行叠置产生一个新的数据层面,叠置的结果综合了原来来年各个或多个层面要素所具有的属性,不仅生成了新的空间关系,还将输入的多个数据层的属性联系起来产生了新的属性关系。叠置分析前提条件:要素层面必须基于相同坐标系统、同一地带,还必须查验叠加层面之间的基准面是否相同。 15.根据操作要素的不同,叠置分析可以分为:点与多边形叠加;线与多边形叠加;多边形

使用距离函数的逐像素位移映射

第8章 使用距离函数的逐像素位移映射 William Donnelly Waterloo大学 在本章中,我们将演示distance map,一种在pixel shader中给对象添加小尺度位移映射的技术。我们把位移映射看作一个光线跟踪的问题,从基表面的纹理坐标开始,计算观察射线相交于偏移后表面的纹理坐标。为了达到这个目的,我们预计算一个三维distance map,它给定了空间中的点和偏移后表面距离的量度。这个distance map给了我们所有快速计算射线和表面相交所需的信息。我们的算法在保持了实时性能的情况下显著地提升了场景可以感知的几何复杂度。 8.1 导读 Cook(1984)引入了位移映射的方法,可以用于向表面增加小尺度细节。不像凹凸映射只影响表面的着色,位移映射调整了表面元素的位置。这可以形成凹凸映射不可能做到的效果,比如互相遮挡的表面特征和非多边形轮廓。图8-1演示了石头表面的渲染,互相的遮挡的特征造成了纵深的感觉。 位移映射一般的实现方式是反复镶嵌一个基表面,沿着基表面的法线挤出顶点,然后重复这个过程,直到产生的多边形接近于一个像素的大小。Michael Bunnell在本书的第7章演示了这个方法,“用于位移映射的自适应表面细分镶嵌”。

Figure 8-1. 一个偏移的石头表面 位移映射(上)给了凹凸映射(下)不可能做到的纵深的感觉。 虽然在流水线的顶点着色阶段实现位移映射看起来更合逻辑,但因为vertex shader没有产生新顶点的能力,所以那是不可能的。这时我们考虑基于pixel shader的技术。使用pixel shader有几个好处: ?目前GPU的像素处理能力远大于顶点处理能力。例如,GeForce 6800 Ultra有16个像素着色流水线,但只有6个顶点着色流水线。另外,一 个像素流水线在每个时钟周期经常可以比顶点流水线进行更多的操作。 ?pixel shader访问纹理更容易。很多GPU不允许在vertex shader中访问纹理;而对于允许的,访问模式也有限制,而且会有性能开销。 ?像素处理的数量随着距离而放缩。vertex shader总是对模型中的每个顶点执行一次,但pixel shader只对屏幕上的每个像素执行一次。这意味 着工作集中于临近的对象,当然它们也是最重要的。 使用pixel shader的缺点是我们不能在shader中改变像素的屏幕坐标。这意味着不像基于镶嵌的方法,使用pixel shader的方法不能用于任意大小的距离。但是,这也不是一个严格的限制,因为在实际中位移几乎总是有限的。 8.2 以前的工作 Parallax mapping是扩展bump mapping的一种简单方式,可以包含视差效果(Kaneko等,2001)。Parallax mapping使用关于表面一个点的高度信息来使纹理坐标靠近或远离观察者。虽然这是一个仅对平滑变化的高度场有效的粗糙近似,但它给出了惊人的好结果,特别是它只需增加三条shader指令就能实现。不幸的是,由于这项技术的内在问题,它不能处理大位移、高频特征或特征间的遮挡。

ArcGIS教程:半变异函数与协方差函数

半变异函数和协方差函数将邻近事物比远处事物更相似这一假设加以量化。半变异函数和协方差都将统计相关性的强度作为距离函数来测量。 对半变异函数和协方差函数建模的过程将半变异函数或协方差曲线与经验数据拟合。目标是达到最佳拟合,并将对现象的认知纳入模型。之后模型便可用于预测。 在拟合模型时,浏览数据中的方向自相关。基台、变程和块金是模型的重要特征。如果数据中有测量误差,请使用测量值误差模型。跟踪这一链接来了解如何将模型与经验半变异函数拟合。 半变异函数 半变异函数定义为 γ(si,sj) = ? var(Z(si) - Z(sj)), 其中 var 是方差。 如果两个位置 si 和 sj,在 d(si, sj) 的距离测量上彼此相近,那么您会希望这两个位置相似,以便缩小两个位置的差值 Z(si) - Z(sj) 的大小。当 si 和 sj 距离逐渐增大时,它们变得越来越不相似,它们的值 Z(si) - Z(sj) 的差异也会增大。在下图中可以看到这一情况,其中显示了典型半变异函数的解析图。 请注意,差值的方差会随距离的增大而增加,因此可以将半变异函数视为相异度函数。与这一函数经常关联的术语也可用在 Geostatistical Analyst 中。半变异函数在其呈平稳状态时所达到的高度称为基台。它通常由两部分组成:原点处不连续(称为块金效应)和偏基台;二者一起形成基台。块金效应可以细分为测量误差和微刻度变化。块金效应就是测量误差和微尺度变化的和,由于任一组件都可为零,因此块金效应可以完全由一个组件或另一个组件形成。变程是半变异函数达到平稳基台处的距离。

协方差函数 协方差函数定义为 C(si, sj) = cov(Z(si), Z(sj)), 其中 cov 是协方差。 协方差是相关性的缩放版。因此当两个位置,si 和 sj 彼此相近时,您会希望这两个位置相似,而他们的协方差(相关性)会变大。当 si 和 sj 距离逐渐增大时,它们变得越来越不相似,并且它们的协方差会变为零。在下图中可以看到这一情况,下图显示典型协方差函数的解析图。 请注意,协方差函数随距离的增大而减小,因此可将其视为一种相似度函数。 半变异函数和协方差函数之间的关系 在半变异函数和协方差函数之间存在以下关系: γ(si, sj) = sill - C(si, sj), 从图中可看出该关系。由于这一相等关系,您可以在 Geostatistical Analyst 中使用两种函数中的任一种来执行预测。(Geostatistical Analyst 中所有半变异函数都拥有基台。) 半变异函数和协方差不是任意函数皆可。为使预测具有非负的克里金标准误差,只有部分函数可以用作半变异函数和协方差。Geostatistical Analyst 提供了多种可接受的选项,您可以为数据尝试不同的选项。您也可以通过同时添加多个模型的方式获得模型 - 此构造提供有效的模型,可以在 Geostatistical Analyst 中添加其中的最多四个模型。有一些当半变异函数存在时,协方差函数却不存在的实例。例如,有一个线性半变异函数,但它没有基台,并且没有相对应的协方差函数。Geostatistical Analyst 中仅使

导数-偏导数-方向导数-梯度及其关系

导数:()()()00' 00 0lim lim x x f x x f x y f x x x ? →?→+?-?==??,导数的意义为函数的变化率。由定义可知,导数是对应一元函数的。 偏 导 数 : ()()() 0000000 ,,,lim x x f x x y f x y f x y x ?→+?-=? ()()() 0000000 ,,,lim y y f x y y f x y f x y y ?→+?-=? 偏导数是对应于多元函数的。其意义是:偏 导数反应的是函数沿坐标轴方向的变化率。 方向导数:设l 为xOy 平面上以()000,P x y 为始发点的一条射线,()cos ,cos l αβ=e 是与l 同方向的单位向量。则该射线的参数方程为: 00cos cos x x t y y t αβ =+=+ ,那么,函数(,)f x y ,在 ()000,P x y 沿l 方向的方向导数为: () ()() 0000000 ,cos ,y cos ,lim t x y f x t t f x y f l t αβ+ →++-?=?。 从方向导数的定义可知,方向导数 () 00,x y f l ??就是函数(,)f x y 在点()000,P x y 沿方向l 的变 化率。方向导数也是对应于多元函数的。方向导数是一个标量值。 方向导数与偏导数的关系:如果函数(,)f x y 在点()000,P x y 可微分,那么函数在改点沿任意方向l 的方向导数存在,且有 () ()()000000,,cos ,cos x y x y f f x y f x y l αβ?=+? ,其中 ()cos ,cos l e αβ=为方向l 的方向余弦。 (若方向()1,0l =e 也就是x 轴方向,则() 0000,(,)x x y f f x y l ?=?,若方向()0,1l =e 也就是y 轴方向,则 () 0000,(,)y x y f f x y l ?=?) 梯度:设函数(,)f x y 在平面区域D 内有一阶连续偏导数,则对于每一个点()000,P x y D ∈都可以定出一个向量()()0000,,x y f x y f x y +i j ,这向量称为函数(,)f x y 在点()000,P x y 的梯度,即()()()000000 ,,,x y f x y f x y f x y =+grad i j 。

4.7—方向导数与梯度

第七节 方向导数与梯度 教学目的:(1) 理解方向导数和梯度的概念; (2) 掌握方向导数和梯度的计算方法。 教学重点:方向导数和梯度的计算 教学难点:方向导数和梯度的概念 教学方法:讲练结合 教学时数:2课时 一、问题的提出 实例:一块长方形的金属板,四个顶点的坐标是(1,1),(5,1),(1,3),(5,3).在坐标原点处有一个火焰,它使金属板受热.假定板上任意一点处的温度与该点到原点的距离成反比.在(3,2)处有一个蚂蚁,问这只蚂蚁应沿什么方向爬行才能最快到达较凉快的地点? 问题的实质:应沿由热变冷变化最骤烈的方向(即梯度方向)爬行. 二、方向导数 由偏导数的定义知,z x ??是),(y x f z =沿x 轴方向的变化率, z y ??是),(y x f z =沿y 轴 方向的变化率. 问题: 讨论函数),(y x f z =在一点0P 沿其他方向的变化率. 1、定义:设函数),(y x f z =在点0P 的某一邻域)(P U 内有定义,l 为非零向量,其方向 角为αβ和. 若极限 00000(cos ,cos )(,) lim f x y f x y ρραρβρ + →++-存在,则称 这极限为函数),(y x f z =在点0P 沿方向l 的方向导数。 记为 00000 (cos ,cos )(,) lim P f x y f x y f l ρραρβρ + →++-?=? 说明: (1). 方向导数的几何意义:方向导数 P f l ??就是函数),(y x f z =在点0P 沿方向l 的变化率. (2).依定义,函数),(y x f z =在点0P 处的偏导数存在时,则函数),(y x f 在点0P 沿着x 轴正向}0,1{1=e 、y 轴正向}1,0{2=e 的方向导数分别为y x f f ,;沿着x 轴负向、y 轴负向 o y x l (,) P x y ?? 000(,) P x y l e

相关文档
最新文档