数学建模综合题影院座位设计问题

1、数学模型张峰华材料学院材料成型及控制工程04 班20123631刘泽材料学院材料成型及控制工程04 班20123627杨海鹏材料学院冶金工程 03 班20123203精选文库一、问题重述影院座位的满意程度主要取决于视角 和仰角 ，视角是观众眼睛到屏幕上下边缘的视线的夹角，越大越好；仰角是观众眼睛到屏幕上边缘视线与水平线的夹角，太大使人的头部过分上仰，引起不适，一般要求仰角不超过 300 ；记影院的屏幕高为 h ，上边缘距离地面高为 H , 影院的地板线通常与水平线有一个倾角，第一排和最后一排与屏幕水平距离分别为d, D ，观众的平均座高为c（指眼睛到地面的距离） ，已知参数 h =1.8.H

2、=5， d4.5, D19 ， c =1.1( 单位 m)。求解以下问题：(1) 地板线的倾角100 时，求最佳座位的所在位置。(2) 地板线的倾角 一般超过 200 ，求使所有观众的平均满意程度最大时的地板线倾角。二、问题的分析电影院座位的设计应满足什么要求，是一个非常现实的问题。根据题意观众对座位的满意程度主要取决于观看时的视角 和仰角 ， 越大越好，而 越小越好，最佳位置就是要在这两者之间找到一个契合点，使观众对两者的综合满意程度达到最大。本文通过对水平视角 和仰角 取权重，建立适当的坐标系，从而建立一个线形型满意度函数。针对问题一，已知地板线倾角，求最佳座位所在，即将问题转化求综合满意

3、度函数的最大值，建立离散加权的函数模型并利用 Matlab 数学软件运算求解；针对问题二，将所有观众视为离散的点，要使所有观众的平均满意程度达到最大，即将问题转化求满意度函数平均值的最大值。对此利用问题一所建立的满意度函数，将自变量转化为地板线倾角；在问题二的基础上对地板线形状进行优化设计， 使观众的平均满意程度可以进一步提高。本文在满意度呈线性的基础上来建立模型的，为使模型简化，更好地说明问题，文中将作以下假设。三、模型假设1. 忽略因视力或其他方面因素影响观众的满意度；2. 观众对座位的仰角的满意程度呈线性；3. 观众对座位的水平视角的满意程度呈线性；4. 最后排座位的最高点不超过屏幕的上

4、边缘；5. 相邻两排座位间的间距相等，取为0.8 m ；6. 对于同一排座位，观众的满意程度相同；7. 所有观众的座位等高为平均座高；8. 影院的的地板成阶梯状。-1精选文库四、符号说明水平视角视高差，即从眼睛到头顶的竖直距离仰角S观众对水平视角为的满意程度地板线与水平线的倾角S观众对仰角为的满意程度d第一排离屏幕水平距离S平均满意程度D最后一排离屏幕水平距离c , c视角 、仰角在综合满意度 Si 中的权重h屏幕的高度l相邻两排座位间沿地板线方向的间距H 屏幕上边缘离地面的高度五、模型的建立与求解5.1 问题一每一个到影院看电影的观众都想坐在最佳位置， 而对座位的满意程度主要取决于两个因素：

