全国大学生建模D题答案

2015-2016学年第一学期《数学建模算法与实践》课程报告姓名姚志力学号20140488专业电气工程及其自动化成绩众筹筑屋规划方案设计众筹筑屋是互联网时代一种新型的房地产形式。现有占地面积为102077.6平方米的众筹筑屋项目。项目推出后，有上万户购房者登记参筹。为了信息公开及民主决策，需要将这个众筹筑屋项目原方案（称作方案Ⅰ）的成本与收益、容积率和增值税等信息进行公布，根据国家相关政策，不同房型的容积率、开发成本、开发费用等在核算上的要求，以及相关条例与政策，对众筹筑屋项目原方案（称作方案Ⅰ）的成本与收益、容积率和增值税等进行核算。我们需考虑诸多因素，如容积率、开发成本、税率、预期收益等，来判断该方案计划是否正确并结合参筹者满意度对方案进行改进来一、模型建立对原方案（称作方案Ⅰ）的成本与收益、容积率和增值税四项进行核算。首先核算容积率，我们设容积率为r,mi表示房型i的房型面积，ti表示房型i的建房套数，我们知道方案Ⅰ的总面积为102077.6平方米,由附件1-1可知，建筑面积参与容积率计算的为前则容积率为r为：r=由附件1-1知，国家规定的最大容积率要求为2.28，r应该小于2.28才是符合国家规定的。对成本进行的核算，我们设成本为C，本项目的成本包括两部分，分别为取得土地使用权所支付的金额，由附件1-2可知，为777179627；另一部分为房地产开发成本，我们设开发成本为c,mi表示房型i的房型面积，ti表示房型i的建房套数,ci表示房型则开发成本c为：c=则成本C为：C=对收益进行的核算，本项目的收入为卖房所得，我们设收入为S，mi表示房型i的房型面积，ti表示房型i的建房套数,si则收入S为：S=除了以上数据外，增值额也需要根据国家政策重新分类核算，我国增值额为以扣除项目金额Q的累进计算，我们设本项目扣除项目金额Q为，扣除项目金额Q主要包括两个部分，一个为取得土地使用权所支付的金额，由附件1-2可知，为777179627；另一个是开发土地的成本，由附加1-1知，有3种房型在计算增值额时不允许扣除成本，属于“其他”类型的住宅的开发成本需要按照已有普桐宅和非普通宅建筑面积比分摊计算。针对普通宅、非普通宅，我们设普通宅、非普通宅的开发成本C1,由附件1-1内要求C针对其他类住宅，我们设其他类住宅的开发成本C2,有要求我们需要先对其他类住宅建筑面积进行分摊普通宅的建筑面积M1为M非普通宅的建筑面积M2为M由分摊可计算其他类住宅的建筑面积M3为房型9普通宅分摊：M房型10普通宅分摊：M由此，我们核算其他类住宅的开发成本C2,为C那么，扣除项目金额Q为：Q=增值额Z为：Z=S按照附件2.表1.2的累进税率表计算增值税我们设增值税U，按照累进税率表我们把增值额分类进行不同公式的核算，分类如下：增值额未超过扣除项目金额50%部分即ZU=Z增值额超过扣除项目金额50%部分，未超过100%部分即0.5×Q<U=Z增值额超过扣除项目金额100%部分，未超过200%部分即1×Q<U=Z增值额超过扣除项目金额20%部分即ZU=Z又由附件2.1.3土地增值税的税收优惠可知，判断本方案是否有优惠情况：建筑普通标准住宅的税收优惠。先判断本情况是否符合优惠条件：纳税人建造普通住宅出售，增值额未超过扣除项目金额20%，免征土地增值税。所以，我们需要对普通住宅出售增值额Z普进行核算，并将普通住宅出售增值额Z普与扣除项目金额普通住宅出售增值额Z普Z由上面收益核算可得S普为S由上面开放成本核算可得Q普为Q如果Z则Z优惠为Z增值税U优惠计算方法同U投资回报率f为：f=现需结合参筹者满意度和其他约束设计最佳方案，我们分析有4个限制，分别为套数限制、满意度限制、容积率限制、投资回报率限制。先找出符合套数限制和容积率限制的所有情况：即t得出不同的结果，并分别计算不同结果的投资回报率即：f=选取投资回报率最高的前150种方案进行满意度分析。对每一个房型生成1024个0或1的随机数组，得到A=1024×11的矩阵，并对第i行进行求和作为第12列，则满意房型i的套数aka再判断mydk是否符合附件1-4myd对于满足的方案，投资回报率越大，容积率越小，为最优解。二、模型求解下面进行模型求解，C语言编程如下：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];/开发成本/income+=ft[i]*fm[i]*fs[i];/销售收入/ } cost+=777179625; for(i=0;i<7;i++)/扣除项目计算/ 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)/判断是否可以优惠/ { flag=1; z=zp; } if(z<=0.5*zcost)/增值税累进计算/ 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,"有增值优惠\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,"投资回报率：%lf\n",hbl); fclose(fp); return0;}运行结果为：有增值优惠总面积277025.000000容积率2.275230总成本2091566725.000000总收入3246720000.000000增值额1376357420.113015增值税457024839.050857净收益919332581.062158投资回报率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]++) /循环判断是否符合/ 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