2016五一联赛A题要点

1、2016年第十三届五一数学建模联赛编号专用页竞赛评阅编号（由竞赛评委会评阅前进行编号）：评阅记录评阅人评分备注裁剪线 裁剪线 裁剪线竞赛评阅编号（由竞赛评委会评阅前进行编号）：参赛队伍的参赛号码：（请各参赛队提前填写好）：2016年第十三届五一数学建模联赛A题：购房中的数学问题影响消费者选购住房的有平均日照时间、价格、交通、环境和噪音等多方面因素， 其中最重要的自然因素莫过于日均采光时间，我们通过建立影子随时间变化的物理模 型，给出了一种基于离散思想的日均采光时间的计算方法，并以东经117.17°,北纬34.18°处一高层建筑小区的14-2-802房间（客厅）为例，分别求得

2、其冬至日9:00-16:00间可以享受日照的时间区间为：9:00-10 : 20和12: 40-1:32 和15:36-16 : 00,全年365天每一天可以享受日照的累计时间为：208110分钟，全年享 受日照时间超过6小时的天数和日期为：2月25号到10月16共计235天，并在仅考虑 采光影响的条件下给出最优选房方案为：基于以上采光时间的计算结果并结合其他多方面因素，我们建立了个性化选房模型。针对问题一：由于题设前提不考虑天气等影响日照的因素，因此临近的高层建筑的 遮挡是唯一影响采光时间的因素，因此我们建立了障碍物影长随时间变化的物理模型， 为避免公式推导太阳及影子变化轨迹引起计算繁琐的情

3、况，我们考虑将太阳位置离散化处理，在保证精度的前提下大幅度简化了运算，据此计算出一日内任意时刻的太阳高度 角及太阳方位角，继而根据障碍物高度，楼宇间距等数据计算出其冬至日9:00-16:00间可以享受日照的时间区间为：9:00-10 : 20和12: 40-1:32和15:36-16 : 00针对问题二：对同一房间，一年内可能遮挡它的建筑物是固定的，而障碍物的影长 随太阳高度角，赤纬角的变化而变化；因此我们根据赤纬角随日期的变化公式，计算了 一年内每一天的太阳高度角和赤纬角，带入问题一的模型并求和，即可描述出全年365天每一天可以享受日照的累计时间为：2月25到10月16共计235天针对问题三

4、：跟据问题仅考虑采光影响前提，我们通过提取小区住宅楼的布局关系， 建立了各楼层互相遮挡的关系矩阵，代入问题二的模型得到小区内18栋楼每一层的年日照累计和，结合附件2的已售楼房信息，我们可以确定一批最优选房方案。针对问题四：我们根据小区平面图对小区各个影响买房的方面做了评估，利用反复取交集的方法为客户选择房子。最终为高层偏好的人选择了9号楼33层，中层偏好的人选择9号楼13层，低层偏好的人选择了 5号楼1层。针对问题五：通过以两个电梯出口为原点建立空间坐标系，计算车位到电梯口的距 离和此车位住户到电梯口的距离，发现并不符合最长和最短距离的结合原则，则车位安 排不合理。对于安排车位，我们将充分考虑

5、各种因素，包括低层住户长时间等待时是否 乘坐电梯，多种电梯的运行模式的选择，出电梯口时距离电梯口的最近车位，车位与住 户楼层的最短路距离，得出最优车位。关键字：影子变化物理模型离散化处理个性化选房模型电梯控制策略最短路问题重述当我们买房时，房子的选择是一件很头疼的一件事，地理位置、周边环境、交通便 利性、住房户型、住房价格、采光、噪音污染、空气污染等因素影响我们的选择，同时 开发商总会追求利润的最大化，所以现在高楼林立，面对现在的高层建筑，选房的问题 变得更加复杂，针对东经117.17°,北纬34.18°地理位置的高层建筑，不考虑天气的影 响，正确解决购房问题。问题一，解决

