数理经济学教学大纲 内容架构 资料来源 : 集合与映射 (JR,pp ) 1. 数学 : 概念 定理 证明, 重在直觉, 多用图形 ; 2. 举例 : 数学例子 ; 3. 应用 : 经济学实例 ;

Similar documents
Fundamentals of Heat and Mass Transfer, 6 th edition

Integrated Algebra. Simplified Chinese. Problem Solving

Chapter 2 THE MATHEMATICS OF OPTIMIZATION. Copyright 2005 by South-Western, a division of Thomson Learning. All rights reserved.

d) There is a Web page that includes links to both Web page A and Web page B.

Lecture English Term Definition Chinese: Public

DYNAMIC LECTURE 5: DISCRETE TIME INTERTEMPORAL OPTIMIZATION

2012 AP Calculus BC 模拟试卷

Chapter 2 Bayesian Decision Theory. Pattern Recognition Soochow, Fall Semester 1

2. The lattice Boltzmann for porous flow and transport

Source mechanism solution

QTM - QUALITY TOOLS' MANUAL.

= lim(x + 1) lim x 1 x 1 (x 2 + 1) 2 (for the latter let y = x2 + 1) lim

上海激光电子伽玛源 (SLEGS) 样机的实验介绍

Chapter 2 the z-transform. 2.1 definition 2.2 properties of ROC 2.3 the inverse z-transform 2.4 z-transform properties

THE INVERSE DERIVATIVE

Rigorous back analysis of shear strength parameters of landslide slip

Riemann s Hypothesis and Conjecture of Birch and Swinnerton-Dyer are False

Lecture Note on Linear Algebra 14. Linear Independence, Bases and Coordinates

Principia and Design of Heat Exchanger Device 热交换器原理与设计

The general programming problem is the nonlinear programming problem where a given function is maximized subject to a set of inequality constraints.

Galileo Galilei ( ) Title page of Galileo's Dialogue concerning the two chief world systems, published in Florence in February 1632.

ECON 582: Dynamic Programming (Chapter 6, Acemoglu) Instructor: Dmytro Hryshko

Numerical Analysis in Geotechnical Engineering

Design, Development and Application of Northeast Asia Resources and Environment Scientific Expedition Data Platform

通量数据质量控制的理论与方法 理加联合科技有限公司

EC /11. Math for Microeconomics September Course, Part II Problem Set 1 with Solutions. a11 a 12. x 2

The Lagrange Mean Value Theorem Of Functions of n Variables

Lecture Note on Linear Algebra 16. Eigenvalues and Eigenvectors

An IntRoduction to grey methods by using R

Linear Algebra 18 Orthogonality

Concurrent Engineering Pdf Ebook Download >>> DOWNLOAD

The dynamic N1-methyladenosine methylome in eukaryotic messenger RNA 报告人 : 沈胤

Tsinghua-Berkeley Shenzhen Institute (TBSI) PhD Program Design

Mathematical Economics. Lecture Notes (in extracts)

EC /11. Math for Microeconomics September Course, Part II Lecture Notes. Course Outline

New Newtonian Mechanics and New Laws of Motion

Conditional expectation and prediction

Math Camp Notes: Everything Else

Desertification Monitoring in China. Wang Junhou China National Desertification Monitoring Center July, 2012

On the Quark model based on virtual spacetime and the origin of fractional charge

Brainwashed Tom Burrell Pdf Download >>> DOWNLOAD

Bending Deflection.

A Tutorial on Variational Bayes

Week 4: Calculus and Optimization (Jehle and Reny, Chapter A2)

GRE 精确 完整 数学预测机经 发布适用 2015 年 10 月考试

Torsion.

Atomic & Molecular Clusters / 原子分子团簇 /

There are only 92 stable elements in nature

Dynamic Optimization Problem. April 2, Graduate School of Economics, University of Tokyo. Math Camp Day 4. Daiki Kishishita.

Modeling effects of changes in diffuse radiation on light use efficiency in forest ecosystem. Wei Nan

能源化学工程专业培养方案. Undergraduate Program for Specialty in Energy Chemical Engineering 专业负责人 : 何平分管院长 : 廖其龙院学术委员会主任 : 李玉香

A new operational medium-range numerical weather forecast system of CHINA. NWPD/NMC/CMA (Beijing,CHINA)

Adrien-Marie Legendre

Lecture 2: Introduction to Probability

USTC SNST 2014 Autumn Semester Lecture Series

沙强 / 讲师 随时欢迎对有机化学感兴趣的同学与我交流! 理学院化学系 从事专业 有机化学. 办公室 逸夫楼 6072 实验室 逸夫楼 6081 毕业院校 南京理工大学 电子邮箱 研 究 方 向 催化不对称合成 杂环骨架构建 卡宾化学 生物活性分子设计

Microeconomics I. September, c Leopold Sögner

CHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS. W. Erwin Diewert January 31, 2008.

