|
第十届华为杯全国研究生数学建模竞赛--第四轮模型建立 - f( v3 S: \2 c: u$ S6 {; H
" u7 Q1 M! T8 Z2 u( h
8 I6 v2 g$ g) e4 p第十届华为杯全国研究生数学建模竞赛 --第四轮解析之模型建立与求解
1 I G5 D" {9 T( e
6 K1 }8 u5 S+ S3.1 问题一模型建立与求解 $ L1 n; H. R0 ?& |. q5 B! F& G+ A$ `
3.1.1 新洛伦兹模型构建
7 Y& B8 o4 z l- p假设基于统计分布表示的收入分配的密度函数是 f ( x),其中 x 表示收入,对应的分布函数为 F ( x),则 p = F ( x)表示收入低于或等于x的人口比例。记收入低于或等于 x 的人口群体拥有收入占总收入的比例为L(p),则应有 ! V h, ?4 S# u% P0 B- M
1 x L ( p )=mò0 tf (t ) dt , p = F ( x) file:///C:\Users\ADMINI~1\AppData\Local\Temp\ksohtml9080\wps1.png 7 S8 [7 A4 x8 |. t( a, O9 E
L ( p)称之为收入分配的洛伦兹曲线。收入分配的洛伦兹曲线 L ( p)在收入分配分析中具有重要地位,它表示人口份额等于 p 的低收入端拥有的总收入份额,因此 L ( p)是定义于[0,1]区间上的函数。按经济意义,它应满足如下条件:
1 f2 W( ]% u- jL(0,t)=0, L(1,t)=1, L¢( p,t)30, L¢¢( p,t)30 * X$ G7 R$ B/ \- `
即 L( p,t)在[0,1]上是凸增函数。
8 T8 J0 F$ h7 S7 `/ Q在分析与测算洛伦兹曲线的实际工作中,在只有分组数据可用的条件下,可以先估计收入分配的密度函数,从而得到相应的洛伦兹曲线,或直接估算洛伦兹曲线。国内外学者所提出的模型可以概括为三大类:几何计算法、分布函数法和曲线拟合法。几何分析法,是根据分组数据刻画洛伦兹曲线,利用这一方法不能得到洛伦兹曲线的表达式,只能用来计算基尼系数,但由于在计算分块面积时用直线近似地代替曲线,所估计的基尼系数要小于实际值,尤其在数据点较少时,误差较大;分布函数法,是基于对指标的概率密度函数或概率分布函数的假设,估计其分布参数,然后对洛伦兹曲线进行估计,这类方法较为复杂,同时由于计算收入分配的概率密度的复杂性,很难提出合适的概率函数;曲线拟合法首先假设收入分配服从某一特殊的统计分布函数形式,如对数正态分布、伽马分布、贝塔分布、帕累托分布、威布尔分布等,再据此为洛伦兹曲线选择适当的参数方程直接进行拟合、确定参数,由此估计出洛伦兹曲线。 1 n g& M: C, c; `8 `' o
为了更准确地描述洛伦兹曲线和精确地估计基尼系数,本文通过分析洛伦兹曲线的特性,构建新的洛伦兹曲线模型,对洛伦兹曲线直接进行拟合。根据洛伦兹曲线应满足的性质入手,在查找相关文献的基础上,总结前人所做研究,构建 出满足洛伦兹曲线的新模型 L ( p,t ) 。 Chotikapanich(1993)提出较早的洛伦兹模型:
' |) f# i: G( D' L, B Q6 2 h( {) F' @2 t
, l5 d9 [# ~8 c$ r8 \2 B: c( L# m5 i7 [5 V% X( A) f
| | |
1 v3 ~$ W, \" _" w( d
| | ( p6 {" L+ e J& y4 I/ \
|
, }0 \1 I- ^* P( N7 `: O b" q
| | |
/ O( k4 p: X! h9 |' a1 S
|
$ {! A, P5 p6 q. u& WSarabia 等提出了基于经典帕累托分布的洛伦兹模型,其基本形式为:
' u) j n9 z; D& U: D& hL ( p ) [size=10.5000pt]= 1 [size=10.5000pt]- (1 [size=10.5000pt]- p) b , b [size=10.5000pt]Î[0,1] | |
* V0 J; A# x& B6 W& V- ]
| & {: m. \2 A: i4 A8 P" W9 P3 J
| : }( E9 U# D( _. h
| ) i9 f) v# _( ]7 b" a! p
从(1)式出发,他们有提出了一族广义帕累托族的洛伦兹曲线模型,具体形式为: B1 g4 y9 A, s% U* o" b+ x& t
L ( p )= pa(1-(1- p)b),a30,bÎ(0,1] | 1 n/ E$ Z; n+ `. x n* B" ]
| ' @" z4 B( s& f3 e# x
| |
A8 s% o3 k. O6 }8 W9 V# i; I
|
2 q" R: \4 C) o, v
|
+ x/ |2 ]' z% Q! @* k
|
3 d4 \+ s2 P# |$ i+ j2 I
|
5 \ Q; R" j3 q$ B/ O" v. [
| ) N+ D' T! P# O, P% [
| / Z0 Y* W! L+ i- p3 G% w
| & y0 j" T! G3 g3 z& Y+ s
|
9 s+ ~4 F- ]5 u. R" |2 v6 T
|
: H I& q% @5 b
| 2 M0 @7 g& _7 `
| & h% Q: Q1 v& a _$ F" x( }' j
| ) L/ Y" L6 d- y- l* q( g; V, m2 d
| L ( p )= pa(1-(1- p)b)n,bÎ(0,1] |
7 C: t0 n* R6 m+ { V
|
9 F: H# @9 X/ ^9 M3 z8 ]
| : G) H4 l" D9 k/ N- X
| |
% x. S5 A$ B# t# h# l/ Y! l
| L ( p )=(1-(1- p)b)n,bÎ(0,1] | % }& M" ^: W: @! c
|
, h8 X$ b# G9 R% A
| i' B8 y8 R& H" o& ?( E2 H) W
| | 3 o* s9 r) v) u( y9 X. v
| * l. v- q8 L" W7 C9 [& |: {% n
| 3 C( F+ }! \0 T' z7 [# D5 i
|
4 n9 I* D5 y, w7 b
|
; V% u& r7 Y1 d, j$ W' u
|
2 D& M7 S& y w+ d0 Y- @, m) R: v
|
( U* J: `/ D0 a; U& l; K, V# L
| 7 h) h. P8 S( }, h0 z
| % R; ^9 K# L; P h3 x
|
5 F5 x7 n1 l F* e
|
' c5 P) m" o$ S7 x" H- h# Q
|
8 B+ H" ]5 D5 O1 ~2 q
| . |: v" s X3 G+ T# a7 z* N
| | ! L$ G$ ^% Y5 W, W, P& i( @. j# Q
| . T7 ~0 p }% F0 X9 H
|
; Y0 f4 a4 [8 l# w# s$ @
| * @; E4 C/ ^3 Z. M& q7 m) V: s
| 7 z! m# S( w4 L
| : U: |% q) Q; y' V8 L
| | | | , b Î (0,1], l Î ( -¥, 0) è (0, ln b | | 5 D' ^% \& c4 O% t) t3 P& s' j+ K
| 4 J0 [& b% G' v1 G2 m( j& n$ V
|
7 e4 T( n' A2 B% C/ m
|
M7 a8 N/ f7 W+ P6 w! V! f
|
9 T2 s" A' m4 `) Z2 T0 A
| |
) h4 G8 ^& ~/ I; H$ a) r
| |
* p$ F' n3 `3 T6 f, F- ]% _
|
. D4 c7 x+ ?& s
| 2 `% q- B2 |0 {" V. h5 o
| + f" [& q( E: Q8 q2 E
| + n4 ]# H; o T/ h" u& u7 s3 D ?
| ) o$ J% U$ E# ^( h
|
* ? M5 e! y: ]
| | ' }/ I( x% Z8 f& e
|
( F. d T" l6 K# y! \5 X* M
| ( d Y4 |% D7 W! k$ I2 |' p
| | |
, b+ M) {, v0 V+ P9 w
| 1 j& ]& z( g3 D) R5 R0 e, s D
|
2 @ v) x- E5 G5 I7 u
| |
# r2 p0 v- I& r3 x7 \
|
4 J+ w9 A- m# x; p. z9 T1 b
| % D/ {- O ?# o# k5 Y d
| & Q7 U$ u& c) m
|
; H3 y2 `0 y- [: z
|
6 j( C z3 D: w2 v7 U, v. l
| " O7 y9 Y( a' {/ D3 k. p
| | & C1 ]3 c3 {/ t2 d3 R! C' X
| 6 s( C, J, G# I$ N- d# ^
|
7 C; W; j9 w f/ S
| | | | |
6 i& c5 [4 f# x& f& ]
|
" I4 l. t! J( U- V6 M$ D
| | . G8 h1 h7 J0 Y. g
|
+ g2 ]& x6 ]5 d: o
|
+ I$ h3 _6 z7 `. `* I
|
# N! j5 B- \* N" F
|
0 y- E; h- z& k; J; y
|
* L6 S ~5 @1 a5 M3 r3 z. T9 Z
| | % H- a5 ^3 H4 p7 {* I" D7 k4 F: C
| |
( V3 v$ }- _. T0 B& h1 H3 e. n* i
|
7 U# F) Z! {* J- C
| % n6 h! g. B# P2 }
| |
8 q: n8 m) W! k5 l5 C) x
| o5 f- F" W ~4 o
Ogwang 和 Rao(2000)曾提出用两种混合的方法建立洛伦兹模型:加权积(weighted product)及洛伦兹模型的凸组合(convex combination)。并且得到凸组合模型(8)满足洛伦兹曲线性质:
8 W3 [# N( J1 a: F* C | 3 R) r0 U* Y. C
|
* @9 ^/ r; l+ D3 @
| |
% y- s& x0 O+ h+ a* z* }3 _4 w
| |
1 b) k' q$ y: U% I# e- n. _9 c$ m
| | * @8 c7 R3 p$ @& F
| | | | | ) j3 C4 @$ C- ]
| | ' k) L3 C; r7 B; [: _* k
| | | . M& D- b& F6 S: V8 c. g8 A
| |
' C$ @4 ^( W9 W. X; T
|
更一般化的经典洛伦兹曲线模型为王祖祥(2007)提出的如下的二元参数模 6 N, X' l. |9 p( z7 g: l& ^
型:
1 x2 O+ T' Y& t* X. @' H4 W | ( p ) [size=12.5000pt]= 1 [size=12.5000pt]- (1 [size=12.5000pt]- p ) b e[size=14.0000pt]-g p | | | . U- o7 ^. W( ]$ m
(9)式所表示的二元参数模型作为一个洛伦兹曲线的参数估计模型使用时
8 O. u/ K# H% Y* z0 o7 d1 L4 k" J比基于帕累托的广义洛伦兹模型具有更好的性质。
* d; n5 X$ `; n: m4 ?现在的研究成果已经证明:假定 L ( p) 为洛伦兹曲线,则对于任意的a 3 0 和 V# D3 E. _# g% ~; U& x
| |
0 f) x- b$ w; E1 F' B: L
| ! H& Q% ]. Z. T6 B
| |
2 X1 W7 W. r Z+ T+ P3 g$ [# \2 K
|
! @% g& } @4 w: H4 b
| | | | ) q& K0 |! P- f, w! [; b( F* U
|
, F# `1 O) C0 L9 F
| | - j7 b# D& j+ z, Q- {( u+ [& g
|
7 ~ k6 H; N! j9 p0 n3 Y
| 都有 L' ' '( p)30,则当a30,h31 / 2 | |
/ d, w- w- F+ G
|
6 [9 J0 U- |/ I# m4 g* j
| / b( F! {* R; {1 J/ g
| 且a + h 31时,L ( p) 也是洛伦兹曲线。 | - L, t) g& j2 j% B
| 6 w7 W$ M+ R1 i# H
国内外的专家学者在研究洛伦兹模型方面做了大量的工作,王祖祥、 $ B `+ V' C- M, [/ A% z
Sarabia 等提出了一系列的洛伦兹曲线模型,除上面提到的公式(9),如:
/ x1 J* J$ m+ l8 q5 V4 Y. `" M- ^ P
/ q/ ~; n5 F2 ]5 \) o, o0 c" p
7 8 q' w4 F/ r U
' a( V/ b2 b7 y7 d* }L ( p ) [size=11.5000pt]= pa | 3 \, U! B* i6 }2 |5 s' v# c
| | : v0 j* ^7 B+ a) H( k4 K
| | 3 i! M0 \& u5 B, X2 S2 Z! u* r
|
/ b7 l K- u+ P
|
8 B6 S' X( @; m) r$ I E% M* k1 u
| @3 m% @# k: q5 Q* r" R/ |
| # M# y6 k1 D5 ?
| $ \1 o: ?" v( Y3 }% j, Z" G
| . p4 f( d3 ~' n' F) [ Q
|
2 |( E2 ?/ Z# ?* C
|
9 K# Q8 ?0 n1 G. C' |# K
|
' W5 X3 `6 [0 q5 \$ @
| , f$ E, i- o: X+ Y7 A8 @# K! n
| |
5 J5 x) O* T% B+ z! b- i8 `/ ?
| | 8 r/ }# |% r& c) c- I2 D3 Q
|
% f& e, f- ?/ m6 h& i
|
S6 a+ C9 o- F* b+ F
|
# D" A0 [$ ~! Q5 O) M' E
| ) r& N& Y3 P+ _* ^/ _
| w; k3 \) V' U |& P0 D4 W$ z) r2 f
|
9 \: b( l/ o% @3 n7 i
| 2 d. H: W0 ?+ X" w( O# w
|
) [5 I: s$ A# u! g0 V
| | . l8 d" }$ v. { `* R) U
| | ( I3 u' ]& ]6 Y& f" g" U4 O
| 4 }( j! D8 a7 C) e5 j* u
|
9 n- h; c0 v3 [; f/ O
|
9 P- Q5 T# P6 j
| 1 V" E3 }' t( E; p1 e
|
' Z0 S: |7 i# j6 O7 y6 q
|
8 C0 X8 c" E3 u( n2 G
| . G- Y- f* ]$ }8 c
|
* I( N/ H$ q1 R
| ) z/ x( _3 o& q: p5 G' P
| | / j5 X' V4 W1 ~ B! l* Z
| 2 ]$ E& ~9 o: w. y* s' a y
| . {5 d- e; l5 _' N- S, U
| |
& r3 s0 H' e4 x# q3 c- N1 d8 m1 M
| ( V" G8 n8 T C; Y6 q& K
| $ b( S9 Y2 ]% z1 H" L4 Q% b
| |
y: ~+ C" R3 y
| : I- D) C8 g2 c; H* t
| |
$ O, L1 S. a6 d$ ~) G/ o! Z x& w
|
6 S( b( Y* I0 t( g5 c8 o7 w
|
. w, v0 |3 H% B; T+ ~: V( w! y K
|
/ i& S, {: S6 U% o
| ( _6 M& k, n7 |" b9 E5 H1 V
| |
2 e! _5 _4 a: c" T9 T5 Y0 ^3 s
| | | 5 H4 z* k0 z9 A2 k8 F1 r
| 8 D% q5 C# f& j/ `
| . ?& b: ?9 A0 I: P) _
| |
. [0 `5 M; R, S
|
2 b) B0 v' u( t1 V6 X7 _2 j
| - G* }5 I& \, ^* q4 ?- D
|
! W5 `0 ~- \# n6 b5 b: E) w' G
| | |
% i* N! ?3 @2 k0 C! E7 c
|
7 w2 ], Q; c- a. b6 ?( m- K+ `
| 9 F/ k. D! L5 r5 Q, W+ O2 i9 A! O
| | " c0 M; D0 T C! W g4 R& i/ t
| | 0 m9 q4 h5 e- Q, h/ Y
| # P) |* s. T% f [) E
| $ \* h. _9 }' F w, a
|
& l9 p, ~+ R: n& _! ?0 K
| ( H' i" O6 S& z: T' l
| " q9 m3 N. O4 U; { h" A* G
|
V8 }6 t5 E `- F
|
; q" Q1 f9 q7 e3 w
| |
7 q! B4 U5 Z8 a. y
|
1 c& `& X" W9 z0 J- y
| ) t/ e5 h- S+ V+ ? U! l
|
) h$ i2 W3 \# I* n0 d: i
|
$ W5 f/ A* B; `' g6 f+ y4 \; u
|
0 |. W6 c" ^7 P* U8 S' L1 ~
| |
# o, P, w2 s. p( @# t. X
| L ( p ) [size=12.0000pt]= [1 [size=12.0000pt]- (1 [size=12.0000pt]- p ) b ]a [1 [size=12.0000pt]- (1 [size=12.0000pt]- p)h ]l | $ J+ _( n& ?/ A* s
| $ V* V, G ?, w2 B8 b
| " }( J4 a# c5 E/ ?
| 3 H0 m$ \, d( X9 Y: M# k& U4 w
| " O: g: s3 ?3 s6 m9 s; O, Q$ R3 s
| |
/ w: P3 n4 g* M. ^# |
| L ( p ) [size=12.0000pt]= p | | | | | | 2 P- R% s' {: y; p+ w0 j! g% f+ f
| | | | | |
. _) s$ C5 w3 b: p
| 7 R, u, e6 [- R
|
* p0 D* L: `6 F+ L' `5 x$ t5 q
|
) P! k" x% c$ L/ V: m& F
| % P/ ]" O% L' w% |9 A
| | 2 N9 G7 |+ F) _2 ]" h2 }
|
1 F- H8 N; [5 j8 S2 d0 [
| | | 4 b) }6 v& p4 {+ Y
|
9 N0 H# Z! g% r, B5 i
|
在研究现有文献基础上,按照洛伦兹曲线应满足的性质以及推导定理,本文构建了基于指数成份的洛伦兹模型 L( p,t ) ,形式如下: 1 x8 U( o/ @# G% x h
L ( p )=[1-(1- p )a e -bp][1-(1- p )g e-hp] | |
1 [+ X6 `6 I( U
|
% G9 A0 @' W: C4 Y( x, A
| ( B2 Z+ U1 R8 G7 o" H) q0 c5 v3 v
| 其中, 0 £ a + b £1, 0 £ g + h £1。 |
7 Z$ [! z7 I; L* G
| 7 k3 e8 z2 c& i9 A; @
| |
/ p3 \8 Q. R' Q3 Z: }: O9 Q$ v
| $ V, ^, g" j+ Z( ?0 }: w
| | 4 E$ _9 f* ^" t* I* j; D8 w
|
6 v4 ~( [# Y9 q' y( P" g B
| |
e8 z; X7 i5 F( `( {
|
+ z5 R: D6 F! J+ p& U
|
/ E. q6 y4 k$ v8 |L'( p )=[a(1- p )a-1 e -bp][1-(1- p )g e-hp]+ ( o7 \& H4 s2 F8 O6 Y' ]
[1 - (1 - p )a e - b p ][g (1 - p )g -1 e -h p + h p (1 - p )g e-h p ] ; a7 O9 R+ Q; t1 g% J% T
在满足条件 a + b 3 0 , g + h 3 0 的情况下,当 p Î[0,1] 时,满足
% p- a! C. ^' C! F) hL'( p )30。
: P- ` @! |0 z. i, j$ s1 RL' '( p )=[a(1-a)(1- p )a-1 e -bp-ab p (1- p )a-1 e -bp+b(1- p )a e-bp
" K8 [8 K9 s4 ^% y$ c- ab p (1- p )a-1 e -bp-b2 p 2(1- p )a e -bp][1-(1- p )g e-hp] / `+ C! r$ ^+ C0 e3 l
+ 2[a (1 - p )a -1 e - b p + b p (1 - p )a e - b p ][g (1 - p ) g -1 e-h p 8 \) C9 P7 r+ h* `! Y% W
+ h p (1- p )g e -hp][g(1-g)(1- p )g-1 e-hp 1 N5 l a1 C- ` h- j h0 Z: u. L% A% M
- gh p (1- p )g-1 e -hp+h(1- p )g e-hp ! f5 P0 M# _; x2 {" v/ X9 R
- gh p (1- p )g-1 e -hp-h2 p 2(1- p )g e -hp][1-(1- p )a e-bp]
6 e6 f; o) C7 ]- R+ z. k. ?) l) q' h# y, a= [(a + b ) - (a + b ) 2 - b 2 p 2 (1 - p )](1 - p )a -1 e-b p " `9 E& T' ^0 }+ Y$ y/ B) _
+ [(g + h ) - (g + h ) 2 - h2 p 2 (1 - p )](1 - p ) g -1 e-h p
* V5 q& A) |3 ]: k当满足条件a + b £1,g + h £1的情况下,当 p Î[0,1]时,L'' ( p) 3 0 。
: l3 B4 Z1 n. d/ l8 k综上所述,新构建的洛伦兹模型模型(14)式满足洛伦兹曲线的定义及性质,
* y7 Y1 T3 ]8 N( K" I" `8 d8
2 p( {1 }! j$ I$ i6 |
i! b6 Q3 R9 t6 z6 N; x: u7 \" O7 }
可以用来拟合题目中给出的数据,并且与现有经典洛伦兹模型做出比较。 1 \& M$ U: U3 a9 e/ H) e6 E& p
3.1.2 模型的计算与比较 1 L& f# k1 j Z5 ?
有关收入与人口的数据一般情况下可以得到收入人口分布的分组数据,这种数据的完整形式为(pi ,xi file:///C:\Users\ADMINI~1\AppData\Local\Temp\ksohtml9080\wps2.pngm),i=1,2, ,n或者(pi,Li),i=1,2, ,n,其中xi是收 9 G- e( F# z: N9 L( n
入区间点,满足 0 £ x1 < x2 <
, a8 h4 z( P$ W! Z5 Z# K |