6、A小区14-2-802房间（客厅）在冬至时特定时间段可以享受到日照的时 间段结合附件1和4得出正确结果。问题二，描绘A小区14-2-802房间（客厅）全年日照累计时间，并得出该房间全年日 照时间超过六小时的天数和具体日期。问题三，在部分房间已售出的情况下，仅考虑光照条件，给消费者提供最佳选房方 案。问题四，在考虑日照条件下，结合价格、交通、环境和噪音的影响，给出此时消费 者C最优选房方案。楼价在不同楼层相应有所不同，周围环境，北侧有一条河流，并有 若干配套设施：地铁，铁路，国道，发电烟囱。综合各种因素得出最佳选房方案。问题五，该小区已建成地下停车场，建立合理模型验证该停车场是否合理，如果不 合

7、理重新设置该停车位。附件三中方格第一行是车位号，第二行是对应的房间号问题分析问题一的分析在题设前提条件（不考虑天气因素）下，14-2-8-2房间可以享受日照的时间区间仅受 前方障碍物遮蔽的影响，根据附件住宅楼分布情况分析，#7和#8楼有可能遮挡#14楼障碍物的阴影位置与障碍物的尺寸，楼宇密度（前后楼间距） ，太阳方位，太阳高度等 因素有关，其中障碍物尺寸和建筑密度可从附件查得，太阳方位角和太阳高度角可由相 关公式推导得到，进而可以建立楼宇遮光时间的物理模型，求得可以获得日照的时间区 问。问题二的分析相比于问题一，问题二需要考虑日期的变化；根据相关天文学知识，赤纬角随每年积日数近似的周期变化，进

8、而影响日出日落时间和太阳高度角。因此我们考虑将365天每天的赤纬角值代入第一问的模型并求和，即可较精确的得出一年内享受日照的累计时间。随着赤纬角每年的连续性变化，每日日照时间必定是连续变化，且在夏至日最长（6h）,冬至日最短（6h）,因此在夏至与冬至日间使用二分搜索算法，即可快速确定6小时日照时间的分解日期，求得每日日照时间超过 6小时的日期及天数。问题三的分析相对于问题二，我们需要考虑小区内每一栋楼的每一层的光照时间，因此我们分析 并简化了小区住宅楼的相对位置关系，将楼层（1-34）也作为变量代入问题二的模型, 即可得到任一栋楼的任一层的年受光累计时间；根据采光因素优先的前提，结合附件二 的

9、销售情况，我们即可给出消费者选房的最优方案。问题四的分析问题四要求我们考虑价格，交通，环境与噪音的影响，给出消费者c的最优选房方案。其中价格主要与楼层和是否是河景房有关，交通主要与离车库远近，离地铁远近有 关，环境主要与离发电厂和河流远近有关，噪音主要与离马路和高架远近有关。由于权 重带有很大主关性，我们买房时总是想所有条件都满足，如果不行的话就退而求其次， 在从其中选择最好的。有一些论文可知，用户楼层的选择也反应了用户的偏好。问题五的分析针对车位安排是否安排合理，首先建立空间直角坐标系，并以出楼梯口为原点，通 过将车位和各住户的高度确定在同一坐标系上， 计算图中所给车位到坐标系的距离和楼 层

10、的距离，则最佳车位位同事去的最小，将车位的距离做升序排列楼层的高度作降序排 列，则最佳排序为两者之和。针对合理安排车位，我们将充分考虑各种因素，包括低层 住户长时间等待时是否乘坐电梯，电梯的运行模式的选择，出电梯口时距离电梯口的最 近车位，车位与住户楼层的距离，得出最优车位。三、模型假设1对各个房子的影响因素的评价是正确的。2上班高峰期时各楼层乘客充足。3楼房的倾斜形状对太阳影子没有影响 。4楼房的厚度对太阳影子没有影响。5地球是球形的。6天气对太阳没有影响。7太阳高度角是15度时，光线才足够量。8大气层对太阳光线影响忽略。四、符号说明符号含义H太阳局度角5太阳赤纬<P地球纬度t0时角A