第五届控制科学与工程前沿论坛 高志强. Center for Advanced Control Technologies

Type and Propositions

碳酸盐岩缝洞型油藏开发关键技术. Key Petroleum Development Technologies in Fractured & Caverned Carbonate Reservoirs Illustrated by the Example of Tahe Oilfield

Cooling rate of water

Molecular weights and Sizes

国际数值预报现状和发展 沈学顺 中国气象局数值预报中心

大原利明 算法点竄指南 点竄術 算 額 絵馬堂

个人简历 冯新龙 : 男, 中共党员, 理学博士, 教授, 博士生导师通信地址 : 新疆大学数学与系统科学学院邮编 :830046

三类调度问题的复合派遣算法及其在医疗运营管理中的应用

tan θ(t) = 5 [3 points] And, we are given that d [1 points] Therefore, the velocity of the plane is dx [4 points] (km/min.) [2 points] (The other way)

Final Exam - Math Camp August 27, 2014

ILC Group Annual Report 2018

Synthesis of anisole by vapor phase methylation of phenol with methanol over catalysts supported on activated alumina

Effect of lengthening alkyl spacer on hydroformylation performance of tethered phosphine modified Rh/SiO2 catalyst

( 选出不同类别的单词 ) ( 照样子完成填空 ) e.g. one three

A new approach to inducing Ti 3+ in anatase TiO2 for efficient photocatalytic hydrogen production

An Introduction to Generalized Method of Moments. Chen,Rong aronge.net

Numerical Heat Transfer ( 数值传热学 ) Chapter 4 Numerical Solution of Diffusion Equation and its Applications(2)

Lecture 13 Metabolic Diversity 微生物代谢的多样性

Optimization. A first course on mathematics for economists

Service Bulletin-04 真空电容的外形尺寸

Workshop on Numerical Partial Differential Equations and Scientific Computing

Mathematical Economics: Lecture 16

Explainable Recommendation: Theory and Applications

Stress State.

Study on property and deterioration of whole boned bolt in heat-harm tunnel

强子谱和强子结构研究前沿简介 邹冰松 中国科学院高能物理研究所 中国科学院大科学装置理论物理研究中心

Bending Internal Forces.

minimize x subject to (x 2)(x 4) u,

Global and Regional Precipitation Measurement and Applications

Geomechanical Issues of CO2 Storage in Deep Saline Aquifers 二氧化碳咸水层封存的力学问题

铁硒超导薄膜的原位制 备和电子结构测量 谢斌平 复旦大学 杭州 2014 全国表面分析应用技术学术交流会及赛默飞世尔科技 XPS 应用及设备维护技术培训班

La pietra da altre colline può levigare la giada di questa qui Il Classico dei Versi 可 以 攻 玉

Influence of surface strain on activity and selectivity of Pd based catalysts for the hydrogenation of acetylene: A DFT study

Algorithms and Complexity

系统生物学. (Systems Biology) 马彬广

Key Topic. Body Composition Analysis (BCA) on lab animals with NMR 采用核磁共振分析实验鼠的体内组分. TD-NMR and Body Composition Analysis for Lab Animals

Week 3: Sets and Mappings/Calculus and Optimization (Jehle September and Reny, 20, 2015 Chapter A1/A2)

Photo induced self formation of dual cocatalysts on semiconductor surface

复合功能 VZθ 执行器篇 Multi-Functional VZθActuator

Transcription:

数理经济学教学大纲 内容架构 1. 数学 : 概念 定理 证明, 重在直觉, 多用图形 ; 2. 举例 : 数学例子 ; 3. 应用 : 经济学实例 ; 资料来源 : 杰里和瑞尼 : 高级微观经济理论, 上海财大出版社和培生集团,2002 年 (JR) 马斯 - 科莱尔 温斯顿和格林 : 微观经济学, 中国社会科学出版社,2001 年 (MWG) 巴罗和萨拉伊马丁 : 经济增长, 中国社会科学出版社,2000 年 (BS) 迪克西特 : 经济理论中的最优化方法, 上海三联和上海人民出版社,2006 年 (KD) 德雷泽 : 宏观经济学中的政治经济学, 经济科学出版社,2003 年 (AD) 斯托基和卢卡斯 : 经济动态的递归方法, 中国社会科学出版社,1999 年 (SL) 集合与映射 (JR,pp.387-435) 凸集开球 闭球开集 闭集分离超平面 (MWG p.1346) 有界集紧集 连续性 ( 柯西 ) 连续性 D 中的开集 D 中的闭集连续性及其逆象 数列收敛的数列有界的数列子数列 威尔斯拉斯定理布劳威不动点定理

