3广东的数据及拟合重构函数如下： 函数代码：：
x=[316 1067 1691 2118 2209 2662 2710 3076 3197 3215 3296 3496 3658 3842 3849 3996 4199 4378 4514 4544 4850 5771 5919 5996 6039 6327 6565 7888 8532 8686 9140 10321];
y=[26450 14486 94410 11023 27297 57037 32214 63786 27359 24352 127965 43888 125220 44679 21526 27620 54836 64010 13704 44708 64807 50412 37608 109281 40081 64867 46543 123986 42011 88501 97935 94805]; nh3=polyfit(x,y,3); subplot(2,2,1);
plot(x,y,'*',x,polyval(nh3,x),'b-'); legend('实验数据','拟合曲线'); title('广东'); subplot(2,2,3); x1=[100:100:8500];
y1=2.3375*10^(-7)*x1.^3-0.003*x1.^2+16.7*x1+1917.6; plot(x1,y1)
title('重构曲线 短期总成本：Ystc'); subplot(2,2,2);
y2=2.3375*10^(-7)*x1.^2-0.003*x1+16.7+1917.6./x1; plot(x1,y2);
title('短期平均成本：Ysac'); subplot(2,2,4);
y3=2.3375*3*10^(-7)*x1.^2-0.003*2*x1+16.7; plot(x1,y3);
title('短期边际成本：Ysmc'); 相应的函数图形如下：
14121086420x 104广东40实验数据拟合曲线3530252015105短期平均成本：Ysac05000100001500002000400060008000100008765432104重构曲线 短期总成本：Ystcx 10短期边际成本：Ysmc1816141210864020004000600080001000020200040006000800010000
4广西的数据及函数拟合重构如下： 函数代码：
x=[316 1067 1691 2118 2209 2662 2710 3076 3197 3215 3296 3496 3658 3842 3849 3996 4199 4378 4514 4544 4850 5771 5919 5996 6039 6327 6565 7888 8532 8686 9140 10321];
y=[26450 14486 94410 11023 27297 57037 32214 63786
27359 24352 127965 43888 125220 44679 21526 27620 54836 64010 13704 44708 64807 50412 37608 109281 40081 64867 46543 123986 42011 88501 97935 94805]; 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.3375*10^(-7)*x1.^3-0.003*x1.^2+16.7*x1+1917.6; subplot(2,2,3); plot(x1,y1);
title('重构曲线 短期总成本：Ystc');
y2=2.3375*10^(-7)*x1.^2-0.003*x1+16.7+1917.6./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'); 相应的函数图形如下：
14121086420x 104广西40实验数据拟合曲线3530252015105短期平均成本：Ysac05000100001500002000400060008000100008765432104重构曲线 短期总成本：Ystcx 10短期边际成本：Ysmc1816141210864020004000600080001000020200040006000800010000
5贵州的数据及数据拟合函数重构如下：
函数代码：
x=[362 867 882 1266 1553 2156 2271 2468 3007 3367]; y=[17527 14259 17646 25004 15329 15052 21951 16663 21631 21544];
nh3=polyfit(x,y,3); subplot(2,2,1);
plot(x,y,'*',x,polyval(nh3,x),'b-'); legend('实验数据','拟合曲线'); title('贵州');
x1=[100:100:4000];
y1=8.2142*10^(-7)*x1.^3-0.0040*x1.^2+6.367*x1+14953; subplot(2,2,3); plot(x1,y1);
title('重构曲线 短期总成本：Ystc');
y2=8.2142*10^(-7)*x1.^2-0.0040*x1+6.367+14953./x1; subplot(2,2,2);
plot(x1,y2);
title('短期平均成本：Ysac');
y3=8.2142*3*10^(-7)*x1.^2-0.0040*2*x1+6.367; subplot(2,2,4); plot(x1,y3);
title('短期边际成本：Ysmc'); 相应的函数图形如下：
x 104贵州160实验数据拟合曲线140120短期平均成本：Ysac2.62.42.221.81.61.4100806040200100020003000400000100020003000400034x 10重构曲线 短期总成本：Ystc短期边际成本：Ysmc1412102.58624201.501000200030004000-201000200030004000
6河北的数据及数据拟合和函数重构如下： 函数代码：
x=[788 1754 1912 2843 2892 3222 3342 4438 4844 4861 5266 5362 6396 6452 7147 7575 7758 7926 8128 8935 9550];
y=[18790 21074 27063 25989 30605 14220 29304 51701
22415 38926 38449 32561 38915 43431 69814 32117 48696 69844 84309 97443 74039]; nh3=polyfit(x,y,3); subplot(2,2,1);
plot(x,y,'*',x,polyval(nh3,x),'b-'); legend('实验数据','拟合曲线');
title('河北');
x1=[100:100:8500];
y1=3.1151*10^(-7)*x1.^3-0.0037*x1.^2+17.94*x1-600; subplot(2,2,3); plot(x1,y1);
title('重构曲线 短期总成本：Ystc');
y2=3.1151*10^(-7)*x1.^2-0.0037*x1+17.94-600./x1; subplot(2,2,2); plot(x1,y2);
title('短期平均成本：Ysac');
y3=3.1151*3*10^(-7)*x1.^2-0.0037*2*x1+17.94; subplot(2,2,4); plot(x1,y3);
title('短期边际成本：Ysmc') 相应的函数图形如下：
10x 104河北16实验数据拟合曲线短期平均成本：Ysac8146124102800200040006000800010000602000400060008000100008765432104重构曲线 短期总成本：Ystcx 10短期边际成本：Ysmc252015105020004000600080001000000200040006000800010000
7黑龙江的数据及数据拟合如下：
函数代码：
x=[252 323 884 1004 1008 1204 1474 1521 2210 2330 3066 3073 3077 3342 3554 3877 3916 4429 5184 5586 7560];
y=[16878 19431 28568 16336 18560 16909 10238 11392
28600 3258 37979 23735 27827 19951 34217 54775 54838 24137 37247 32167 94174]; nh3=polyfit(x,y,3); subplot(2,2,1);
plot(x,y,'*',x,polyval(nh3,x),'b-');