前言:

柯布道格拉斯生产函数(Cobb-Douglas Production Function)是微观经济学中最适合用Python进行可视化演示的知识之一。本文分别作出三种参数情况下的产量曲线,另外还有产量曲面与等产量线的绘制。该生产函数基本形式如下:

代码实现:

import matplotlib.pyplot as plt import numpy as np n = 100 #劳动与资本投入量 Labor = np.linspace(0,10,n) Kapital = np.linspace(0,10,n) #技术水平参数 Tech = 3 Output1 = [] Output2 = [] Output3 = [] OutputN = [] #生产函数参数设置 for i in range(len(Labor)): labor = Labor[i] kapital = Kapital[i] #三组指数 a1 = 0.5 b1 = 0.3 a2 = 0.5 b2 = 0.5 a3 = 0.7 b3 = 0.8 out1 = Tech * (labor**a1) * (kapital**b1) out2 = Tech * (labor ** a2) * (kapital ** b2) out3 = Tech * (labor ** a3) * (kapital ** b3) # 对照组曲线:生产要素增长 outn = Tech * labor Output1.append(out1) Output2.append(out2) Output3.append(out3) OutputN.append(outn) fig = plt.figure(1,figsize=(8,8)) #图形1:生产要素指数之和小于1 ax1 = plt.subplot(1,3,1) ax1.plot(Labor,Output1,'red',label='a+b<1') ax1.plot(Labor,OutputN,'black') ax1.set_ylabel('Output(Q)') ax1.set_xlabel('Input(L,K)') ax1.spines['top'].set_visible(False) ax1.spines['right'].set_visible(False) ax1.spines['bottom'].set_position(('data', 0)) ax1.spines['left'].set_position(('data',0)) plt.legend(['a+b<1']) #图形2:生产要素指数之和等于1 ax2 = plt.subplot(1,3,2) ax2.plot(Labor,Output2,'green',label='a+b=1') ax2.plot(Labor,OutputN,'black') ax2.set_ylabel('Output(Q)') ax2.set_xlabel('Input(L,K)') ax2.spines['top'].set_visible(False) ax2.spines['right'].set_visible(False) ax2.spines['bottom'].set_position(('data', 0)) ax2.spines['left'].set_position(('data',0)) plt.legend(['a+b=1']) #图形3:生产要素指数之和大于1 ax3 = plt.subplot(1,3,3) ax3.plot(Labor,Output3,'blue',label='a+b>1') ax3.plot(Labor,OutputN,'black') ax3.set_ylabel('Output(Q)') ax3.set_xlabel('Input(L,K)') ax3.spines['top'].set_visible(False) ax3.spines['right'].set_visible(False) ax3.spines['bottom'].set_position(('data', 0)) ax3.spines['left'].set_position(('data',0)) plt.legend(['a+b>1']) fig.suptitle('Cobb-Douglas Function',fontsize=14,weight='bold') # plt.savefig('Cobb-Douglas.png') #图形4:生产函数3D绘图,

叠加等产量曲线 fig = plt.figure(2) ax4 = plt.axes(projection='3d') #另外一组生产要素指数参数 alpha = 0.8 beta = 0.7 X,Y = np.meshgrid(Labor,Kapital) Z = Tech * (X**alpha) * (Y**beta) Production = ax4.plot_surface(X,Y,Z,cmap='rainbow',alpha=0.8) ax4.contour(X,Y,Z,zdir='z',offset=0,cmap='coolwarm') ax4.set_xlabel('Labor') ax4.set_ylabel('Kapital') ax4.set_zlabel('Production') fig.colorbar(Production, shrink=0.5, aspect=20) #图形5:单独绘制等产量曲线 fig = plt.figure(3) ax5 = plt.subplot(1,1,1) Output = ax5.contour(X,Y,Z,cmap='coolwarm') ax5.set_xlabel('Labor') ax5.set_ylabel('Kapital') plt.clabel(Output,inline=True,fontsize=12) plt.show()

图形效果:

三种指数条件下的产量曲线

产量曲面

等产量曲线,线内数值为产量规模

