众筹筑屋规划方案设计

赛区评阅编号（由赛区组委会填写）：2015高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号（从A/B/C/D中选择一项填写）：D 我们的报名参赛队号（12位数字全国统一编号）：201523036005参赛学校（完整的学校全称，不含院系名）：宜宾学院参赛队员(打印并签名)：1.邓银梅2.刘明茂3.王富文指导教师或指导教师组负责人(打印并签名)：全靖 日期：2015年9月14日赛区评阅编号（由赛区组委会填写）：2015高教社杯全国大学生数学建模竞赛编号专用页赛区评阅记录（可供赛区评阅时使用）：评阅人备注送全国评奖统一编号（由赛区组委会填写）：全国评阅统一编号（由全国组委会填写）：II是可行的。5.3.3.6约束条件对投资回报率的影响（1）扣除比例对投资回报率的影响在计算房地产开发费用时，假设凡不能按转让房地产项目计算分摊利息支出或不能提供金融机构证明的，房地产开发费用按取得土地的使用权所支付的金额和房地产开发成本规定计算的金额之和的5%计算扣除；在实际生活中，房地产开发费用可能不扣除或按10%内计算扣除（没有固定扣除比例）。依据上述说明，应用MATLAB软件可画出投资回报率与扣除比例的关系图5-3-3-2（详见附录中程序5）：图5-3-3-2（“红线”表示0-10%所对应的投资回报率，“黑点”表示方案Ⅱ的投资回报率）由图5-3-3-2可知，无论开发费用不扣除或按10%内计算扣除，其投资回报率均达到25%以上；同时得到：扣除比例越大，其对应的投资回报率越大。开发成本对投资回报率的影响在建立非线性规划模型时，将方案Ⅱ的开发成本控制在方案Ⅰ的开发成本内，没有设置其它开发成本发生的情况，从而得出的投资回报率达到25%以上。若将未设置的情况考虑其中，会发生以下变化：根据逐步逼近法（详见附录中程序6）得出，开发成本与投资回报率的关系图（如图5-3-3-3）（详见附录中程序7）：图5-3-3-3（“·”表示投资回报率为25%所对应的点，“·”表示方案Ⅱ的投资回报率和对应的开发成本）由图5-3-3-3可知：当开发成本低于投资回报率为25%所对应的开发成本时，其对应方案的投资回报率达不到25%，因此，该方案不被成功执行；反之，则被成功执行。对于方案Ⅱ,它的投资成本高于投资回报率为25%所对应的开发成本，因此，该方案被成功执行。容积率对投资回报率的影响容积率是总建筑面积与总用地面积的比率，在总用地面不变的情况下，容积率与总建筑面积成正比，所以开发商希望越大越好，这样出售的面积大，赚的钱就多；住户希望越小越好，这样小区环境就好，绿化公共设施相对就多。方案Ⅱ是依据参筹户的购买意愿得出的容积率为1.99，忽略了开发商的收益，在成熟的企业会在尽量满足群众购买意愿的基础上考虑其收益，因此会适量地增大容积率。运用逐步逼近法（详见附录中程序8），得出容积率与投资回报率的关系图（如图5-3-3-4）（详见附录中程序9）：图5-3-3-4（“·”表示方案Ⅱ的容积率所对应的投资回报率）由图5-3-3-4可知，在约束范围内内的条件下容积率越大，其对应的投资回报率越高。六．模型的评价与推广6.1模型的评价众筹筑屋使购房人有机会在土地阶段就参与房产项目众筹，获得定制化服务与更大的购房优惠空间。同时，房产开发商可以通过众筹大幅降低融资成本，并提前锁定了购房客户，降低项目的销售风险与销售成本。本文建立的非线性规划模型，以满足群众的购买意愿为前提，得到新的建设规划方案。但在确立约束条件时，没有深入理解，得到的方案不佳，考虑的要素也不全面。建立的相关分析模型，涉及因素太少，具有局限性，各种因素之间分析的不够全面，若结合层次分析法对所涉及的因素进行权重分析，得到的效果将会更好。6.2模型的推广本文建立的模型，既满足了参筹者的购买意愿，又考虑了开发商的收益，并适用于为人民服务的公共设施的建造设计参数。参考文献[1]韩卫涛、赵硕，土地增值税的计算方法，乡镇企业科技，05期：31-32，1997年；[2]司守奎、孙玺菁，数学建模算法与应用，北京：国防工业出版社，2012年；[3]张友生、王勇，系统分析师考试全程指导，北京：清华大学出版社，2009年；[4]韩中庚，数学建模算法与应用（第2版），北京：高等教育出版社，2005年附录程序1（%运用MATLAB7.1软件求解普通宅和非普通宅的增值税）clearall;clc;l=load('m.txt');m1=0;cc1=0;b1=0;m2=0;cc2=0;b2=0;q=777179627;%土地使用权所支付的金额cc=l(3,1)*l(3,2)*l(3,3)+l(8,1)*l(8,2)*l(8,3)+l(11,1)*l(11,2)*l(11,3);fori=1:3m12=l(i,1)*l(i,2);b12=l(i,1)*l(i,2)*l(i,4);m1=m1+m12;b1=b1+b12;endm1;%普通宅的建筑面积b1;%普通宅的转让房地产收入fori=[4:8,11]m21=l(i,1)*l(i,2);b21=l(i,1)*l(i,2)*l(i,4);m2=m2+m21;b2=b2+b21;endm2;%非普通宅的建筑面积b2;%非普通宅的转让房地产收入m9=l(9,1)*l(9,2);m10=l(10,1)*l(10,2);%房型9、房型10的建筑面积m91=(m1/(m1+m2))*m9;m92=m2/(m1+m2)*m9;%房型9普通宅和非普通宅的建筑面积m101=(m1/(m1+m2))*m10;m102=m2/(m1+m2)*m10;%房型10普通宅和非普通宅的建筑面积q1=(m1/(m1+m2))*q;q2=(m2/(m1+m2))*q;%把土地使用权所支付的金额按普通宅和非普通宅建筑面积比分摊b91=m91*l(9,4);b92=m92*l(9,4);%房型9普通宅和非普通宅的转让房地产收入b101=m101*l(10,4);b102=m102*l(10,4);%房型10普通宅和非普通宅的转让房地产收入cc91=m91*l(9,3);cc92=m92*l(9,3);%房型9普通宅和非普通宅的开发成本cc101=m101*l(10,3);cc102=m102*l(10,3);%房型10普通宅和非普通宅的开发成本fori=[12]cc12=l(i,1)*l(i,2)*l(i,3);cc1=cc1+cc12;endcc1;cc1=cc1+cc91+cc101;%普通宅的开发成本fori=4:7cc21=l(i,1)*l(i,2)*l(i,3);cc2=cc2+cc21;endcc2;cc2=cc2+cc92+cc102;%非普通宅的开发成本b1=b1+b91+b101;b2=b2+b92+b102;%普通宅和非普通宅的转让房地产总收入kj1=b1*0.0565+(cc1+q1)*(1+0.05+0.20);%求解普通宅扣除项目金额ze1=b1-kj1;%求解普通宅增值额g1=ze1/kj1kj2=b2*0.0565+(cc2+q2)*(1+0.05+0.20);%求解非普通宅扣除项目金额ze2=b2-kj2;%求解非普通宅增值额g2=ze2/kj2;%%%核算模式1求解增值税ifg1<=0.5zs1=ze1*0.3;elseifg1>0.5&&g1<=1zs1=kj1*0.5*0.3+(ze1-kj1*0.5)*0.4;elseifg1>1&&g1<=2zs1=kj1*0.5*0.3+(kj1*1-kj1*0.5)*0.4+(ze1-kj1*1)*0.5;elseifg1>2zs1=kj1*0.5*0.3+(kj1*1-kj1*0.5)*0.4+(ze1-kj1*1)*0.5+(ze1-kj1*2)*0.6;endzs1ifg2<=0.5zs2=ze2*0.3;elseifg2>0.5&&g2<=1zs2=kj2*0.5*0.3+(ze2-kj2*0.5)*0.4;elseifg2>1&&g2<=2zs2=kj2*0.5*0.3+(kj2*1-kj2*0.5)*0.4+(ze2-kj2*1)*0.5;elseifg2>2zs2=kj2*0.5*0.3+(kj2*1-kj2*0.5)*0.4+(ze2-kj2*1)*0.5+(ze2-kj2*2)*0.6;endzs2zs=zs1+zs2%增值税b=b1+b2;%转让房地产总收入c=cc+cc1+cc2;%开发成本C=c+q+b*0.0565+zs%开发总成本sy=b-C%求解其收益hui=sy/C%投资回报率%%%核算模式2求解普通宅增值税ifg1<=0.5zs1=ze1*0.3;elseifg1>0.5&&g1<=1zs1=ze1*0.4-kj1*0.05;elseifg1>1&&g1<=2zs1=ze1*0.5-kj1*0.15;elseifg1>2zs1=ze1*0.6-kj1*0.35;endzs1;ifg2<=0.5zs2=ze2*0.3;elseifg2>0.5&&g2<=1zs2=ze2*0.4-kj2*0.05;elseifg2>1&&g2<=2zs2=ze2*0.5-kj2*0.15;elseifg2>2zs2=ze2*0.6-kj2*0.35;endzs2;zs=zs1+zs2%增值税b=b1+b2;%转让房地产总收入c=cc+cc1+cc2;%开发成本C=c+q+b*0.0565+zs%开发总成本sy=b-C%求解其收益hui=sy/C%投资回报率程序2（%运用LINGO14.0软件求解最大满意率情况下的建房套数）max=(0.4*x1+0.6*x2+0.5*x3+0.6*x4+0.7*x5+0.8*x6+0.9*x7+0.6*x8+0.6*x9+0.3*x10+0.4*x11)/(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11);(x1*77+x2*98+x3*117+x4*145+x5*156+x6*167+x7*178+x8*126)/102077<2.28;x1*77*0.4263+x2*98*0.4323+x3*117*0.4532+x4*145*0.5288+x5*156*0.5268+x6*167*0.5533+x7*178*0.5685+x8*126*0.4323+x9*103*0.2663+x10*129*0.2791+x11*133*0.2982<131438.71;@bnd(50,x1,450);@bnd(50,x2,500);@bnd(50,x3,300);@bnd(150,x4,500);@bnd(100,x5,550);@bnd(150,x6,350);@bnd(50,x7,450);@bnd(100,x8,250);@bnd(50,x9,350);@bnd(50,x10,400);@bnd(50,x11,250);@gin(x1);@gin(x2);@gin(x3);@gin(x4);@gin(x5);@gin(x6);@gin(x7);@gin(x8);@gin(x9);@gin(x10);@gin(x11);程序3（%运用MATLAB7.1软件求解方案Ⅱ的投资回报率）clearall;clc;l=load('m1.txt');m1=0;cc1=0;b1=0;m2=0;cc2=0;b2=0;q=777179627;%土地使用权所支付的金额cc=l(3,1)*l(3,2)*l(3,3)+l(8,1)*l(8,2)*l(8,3)+l(11,1)*l(11,2)*l(11,3);fori=1:3m12=l(i,1)*l(i,2);b12=l(i,1)*l(i,2)*l(i,4);m1=m1+m12;b1=b1+b12;endm1;%普通宅的建筑面积b1;%普通宅的转让房地产收入fori=[4:8,11]m21=l(i,1)*l(i,2);b21=l(i,1)*l(i,2)*l(i,4);m2=m2+m21;b2=b2+b21;endm2;%非普通宅的建筑面积b2;%非普通宅的转让房地产收入m9=l(9,1)*l(9,2);m10=l(10,1)*l(10,2);%房型9、房型10的建筑面积m91=(m1/(m1+m2))*m9;m92=m2/(m1+m2)*m9;%房型9普通宅和非普通宅的建筑面积m101=(m1/(m1+m2))*m10;m102=m2/(m1+m2)*m10;%房型10普通宅和非普通宅的建筑面积q1=(m1/(m1+m2))*q;q2=(m2/(m1+m2))*q;%把土地使用权所支付的金额按普通宅和非普通宅建筑面积比分摊b91=m91*l(9,4);b92=m92*l(9,4);%房型9普通宅和非普通宅的转让房地产收入b101=m101*l(10,4);b102=m102*l(10,4);%房型10普通宅和非普通宅的转让房地产收入cc91=m91*l(9,3);cc92=m92*l(9,3);%房型9普通宅和非普通宅的开发成本cc101=m101*l(10,3);cc102=m102*l(10,3);%房型10普通宅和非普通宅的开发成本fori=[12]cc12=l(i,1)*l(i,2)*l(i,3);cc1=cc1+cc12;endcc1;cc1=cc1+cc91+cc101;%普通宅的开发成本fori=4:7cc21=l(i,1)*l(i,2)*l(i,3);cc2=cc2+cc21;endcc2;cc2=cc2+cc92+cc102;%非普通宅的开发成本b1=b1+b91+b101;b2=b2+b92+b102;%普通宅和非普通宅的转让房地产总收入kj1=b1*0.0565+(cc1+q1)*(1+0.05+0.20);%求解普通宅扣除项目金额ze1=b1-kj1;%求解普通宅增值额g1=ze1/kj1kj2=b2*0.0565+(cc2+q2)*(1+0.05+0.20);%求解非普通宅扣除项目金额ze2=b2-kj2;%求解非普通宅增值额g2=ze2/kj2;%%%核算模式1求解增值税ifg1<=0.5zs1=ze1*0.3;elseifg1>0.5&&g1<=1zs1=kj1*0.5*0.3+(ze1-kj1*0.5)*0.4;elseifg1>1&&g1<=2zs1=kj1*0.5*0.3+(kj1*1-kj1*0.5)*0.4+(ze1-kj1*1)*0.5;elseifg1>2zs1=kj1*0.5*0.3+(kj1*1-kj1*0.5)*0.4+(ze1-kj1*1)*0.5+(ze1-kj1*2)*0.6;endzs1ifg2<=0.5zs2=ze2*0.3;elseifg2>0.5&&g2<=1zs2=kj2*0.5*0.3+(ze2-kj2*0.5)*0.4;elseifg2>1&&g2<=2zs2=kj2*0.5*0.3+(kj2*1-kj2*0.5)*0.4+(ze2-kj2*1)*0.5;elseifg2>2zs2=kj2*0.5*0.3+(kj2*1-kj2*0.5)*0.4+(ze2-kj2*1)*0.5+(ze2-kj2*2)*0.6;endzs2zs=zs1+zs2%增值税b=b1+b2;%转让房地产总收入c=cc+cc1+cc2;%开发成本C=c+q+b*0.0565+zs%开发总成本sy=b-C%求解其收益hui=sy/C%投资回报率程序4（%运用MATLAB7.1软件计算相关系数）x=[1429476711450118000666619497121872162738250013419200003.723.723.723.72826.228138.8102569440];r=[];fori=1:4r1=(11*x(5,i)-x(1,i)*x(2,i))/(sqrt(11*x(3,i)-x(1,i)*x(1,i))*sqrt(11*x(4,i)-x(2,i)*x(2,i)));r=[rr1];endr程序5（%运用MATLAB7.1软件绘制投资回报率与扣除比例关系图）clearall;clc;l=load('m1.txt');%导入方案2的相关数据i=length(l);b=0;sy=0;Hui=[];c=1146091100;%开发成本q=777179627;%土地支付的金额fork=1:ib1=l(k,1)*l(k,2)*l(k,4);b=b+b1;endb;%转让房地产总收入cc=l(3,1)*l(3,2)*l(3,3)+l(8,1)*l(8,2)*l(8,3)+l(11,1)*l(11,2)*l(11,3);ct=c-cc;forj=0:0.001:0.1kj=(ct+q)*(1+j+0.2)+b*0.0565;%求解扣除项目金额ze=b-kj;%求解增值额g=ze/kj;%%%核算模式2求解增值税ifg<=0.5zs=ze*0.3;elseifg>0.5&&g<=1zs=ze*0.4-kj*0.05;elseifg>1&&g<=2zs=ze*0.5-kj*0.15;elseifg>2zs=ze*0.6-kj*0.35;endC=c+q+b*0.0565+zs;%开发总成本sy=b-C;%求解其收益hui=sy/C;Hui=[Huihui];endHuij=0:0.001:0.1;plot(j,Hui,'r');xlabel('扣除比例');ylabel('投资回报率');holdonplot(0.05,0.282,'.','MarkerEdgeColor','k','MarkerSize',20);程序6（%运用LINGO14.0软件得出开发成本与投资回报率的结果）max=(0.4*x1+0.6*x2+0.5*x3+0.6*x4+0.7*x5+0.8*x6+0.9*x7+0.6*x8+0.6*x9+0.3*x10+0.4*x11)/(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11);(x1*77+x2*98+x3*117+x4*145+x5*156+x6*167+x7*178+x8*126)/102077<2.28;x1*77*0.4263+x2*98*0.4323+x3*117*0.4532+x4*145*0.5288+x5*156*0.5268+x6*167*0.5533+x7*178*0.5685+x8*126*0.4323+x9*103*0.2663+x10*129*0.2791+x11*133*0.2982>90000;！90000表示从90000开始，以5000为步长，到150000结束的所有数x1*77*0.4263+x2*98*0.4323+x3*117*0.4532+x4*145*0.5288+x5*156*0.5268+x6*167*0.5533+x7*178*0.5685+x8*126*0.4323+x9*103*0.2663+x10*129*0.2791+x11*133*0.2982<95000;！95000表示从90000开始，以5000为步长，到155000结束的所有数c=x1*77*0.4263+x2*98*0.4323+x3*117*0.4532+x4*145*0.5288+x5*156*0.5268+x6*167*0.5533+x7*178*0.5685+x8*126*0.4323+x9*103*0.2663+x10*129*0.2791+x11*133*0.2982;b=x1*77*1.2000+x2*98*1.0800+x3*117*1.1200+x4*145*1.2800+x5*156*1.2800+x6*167*1.3600+x7*178*1.4000+x8*126*1.0400+x9*103*0.6400+x10*129*0.6800+x11*133*0.7200;cc=x3*117*0.4532+x8*126*0.4323+x11*133*0.2982;ct=c-cc;kj=b*0.0565+(ct+77717.9627)*(1+0.05+0.2);ze=b-kj;g=ze/kj;zs=@if(g#le#0.5,ze*0.3,@if(g#ge#0.5#and#g#le#1,ze*0.4-kj*0.05,@if(g#ge#1#and#g#le#2,ze*0.5-kj*0.15,@if(g#ge#2,ze*0.6-kj*0.35,zs))));sy=b-(c+77717.9627+b*0.0565+zs);hui=sy/(c+77717.9627+b*0.0565+zs);@bnd(50,x1,450);@bnd(50,x2,500);@bnd(50,x3,300);@bnd(150,x4,500);@bnd(100,x5,550);@bnd(150,x6,350);@bnd(50,x7,450);@bnd(100,x8,250);@bnd(50,x9,350);@bnd(50,x10,400);@bnd(50,x11,250);@gin(x1);@gin(x2);@gin(x3);@gin(x4);@gin(x5);@gin(x6);@gin(x7);@gin(x8);@gin(x9);@gin(x10);@gin(x11);程序7（%运用MATLAB7.1软件绘制开发成本与投资回报率的关系图）x=[99.51010.51111.51212.51313.51414.51515.5];y=[0.20311610.22142700.23849310.25451150.26926730.28307060.29561890.30738990.31845430.32846780.33493490.33465540.34123680.3499369];plot(x,y,'r')xlabel('开发成本（亿元）');ylabel('投资回报率');gridonholdonplot(11.47,0.282,'.','MarkerEdgeColor','k','MarkerSize',20);plot(10.35,0.25,'.','MarkerEdgeColor','b','MarkerSize',20);程序8（%运用LINGO14.0软件求解得出容积率与投资回报率)max=(0.4*x1+0.6*x2+0.5*x3+0.6*x4+0.7*x5+0.8*x6+0.9*x7+0.6*x8+0.6*x9+0.3*x10+0.4*x11)/(x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11);(x1*77+x2*98+x3*117+x4*145+x5*156+x6*167+x7*178+x8*126)/102077<2.28;(x1*77+x2*98+x3*117+x4*145+x5*156+x6*167+x7*178+x8*126)/102077>1.98;！1.98表示从1.98开始，以0.02为步长，到2.26结束的所有数c=x1*77*0.4263+x2*98*0.4323+x3*117*0.4532+x4*145*0.5288+x5*156*0.5268+x6*167*0.5533+x7*178*0.5685+x8*126*0.4323+x9*103*0.2663+x10*129*0.2791+x11*133*0.2982;b=x1*77*1.2000+x2*98*1.0800+x3*117*1.1200+x4*145*1.2