11、太阳方位角s太阳位置o地球中心HH该地的水平向EE赤道b楼层图度Ti电梯运行时间x坐标轴2yz lit坐标轴坐标轴车位到坐标轴原点的距离tZi k楼层到坐标轴原点的跑离 电梯乘客数极限五、模型建立与求解5.1问题一的模型建立与求解5.1 . 1问题一的模型建立5.1 . 1. 1太阳高度角H的计算太阳高度角是指太阳光的入射方向和地平面的夹角。我们假设太阳光是不存在折射的平行光1,参考资料得出太阳高度角H的计算公式；sin H = sin :sin 4 -cos :cos、cost式中，太阳赤纬用6表示，地理纬度用中表示（太阳赤纬和地理纬度都是北纬为正， 南纬为负），地方时角用t表示。日出或日末

12、的时刻，可有式取 H =°,而求得其时角to, 即得：cost0 Jsinsin, -tantan、. cos cos、图一太阳高度角示意图又日出或日末似的太阳方位角 %可取H=0,即得：sin 1 0 二 cos sint0求得太阳高度角H ,即得：,H =180： -（90, ： 一、）=90； -（： 一、）图二太阳直射点纬度与高度角关系示意图s太阳包置；o地球中心；邛纬度；H正午太阳图度角；HH该地的水平面；EE赤道 5.1 . 1. 2太阳赤纬角6的计算太阳赤纬及太阳直射点纬度，查阅资料得到赤纬角6的计算方式；、=0.3723+23.2567sin i 0.1149sin

13、21-0.1712sin 31- 0.7580cos10.3656cos 210.0201cos3式中，e称谓日角，即8 = 2C%cx>4。而R由两部分组成，即N=N-N0,其中，N为 365.224积日，即当天日起到当年一月一日的天数； * *N0 =79.6764+0.2422父(N -1985) -INT(N -1985)/4 5.1 . 1. 3时角工的计算时角为OP线在地球赤道平面上的投影与当地时间 12点时地中线与赤道平面上的投影 之间的夹角。t =s。T。- 订-：式中，其中So是当天平时的恒星s0=6h40m+dM3m56s,元旦子夜时的恒星时是6h40m, d是从元旦

14、起算的天数。T0是当时北京时间，AT =120° -九是当地的地理经度与东经120，的差，a是恒星的赤经。5.1 . 1. 4太阳方位角A的计算八 .,cos; sin 、A 二 a sin()cosHA太阳方位角；H太阳高度角；6太阳的赤纬；T太阳的时角5.1 . 1. 5遮光区域占的求解的计算、=arctanH ssi楼间距；H太阳高度角得到太阳方位角和太阳高度角可以在天穹上确定太阳的位置2,则障碍物在地面上的影子可以由投影定理的出。但考虑到物体影子不但有长度变化还有方位变化，如下图 （三-C）正方形物体在地面上的影子范围变化示意图，其大小方位受日期，时间，维度等 因素影响，转化

15、为数学表达很繁杂；其次考虑从观测点视角看障碍物ABCD ,在天穹面上得到阴影部分aPK ,当太阳 处于阴影部分时即被遮住，但阴影部分 部孤 为空间曲面，几何关系复杂，求其与太阳 轨迹方程的解析解较为复杂。不利于求解。我们考虑将太阳位置离散化，以1分钟为步长，如下图（三-D）计算太阳位置，在空 间坐标系中求目标点M与太阳连线的直线方程，将楼宇间距y代入直线方程，可解得 一高度值h ,通过对比h与障碍物的高度H即可快捷判断此时此刻太阳是否被障碍物遮 住3,通过计数被太阳遮盖点的个数即可推算出某日光照时长与光照时间范围其次：对于目标楼层n,以其为原点建立的坐标系中障碍物的高度会相对减少n层,当太阳不