Python+经济学:柯布道格拉斯生产函数相关相关推荐

  1. matlab进行道格拉斯筛选,用MATLAB计算柯布道格拉斯函数,需要计算程序,急!谢谢大神...

    用MATLAB计算柯布道格拉斯函数,需要计算程序,急!谢谢大神0 产值.资金.劳动力关系函数问题具体内容增加生产.发展经济所依靠的主要因素有增加投资.增加劳动力以及技术革新等,在研究国民经济产值与这些 ...

  2. 柯布-道格拉斯生产函数

    柯布-道格拉斯生产函数(Cobb–Douglas production function)最初是美国数学家查尔斯·柯布(Charles Wiggins Cobb)和经济学家保罗·道格拉斯(Paul H ...

  3. 从柯布-道格拉斯生产函数看云计算对传统软件工程的影响

    随着软件技术的不断发展,软件的功能变得越来越强大,软件的结构也变得越来越复杂,需要处理的数据变得越来越庞大,计算的复杂度也越来越高.如何开发能够应对这样日益增长的需求的软件便成为了亟待解决的问题. 2 ...

  4. python经济学函数_有没有python计量经济学的教程?

    多元线性回归模型假设: 假设中国2013年各地区人均现金消费支出与工资性收入.其他收入之间的关系为: Y= β 0 =\beta_0=β0​+β 1 X 1 \beta_1X_1β1​X1​+β 2 ...

  5. python kafka kerberos 验证 消费 生产

    2019独角兽企业重金招聘Python工程师标准>>> [toc] 安装 pykafka github $ pip install pykafka$ conda install -c ...

  6. 【内附PDF资料】Python实现下载图片并生产PDF文件

    转载请注明出处:http://blog.csdn.net/gamer_gyt  博主微博:http://weibo.com/234654758  Github:https://github.com/t ...

  7. python经济学函数_Sargent的量化经济学(8):Python工作原理

    原标题:Sargent的量化经济学(8):Python工作原理 许文立 武汉大学经济与管理学院 xuweny87@163.com 许坤 安徽大学经济学院 kunxu2014@126.com Pytho ...

  8. python经济学函数_Sargent的量化经济学(6):Python基础

    原标题:Sargent的量化经济学(6):Python基础 许文立 武汉大学经济与管理学院 xuweny87@163.com 许坤 安徽大学经济学院 kunxu2014@126.com 本讲呈现了阅读 ...

  9. Python+经济学:资本边际效率与IRR

    前置知识:贴现与贴现率 由于投资不再是外生变量,因此在IS-LM模型的开篇,课本一般先引入资本边际效率(Marginal Efficiency of Capital,MEC)这一概念来探讨投资规模的决 ...

最新文章

  1. js 动态加载select触发事件
  2. JDK 1.6 HashMap 源码分析
  3. 计算机系统存数及取数方式,计算机系统结构总结2
  4. JS对url进行编码和解码(三种方式区别)
  5. 前端学习(3124):react-hello-react之对props进行限制
  6. Jython中文乱码问题
  7. 下班搜包被诉 苹果支付近3000万美元与加州门店员工和解
  8. OPPO正式公布ColorOS 6系统:Reno新机4月10日首发适配
  9. ACR2010_常规医疗环境下TNF拮抗剂对RA骨侵蚀的修复作用
  10. 13. ZooKeeper最佳实践
  11. 如何设计可靠的灰度方案
  12. iphone双重认证关闭不了怎么办_苹果手机双重认证怎么关闭-互盾苹果恢复精灵...
  13. 面授班命令记录(更新中)
  14. linux底层播放器,[转载]媒体播放器三大底层架构:MPC、MPLAYER、VLC
  15. 51ditu Maps API 类参考文档
  16. VMware 配置局域网内访问
  17. Android 切割图片
  18. 大数据项目实践过程笔记
  19. luoguP3518 [POI2011]SEJ-Strongbox
  20. 牛客网Chino with Expectation【思维+位运算】

热门文章

  1. Java中如何实现一个函数返回多个值
  2. 晨风机器人插件编辑器_晨风机器人文本编辑器
  3. LumaQQ2006的安装
  4. 独孤思维:减压还能赚钱?让压力变有利可图的项目
  5. 【Hyper-V】Windows10家庭版安装Hyper-V虚拟机
  6. 最大化 Google Chrome 浏览器
  7. 计算机与汽车智能信息化管理,谈计算机控制系统在汽车行业的应用
  8. Matthew Ball:为什么说AR/VR元宇宙是Meta和苹果的下一个战场?
  9. 百度地图-删除替换标注
  10. 成矿远景叠加分析matlab代码,深水潭金成矿地质特征及成矿远景分析