摘要：針對2015年全國大學(xué)生數(shù)學(xué)建模競賽題D題：眾籌筑屋規(guī)劃方案設(shè)計(jì)問題，在文中建立了規(guī)劃方案的核算與優(yōu)化模型。本文重點(diǎn)說明通過編輯MATLAB程序，對復(fù)雜的核算問題進(jìn)行計(jì)算，從而得出相關(guān)的核算數(shù)據(jù)，使用該程序核算非常方便，可推廣到其他有關(guān)增值稅的核算問題中。

關(guān)鍵詞：眾籌筑屋 方案 核算 程序

中圖分類號(hào)：TU12 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1007-9416（2016）07-0068-02

2015年全國大學(xué)生數(shù)學(xué)建模競賽題D題：眾籌筑屋規(guī)劃方案設(shè)計(jì)。在文中建立了眾籌筑屋規(guī)劃方案的核算與優(yōu)化模型，本文通過編輯MATLAB核算模型程序，對規(guī)劃方案Ⅰ、改進(jìn)方案Ⅱ和優(yōu)化方案Ⅲ進(jìn)行全面核算。

1 眾籌筑屋規(guī)劃方案的MATLAB核算程序

%計(jì)算每種房型的建筑面積d，開發(fā)成本C，售房收入R，各項(xiàng)成本XC，容積率r，增值稅Z，純利潤L，投資回報(bào)率ROI，%b是各種房型的建筑套數(shù)

function [d，R，r，XC，Z，L，ROI]=hesuan（b）

a=[77，98，117，145，156，167，178，126，103，129，133]； %每種房型的面積

c=[4263，4323，4532，5288，5268，5533，5685，4323，2663，2791，2982]；%每種房型的開發(fā)成本

p=100*[120，108，112，128，128，136，140，104，64，68，72]； %每種房型的售價(jià)

d=a.*b；%每種房型的總面積

C=c.*d；%每種房型的開發(fā)總成本

R=p.*d；%每種房型的總收入

r=sum（d（1：8））/102077.6；%容積率

d1=[d（1：8），d（11）]；%提取房型1-8、11的面積

alf=sum（d1（1：3））/sum（d1（4：9））；%普通宅與非普通宅的面積之比

f=[C（9）；C（10）]*[alf/（1+alf），1/（1+alf）]；%將房型9、10的開發(fā)成本按比例分?jǐn)倿槠胀ㄕc非普通宅

g=[R（9）；R（10）]*[alf/（1+alf），1/（1+alf）]；%將房型9、10的收入按比例分?jǐn)倿槠胀ㄕc非普通宅

f=round（f）；g=round（g）；

%計(jì)算每種房型的五種扣除項(xiàng)目金額

CF1=[C（1：2），0，C（4：7），0，f（1，：），f（2，：），0]；%房型3、8、11的開發(fā)成本不允許扣除記為0，房型9、10按分?jǐn)傆?jì)

CF=[C（1：8），f（1，：），f（2，：），C（11）]； %各種房型的開發(fā)成本，房型9、10按分?jǐn)傆?jì)

RF=[R（1：8），g（1，：），g（2，：），R（11）]；%各種房型的收入，房型9、10按分?jǐn)傆?jì)

dF=[d（1：8），d（9），0，d（10），0，d（11）]；

w=777179627；

s=w/sum（d）；%每單位建筑面積分?jǐn)偟耐恋刂Ц督痤~

S1=s*dF；%每種房型分?jǐn)偟耐恋刂Ц督痤~

S2=CF1；%房型3、8、11的開發(fā)成本不允許扣除

S3=（S1+S2）*0.1；%每種房型的開發(fā)費(fèi)用

S4=RF*0.0565；%每種房型與轉(zhuǎn)讓房地產(chǎn)有關(guān)的稅金

S5=（S1+S2）*0.2；%每種房型的其他扣除項(xiàng)目

S=S1+S2+S3+S4+S5；%每種房型總的扣除金額

bet=RF./S；%每種房型的收入與扣除項(xiàng)目金額之比

S=round（S）；

%以下計(jì)算增值稅

X=zeros（1，13）；

for i=1：13

if bet（i）>1.5&bet（i）<=2

X（i）=（0.4*bet（i）-0.45）*S（i）；

elseif bet（i）>2&bet（i）<=3

X（i）=（0.5*bet（i）-0.65）*S（i）；

elseif bet（i）>=3

X（i）=（0.6*bet（i）-0.95）*S（i）；

end

end

for i=1：3

if bet（i）<=1.2