16、处于天穹阴影的边界角内时，太阳即使低于障碍物也不会被遮挡。 1图三阴影部分示意图5.2 数据的处理与分析利用matlab编写程序解得当位于视野角范围内比较 &与78.3米的差。 = 8-78.3当至。时没有遮光，当AW。时遮光，通过matlab得出结果为日照时间为9:00-10: 20和12: 40-1:32和 15:36-16: 005.2问题二的模型建立与求解通过简化图得到高楼A受到前方区域影响的视野角图四高楼遮挡视野角图表一各个楼的遮盖关系表视野角范围轴线距离（m）视野角范围 轴线距离（m）2F3F-73°51°54.26711F4F-24 °-9。

17、170.2573F4F-79°-67 °62.79112F3F22 °38°145.7525F6F-73°-66 °68.45412F4F-3 °14°153.9527F17F-6。41°54.29612F7F-61 °-48 °72.4747F18F-44 °-19 °82.99412F17F-34 °-17 °127.8778F17F56°74°86.11313F7F-28 °14°58.0388F18

18、F-8。40°52.58613F8F-51 °-32°97.779F1F-14 °14°95.9814F7F34°59°86.2239F2F-47°-31 °125.0314F8F-14 °21°74.92810F1F30°47°126.34817F4F57°72°75.99510F2F-14 °14°97.6917F5F-25 °26°48.7910F3F-24°-19 °139.07

19、717F6F45°27°104.5811F2F20°40°118.79418F5F47°72°69.97511F3F-5 °14°135.90218F6F-19 °20°66.1216在第二问的基础上遮光区域6的求解的计算、=arctanH sS楼间距；H太阳高度角通过调整太阳赤纬角6和时角T利用穷举法得出各个高楼的遮挡时间， 如下表得A小区 14-2-802房间（客厅）到全年享受日照累计时间为 208110分钟。14-2-802房间（客厅）全年 享受日照时间超过6小时的天数为235天，日期为2

20、月25号到10月165.2问题二的模型求解图五 各楼层年遮光时间随楼层高度变化图最优解如附表15.4.2 .模型的建立（一） 选房步骤买房时，我们把选房定为两个步骤，第一步，选择区位，第二步，选择楼层。影响区 位的因素主要是出行便利，环境，噪音。影响楼层选择的因素主要是价格。消费者会从 各种不同的因素考虑买房，我们把满足因素一最好的一批房子列为集合M满足因素二最好的一批房子列为集合，犷法满足集合二次好的一批房子.场&。以此类推。1令周言n也,二门工1n一班"二S 02若5上& 则根据价格，采光选择具体楼层。2若5=%则降低标准，令M址n.T&bR.%R且 =

21、6,直到S=O为止，然后执行第二（二）集合分类原则1 .环境：根据附近有无河流，电厂，垃圾房。2 .噪音：根据楼房附近马路数量。3 .比较每栋居民楼与商业楼和街道口的距离，选择最小距离，并排序。a 6F 2+ （或一3）：B =;应/）二区T）：DI£=MIN®B）（三）城市高层住宅楼层偏好1 .高层用户比较偏好环境，噪音，但由于楼层较高，一般楼底的噪音对高层没有影响。2 .中层用户考虑比较综合。3 .底层用户优先考虑采光，并且根据第二问计算结果，中高层用户基本不需要考虑采光1-4流程图5.4.3模型的求解在小区平面图上取坐标，并根据公式计算表二DIS距离表跑离排名DIS

22、（距离）1#4334.22#3411.53#14614#96005#5703.56#28317#17965.78#61074.59#101113.610#181290.311#71428.512#111500.113#121524.914#81692.5根据图片 环境a等集合孤!工=#9, #10, #11, #12, #13, #14环境b等集合 叫广#7, #8, #17, #18 环境C等集合 U"=#1, #5, #6 环境d等集合.%或=#2, #3, #4噪音 a 等集合 =#10, #11, #12, #13, #7, #17噪音 b 等集合 网|=#9, #2, #3