实值函数递增的 严格递增与强递增函数递减的 严格递减与强递减函数水平集上优集与下劣集 凹函数严格凹函数拟凹函数严格拟凹函数凸与拟凸函数凸与严格凸函数 微积分 (JR,pp.436-483) 链式运算法则凹性与一和二阶导数方向导数梯度海赛矩阵杨格定理 齐次函数齐次函数的偏导数欧拉定理 隐函数定理 (BS p.490) 泰勒定理 (BS p.491) 静态最优化 (JR, pp.436-483) 实值函数局部内点最优化的一阶必要条件实值函数局部内点最优化的二阶必要条件海赛矩阵负定与正定的充分条件实值函数局部内点最优化的充分条件 ( 无约束的 ) 局部与全局最优化严格凹性 / 凸性与全局最优化的惟一性惟一全局最优化的充分条件 约束最优化拉格朗日方法

拉格朗日定理 加边海赛矩阵两变量 一约束最优化问题中的一个局部最优化的充分条件等式约束条件下的局部最优化的充分条件不等式约束受非负性条件约束的实值函数最优化的必要条件 非线性规则问题库恩 - 塔克条件受不等式条件约束的实值最优化的 ( 库恩 - 塔克 ) 必要条件 值函数包络定理 微分方程 (BS,pp.438-466) 微分方程一阶线性微分方程二阶线性微分方程微分方程的稳定性一阶线性微分方程的稳定性线性微分方程系统的稳定性非线性微分方程系统的稳定性 动态优化 最优控制原理 (BS, pp.466-487; KD, pp.130-144) 典型问题极大值原理横截性条件现值和当期值汉密尔顿函数多个变量 动态规划导论 (KD,pp.145-159;AD,pp.28-33) 贝尔曼方程

CHAPTER 1 SETS AND MAPPINGS 1. Convex Set S n is a convex set if for all x 2 S, we have tx 1 1 tx 2 S for all t 0, 1. Convex set: the Intersection of Convex Sets is Convex 2.Open and Closed Balls Open Ball with center x 0 and radius 0 is the subset of points in n, B x 0 x n dx 0, x Closed Ball with center x 0 and radius 0 is the subset of points in n, B x 0 x n dx 0, x Open set: S n is an open set if, for all x S, there exists some 0 such that B x S. Close set: S is closed set if its complements S C is an open set. Open sets: The empty set,, is an open set. The entire space, n, is an open set. The union of open sets is an open set. The intersection of any finite number of open sets is an open set. Every Open Set is a Colletion of Open Balls: Let S be an open set. For every x S, choose some x 0 such that B x x S. Then, S x S B x x. Closed Sets: The empty set,, is a closed set. The entire space, n, is a closed set. The union of any finite collection of closed sets is a closed set. The intersection of closed sets is a close set. 3.Separating Hyperplanes Given p n with p 0, and c, the hyperplane generated by p

and c is set H p,c z n p z c. The sets z n p z c and z n p z c are called, respectively, the half-space above and the half-space below the hyperplane H p,c. Separating Hyperplane Theorem: Suppose that B n is convex and closed, and that x B. Then there is p n with p 0, and a value c such that p x c and p y c for every y B. More generally, suppose that the convex sets A, B n are disjoint (i.e., A B ). Then there is p n with p 0, and a value c, such that p x c for every x A and p y c for every y B. That is, there is a hyperplane that separates A and B, leaving A and B on different sides of it. Supporting Hyperplane Theorem: Suppose that B n is convex and that x is not an element of the interior of set B (i.e., x Int B). Then there is p n with p 0 such that p x p y for every y B. 4.Bounded sets Bounded sets: S is bounded if there exists some 0 such that S B x for some x n. Lower bound: any real number l for which l x for all x S is called a lower bound for S. Upper bound: any real number u for which x u for all x S is called an upper bound for S. S is bounded from below if it has a lover bound, and is bounded from above if it has an upper bound. The largest number among the lower bounds is called the greatest lower bound (g. l. b. ) of S. The smallest number among the upper bounds is called the least upper bound (l. u. b. ) of S. Upper and Lower Bounds in Subsets of Real Numbers: 1. Let S be a bounded open set in and let a be the g. l. b. of S and b be the l. u. b. of S. Then a S and b S. 2. Let S be a bounded closed set in and let a be the g. l. b. of S and S and b be the l. u. b. of S. Then a S and b S. (Heine-Borel) Compact sets: a set S in n and bounded. is called compact if it is closed

5. (Cauchy) Continuity Let D be a subset of n, and let f : D n. The function f is continuous at the point x 0 D if for every 0, there is a 0 such that fb x 0 D B fx 0. If f is continuous at every point x D, then it is called a continuous function. Open sets in D : Let D be a subset of n. Then a subset S of D is open in D if for every x S there is an 0 such that B x D S. Closed sets in D : Let D be a subset of m. A subset S of D is closed in D if its complement in D is open in D. Continuity and Inverse Images: Let D be a subset of m. The following conditions are equivalent: 1. f : D n is continuous. 2. For every open ball B in n, f 1 B is open in D. 3. For every open set S in n, f 1 S is open in D. Theorem The Continuous Image of a Compact Set is a Compact Set Let D be a subset of m and let f : D n be continuous function. If S D is compact in D, then its image fs n is compact in n. 6. Sequence Sequences in n : A sequence of n is function mapping some infinite subset I of positive integers into n. We shall denote a sequence by x k } k I, where x k n for every k I. Convergent sequences: The sequence x k k I converges to x n if for every 0, there is a k such that x k B x for all k I exceeding k. Bounded sequences: A sequence x k } k I in n is bounded if for some M, x k M for all k I. Subsequences: x k } k J is a subsequence of the sequence x k } k I in n, if J is an infinite subset of I. Bounded sequences: Every bounded sequence in n has a convergent subsequence.

