




已阅读5页,还剩14页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1(一)函数1 凹(凸)函数1.1 凸集凸集:对于任意两点 和 ,且对于每一个 ,当且仅当uSv0,1为真时,集合 为凸集。(1)wuvnR凸集要求集合内两点之间的连线必须也在集合内,即该集合不存在任何孔,它的边缘也不能有缩进。例如,平面中,一条线段就是一个凸集,而一个圆圈则不是。1.2 凹(凸)函数介绍凸集是为了引入凹(凸)函数:不管是凹函数还是凸函数都要求其定义域是凸集。我们可以先举个例子直观感受下凹(凸)函数的特征,比如函数 就是一24yx个凹函数,它在定义域内呈现出峰形;函数 就是一个凸函数,它在定义域24yx内呈现谷底。现在具体给出凹(凸)函数的定义:对于函数 ,其定义域内任意两个不同的点 和 ,当且仅当:fDR12x121(x)(x)()(0)ttffttt时,函数 f 为凹函数。对于函数 ,其定义域内任意两个不同的点 和 ,当且仅当: 1x2121(x)(x)()(0)tftffttt时,函数 f 为凸函数。若将不等号“ ” 和“ ”分别变换成严格不等号“ ”和“ ”,上述定义便成了严格凹函数和严格凸函数的定义。因为凹函数的定义域为凸集,因此点 也一定在函数的定义域内。12x()tt我们可以利用凹(凸)函数和严格凹(凸)函数判断函数极值的情况。凹函数一定存在绝对极大值,但绝对极大值可能不是唯一的,因为如果山峰包含一个平顶,则可能存在多重绝对极大值。仅当我们限定它为严格凹形函数时,绝对值才可能是唯一的。1.3 凹(凸)函数与凸集的关系首先我们必须区别凸集与凸函数的概念。根据定义,可知当“凸的”在描述集合时,它要求该集合不能出现任何孔,边缘也不能有缩进。这不同于之前的凹(凸)函数:当“凸的”在描述函数时,它确定的是一条曲线或曲面是如何弯曲的。但凹(凸)函数确实与凸集有关。除了定义域都要求是凸集之外,它们都可以引致一个凸集。定理是凹函数 是凸集;(x)f(x),(x)AyDfy,Comment z1: 这里有个问题,我的概念有些模糊,二元的效用函数下,介无差异曲线凸向原点也保证其上等值集为凸集,感觉那时上等值集是平面上可以画出来的。但现在这里的上等值集为凸集应该是三维的,两个凸集有关系吗?2是凸函数 是凸集。(x)f(x),()AyDfy,即,由函数上的点以及函数曲线(曲面)之下的点组成的集合若是凸集 该函数为凹函数;由函数上的点以及函数曲线(曲面)之上的点组成的集合若是凸集 该函数为凸函数。2拟凹(拟凸)函数不管是凹(凸)函数还是严格凹(凸)函数,它们对函数都有比较强的设定。但是通常,理论研究的工作之一是为保证获得结果,识别出我们需要对函数进行的最弱的可行设定。拟凹(拟凸)函数则是一个相对而言更弱的条件。拟凹(拟凸)函数的定义如下:对于函数 ,其定义域内任意两个不同的点 和 ,当且仅当:fDR1x2121min(x),()(0,)fftforalt时,函数 f 为拟凹函数。对于函数 ,其定义域内任意两个不同的点 和 ,当且仅当: 12121a(),(x)(,)fftfrlt时,函数 f 为拟凸函数。若将不等号“ ” 和“ ”分别变换成严格不等号“ ”和“ ”,上述定义便适用于严格拟凹函数和严格拟凸函数的定义。我们也可以通过更直观的方法检验函数的拟凹性和拟凸性。设 为函数 在 水平上的上等值集,00()x,()SyDfy(x)f0y为函数 在 水平上的下等值集。I (x)f0定理对于值域内的所有 y值, 都是凸集 是拟凹函数S:fDR对于值域内的所有 y值, 都是凸集 是拟凸函数()I经济学中常假设拟凹的效用函数。根据定理,拟凹的效用函数保证了其上等值集为凸集。3函数间关系(1) 是(严格)凹函数 是(严格)凸函数;(x)f (x)f(2) 是(严格)拟凹函数 是(严格)拟凸函数;(3) 是(严格)凹函数 是(严格)拟凹函数(反之不成立) ;()f()f(4) 是(严格)凸函数 是(严格)拟凸函数(反之不成立) ;xx(5)单调函数既是拟凹函数也是拟凸函数Comment z2: 稳态值是指选择变量的最优解 还是指函数的*12(,.)nx最优解?产生疑问是f因为蒋中一那本书里提到的是稳定值的概念,用的是后一种表述。前一种表述是高微笔记上记的。3(6)凹(凸)函数相加仍为凹(凸)函数,拟凹和拟凸函数则没有类似关系。(二)无约束的最优化问题1一元函数的无约束极值本讲义将讨论的函数范围限定在二次连续可微函数的范围里。给定一个二次连续可微的一元函数, 。易知,它在 处取得极值的一阶()yfx0x必要条件为: 。而该极值究竟是极大值还是极小值得看 的符号:若()0fx()f,则 为唯一的绝对极大值;若 ,则 为唯一的绝对极小值。()f()f0利用上述极值的导数条件,我们可以推导出极值的微分条件,即:极值的一阶必要条件:对于任意非零 ,函数的一阶全微分为()0dyfxdx零;对于任意非零 ,我们也可以通过222()()ffxdfx计算函数的二阶全微分来判断极值的情况。综上,当函数为二次连续可微时,它取得极值的必要条件为:(1)函数在 取得绝对极大值 ,对于任意非零 都成立; *x*22()0yfdxdx(2)函数在 取得绝对极小值 ,对于任意非零 都成立。*22()f在满足必要条件的前提下,函数取得唯一的绝对极值时充分条件为,对于任意非零 都成立 函数在 取得唯一绝对极大值;2*2()0dyfxdx*x,对于任意非零 都成立 函数在 取得唯一绝对极小值。只要将 改为一阶微分向量 ,以上极值的微分条件能直接从单变量的情况推广至两个甚至多个变量的情况。2多元函数的最优化问题2.1一阶条件稳态值: 上的函数 的稳态值 ,在该点处,nR12(,.)nyfx*12(,.)xTnx下面几个等式同时成立: *21*,.0,().,nnfxComment z3: 最优解能不能这么表示?在 Reny的附录里有表达式*12(,.)(xnfxf4定理如果在点 ,我们可能得到局部最大(小)值,即对于*12(,.)xTnx一个尽可能小的邻域内,所有点 都有 *12(,.)Tn 12(,.)xTnx,那么稳态条件必然满足。*12(,.)nff2.2二阶条件直觉上,多元函数与一元函数一样,在稳态值取得最大值还是最小值与 的符号有关。2dy我们先对 进行微分,可得:dy2*12 *11212 2*1122(x)().(x)()().()xnnnnnnnfdfdffxdffxfd*12*12()()nTnnfdfff其中, 为海塞矩阵。根据杨格定理:11*2*12.x().x().()nnHfff,因此海塞矩阵为对称矩阵。ijif在判断 的符号之前,我们先正(负)定矩阵及其判定方法。2dy定义若对于所有的 , 始终成立,则称 正定,A 为正定矩阵;x0()=x0TqA(x)q若对于所有的 , 始终成立,则称 负定,A 为负定矩阵;若对于所有的 , 始终成立,则称 半正定,A 为半正定矩阵;()T()若对于所有的 , 始终成立,则称 半负定,A 为半负定矩阵。x0=x0qAxq根据以上定义,若要判断 的符号,我们只需判定与其对应的海塞矩阵的正(负)定。2dyComment z4: 高微笔记里是前主子阵。看以后感觉那里记的概念不准确,就用了金老师上课用的解释,但只有顺序主子式的介绍,不知道适不适合前主子阵。5其实,通过判定海塞矩阵的正(负)定,我们也可以判定函数的凹(凸)性,即对于二次连续可微函数 , 12(,.)nyfx(1)其海塞矩阵 负定 函数为严格凹函数 存在唯一绝对极大值;H(2)其海塞矩阵 正定 函数为严格凸函数 存在唯一绝对极小值。()接下来介绍正负定的判定方法。定义主子阵:对 矩阵 A,由 A 的 k个主对角线元素及其对应的非对角线元素来得到n的矩阵,称为 A的 k阶主子阵;由 A 的 前 k个主对角线元素及其对应的非对角线元素来得到的矩阵,为 k阶前主子阵。主子阵的行列式为主子式;前主子阵的行列式为顺序主子式。我们用 表示 的 k阶顺序主子式(其中 ) ,如:kD(x)H1,23.n,1()f,122()xf。121212().(x).(x).()nknnffDf定理对于二次连续可微函数, 12,.yx(1) 海塞矩阵正定;0(,.)k(2) 海塞矩阵负定。),Dn用 表示海塞矩阵 H的指标(1,2,3,n)的任意排序, 为 的 k阶顺序主 kDH子式,则(3) 海塞矩阵半正定;0,(1,2.)k(4) 海塞矩阵半负定。),从而,我们给出极值的充分条件:已知二次连续可微函数 , 12(,.)(xnyfxf*12(,.)0(,.)infxi6(1)其海塞矩阵 负定 严格凹函数 为函数的唯一绝对极大值;(x)H*(x)f(2)其海塞矩阵 正定 严格凸函数 为函数的唯一绝对极小值。3 举例:二元函数的无约束极值问题有一个二次连续可微函数 ,12(,)yfx可知其海塞矩阵为 ,则21Hf, ,1Df212f,12,f 1212,ff根据之前的判定规则, (1) , 为严格凹函数;10Df22110ff12(,)yfx(2) , 为严格凸函数;2(3) , 为凹函数;12,f211ff12(,)yfx(4) , 为凸函数;0D220若 ,我们就可以根据函数的凹凸性来判定函数在点*121(,)(,)fxfx取得的是绝对极大值还是绝对极小值。*2,7(三)具有约束条件的最优化问题之前的部分只是考虑了无约束条件的最优化问题,这即是说在求极值的过程中,我们没有对选择变量的值进行约束,从而求得的解可能是负值,也可能很大。然而考虑到经济学是建立在稀缺的资源如何配置的问题上的,因而在经济学的最优化求解过程中,我们通常不得不面临资源的稀缺性即对选择变量的值加上约束条件。约束条件大致分三类:等式约束、非负约束以及更普遍的,其它形式的不等式约束。我们将依次介绍对应的求解方法。从现在开始,讨论将以最大化问题为主,在解决最大化问题后会稍微提及解决最小化问题的方法。1 等式约束关于解决等式约束的方法,其实我们已经学过了,就是利用拉格朗日方法求解的过程。现在简要回顾拉格朗日函数。1.1 二元目标函数、一个等式约束的约束最优化条件考虑二元函数下,具有约束条件的最优化问题 12, 12max(,)yfx.stgc其中 c 是一个常数,z 和 g 都是二次连续可微函数。该问题的拉格朗日函数为: 1212(,)(,)Lfxcx一阶条件要求: *12121(,)(,)0fgxx*12122(,)(,)fL*12(,)0cgx求出上述一阶条件,可得 。*12(,)x二阶条件:将拉格朗日乘子也看作是变量,则最大化拉格朗日函数的过程可视为无约束最优化过程。这也就是说,如果解 满足 L 的无约束极值中极大值的二阶条件,我们便可*12(,)x确定 是我们约束最优化问题的解。事实上在二阶条件求导过程中,这里与无约*12(,)x束最优化关键区别在于, 与 的取值不再是任意非零即可,等式约束中 与 的1dx2 1dx2取值有关。对等式约束 两边求一阶全微分,可得:12(,)gc8, 11220gdgxdxd因此等式约束要求 。12x对函数 进行二阶全微分可得:12(,)yf2122 221212121() () dydxddxfffxfdxffdxxd对 进行二阶全微分,化简可得:12(,)gc2211222 1() 0( )gxgxdgxdd将上式代入 ,可得:2y2 21121()dxLgLg而 。2211210xyyxyLggL定义 为加边海塞矩阵, 它是由海塞矩阵和一阶导数(边)构成0xyyxyHgL的矩阵,用 表示,H 上面的- 表示边。 。则0gxyyxyDL221()dg综上,我们可以得到目标为二元函数、仅包含一个等式约束的最优化条件:当 满足拉格朗日函数的一阶条件时, 为约束极小(大)*12(,)x *12()0,)Dx值。91.1.1 严格拟凹(拟凸)函数与约束极值的关系当函数 y 是二次连续可微时,我们还可以用函数的一阶导数和二阶导数(整理成加边行列式)的方法来检验:设 1220fBff(1) z 为严格拟凹函数;0B(2) z 为严格拟凸函数。将 B 与之前的加边海塞矩阵 进行比较,可以发现两个不同之处:H一为 B 中的加边元素是函数 f 而非 g 的一阶偏导数,二为 B 中的其余元素是 f 而非拉格朗日函数 L 的二阶偏导数。然而,在线性等式约束 的特定情1212(,)xxc况下(这类等式约束在经济学中经常遇到) , 可简化为 ,即 。ijLijfijijLf从而,拉格朗日函数为 1212()()Lfxcx从而且 。0iiifijijLf回到“边” ,我们注意到线性约束函数产生一阶导数 ,因而一阶条件可写为iig。因此 B 中的边只不过是 的边被正的标量 乘。通过顺序提取 的横边和纵iifgHH边的公因子,得到 1212200gfLDff结果,在线性约束情况下, 与 总有相同的符号。BD由此可知,在线性约束的条件下,我们可以通过直接判断目标函数的严格拟凹(凸)性去判断约束极值的情况:(1)目标函数为严格拟凹函数 函数在稳态值取得唯一的约束绝对极大值;(2)目标函数为严格拟凸函数 函数在稳态值取得唯一的约束绝对极小值。1.1.2 拟凹(凸)函数与凹(凸)函数的关系平滑、递增、拟凹的效用函数 上等值集为凸集 凸的向下倾斜的无差异曲线。因为等产量曲线的概念几乎与无差异曲线是一致的,我们可以类推:平滑、递增、拟凹的生产函数 上等值集为凸集 凸的向下倾斜的等产量曲线。Comment z5: 自己的总结,逻辑对么?101.2多元目标函数、m 个等式约束的约束最优化条件现在将拉格朗日方法应用于多元函数。面临的最优化问题为: 12,. 12ax(,.)n nyfx122,(,.)nmmstgcx拉格朗日函数为: 1(,)+(jjLfcg一阶条件: *1(x)()0,2.)jmjii if nx*(),.jjLcgj二阶条件:此时加边海塞矩阵为: 1110nmnnngHL定义 。1110,1,.()kmk kkkgDnmLg利用 我们直接给出多元目标函数、 m个等式约束的约束最优化条件:为满足一阶必要条件的解,则*(x,)(1) 函数在点 取得唯一的约束绝对极大值;10(1,2.,)kDn*x11(2) 函数在点 取得唯一的约束绝对极小值。0(1,2.,)kDmn*(x)2 非负约束考虑一元可微函数: a()fx.0st由于约束条件 ,因此可能会出现三种情况:0x(1)在 x 大于零时 取得绝对极大值。此时我们得到了一个内点解。在这种情况下,一阶条件是 ,和经典问题一样。()df(2)x 等于零时 取得绝对极大值。此时我们得到了一个边界解,但仍然成立。()0f(3)x 小于零时 取得绝对极大值。此时我们也得到了边界解,但因为作为非线性约束问题中的一个局部极大值,候选点必须必可行域中的邻近点高,从而要求。()df综上,为了在 内找到 的极大值, 必须满足以下三个条件中的一个:0,)x(1) 且 ; ()f(2) 且 ;x(3) 且 ;()0f将上述三个条件合成一个论述:, 且 ,()0fx()0xf其中,第三个等式表达了三个条件的一个共同特点,即 x 和 至少有一个是零,因()f此两者的乘积一定是零。这个特点是指 x 与 互补松弛。()f当问题包含 n 个选择变量时: 12ma(,.)nfx.0jstx解决的思路与一元函数相同,这里我们直接给出该约束最优化的必要条件:(1)给定非负约束,多元函数 在稳态值处取得约束极大值,则12(,.)nf满足*12(,.)nfx12*(x)()0,12,.,.()0,12.iiifi nfixn(1)给定非负约束,多元函数 在稳态值处取得约束极小值,则12,.)nf满足*12(,.)nfx*(x)()0,12,.,.()0,12.iiifi nfixn3 其它形式的不等式约束现在我们在非负约束的基础上,再引入不等式约束。为简化,我们先处理两个选择变量和一个不等式约束条件的问题: max()fx.0stgr在虚拟变量 s 的帮助下,我们可以将上述问题变换为: ax()fx.,0stgsr若没有非负约束,则我们可利用拉格朗日函数求解: ()()Lfxrxs一阶条件为: 0Ls但由于 x 与 s 必须是非负的,因此 1.1 部分的思路,上述一阶条件应改为:且 ;0,Lx且 ;,s0Ls。13注意 仍然成立,为何?因为 一定成立 。0L()0Lrgxs进而,将 以及 代入第二个条件,则第二个条件与第三个条件可变为:s()srgx且12(,)0,()0rgx从而我们可以用没有虚拟变量的等价形式来表达的一阶条件(这时 )()()Lfxrgx:且 ;0,xLfg0x且 。(),rL上述讨论可以直接的方式应用于 n 个选择变量和个约束条件的问题。拉格朗日函数 L 的形式为: 12 121(,.)(,.)mini nLfxrgx则该非线性约束问题的库恩-塔克条件为(极大化):且 ;0,jjjjfgxx0(,.)jLjx且 。12(,.),iijniLrii如果问题是求极小值,那么可以将它转化为极大化问题,然后应用以上条件求解。此处仅介绍用库恩-塔克条件求解非线性规划中的极大值问题。按水平的方向解读上式,我们可以看见库恩-塔克条件在极大化问题中包括了一组与选择变量与拉格朗日乘数的条件。从垂直的方向来解读,对于每一个选择变量和拉格朗日函数,都有一个边际条件(第一列) 、一个非负约束(第二列)和一个互补松弛条件(第三列)。在任一个给定条件下,与选择变量相关的一组边际条件与拉格朗日成熟的一组边际条件在不等号方向上是不同的。若满足约束规范,库恩-塔克条件极大(小)化条件可以作为总体极大(小)值的必要条件。1.2 举例我们现在就将库恩-塔克条件应用于效用最大化问题: maxUy.104,stxy拉格朗日函数为:1412(0)(40)Lxyxyx该问题的约束条件为线性的,因而它一定满足约束规范(之后解释) ,故可利用库恩-塔克条件求解,可得:且 ;120,yxx0L且 ;1,Ly且 ;110,0x10L且 。224,L2写出库恩-塔克条件后,典型的方法是通过试错法来求解。步骤如下:(1)首先给选择变量赋值为零。通过消除某些项来使条件简化。如果适当的非拉格朗日乘数可以满足所有边际不等式,那么零解将是最优的。对于当前这个例子,当 或 时, 没有意义,因此该问题中 x 和 y 都是0xy0U正数。此步骤跳过。(2)如果零解违反一些不等式,那么可以尝试让一个或更多选择变量为正数。对于每个正的选择变量,我们可以通过互补松弛条件使不等式边际条件转换为严格等式边际条件。应用于当前例子,则:;1212210=0,0yxLyyxx 由于 ,其它条件可变为:1x=0;10Lxy这时我们仍然不能简化条件 且 。因此我们需要步骤(3)。224,20L(3)假设函数对拉格朗日乘数的偏导取不等号;若这个假定导致矛盾,那我们应该将该偏导等于零进行测算。15因此,我们先假设 则 ,那么有 。 但这一解违反240,Lx20yx50xy了约束 ,故舍去。40x那么,我们假设 则 。2,x4,6y1.3 约束规范之前多次强调,库恩-塔克条件只有在满足规范约束时,才是极值的必要条件。那么约束规范具体是什么?先介绍几个概念。令 是可行区域边界上的一个可能的解点,并令*12(,.)nxx表示由所提到的边界点移动的特定方向。dd测试向量:若某一向量满足条件(i)如果第 j 个选择变量在点 x*处取得零值,那么只允许在 xj 轴上有非负变化,即:若,那么 ;*0jx0jx(ii)如果在点 x*处恰好满足第 i 个约束条件的等式约束,那么将只允许 的12,.,ndx取值使得约束函数值 不增加(对极大值问题) ,即:()ig若 ,则 必须成立;(*)iigxr12*.0iiiindxgdx则该向量为测试向量。规范弧:若某一弧段满足条件(i)从点 x*处出发;(ii)整个包含在可行区域内;(iii )与已知测试向量相切;则该弧段我们成为该测试向量的规范弧。有了这些预备知识后,约束规范可简单地表达为:如果对可行区域边界上的任意点 x*,对每一测试向量 dx,存在一规范弧,那么,就满足约束规范。约束规范的定义看上去有些抽象,有兴趣的同学可以借助参考书中的例子求解理解掌握这一概念。如果实际求解中,我们遇到可行区域是仅由线性约束形成的凸集,那么约束规范总是满足,且库恩-塔克条件在最优解处总成立。这时就免去了检验约束规范是否满足的步骤。16(四)最优化的其他主题接下来我们将回到经典的等式约束最优化领域来讨论包络定理。1 极大值函数极大值函数是当选择变量都是最优值时候的目标函数。这些选择变量的最优值是外生变量和参数的函数。一旦选择变量的最优值代入原目标函数中,那么目标函数就间接地称为参数的函数。因此,极大值函数也称间接目标函数。它是当参数发生变化的时候,目标函数极大(小)值变化的轨迹。举例:通过求解 max(,)Uuxy.,0stpB我们可得 。将最优值代入 U=u(x,y),则可得间接*(,),*()xyxypB效用函数: 。V2 包络定理包络定理:即使在外生变量可能作为内生选择变量的解的一部分间接进入极大值函数的情况下,也只有外生变量参数变化的直接效应才需要考虑。为了阐释这一概念,考
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 物理基础考试试题及答案
- java路由配置方法面试题及答案
- 赌博知识考试题及答案
- 2025年河南省中考语文真题(含答案)
- 石油测井考试题及答案
- 长春燃气考试题及答案
- 函数极限试题及答案
- 货币知识测试题及答案
- 手卫生院感考试试题及答案
- 校园业务知识培训内容课件
- 部编小学语文单元作业设计五年级上册第二单元
- 企业社会责任报告模板
- 25题后期-剪辑-特效岗位常见面试问题含HR问题考察点及参考回答
- 2024年辽宁交投集团招聘笔试参考题库附带答案详解
- 高考英语必背1500个真题高频词汇- 高考英语一轮复习
- 人体足解剖学
- 机械基础 第三版 课件 (郁志纯)模块三 机械零件的精度
- 环境监测仪器设备采购投标方案(技术标)
- 【食品零食】桂格燕麦食品抖音账号运营方案
- 食材供应服务投标方案(完整技术标)
- 诺丁山-Notting-Hill-中英文剧本
评论
0/150
提交评论