我们认为该高校需紧缩的招生数应该控制在20000~25000之间。
综上检验,我们发现我们的模型在一定的程度上可以很好的达到我们预期的效果。当然我们的模型还是存在一定的缺陷的,我们将会继续完善。
8、参考文献及网站
[1]张珠宝,数学建模与数学实验[M],高等教育出版社,2005
[2]高鸿业,西方经济学(微观部分)第四版[M],中国人民大学出版社,2007.3 [3]张圣勤,MATLAB7.0实用教程[M],机械工业出版社,2008.1 [4]朱旭等,MATLAB与基础数学实验[M],西安交通大学出版社,2008 [5]罗万成等,大学生数学建模案例精选[M],西南交通大学出版社,2007
9、附录
附录一:短期成本函数、短期生均成本和短期边际成本函数的求解
1.安徽的数据分析: 数据分布及拟合曲线: MATLAB代码:
x=[1088 1098 1594 1748 1765 1770 1883 2013 2532 2770 3061 3118 3217 3523 3833 3919 3986 4258 4969 5532 5632 5709 5771 6119 6180 6925 6998 8364]; y=[9053 11152 21336 15462 23805 23504 13536 12291 17619 25187 26427 21591 22706 32560 16259 32051 29438 19933 34347 33149 24545 38998 42994 30650 36124 44501 56392 70916]; nh3=polyfit(x,y,3);
plot(x,y,'*',x,polyval(nh3,x),'b-'); legend('实验数据','拟合曲线'); title('安徽');
8765432101000x 104安徽实验数据拟合曲线20003000400050006000700080009000
在上面选取4个点,带入三次函数方程,使用MATLAB解线性方程,得出相应的系数:
MATLAB代码如下:
A=sym('[1748^3 1748^2 1748 1;3217^3 3217^2 3217 1;5532^3 5532^2 5532
1;6925^3 6925^2 6925 1]'); B=[17300;24040;33400;44500]; x=A\\B
根据相应的系数重构三次函数(短期成本函数)如下: MATLAB代码: x=[100:100:8500];
y=2.323*10^(-7)*x.^3-0.0026*x.^2+13*x+1243.7; plot(x,y)
7x 10465432100100020003000400050006000700080009000
相应的短期平均成本MATLAB代码: x=[1000:100:8500];
y=2.323*10^(-7)*x.^2-0.0026*x+13+1243.7./x; plot(x,y)
相应的短期边际成本MATLAB代码:
x=[100:100:8500];
y=2.323*3*10^(-7)*x.^2-0.0026*2*x+13; plot(x,y)
综合绘出所有的函数曲线如下:
x=[1088 1098 1594 1748 1765 1770 1883 2013 2532 2770 3061 3118 3217 3523 3833 3919 3986 4258 4969 5532 5632 5709 5771 6119 6180 6925 6998 8364]; y=[9053 11152 21336 15462 23805 23504 13536 12291 17619 25187 26427 21591 22706 32560 16259 32051 29438 19933 34347 33149 24545 38998 42994 30650
36124 44501 56392 70916]; nh3=polyfit(x,y,3); subplot(2,2,1);
plot(x,y,'*',x,polyval(nh3,x),'b-'); legend('实验数据','拟合曲线'); title('安徽');
x1=[100:100:8500];
y1=2.323*10^(-7)*x1.^3-0.0026*x1.^2+13*x1+1243.7; subplot(2,2,3); plot(x1,y1);
title('重构曲线 短期总成本:Ystc');
y2=2.323*10^(-7)*x1.^2-0.0026*x1+13+1243.7./x1; subplot(2,2,2); plot(x1,y2);
title('短期平均成本: Ysac');
y3=2.3375*3*10^(-7)*x1.^2-0.003*2*x1+16.7; subplot(2,2,4); plot(x1,y3);
title('短期边际成本:Ysmc'); 函数曲线:
8x 104安徽30实验数据拟合曲线2520短期平均成本: Ysac6415210500200040006000800010000020004000600080001000076544重构曲线 短期总成本:Ystcx 10短期边际成本:Ysmc20151032100200040006000800010000002000400060008000100005
2重庆的数据及数据拟合及重构函数如下:
x=[373 807 995 1370 1895 2058 4561 4836 5257 5343 5351];
y=[5281 17456 23739 22595 19486 28997 26920 18060 28329 37739 46553]; nh3=polyfit(x,y,3); subplot(2,2,1);
plot(x,y,'*',x,polyval(nh3,x),'b-'); legend('实验数据','拟合曲线'); title('重庆');
x1=[373:100:8500];
y1=1.1187*10^(-6)*x1.^3-0.0097*x1.^2+27*x1+1328; subplot(2,2,3); plot(x1,y1);
title('重构曲线 短期总成本:Ystc');
y2=2.323*10^(-7)*x1.^2-0.0026*x1+13+1243.7./x1; subplot(2,2,2); plot(x1,y2);
title('短期平均成本:Ysac');
y3=1.1187*3*10^(-6)*x1.^2-0.0097*2*x1+27; subplot(2,2,4); plot(x1,y3);
title('短期边际成本:Ysmc');
相应的函数图形:
54310281064x 104重庆16实验数据拟合曲线1412短期平均成本:Ysac020004000600005000100002.521.510.505重构曲线 短期总成本:Ystcx 10短期边际成本:Ysmc1201008060402000500010000-200500010000