23、, #4, #5, #14噪音C等集合 此产#1, #6 出行便利a等集合义工=#4, #3 , #1, #9 出行便利a等集合 叫产#5, #2, #17, #6 出行便利a等集合 洱!产#10, #18, #7, #11 出行便利a等集合 ”#12, #8, #13, #14（一）对高层用户首先考虑环境，其次生活便利。小口=用9所以选择9号楼，由于33层环境较好，所以选择9号楼33层。（二）对中层用户，考虑采光以外所有因素-%# c 2%占'占=g'j根据第二题数据，9号楼采光差别不大，选择价格最低 13层，对低层用户，采光首先 不能太差。采光a等集合周二幽一科峋采光b等集

24、合ML9, m 811，-二, -'Gt 介2口二%d 二0端物c口，24 = 1#9如果将楼层限定在6层以下。%咚咚-“蹶,#6)由于5号楼与6号楼其他条件相同的同时，噪音比较小，生活更便利。由于 5号楼 1-6层采光没有明显差距，选择价格最低的 5号楼1层。1.如果用户对价格考虑较多，那么只考虑第一层与顶层，5号第一层4250元/峭，9号第一层3的3元/杆。5.5问题三的模型建立与求解5.5 . 1双坐标系的建立将附件一所给的图像抽象到 Excel表格上，并赋予每个地下车位双坐(x：y：)(x2py2p)坐标，以两个电梯出口为原点(x0y0 ),(x(py0P ),建立双坐标系，并

25、以楼房的高度为分别为*z0P,我们将通过地下车位唯一的确定下来。车位到(x0 y0 ),(x°p y°p )的距离可以简化为，可以简单的求1；=年7；和1=x+yp ,求得 =1： -1的最小值则能唯一确定确定车位的位置， 通过图上图形的分 布很容易确定一点到(x0 y0),(x°p y )的最近距离，由图形的分布特点，位于两电梯口 的中位线左侧的用左坐标系，反之用右坐标系，则简单的计算空间距离。要求1： =x：+y：和1=xiP+yiP分别按升序排列，1t ：1；,口附：1P 11"吠zt,ZiP按降序排列,Z:*如反P P P PZ ： Z1 , Z

26、2 Z则两者相互组合形成最佳车位。图中不符合此最佳排列，所以不合理。常见的几种 电梯运行方案的比较。为简化描述时又不失一般性，我们假设两台电梯同时独立运行。电梯方案的比较有 多重标准，在此我们考虑如何在上下班高峰期时，将所有等待的乘客快速送到目的地， 一次我们采用侧重乘客等待时间的优化，时间最大、最小群控方法。即用乘客等待时间 和被运送时间和最小，最为评价依据，并依据电梯运行周期与时间之比等于电梯在一个 周期内运送的乘客数于总乘客数之比的比例原则，对常见的运行模式进行描述，常见电10梯运行方案的描述。假设每层楼工作人员的人数(即电梯需要向每层楼运送乘客的人数)均为 m电梯最 多能运送k人，假设