Sequence, Sets, and Continuous Functions: Let D be subset of n, and let f : D m. Then 1. D is open if and only if for each x D, if x k k1 converges to x, then for some k, x k D for all k k. 2. D is closed if and only if for every x k k1 of points in D converging to some x n, it is also the case that x D. 3. f is continuous if and only if whenever x k k1 in D converges to x D, fx k k1 converges to fx. 7. Some Existing Theorems (Weierstrass) Existence of Extreme Values: Let f : S be a continuous real-valued mapping where S is a nonempty compact subset of n. Then these exists a vector x S such that f x fx fx for all x S. The Brouwer Fixed-Point Theorem: Let S n be a nonempty compact and convex set. Let f : S S be continuous mapping. Then there exists at least one fixed point of f in S. That is, thers exists at least one x S such that x fx. 8. Real-Valued Functions Real-valued function: f : D R is a real-valued function fi D is any set and R. Increasing, Strictly Increasing and Strongly Increasing Functions: Let f : D, where D is subset of n. f is increasing if fx 0 fx 1 whenever x 0 x 1. f is strictly increasing if fx 0 fx 1 whenever x 0 x 1. f is strongly increasing if fx 0 fx 1 whenever x 0 x 1 and x 0 x 1. Decreasing, Strictly Decreasing and Strongly Decreasing Functions: Let f : D, where D is subset of n. f is decreasing if fx 0 fx 1 whenever x 0 x 1. f is strictly decreasing if fx 0 fx 1 whenever x 0 x 1. f is strongly decreasing if fx 0 fx 1 whenever x 0 x 1 and x 0 x 1. 9.Related sets Level sets: Ly 0 is a level set of the real-valued function f : D R iff Ly 0 x x D, fx y 0, where y 0 R.

Level sets relative to a point: Lx 0 is a level set relative to x 0 Lx 0 x x D, fx fx 0. if Superior set: Sy 0 x x D, fx y 0 is called the superior set for level y 0. Inferior set: Iy 0 x x D, fx y 0 is called the inferior set for level y 0. Strictly superior set: S y 0 x x D, fx y 0 is called the strictly superior set for level y 0. Strictly inferior set: I y 0 x x D, fx y 0 is called the strictly inferior set for level y 0. Superior, Inferior, and Level Sets: For any f : D R and y 0 R : 1. Ly 0 Sy 0. 2. Ly 0 Iy 0. 3. Ly 0 Sy 0 Iy 0. 4. S y 0 Sy 0. 5. I y 0 Ly 0. 6. S y 0 Ly 0. 7. I y 0 Ly 0. 8. S y 0 I y 0. 10. Concave Functions Assumption: (1) D n is a convex set. (2) When x 1 D and x 2 D, x t tx 1 1 tx 2, for t 0, 1, denote the convex combination of x 1 and x 2. Concave function: f : D R is a concave function if for all x 1, x 2 D, fx t tfx 1 1 tfx 2 t 0, 1 Theorem Points on and below the graph of a concave function is always form a convex set Let A x, y x D, fx y be the set of points "on and below" the graph of f : D R, where D, Then, f is a concave function A is a convex set. Strictly concave function: f : D R is a strictly concave function iff, for all x 1 x 2 in D, fx t tfx 1 1 tfx 2 for all t 0, 1

Quasiconcave function: f : D R is quasiconcave iff, for all x 1 and x 2 in D, fx t minfx 1, fx 2 for all t 0, 1 Quasiconcavity and the superior sets f : D R is a quasiconcave function iff Sy is a convex set for all y. Strictly quasiconcave function: A function f : D R is strictly quasiconcave iff, for all x 1 x 2 in D, fx t minfx 1, fx 2 for all t 0, 1. Concavity implies quasiconcavity: A concave function is always quasiconcave. A strictly concave functions is always strictly quasiconcave. 11. Convex and quasiconvex functions Convex function: f : D R is a convex function iff, for all x 1, x 2 in D, fx t tfx 1 1 tfx 2 for all t 0, 1. Strictly convex function: f : D R is a strictly convex function iff, for all x 1 x 2 in D, fx t tfx 1 1 tfx 2 for all t 0, 1. Concave and convex functions: fx is a (strictly) concave function fx is a (strictly) convex function. Points on and above the graph of a convex function always form a convex set: Let A x, y x D, fx y be the set of points "on and above" the graph of f : D R, where D n is a convex set and R. Then f is a convex function A is a convex set. Quasiconvex function: f : D R is quasiconvex x 1, x 2 D, fx t maxfx 1, fx 2 Strictly quasiconvex function: f : D R is strictly quasiconvex x 1, x 2 D, fx t maxfx 1, fx 2