5、水平视角 和仰角 ，且视角是观众眼睛到屏幕上下边缘的视线的夹角，越大越好，仰角是观众眼睛到屏幕上边缘视线与水平线的夹角，太大使人的头部过分上仰，引起不适，要求不超过 300 。模型的建立：仰角在满足条件的范围内, 观众满意度只取决于视角以第一排观众的眼睛为原点，建立平面直角坐标系，如图1 所示：屏yh 幕BEH c h x tan视觉线SH hP地板cOcx tanx地面MxNdD d图 1影院座位设计的剖面图其中， AB 为屏幕， MS 为地板线， OE 为所有的观众的眼睛所在的直线。则由图可设视觉线 OE 上任意一点 P 的坐标为 ( x, x tan) ，屏幕上下点的坐标分别为A( d,

6、 H c) ，-2B( d, Hh c) ， AP 的斜率记为 kAP , BP 的斜率记为 kBP 。由斜率公式得：k APtanx tanH c ， k BPtan()x tanHh cx( d )x (d )则直线 AP 和 BP 的斜率与夹角满足如下关系：tank BPk APd ) 2h(xd )1 kBP kAP( x(x tanHc)( x tanHh c)仰角满足条件： 0,30所以：0tan3 3x tanHc3 30( d )xH c3 3 dxHc3 3tantan由公式 (1.1) (1.2)得到模型为：maxarctanh( xd )( xd) 2(x tanHc)(

7、 x tanHhc)0xDds.t. Hc33dxHc3 3tantan精选文库(1.1)(1.2)(1.3)模型的求解当10 时，用 Matlab 软件运算求解（程序见附录1），得最大视角为13.9522 ，仰角为30 ， x1.7274米。即 P 点的坐标为 (1.7274,0.3046) 为最佳位置。离屏幕的水平距离为 4.51.72746.2274米 。模型的建立：离散加权模型在地板线上的座位可视为是离散的点，设两排座位在地板线方向上的前后间距为 l （查阅相关资料间距一般取 0.8 米），则在水平方向的间距为 l cos ，考虑仰角和视角对观众的满意度为主要因素。对模型进行修正，将座

8、位连续情况进行离散化可以得到：x tanHc(k 1)l cos tanHc(2.1)tanx ( d)(k1)l cos( d)tanh( k1)l cosd )1)l cosd ) 2(k1)l cos tanHc)( k1)l cos tanH h c)( k(2.2)其中， k1,2,3, n ， n 为地板线上的座位的总排数，且n 14.5 1 19。l cos一般说来，人们的心理变化是一个模糊的概念。本文中观众对某个座位是否满意的-3Hxtan精选文库看法就是一个典型的模糊概念。由模糊数学隶属度的概念和心理学的相关知识，根据题意，在假设条件下，对于第k 排座位，建立观众对视角、仰角

9、的满意度函数 1 如下：tanS ktankmaxtantanmin(2.3)minS k1tan tankmaxtantanmin(2.4)min式中k , k 为第 k 排座位上观众视角和仰角， max ,max 表示在 给定的情况下最优满意度，min , min 表示在 给定的情况下最差满意度。视角、仰角 在综合满意度 Sk 中的权重分别为 c, c ，建立第 k 排座位综合满意度函数如下：Skc S kc Scc根据地板线倾角k(2.5)10 ，通过计算可以得出 5.421015.8975 ，4.045140.9149 ，主观给定权重 C0.6,C0.4 ，根据模型的建立，可以得出：S

10、kc S kcS k0.6S k0.4S k3.1596 tank 0.5025 tan k 0.1357(2.6)cc0.60.4将式 (2.1)和式 (2.2)带入公式 (2.6)得到优化模型为：max Sk(k1)l cos0.5025( k1)l cos tan(k1)l cos0xDdHc3 3 ds.t.3 3tanx(k1)l cos3.1596* h(k1)l cos d )d ) 2( k 1)l cos tanH c)( k 1)l cos tanH h c)Hc( d )0.1357c, k1,2,3,195.1.4 模型的求解用 Matlab 软件运算求解（程序见附录2

11、）可得： x2.3635米， k 4 排，最大满意度为 S4 0.6176 , 最大视角为13.1282，仰角为26.9084 ,最佳位置离屏幕的水平距离为 4.5 2.3635 6.8635米 。5.2 问题二5.2.1 模型的建立要使所有观众的平均满意程度达到最大，即需求S的最大值。由模型可知，第 kn排观众的满意度为 S, 则观众平均满意程度函数为： SSk n ，平均满意度 S 的大小由每一排的满意度所决定，而又是由仰角和视角k 1所决定。所以，要使观众的满意程度达到最大，取决于两个方面： (1) 仰角不超过条件的座位所占的比例越大，观众的平均满意程度就越大； (2) 所有座位的视角的

12、均值越大，观众的平均满意程度就越大。由式 (1.1) 可知，地板线倾角的改变将同时使所有座位的仰角和视角的大小发生-4精选文库改变，且在某一座位 ( 即 x 取某一定值 ) ，在逐渐增大的过程中仰角逐渐减小，视角逐渐增大，见图 2 所示。仰角不超过条件的区域扩大，即地板线倾角 越大 , 仰角不超过条件的座位所占的比例越大。7.15随 的变化曲线16随的变化曲线7.114127.05化 7化 10变变 8角 6.95角66.946.8526.800246810121416182002468101214161820角变化角变化图 2视角和仰角随变化的变化曲线第一排观众的仰角为40.9149，不满足

13、仰角的条件，由模型可知第k 排座位所对应的仰角的正切值：(k 1)l cos tanHc1,2,3,ntan k( d ), k( k 1)l cos其中 n 为地板线上的座位的总排数： n14.5 1，随着地板线倾角的变化，相l cos邻两排座位间的间距 l 不变，但相邻两排座位间的水平间距会发生改变。由于地板线倾角 不超过 20 ，所以 19 n 20 ，并限制最后一排观众的视高不要超过屏幕的上边缘，即15.0543 。由模型可求出第k 排座位所对应的水平视角的正切值为：tanh( k1)l cosd )d ) 2(k 1)l costanHc)( k 1)l costanH h c)(

14、k 1)l cos5.2.2 模型的求解让地板线倾角在 0,20 内逐一取值 , 步长为 0.01 ；让 x 在 0,14.5 内逐一取值 , 步长为 0.01 。对一个取定的，判断 x 所在的位置仰角是否超过30 ，若超过，则该座位的综合满意度必须同时考虑仰角和视角的取值；否则，只需要考虑视角的取值，把所有座位的综合满意度相加， 并求出观众的平均综合满意度， 判断此时的平均满意度是否最大，最后一排的高度是否超过屏幕的上边缘，并记下最大值时的取值。当取地板线倾角为变化时，通过计算可以得出5.114315.8975 ，040.9149 。-5精选文库由模型的 (2.5)c S kcS k0.6S