27、该办公楼总共b+1层。假设每层楼之间电梯的平均运行时间是 i秒, 则平均停留时间为t2秒。5.5 . 1. 1分层法该方案允许电梯可以在任意层停靠，随机运行，两台电梯平均运行周期均为(2kb*ti+bt)秒，共运送乘客2*女人，运送所有乘客共b*m人,所用时间为T ,并依据“电 梯运行周期与运行总时间之比等于电梯在一个周期内运送的乘客数与乘客总数之比”的“比例”原则得：2k 2bti bt2 bm T 解得；b2m(2ti t2) I 二2k5.5 . 1. 2奇偶层运行方案该方案规定两台中的一台停靠在奇数层，另一台停放在偶数层。因为一共b + 1层，k *,所以停靠在奇数层的电梯运行时间为(

28、2 bt1+bt2%)秒，完成运送乘客时间T1,并依据“电梯运行周期与运行总时间之比等于电梯在一个周期内运送的乘客数与乘客总数之 比”的“比例”原则：kkm2bti 当2.2,、b m(4t1 t2)T 二4k而停靠偶数层的电梯的运行周期为(2kb*-1)t+1bt/)秒，共运送乘客20人,运送所 有乘客完成运送至偶数层的乘客所用的时间2T,并依据“电梯运行周期与运行总时间之比等于电梯在一个周期内运送的乘客数与乘客总数之比”的“比例”原则：一 2 一一丁 b m(4bt1 bt2 t2)T 二4k5.5.1.3分段运行方案此方案以b*n十1层为界分为上下两段，一台运行在1- b*n+1层，运送

29、所有乘客用 时工，并依据“电梯运行周期与运行总时间之比等于电梯在一个周期内运送的乘客数与 乘客总数之比”的“比例”原则：k 2bti (b - nb)t2(1 - n)bmT1k _ 2nbt1 nbt2nbm11整理得:b2m(1 - n)(2tl t2 - nt2)b2n2m(2t2 t)T =max(U)最优方案为 2*, 2 , *、 2、丁 b m(1 m )(2t1t2nt2) b (n ) m(2t1t2)T 二二5.5.1.4分层次法该万案同样将以b n+1层为界分为上下两段，第一台运行第1层，第bn + 1层至第 b+1层，而另一台电梯则可以停在所有楼层，在等可能几率乘坐的情

30、况下，故所需的总时间为运送目的地为第2层至nb*十1层的乘客所用时间，第一台电梯完成运送所有乘客所用时间1T,并依据“电梯运行周期与运行总时间之比等于电梯在一个周期内运送的乘 客数与乘客总数之比”的“比例”原则：k , nbbbnm2bti bt2Tbm（2Lt2）I -k由于总时间为运送目的地为第 2-层的乘客所用时间4,而这些楼层之间仅有一台电 梯按平均原则运送乘客，故当b*n 1 ,即b*n22时无论n为何值，该电梯每运行一次 可向第2层-（b*n+1）层运行（必）名乘客，运行总数为（m*%）,运行一次所用时间（2 bt1 +bt1）,故即与n无关 5.5.15逐层运输法两台电梯运送完一

31、层运送另一层。故时间：T 二T5.5.2 常见电梯运行方案比较（1）、（2）、（3）、（4）、（5）式。我们容易得到：综上考虑电梯的运行效率可得：逐层运输方案 分段运行方案 奇偶层运行方案 分层次与分段 相结合 分层次运行方案的方案。不过考虑到逐层运输不符合实际。因此我们得出结论：分段运行方案是最及时的将所有等待的乘客快速运至目的地，尽快地疏散等候区的乘客的最优调度方案5.5.3 实际问题的求解假定b =34,t1 =3秒，t2=10秒代入公式的*（2t1 t2）2 /-& t2）n 二 n =23n = 0.535则得出层数为b*+1=18.19,取整的19或者18带入得12.2_*

32、 、，-、, 2 , * 、 2、b m(1 - m )(2ti t2 - nt2) b (n ) m(2ti t?)I 二二当取19时，n计算T19 =288.8 ,当取18时，n =18计算18 = 259.2,则应选183434层为分界点。则可以得到各个楼层到达地面时间，ti =3*( i 1)十10楼梯口到达车位所用时间t；,则1 =七+t；按照ti的升序和t：的降序结合获得最优组合。考虑到较低层住户会选 择不乘坐电梯，结合楼梯周围车位数为16,则1-4楼住户走楼梯，并且参照正常爬楼梯 时所用时间，则四楼住户用时最长为t4=80,远大于需要等待时间206秒，5-18楼住户 乘坐一电梯，