Quasiconvexity and the inferior sets: f : D R is quasiconvex function Iy is a convex set for all y. Quasiconcave and quasiconvex functions: fx is a (strictly) quasiconcave function function. fx is a (strictly) quasiconvex Remark f is concave the set of points beneath the graph is convex f is convex the set of points above the graph is convex f quasiconcave superior sets are convex sets f quasiconvex inferior sets are convex sets f concave f quasiconcave f convex f quasiconvex f (strictly) concave f (strictly) convex f (strictly) quasiconcave f (strictly) quasiconvex

CHAPTER 2 CALCULUS Functions of a Single Variable Chain Rule: d dx fgx f gxg x Concavity and first and second derivatives: Let D be a nondegenerate interval of real numbers on which f is twice continuously differentiable. The following statements 1 to 3 are equivalent: 1. f is concave. 2. f x 0 x D. 3. x 0 D : fx fx 0 f x 0 x x 0 x D. Moreover, 4. f x 0 x D f is strictly concave. Functions of Several Variables Directional derivative: x, z, gt fx tz g 0 n i1 f i xz i The term on the right-hand side is known as the directional derivative of f at x in the direction z. Gradient: fx f 1 x,, f n x is called the gradient of f at x. Remark g 0 fxz. Second-order partial derivative: x i fx x 1 2 fx x i x 1 f 1i x Second-order gradient vector: the gradient of the partial with respect to x 1, f 1 x. f 1 x 2 fx x 1 x 1 Hessian matrix:,, 2 fx x n x 1 f 11 x,, f 1n x f 11 x f 12 x f 1n x Hx f 21 x f 22 x f 2n x f n1 x f n2 x f nn x

Young s Theorem: For any twice continuously differentiable function fx, 2 fx x i x j 2 fx x j x i i and j. Single-Variable and Multivariable Concavity: f is a real-value function defined on the convex subset D of n. f is (strictly) concave x D, nonzero z n, gt fx tz is (strictly) concave on t x tz D Negative semidefinite matrix: A is negative semidefinite z n, z T Az 0. A is negative definite z n, z T Az 0. A is positive semidefinite z n, z T Az 0. A is positive definite z n, z T Az 0. Slope,curvature, and concavity in many variables Let D be a convex subset of n with a nonempty interior on which f is twice continuously differentiable. The following statements 1 to 3 are equivalent: 1. f is concave. 2. Hx is negative semidefinite for all x in D. 3. For x 0 D : fx fx 0 fx 0 x x 0 x D. Moreover, 4. If Hx is negative definite for all x in D, then f is strictly concave. Concavity, convexity, and second-order own partial derivatives Let f : D R be a twice diffentiable function. 1. f is concave x, f ii x 0, i 1,, n. 2. f is convex x, f ii x 0, i 1,, n. Homogeneous Functions Homogeneous function: A real-valued function fx is homogeneous of degree k ftx t k fx, t 0. fx is homogeneous of degree 1, or linear homogeneous ftx tfx, t 0. fx is homogeneous of degree zero ftx fx, t 0.

Partial derivatives of homogeneous functions If fx is homogeneous of degree k, its partial deriatives are homogeneous of degree k 1. Euler s theorem (or adding-up theorem) n fx is homogeneous of degree k x, kfx i1 fx x i x i Some Useful Results in Calculus (BS) Implicit function theorem Taylor s theorem

CHAPTER 3 STATIC OPTIMIZATION Unconstrained Optimization Necessary conditions for local interior optima in the single-variable case Let fx be a twice continuously differentiable function of one variable. Then fx reaches a local interior 1. maximum at x f x 0 (FONC), f x 0 (SONC). 2. minimum at x f x 0 (FONC) f x 0 (SONC). FONC for local interior optima of real-valued functions: fx reaches a local interior maximum or minimum at x fx 0 SONC for local interior optima of real-valued functions: Let fx be a twice continuously differentiable. fx reaches a local interior maximum at x Hx is negative semidefinite fx reaches a local interior maximum at x H x is positive semidefinite Sufficient condition for negative and positive definiteness of the Hessian: Let fx be twice continuously differentiable, and let D i x be the ith-order principal minor of the Hessian matrix Hx. 1. 1 i D i x 0, i 1,, n Hx is negative definite. 2. D i x 0, i 1,, n Hx is positive definite. If condition 1 holds for all x in the domain, then f is strictly concave. If condition 2 holds for all x in the domain, then f is strictly convex. Sufficient conditions for local interior optima of real-valued functions Let fx be twice continuously differentiable. 1. f i x 0 and 1 n D i x 0, i 1,, n fx reaches a local maximum at x. 2. f i x 0 and D i x minimum at x. 0, i 1,, n fx reaches a local (Unconstrained) local-global theorem: Let f be a twice continuously differentiable real-valued concave function on D. The following statements are equivalent, where x is an interior point of D : 1. fx 0.