X（i）=0； %普通宅當(dāng)bet《=1.2時(shí)免增值稅

elseif bet（i）>1.2&bet（i）<=1.5

X（i）=（0.3*bet（i）-0.3）*S（i）；

end

end

for i=4：13

if bet（i）<=1.5

X（i）=（0.3*bet（i）-0.3）*S（i）；

elseif bet（9）<=1.2&bet（11）<=1.2

X（9）=0；X（11）=0； %普通宅當(dāng)bet《=1.2時(shí)免增值稅

end

end

Y=round（X）；

Z=[Y（1：8），Y（9）+Y（10），Y（11）+Y（12），Y（13）]；

%計(jì)算投資回報(bào)率

CB=round（S1+CF+S3+S4）； %成本，其中w=sum（S1），sum（CF）=sum（C）

XC=[CB（1：8），CB（9）+CB（10），CB（11）+CB（12），CB（13）]； %計(jì)算房型1-11的各項(xiàng)總成本

L=R-XC-Z；

CC=sum（CB）；

SY=sum（R）； %收益，sum（RF）=sum（R）

ROI=（SY-CC）/CC； %投資回報(bào)率

2 眾籌筑屋方案Ⅰ的程序核算數(shù)據(jù)

在MATLAB命令窗口輸入方案Ⅰ的建設(shè)套數(shù)：

b=[250，250，150，250，250，250，250，75，150，150，75]；

[d，R，r，XC，Z，L，ROI]=hesuan（b）

回車，即可得到方案Ⅰ全面核算數(shù)據(jù)見表1所示。

3 改進(jìn)的方案Ⅱ的程序核算數(shù)據(jù)

通過對參籌者進(jìn)行抽樣調(diào)查，得到了參籌者對11種房型購買意愿的比例，為了盡量滿足參籌者的購買意愿，需重新設(shè)計(jì)建設(shè)規(guī)劃方案（稱為方案Ⅱ），并對方案II進(jìn)行核算。

在MATLAB命令窗口輸入方案II的建設(shè)套數(shù)：

b=[55，355，300，150，100，150，310，170，50，50，60]；

[d，R，r，XC，Z，L，ROI]=hesuan（b）

回車，得到方案II全面核算數(shù)據(jù)見表2所示。

4 優(yōu)化后的方案Ⅲ的程序核算數(shù)據(jù)

一般而言，投資回報(bào)率達(dá)到25%以上的眾籌項(xiàng)目才會(huì)被成功執(zhí)行，而方案Ⅱ的投資回報(bào)率ROI=24.61% ，需要對方案Ⅱ進(jìn)行優(yōu)化改進(jìn)，使其投資回報(bào)率達(dá)到25%以上。

在MATLAB命令窗口輸入優(yōu)化后的方案Ⅲ的建設(shè)套數(shù)：

b=[55，355，300，150，100，150，420，170，50，50，60]；

[d，R，r，XC，Z，L，ROI]=hesuan（b）

回車，得到方案Ⅲ的核算數(shù)據(jù)見表3所示。

由表1、表2可以看出，優(yōu)化后的方案Ⅲ比方案Ⅱ投資回報(bào)率提高7.24%，建房套數(shù)增加300套，純利潤增加136396967元，優(yōu)化效果顯著。方案Ⅲ既考慮投資者的購買意愿，又兼顧了開放商的利益，該方案可成功實(shí)施。

5 結(jié)語

通過編輯MATLAB核算函數(shù)程序，只要輸入規(guī)劃方案的建設(shè)套數(shù)，就能很方便地核算出重要的相關(guān)數(shù)據(jù)，減少了大量的重復(fù)計(jì)算。若核算數(shù)據(jù)未達(dá)到要求，可建立優(yōu)化模型進(jìn)行改進(jìn)使其達(dá)到預(yù)設(shè)目標(biāo)。該方法和程序運(yùn)用到其他行業(yè)的核算中非常簡捷方便。

參考文獻(xiàn)

[1]吉耀武.眾籌筑屋規(guī)劃設(shè)計(jì)的核算與優(yōu)化[J].中小企業(yè)管理與科技，2016.2.

[2]中華人民共和國土地增值稅暫行條例.中華人民共和國國務(wù)院[第538號(hào)令].

[3]胡良劍，孫曉君.MATLAB數(shù)學(xué)實(shí)驗(yàn)[M].北京：高等教育出版社，2006.06.

[4]肖華勇.實(shí)用數(shù)學(xué)建模與軟件應(yīng)用[M].西安：西北工業(yè)大學(xué)出版社，2008.11.