33、19-34乘坐另一电梯，根据到达时间长短安排车位的远近。表三优化后车位安排表单元西户二单元东户一单元西户一单元东户341-1161341-1160 341-1566r 341-1567331-1159331-1334331-1588331-1631321-1335321-1336321-1630321-1629311-1337311-1424311-1412311-1411301-1423301-1422301-1569301-1570291-1554291-1555291-1628291-1627281-1421281-1556281-1626281-1571271-1420271-1557

34、271-1410271-1625261-1558261-1559261-1572261-1409251-1419251-1559251-1624251-1573241-1643241-1642241-1408241-1623231-1817231-1818231-1574231-1407221-1641221-1819221-1622221-1575211-1640211-1639211-1406211-1621201-1820201-1821201-1575201-1405191-1822191-1418191-1620191-1577楼181-1417181-1560181404181-1

35、619层171-1561171-1416171-1578171-1403161-1562161-1637161-1618161-1579151-1636151-1635151-1402151-1617141-1823141-1415141-1580141-1401131-1563131-1634131-1824131-1862121-1414121-1564121-1863121-1864111-1633111-2101111-1865111-1866101-1565101-1632101-1867102-18891-187091-186992-18792-18681-186882-17382

36、-131282-131372-17472-17572-133172-133262-17662-17862-149162-133352-17952-18052-149052-133442-18142-18242-148942-133532-18932-18432-148832-133622-18522-18222-148722-133712-18512-18412-148612-133813六、模型评估与改进6.1优点:模型一，对时间及太阳位置进行离散化处理，大幅度简化了计算，避免了推导复杂 的太阳轨迹方程的同时通过提高离散密度保证了精度模型二，基于模型一，拓展日期，简单易用模型三，精确计算了每

37、一栋楼每一层每一年的的累计日照时间，精度高可靠性强模型四，根据客户的楼层推测客户的偏好，给出个性化的推荐。模型五，充分考虑低层住户不会长时间等待电梯而是通过楼梯下楼。考虑不同的电 梯运行方式得到运行最省时的方式。6.2缺点：模型一，离散化太阳位置求得的阴影变化与现实有着微小的误差；计算障碍物遮挡 范围时将立体的楼简化为一竖直平面，在轴线倾角较大时会产生误差；没有考虑障碍物 的实际形状未考虑太阳折射率（日最大折射误差为 35”左右）模型二，计算以日期推算得到的赤纬；角与实际观测值有一定误差，积累了模型一 的误差模型三，算法时间复杂度较高，对大规模问题拓展性差；楼宇位置关系矩阵是基于 图片推算，有

38、一定的计算误差；模型四，根据售出房间的楼层，环境评价，生活便利程度分类，根据分类结果推断 客户潜在需求。模型五，统计住户用车频率和用车高峰时期，同时结合住户楼层高度，进行车位的 合理分配。没有统计住户用车频率和用车高峰时期，没有结合住户楼层高度，进行车位 的分配。七、参考文献1卓金武，MATLAB数学建模中的应用第2版，北京，北京航空航天大学出版社，2014.9。2韩中庚，数学建模方法及其应用，高等教育出版社，2009.9。3邓维斌、唐兴艳、胡大权、周玉敏，SPSS1甑计分析实用教程中文版，电子工业出 版社 2012.3。4周品、赵新芬，MATLAB数学建模与仿真，国防工业出版社，2011.7

39、。14代码一：（计算日出日落方位角）function arealeft arearight=area(weidu,chiwei)syms t0;weidu=weidu.*pi()/180;chiwei=chiwei.*pi()/180;bianjie=eval(solve(sin(weidu).*sin(chiwei)+cos(weidu).*cos(chiwe i).*cos(t0),t0)*180/pi();if bianjie>90bianjie=90;endarealeft=-bianjie+45;arearight=bianjie-30;代码二(计算某天被遮光的时间)function