2. f achieves a local maximum at x. 3. f achieves a global maximum at x. Stict concavity/convexity and the uniqueness of global optima: 1. x maximizes the strictly concave function f x is the unique global maximizer, i.e., fx fx x D, x x. 2. x minimizes the strictly convex function f x is the unique global minimizer, i.e., f x fx x D, x x. Sufficient condition for unique global optima: Let fx be twice continuously differentiable. 1. fx is strictly concave and f i x 0, i 1,, n x is the unique global maximizer of fx. 2. fx is strictly convex and f i x 0, i 1,, n x is the unique global maximizer of fx. Constrained Optimization Equality constraints Two-variable, one constraint optimization problem: max fx x 1,x 1, x 2 s. t. gx 1, x 2 0 2 Lagrange s method 1.Lagarangian function: Lx 1, x 2, fx 1, x 2 gx 1, x 2 2.F.O.C L fx 1, x 2 x 1 x gx 1, x 2 1 x 1 0 L fx 1, x 2 x 2 x gx 1, x 2 2 x 2 0 L gx 1, x 2 0. The first-order partials fo the Lagrangian function with respect to the x i were L i f i g i. The second-order partials of L would then be

L 11 f 11 g 11 L 12 f 12 g 12 L 22 f 22 g 22. Bordered Hessian of the Lagrangian function H 0 g 1 g 2 g 1 L 11 L 12 g 2 L 21 L 22 Sufficient condition for a local optimum: D H 0 0 x 1, x 2 is a local maximum (minimum) n-variable, m-constraint optimization problem: Lagrange s theorem: Let fx and g j x, j 1,, m, be continuously differentiable real-valued functions over some domain D n. Let x be an interior point of D and suppose that x is an optimum of f subject to the constraints, g j x 0. If g j x are linearly independent, then there exist m unique numbers j, such that Lx, x i fx x i m j1 j g j x x i 0 Bordered Hessian matrix H 0 0 g 1 1 g n 1 0 0 g 1 m g n m g 1 1 g 1 m L 11 L 1n g n1 g nm L n1 L nn Sufficient conditions for local optima with equality constraints: Let the objective function be fx and the m n constraints be g j x 0, j 1,, m. 1. x is a local maximum of fx s.t. the constraints if the n m principal minors in D H alternate in sign beginning with positive

D m1 0, D m2 0,, when evaluated at x,. 2. x is a local minimum of fx s.t. the constraints if the n m principal minors in D H are all negative D m1 0, D m2 0,, when evaluated at x,. Necessary conditions for optima of real-valued functions subject to nonnegativity constraints Let fx be continuously differentiable. 1. If x maximizes fx subject to x 0, then x satisfies fx i. 0, i 1,, n x i ii. x i fx 0, i 1,, n x i iii. x i 0, i 1,, n. 2. If x minimizes fx subject to x 0, then x satisfies fx i. 0, i 1,, n x i ii. x i fx 0, i 1,, n x i iii. x i 0, i 1,, n. Inequality constraints A two-variable nonlinear programming programming max fx x 1,x 1, x 2 2 s. t. gx 1, x 2 0. Equivalent to the three-variable problem with equality and nonnegativity constraints: max x 1,x 2,z 1, x 2 s. t. gx 1, x 2 z 0 z 0. Lagrangian function: Lx 1, x 2, z, fx 1, x 2 gx 1, x 2 z Kuhn-Tucker conditions: L 1 f 1 g 1 0 L 2 f 2 g 2 0 gx 1, x 2 0 0, gx 1, x 2 0

(Kuhn-Tucker) Necessary conditions for optima of real-valued functions subject to inequality constraints Let fx and g j x, j 1,, m, be continuously diffentiable real-valued functions over some domain D n. Let x be an interior point of D and suppose that x is an optimum of f subject to the constraints, g j x 0, j 1,, m. If g j x associated with all binding constraints are linearly independent, then there exists a unique vector such that x, satisfies the Kuhn-Tucker conditions: Lx, x i fx x i m j1 j g j x x i 0, i 1,, n j g j x 0 g j x 0 j 1,, m. Furthermore, the vector is nonnegative if x is a maximum, and nonpostitive if it is a minimum. Value Functions Envelop Theorem: Consider the optimization problem (P1): max fx, a x s.t. gx, a 0 x 0. suppose the objective function and constraint are continuously differentiable in a, let xa 0 uniquely solve P1 and assume that it is also continuously differentiable in the parameters a. Let Lx, a, be the problem s associated Lagrangian function and let xa, a solve the Kuhn-Tucker conditions. Finally, let Ma be the problem s associated maximum-value function. Then, Ma L j 1,, m. a j aj xa,a where the right-hand side denotes the partial derivative of the Lagrangian function with respect to the parameter a j evaluated at the point xa, a.