15、 k0.4S k（3.1）式得：Skc0.60.4c所以，将式 (2.1)和式 (2.2)带入公式 (3.1)得到平均满意度的优化模型为：nSkmax Sk 1n19n20s.t0x15.0543, k 1,2, , n其中 n取整数0D dx(k1)l cos用 Matlab 软件计算（程序见附录3）可得：最大平均满意度为 S0.6572 ，对应地板线的倾角为15.0543 。5.3 在问题二的基础上对地板线形状进行优化设计，使观众的平均满意程度可以进一步提高。模型的建立与求解由上两问可知 , 观众的满意程度与仰角 , 视角和地板线倾角都有关 , 而每一座位到屏幕的水平距离基本固定不变，考虑

16、观众的满意度，就要考虑仰角，视角随着的变化情况。引理 地板线不管设计成什么形状 , 各排的间距不变， 区别在于各排的高度差如何变化，若竖直方向上的两定点，在与它们相距一定水平距离的竖直方向上有一动点，当该动点位于两定点的垂直平分线上时，动点与两定点形成的视角最大。动点距两定点的垂直平分线越近，动点与两定点形成的视角越大。要使每一个座位所对应的视角取最大值, 对应的 y 值应在直线上 . 设计地板线应考虑以下几个方面： (1) 第 k 排座位所在的位置应高于第 k1排座位所在的高度； (2) 前一排的观众不会挡住后一排观众的视线； (3)视角尽可能大 , 即眼睛的位置应尽可能分布在垂直平分线的附

17、近； (4)仰角的座位所占的比例尽可能大。假设每排座位所在的点构成一条折线，任意相邻两排座位水平间距为l ，第 k 排座位地板线倾角为k ，第 k 排座位与第 k1排座位地板线倾角变化为。从而可得：k0(k 1)，故：nn(k1)l tan kH cl tan(k1) H ctank1k 1k(k1)l (d )(k1)l( d )同理可得：-6精选文库tanh(k 1)l cosd )d )2( k1)l costanHc)( k 1)l costanHhc)(k 1)l cosh(k 1)l cosd )d) 2nn(k 1)l cos(l cosl tan( k1)Hc)(l cosl

18、tan( k1)H h c)k1k1n观众平均满意程度函数为：SSknk 1可算出地板线上的座位的总排数为：n14.5 1，则可计算得当2.5 时，l cosSmax 0.6692。但此时(191) 2.545 ，根据一般习惯，要求地板线倾角20，但此时求得最后一排座位的地板线倾角为45 ，这大大超过观众的心理范围， 因此文中将对此进一步的修改。当 (i1)20时，令 (i1)20 。当20时，即将问题转化为问题二中所建立的模型。由于2.5，则地板线倾角增加到第8 排到达 20，然后保持不变。对于这两种情况，分别代入不同的函数，利用matlab 数学软件求得：满意度函数的最大值 Smax0.6

19、6430.6572 。可以通过利用 Matlab 软件来描点，如图3 所示：2.521.510.50024681012141618图 3从上图可以看出，报告厅座位的前8排呈折线状，以2.5 递增 ,当倾角增加到20 时保持不变，且第一排应抬高1.2米。六、模型的评价与推广6.1模型的评价模型的优点 :模型抓住影响观众满意程度的主要因素 ( 仰角和视角 ), 合理构造满意度函数 , 过程清晰明了 , 结果科学合理。模型具有较好的通用性 , 实用性强 , 对现实有很强的指导意义。模型的不足以及需要改进的地方:模型主观假设同一排座位观众的满意程度相同 , 实际情况并非如此 , 这就使得我们的模型对解

20、决实际问题时有一定的局限性。模型建立的过程中 , 以观众眼睛所在的点为坐高点, 没有考虑前排观众额部对后排-7精选文库观众的遮挡 , 需要进一步的考虑在内。6.2模型的推广本文中所建立模型的方法和思想对其他类似的问题也很适用 , 所建立的模型可用于大型场所的座位的设计与安排，以及彩民对中奖率的满意程度等问题上。同时对于已知剖面来分析物体的形状这一类型问题的处理有很好的参考价值 . 例如 : 运用该模型去解决会议厅、报告厅的布局 , 灯塔高度的设计等相关的问题。因此具有很强的实用性和推广性。八、附录：附录一clearclcH=5;h=1.8;D=19;d=4.5;c=1.1;l=0.8;pi=3

