与中微相伴的一学期终于结束了，现将平时写了代码但没时间整理的内容整合一下。

用python学习中级微观经济学，顾名思义，即用py程序(主要是sympy模块)解题。

写这些小程序的初衷是做题目太烦了，而且题目明明有一个很简单的套路，可以很方便的用程序实现(由于不是计算机专业的，代码写的太水了，完全没有章法和美感Orz)

下面是学习目录：显示偏好弱公理(WARP)

计算消费束(Consumption Bundle)

斯勒茨基方程(Slutsky Equation)

博弈论(Game Theory)

交换(Exchange)

必备书目：《微观经济学：现代观点》第9版 范里安

钟根元《中级微观经济学学习指南》(第4版)

准备步骤：安装 python(实例使用的是python 3.6 32-bit)

安装 sympy 模块(sympy官方文档)

sympy模块是python中最好用的科学符号计算模块。它能够快速提供各种符号、高等数学、线性代数运算等等(可能我用到的功能不到其功力的5%)

安装sympy最简易方式：

pip install sympy

sympy模块介绍

WARP 显示偏好弱公理

WARP 显示偏好弱公理讲的是：如果一个消费者在能够消费得起消费束

时，选择了消费束

,则

，消费束1显示偏好于消费束2。

但是并不是每个消费者的消费行为都会满足这个公理，所以有一类题目就是求判断给定的消费行为是否满足此公理。

废话不多说直接上题：价格为($4, $12) 时，某消费者选择了(9 , 4) 消费束。在价格为($8, $4) 时，他选择了(2, 9) 消费束。则他的行为是否服从显示偏好弱公理？

【解析】价格为($4, $12) 时， 4x9+12x4=84; 4 x2 +12 x9 = 116 。说明在该价格水平下，消费者买不起(2 , 9), 因此才选择了(9, 4) 。

价格为($8, $4) 时， 8x9+4x4=88; 8x2+4x9=52 。说明在该价格水平下， 消费者买不起(9, 4) , 因此才选择了(2, 9) 。所以，消费者的行为没有违反显示偏好弱公理。

【思路】

判断是否满足WARP的方式，实际上是构造消费束矩阵C和价格矩阵P。将价格矩阵转置以后，与消费束矩阵点乘即可得到不同消费束+价格下的收入矩阵M。

根据消费者的实际选择，遍历收入矩阵即可判断该消费者是否满足WARP。

【程序】

import numpy as np

def is_warp():

global mat_star

bool = True

for i in range(info):

for j in range(info):

if mat_star[i,j] == 1 and mat_star[j,i] == 1:

bool = False

break

else:

pass

return bool

try:

info = int(input('请输入数据组数：'))

except:

print('给我一个正常的整数嘛~')

p = np.zeros((info,2))

c = np.zeros((2,info))

for row in range(info):

p[row][0]=float(input('输入第'+str(row+1)+'次价格(用空格或逗号隔开)(p1,p2) = '))

p[row][1]=float(input('输入第'+str(row+1)+'次价格(用空格或逗号隔开)(p1,p2) = '))

c[0][row]=float(input('输入第'+str(row+1)+'次消费(用空格或逗号隔开)(x1,x2) = '))

c[1][row]=float(input('输入第'+str(row+1)+'次消费(用空格或逗号隔开)(x1,x2) = '))

mat_m = p.dot(c)

print('数据表如下')

print(mat_m)

mat_star = np.zeros(np.shape(mat_m))

for i in range(info):

for j in range(info):

if mat_m[i][j]

mat_star[i][j] = 1

print('---------------------------')

print('是否\'*\'如下(1代表打星号)')

if __name__ = '__main__':

if is_warp():

print(mat_star)

print('Yes.')

else:

print(mat_star)

print('Hell no!')运行结果

小case了

下篇文章会比较有意思

【Consumption Bundle】有图有真相！