CHAPTER 4 Differential Equations INTRODUCTION A differential equation is an equation that involves derivatives of variables. If there is only one independent variablce, then it is called an ordinary differential equation (ODE). If the highest derivative is an ODE is of order n, then it is an nth-order ODE. When the functional form of the equation is linear, then it is a linear ODE. Example: A first-order linear ODE a 1 ẏt a 2 yt xt 0 where ẏt dyt dt xt : forcing function If xt a 3, then the equation is called autonomous. If xt 0, then the equation is called homogeneous. Solution methods 1.Graphical Used for nonlinear, as well as linear, differential equations; Used only for autonomous equations. 2.Analytical Used only with a limited set of functional forms. 3.Numerical analysis e.g. Matlab has the subroutines ODE23 and ODE45, and Mathematica has the command NDSOLVE. First-order ODE Graphical solutions 1. Consider an automoumous ODE of the form, ẏt fyt Example 1 - Case 1, a 0 ẏt fyt a yt x

Case 2, a 0 - Example 2 ẏt fyt s yt yt s,, 0; 1 Definition: Stability ẏ - If y y 0, then y is locally unstable, ẏ - If 0, then y is locally stable. y y Analytical Solutions The solution to ẏt a is obviously yt b at, where b is an arbitrary constant.

Equations that involve polynomial functions of time ẏt a 0 a 1 t a 2 t 2...a n t n - has the solution yt b a 0 t a 1 t2 2...a n t n1 n 1 The general solution for linear, first-order ODEs - Linear, first-order differential equations with constant coefficients. ẏt a yt xt 0 1 ẏt a yt xt 2 e at ẏt a ytdt e at xtdt The term e at is called the integrating factor. The reason for multiplying by the integrating factor is that the term inside the left-hand side integral becomes the deriative of e at yt with respect to time: e at ẏt a yt d dt e at yt b 0 - Hence, the term on the left-hand side of Eq.(2) equals e at yt b 0. 3 Compute the integral on the right-hand side of Eq.(2). Call the result INTt b 1 4 Mutiply both sides by e at to get yt : yt e at INTt be at where b b 1 b 0 is an arbitrary constant. Problem Sets Ex.1 Show the general solution to ẏt yt 1 0 [Answer: yt 1 be t ] Ex.2 Suppose k t kt k and k0 is the initial value of kt, show that kt k e t k0 k. (Romer, 2001, p.24) Linear, first-order differential equations with variable coefficients. ẏt at yt xt 0 The integrating factor is now e ad 0, so that the left-hand side becomes the derivative of yt e 0 We find that the solution to the ODE is yt e 0 t ad e 0 t t t ad ad. xt dt b e 0 t - Where b is an arbitrary constant of integration. Systems of Linear ODE A system of linear, first-order ODEs of the form ad,

ẏ 1 t a 11 y 1 t...a 1n y n t x 1 t,... ẏ n t a n1 y 1 t...a nn y n t x n t. In matrix notation, the system is ẏt A yt xt, Where yt ẏt y 1 t y n t ẏ 1 t ẏ n t A is an n n square matrix of constant coefficients xt is a vector of n functions. Solutions Methods Phase diagram Analytical Numerical Phase Diagram Diagonal systems: A simple case ẏ 1 t a 11 y 1 t, ẏ 2 t a 22 y 2 t, where a 11 and a 22 are real numbers. Case 1, a 11 0 and a 22 0 : An unstable system - Step One Step Two

Step Three Step Four, use the boundary conditions to see which one of the many paths depicted in the picture constitutes the exact solution. Case 2, a 11 0 and a 22 0 : A stable system Case 3, a 11 0 and a 22 0 : A saddle-path stability.

A nondiagonal example [see BS, p.476)]\

CHAPTER 5 Dynamic Optimization Dynamic Control Typical Problem max ct V0 0 T vkt, ct, tdt s. t. k t gkt, ct, t (1) k0 k 0 0 (2) kt e rtt 0 (3) where V0 objective function as seen from the initial moment; rt average discount rate T terminal planning date, finite or infinite kt state variable ct control variable v instantaneous felicity functions Eq.(1): transition equation or equation of motion Eq.(2): the initial condition Eq.(3): the final constraint - For finite values of T, this constraint implies kt 0, as long as the discount rate rt is positive and finite. If kt represents a person s net assets and T is the person s lifetime, then the constraint in Eq.(3) precludes dying in debt. - If the planning horizon is infinite, then the condition says that net assets can be negative and grow forever in magnitude, as long as the rate of growth is less than rt. This constraint rules out chain letters or Ponzi schemes for debt. - Example vk, c, t e t uct k gkt, ct, t fkt, t ct kt Procedure to Find the First-Order Conditions 1. Hamiltonian function: H vk, c, t t gk, c, t. 2. Take the derivative of the Hamiltonian with respect to the control variable and set it to 0: H v g 0. c c c 3. Take the derivative of the Hamiltonian with respect to the state

