全國(guó)大學(xué)生建模D題答案

2015-2016學(xué)年第一學(xué)期《數(shù)學(xué)建模算法與實(shí)踐》課程報(bào)告姓名姚志力學(xué)號(hào)20140488專業(yè)電氣工程及其自動(dòng)化成績(jī)眾籌筑屋規(guī)劃方案設(shè)計(jì)眾籌筑屋是互聯(lián)網(wǎng)時(shí)代一種新型的房地產(chǎn)形式?，F(xiàn)有占地面積為102077.6平方米的眾籌筑屋項(xiàng)目。項(xiàng)目推出后，有上萬(wàn)戶購(gòu)房者登記參籌。為了信息公開及民主決策，需要將這個(gè)眾籌筑屋項(xiàng)目原方案（稱作方案Ⅰ）的成本與收益、容積率和增值稅等信息進(jìn)行公布，根據(jù)國(guó)家相關(guān)政策，不同房型的容積率、開發(fā)成本、開發(fā)費(fèi)用等在核算上的要求，以及相關(guān)條例與政策，對(duì)眾籌筑屋項(xiàng)目原方案（稱作方案Ⅰ）的成本與收益、容積率和增值稅等進(jìn)行核算。我們需考慮諸多因素，如容積率、開發(fā)成本、稅率、預(yù)期收益等，來(lái)判斷該方案計(jì)劃是否正確并結(jié)合參籌者滿意度對(duì)方案進(jìn)行改進(jìn)來(lái)一、模型建立對(duì)原方案（稱作方案Ⅰ）的成本與收益、容積率和增值稅四項(xiàng)進(jìn)行核算。首先核算容積率，我們?cè)O(shè)容積率為r,mi表示房型i的房型面積，ti表示房型i的建房套數(shù)，我們知道方案Ⅰ的總面積為102077.6平方米,由附件1-1可知，建筑面積參與容積率計(jì)算的為前則容積率為r為：r=由附件1-1知，國(guó)家規(guī)定的最大容積率要求為2.28，r應(yīng)該小于2.28才是符合國(guó)家規(guī)定的。對(duì)成本進(jìn)行的核算，我們?cè)O(shè)成本為C，本項(xiàng)目的成本包括兩部分，分別為取得土地使用權(quán)所支付的金額，由附件1-2可知，為777179627；另一部分為房地產(chǎn)開發(fā)成本，我們?cè)O(shè)開發(fā)成本為c,mi表示房型i的房型面積，ti表示房型i的建房套數(shù),ci表示房型則開發(fā)成本c為：c=則成本C為：C=對(duì)收益進(jìn)行的核算，本項(xiàng)目的收入為賣房所得，我們?cè)O(shè)收入為S，mi表示房型i的房型面積，ti表示房型i的建房套數(shù),si則收入S為：S=除了以上數(shù)據(jù)外，增值額也需要根據(jù)國(guó)家政策重新分類核算，我國(guó)增值額為以扣除項(xiàng)目金額Q的累進(jìn)計(jì)算，我們?cè)O(shè)本項(xiàng)目扣除項(xiàng)目金額Q為，扣除項(xiàng)目金額Q主要包括兩個(gè)部分，一個(gè)為取得土地使用權(quán)所支付的金額，由附件1-2可知，為777179627；另一個(gè)是開發(fā)土地的成本，由附加1-1知，有3種房型在計(jì)算增值額時(shí)不允許扣除成本，屬于“其他”類型的住宅的開發(fā)成本需要按照已有普桐宅和非普通宅建筑面積比分?jǐn)傆?jì)算。針對(duì)普通宅、非普通宅，我們?cè)O(shè)普通宅、非普通宅的開發(fā)成本C1,由附件1-1內(nèi)要求C針對(duì)其他類住宅，我們?cè)O(shè)其他類住宅的開發(fā)成本C2,有要求我們需要先對(duì)其他類住宅建筑面積進(jìn)行分?jǐn)偲胀ㄕ慕ㄖ娣eM1為M非普通宅的建筑面積M2為M由分?jǐn)偪捎?jì)算其他類住宅的建筑面積M3為房型9普通宅分?jǐn)偅篗房型10普通宅分?jǐn)偅篗由此，我們核算其他類住宅的開發(fā)成本C2,為C那么，扣除項(xiàng)目金額Q為：Q=增值額Z為：Z=S按照附件2.表1.2的累進(jìn)稅率表計(jì)算增值稅我們?cè)O(shè)增值稅U，按照累進(jìn)稅率表我們把增值額分類進(jìn)行不同公式的核算，分類如下：增值額未超過(guò)扣除項(xiàng)目金額50%部分即ZU=Z增值額超過(guò)扣除項(xiàng)目金額50%部分，未超過(guò)100%部分即0.5×Q<U=Z增值額超過(guò)扣除項(xiàng)目金額100%部分，未超過(guò)200%部分即1×Q<U=Z增值額超過(guò)扣除項(xiàng)目金額20%部分即ZU=Z又由附件2.1.3土地增值稅的稅收優(yōu)惠可知，判斷本方案是否有優(yōu)惠情況：建筑普通標(biāo)準(zhǔn)住宅的稅收優(yōu)惠。先判斷本情況是否符合優(yōu)惠條件：納稅人建造普通住宅出售，增值額未超過(guò)扣除項(xiàng)目金額20%，免征土地增值稅。所以，我們需要對(duì)普通住宅出售增值額Z普進(jìn)行核算，并將普通住宅出售增值額Z普與扣除項(xiàng)目金額普通住宅出售增值額Z普Z(yǔ)由上面收益核算可得S普為S由上面開放成本核算可得Q普為Q如果Z則Z優(yōu)惠為Z增值稅U優(yōu)惠計(jì)算方法同U投資回報(bào)率f為：f=現(xiàn)需結(jié)合參籌者滿意度和其他約束設(shè)計(jì)最佳方案，我們分析有4個(gè)限制，分別為套數(shù)限制、滿意度限制、容積率限制、投資回報(bào)率限制。先找出符合套數(shù)限制和容積率限制的所有情況：即t得出不同的結(jié)果，并分別計(jì)算不同結(jié)果的投資回報(bào)率即：f=選取投資回報(bào)率最高的前150種方案進(jìn)行滿意度分析。對(duì)每一個(gè)房型生成1024個(gè)0或1的隨機(jī)數(shù)組，得到A=1024×11的矩陣，并對(duì)第i行進(jìn)行求和作為第12列，則滿意房型i的套數(shù)aka再判斷mydk是否符合附件1-4myd對(duì)于滿足的方案，投資回報(bào)率越大，容積率越小，為最優(yōu)解。二、模型求解下面進(jìn)行模型求解，C語(yǔ)言編程如下：A．#include<stdio.h>intmain(){ FILE*fp; inti,flag=0,ft[11]={250,250,150,250,250,250,250,75,150,150,75}; intfm[11]={77,98,117,145,156,167,178,126,103,129,133};intfc[11]={4263,4323,4532,5288,5268,5533,5685,4323,2663,2791,2982}; intfs[11]={12000,10800,11200,12800,12800,13600,14000,10400,6400,6800,7200}; doublerjl=0,u=0,hbl=0,z=0,m=0,m8=0,cost=0,zcost=0,zcost1=0,zcost2=0,income=0; doubleM1=0,M2=0,M9=0,M10=0,sp=0,cp=0,zp; for(i=0;i<8;i++) m8+=ft[i]*fm[i]; rjl=m8/102077.6; for(i=0;i<11;i++) { m+=ft[i]*fm[i];/總面積/ cost+=ft[i]*fm[i]*fc[i];/開發(fā)成本/income+=ft[i]*fm[i]*fs[i];/銷售收入/ } cost+=777179625; for(i=0;i<7;i++)/扣除項(xiàng)目計(jì)算/ zcost1+=ft[i]*fm[i]*fc[i]; zcost1-=ft[2]*fm[2]*fc[2]; for(i=0;i<3;i++) M1+=ft[i]*fm[i]; for(i=3;i<8;i++) M2+=ft[i]*fm[i]; M2+=ft[10]*fm[10]; M9=ft[8]*fm[8]*M1/(M1+M2); M10=ft[9]*fm[9]*M1/(M1+M2);zcost2=M9*fc[8]+M10*fc[9]; zcost=zcost1+zcost2+777179625; z=income-zcost; for(i=0;i<3;i++) { sp+=ft[i]*fm[i]*fs[i]; cp=ft[i]*fm[i]*fc[i]; } sp+=M9*fs[8]+M10*fs[9]; cp+=M9*fc[8]+M10*fc[9]+(M1+M9+M10)/m*777179625;zp=sp-cp; if(zp<=0.2*cp)/判斷是否可以優(yōu)惠/ { flag=1; z=zp; } if(z<=0.5*zcost)/增值稅累進(jìn)計(jì)算/ u=0.3*z; elseif(z<=zcost) u=0.4*z-0.05*zcost; elseif(z<=2*zcost) u=0.5*z-0.15*zcost; else u=0.6*z-0.35*zcost; hbl=(z-u)/(zcost+u); fp=fopen("C:\\Users\\lenovof\\Desktop\\建模比賽題目\\ft.txt","w"); if(flag=1) fprintf(fp,"有增值優(yōu)惠\n",flag); fprintf(fp,"總面積：%lf\n",m); fprintf(fp,"容積率：%lf\n",rjl); fprintf(fp,"總成本：%lf\n",cost); fprintf(fp,"總收入：%lf\n",income); fprintf(fp,"增值額：%lf\n",z); fprintf(fp,"增值稅：%lf\n",u); fprintf(fp,"凈收益：%lf\n",z-u); fprintf(fp,"投資回報(bào)率：%lf\n",hbl); fclose(fp); return0;}運(yùn)行結(jié)果為：有增值優(yōu)惠總面積277025.000000容積率2.275230總成本2091566725.000000總收入3246720000.000000增值額1376357420.113015增值稅457024839.050857凈收益919332581.062158投資回報(bào)率0.395006B．#include<stdio.h>intmain(){ FILE*fp; inti=0,j,k,s,flag=0,newft[10000][8],ft[8]={50,50,50,150,100,150,50,100},min[11]={50,50,50,150,100,150,50,100},max[8]={450,500,300,500,550},fm[8]={77,98,117,145,156,167,178,126}; doubleS,rjl,bz; bz=2.280000; fp=fopen("C:\\Users\\lenovof\\Desktop\\建模比賽題目\\ft.txt","w"); for(ft[0]=min[0];ft[0]<=max[0];ft[0]++) /循環(huán)判斷是否符合/ for(ft[1]=min[1];ft[1]<=max[1];ft[1]++) for(ft[2]=min[2];ft[2]<=max[2];ft[2]++) for(ft[3]=min[3];ft[3]<=max[3];ft[3]++) for(ft[4]=min[4];ft[4]<=max[4];ft[4]++) for(ft[5]=min[5];ft[5]<=max[5];ft[5]++) for(ft[6]=min[6];ft[6]<=max[6];ft[6]++) for(ft[7]=min[7];ft[7]<=max[7];ft[7]++) { for(j=0,S=0;j<8;j++) S=S+ft[j]*fm[j]; rjl=S/102077.6; if(rjl<=bz)/*problem*/ for(k=0,s=0;k<8;k++) { newft[flag][k]=ft[k]; printf("%d",newft[flag][k]); fprintf(fp,"%d",newft[flag][k]); s=s+newft[flag][k]; } fprintf(fp,"%lf",rjl); fprintf(fp,"%d\n",s); flag++; p