21、.1415926;f=10;for Q=0:0.1:20for l=1:floor(14.5/cos(Q/180*pi)+1)x=(l-1)*cos(Q/180*pi);T=tan(Q/180*pi);A=(d+x)*h/(d+x)2+(H-c-T*x)*(H-h-c-T*x);if f>Af=A;endendendfor Q=0:0.1:20for l=1:floor(14.5/cos(Q/180*pi)+1)x=(l-1)*cos(Q/180*pi);T=tan(Q/180*pi);A=(d+x)*h/(d+x)2+(H-c-T*x)*(H-h-c-T*x);if f=Afprint

22、f('Q is:%dn',Q);fprintf('k is:%dn',l);endendendf附录二clear-8精选文库clcH=5;h=1.8;D=19;d=4.5;c=1.1;l=0.8;pi=3.1415926;t=10;for Q=0:0.1:20for l=1:floor(14.5/cos(Q/180*pi)+1)x=(l-1)*cos(Q/180*pi);T=tan(Q/180*pi);B=(H-c-T*x)/(d+x);if t>Bt=B;endendendfor Q=0:0.1:20for l=1:floor(14.5/cos(Q/18

23、0*pi)+1)x=(l-1)*cos(Q/180*pi);T=tan(Q/180*pi);B=(H-c-T*x)/(d+x);if t=Bfprintf('Q is:%dn',Q);fprintf('k is:%dn',l);endendendt附录三clear;%clc;H=5;h=1.8;D=19;d=4.5;c=1.1;Q=0.1763; %tan(10/180*pi);s=0;for x=2.36353.15143.93924.72715.51496.30287.09067.87858.66639.454210.242011.029811.817712

24、.605513.393414.1812-9精选文库t=3.1596*(h*(x+d)/(x+d)2+(x*Q-H+c)*(x*Q-H+h+c)-0.5025*(-(x*Q-H+c)/(x+d)+0.1357;if s<ts=t;endendfor x=2.36353.15143.93924.72715.51496.30287.09067.87858.66639.454210.242011.029811.817712.605513.393414.1812t=3.1596*(h*(x+d)/(x+d)2+(x*Q-H+c)*(x*Q-H+h+c)-0.5025*(-(x*Q-H+c)/(x+

25、d)+0.1357;if s=tfprintf('nX is:%d',x);fprintf('nk is:%d',x/(0.8*cos(10/180*pi)+1);fprintf('nais:%d',(atan(h*(x+d)/(x+d)2+(x*Q-H+c)*(x*Q-H+h+c)/pi*180); fprintf('nb is:%dn',(atan(-(x*Q-H+c)/(x+d)/pi*180);endends附录四clear;clc;H=5;h=1.8;D=19;d=4.5;c=1.1;l=0.8;pi=3.1415926

26、;t=0;for k=1:20for Q=0:0.01:20z=h*(k-1)*l*cos(Q/180*pi)+d)/(k-1)*l*cos(Q/180*pi)2+(k-1)*l*cos(Q/180 *pi)*tan(Q*pi/180)-H+c)*(k-1)*l*cos(Q*pi/180)*tan(Q/180*pi)-H+h+c);if t<zt=z;endendendfor k=1:20-10精选文库for Q=0:0.01:20z=h*(k-1)*l*cos(Q/180*pi)+d)/(k-1)*l*cos(Q/180*pi)2+(k-1)*l*cos(Q/180 *pi)*tan(

27、Q*pi/180)-H+c)*(k-1)*l*cos(Q*pi/180)*tan(Q/180*pi)-H+h+c);if t=zfprintf('Q为： %dn',Q);fprintf('k为： %dn',k);endendendt附录五clear;clc;H=5;h=1.8;D=19;d=4.5;c=1.1;l=0.8;pi=3.1415926;t=0;for k=1:20for Q=0:0.01:20z=h*(k-1)*l*cos(Q/180*pi)+d)/(k-1)*l*cos(Q/180*pi)2+(k-1)*l*cos(Q/180 *pi)*tan(Q*pi/180)-H+c)*(k-1)*l*cos(Q*pi/180)*tan(Q/180*pi)-H+h+c);if t<zt=z;endendendfor k=1:20for Q=0:0.01:20z=h*(k-1)*l*cos(Q/180*pi)+d)/(k-1)*l*cos(Q/180*pi)2+(k-1)*l