variable and set it to equal the negative of the derivative of the multiplier with respect to time: H v g. k k k 4. Step four (transversality condition): Case 1: Finite horizons. Set the product of the shadow price and the capital stock at the end of the planning horizon to 0: T kt 0 Case 2: Infinite horizons with discounting. lim t t kt 0. Case 3: Infinite horizons without discounting. In this case, we use Michel s condition, lim t Ht 0 Present-Value and Current-Value Hamiltonians Most of the models that we deal with have an objective function of the form, T vkt, ct, t dt 0 T e t 0 ukt, ct dt Constructiong the Hamiltonian H e t uk, c gk, c, t. The shadow price t represents the value of the capital stock at time t in units of time-zero utils. Rewrite the Hamiltonian as H e t uk, c qt gk, c, t, where qt t e t. The variable qt is the current-value shadow price. Define H H e t to be the current-value Hamiltonian: H uk, c qt gk, c, t. The first-order conditions can be expressed as 4 H c 0, 5 H k q q. The transversality condition can be expressed as qt e T kt 0. Multiple Variables A general dynamic problem:

max c 1 t,,c nt T uk1 t,...,k m t; c 1 t,...,c n t; t dt 0 s. t. k 1t g 1 k 1 t,...,k m 1; c 1 t,...,c n t; t k 2t g 2 k 1 t,...,k m 1; c 1 t,...,c n t; t... k mt g m k 1 t,...,k m 1; c 1 t,...,c n t; t k 1 0 0,...,k m 0 0, given k 1 T 0,...,k m T 0, free The Hamiltonian is H uk 1 t,...,k m t; c 1 t,...,c n t; t i1 FONC H 0, i 1,...,n, c i t H i, i 1,...,m, k i t The transversality conditions are i t k i T 0, i 1,...,m. m i g i. Dynamic Programming: An Introduction Example: An optimal growth problem max Ut c t,k t t uc t t0 s. t. c t k t1 fk t Where ut instantaneos utility, which is an increasing, concave function of current consumption t individual s subjective discount factor, t 0, 1 c t current consumption k t1 capital to be carried over to the following period c t k t1 fk t where c t current consumption k t1 capital to be carried over to the following period Value function The value function represents the fact that the maximum present

discounted value of the objective function from a point in time onward can be expressed as a function of the state variables at that date. Since the state variables at a point in time fully determine all other variables both currently and (via the transition equations) at all future dates, including those which enter the objective function, they determine the maximum attainable value of the objective function. Example: k k t t + k 1 t+ 2 ct c t + 1 ( ) t u c uc ( t + 1) - The state variable k t determines k t1, which determines k t2, et cetera. - k t therefore determines the utility-maximizing value of c t, c t1, c t2, et cetera, and this maximum attainable value is simply Vk t. - Therefore Vk t max Ut c t,k t t uc t t0 Bellman equation 1. In any period t, where the planner begins with k t, the choice between c t and k t1 can be represented as maximizing uc t Vk t1 2. Combining two observations, we have Vk t maxufk t k t1 Vk t1 (eq.4) k t1 This is the Bellman equation, whose solution is a function V. What lies behind the Bellman equation? 1. The value function allows a nontrivial dynamic problem to be turned into a simple-looking single-period optimization. 2. Dynamic problems are two-period problems balancing "today" against "the infinite future", but this works only when there is consistency between how one treated the future yesterday and how one treats it today. Methods of solution

1. Method of conjecture. 2. Method of successive approximations. The most useful method via the Bellman equation is not to solve for the value function V at all, but to derive the optimal path without finding V itself. Differentiability of V allows us to do this in many cases. With smooth functions and interior solutions, marginal changes in the state variable k imply marginal changes in attainable welfare in the same direction. With differentiability of V, the maximum problem in (eq.4) yields the following first-order condition for k t1 : u fk t k t1 V k t1 (eq.5) Straightforward interpretation: choose k t1 so that the loss in utility from "one less unit" of consumption today is just equal to the (discounted) gain in future "one more unit" of capital carried over (i.e. "saving") would allow. What is the utility gain from higher k t1? This can be calculated by use of the envelope theorem applied to V to find the derivative V, namely, u fk t k t1 u fk t k t1 f k t (eq.6) Straightforward interpretation: the value of another unit of beginning-of-period capital along an optimal path is the marginal utility value of the extra product the higher capital allows evaluated at the optimal level of consumption. Combining (eq.5) and (eq.6), one obtains u fk t k t1 u fk t1 k t1 f k t1 (eq.7) where k t1 k t2. Hence the time-invariant function giving the optimal level of the control variable as a funtion of the state, sometimes called the policy function, which characterize the optimal path. For many functional forms of u and f, it is easy to solve (eq.7) to obtain a closed form for. More generally, we can analyze this equation to obtain certain characteristics of and hence the optimal path.