发明名称 高速公路长隧道入口段昼间照明优化方法
摘要 本发明公开了一种高速公路长隧道入口段昼间照明优化方法,该方法采集隧道环境照度、驾驶人视觉适应时间、驾驶人瞳孔面积数据,建立驾驶人瞳孔面积与照度模型,照度与隧道纵深模型和驾驶人视觉适应模型;利用软件做出驾驶人瞳孔面积、视觉适应时间、相应环境照度之间的三维曲面,根据曲面所得参数建立函数模型并根据计算最终求得隧道纵深、车速与隧道优化照度之间的关系。本方法不仅能消除昼间隧道环境行车时,因照度差异造成的暗适应问题,减少视觉障碍诱发的交通事故,同时能大幅度降低隧道照明设施电能消耗。
申请公布号 CN103383708A 申请公布日期 2013.11.06
申请号 CN201310279802.6 申请日期 2013.07.04
申请人 长安大学 发明人 刘浩学;赵炜华;谢陈江;林淼
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西安恒泰知识产权代理事务所 61216 代理人 李婷
主权项 一种高速公路长隧道入口段昼间照明优化方法,其特征在于,该方法包括以下步骤:步骤一,采集驾驶人瞳孔面积、隧道环境照度、驾驶人视觉适应时间数据;步骤二,根据步骤一采集的驾驶人瞳孔面积、环境照度和驾驶人视觉适应时间数据分别建立驾驶人瞳孔面积与照度模型,照度与隧道纵深模型和驾驶人视觉适应模型;其中驾驶人瞳孔面积与照度模型为:取驾驶人瞳孔面积为Q、隧道环境照度为l,对昼间高速公路隧道入口段数据样本建立ln(Q×l)与ln(l)的关系图并进行回归分析,驾驶人瞳孔面积与环境照度存在如式(1)的关系:ln(Ql)=8.368‑0.576lnl    (1)将式(1)进行转换,整理得到驾驶人瞳孔面积Q与隧道环境照度l的关系函数,如式(2)所示:Q=e8.368l‑1.576    (2)照度与隧道纵深模型:隧道环境照度l与隧道纵深d是乘幂关系,拟合的关系式为:l=17963d‑0.959    (3)驾驶人视觉适应模型:构建驾驶人瞳孔面积Q、驾驶人在隧道中视觉适应时间t和相应隧道环境照度l的二元四次函数:Q=x0+x1t4+x2t3+x3t2+x4t+x5l4+x6l3+x7l2+x8l+x9tl3+x10tl2+x11tl+x12t2l2+    (4)x13t2l+x14t3l+ε式中:ε~N(0,σ2),表示随机误差的变量;通过对驾驶人瞳孔面积Q与驾驶人在隧道中视觉适应时间t和相应隧道环境照度l多次测量统计,取得i组数值:(ti4,ti3,ti2,ti1,li4,li3,li2,li1,tili3,tili2,tili,ti2li2,ti2li,ti3li,yi),i=1,2,…,n.(5)构成方程组如下: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <msup> <msub> <mi>t</mi> <mn>1</mn> </msub> <mn>4</mn> </msup> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>+</mo> <msub> <mi>x</mi> <mn>14</mn> </msub> <msup> <msub> <mi>t</mi> <mn>1</mn> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>&epsiv;</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <msup> <msub> <mi>t</mi> <mn>2</mn> </msub> <mn>4</mn> </msup> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>+</mo> <msub> <mi>x</mi> <mn>14</mn> </msub> <msup> <msub> <mi>t</mi> <mn>2</mn> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>&epsiv;</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>n</mi> </msub> <mo>=</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <msup> <msub> <mi>t</mi> <mi>n</mi> </msub> <mn>4</mn> </msup> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>+</mo> <msub> <mi>x</mi> <mn>14</mn> </msub> <msup> <msub> <mi>t</mi> <mi>n</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>n</mi> </msub> <mo>+</mo> <msub> <mi>&epsiv;</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>其中:ε1,ε2,…,εn相互独立,将方程组转化成矩阵形式如式(7):Y=λX+ε(7)其中: <mrow> <mi>Y</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>&lambda;</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mn>1</mn> </msub> <mn>4</mn> </msup> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mn>1</mn> </msub> <mn>3</mn> </msup> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mn>1</mn> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mn>2</mn> </msub> <mn>4</mn> </msup> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mn>2</mn> </msub> <mn>3</mn> </msup> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mn>2</mn> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mi>n</mi> </msub> <mn>4</mn> </msup> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mi>n</mi> </msub> <mn>3</mn> </msup> </mtd> <mtd> <mo>.</mo> <mo>.</mo> <mo>.</mo> </mtd> <mtd> <msup> <msub> <mi>t</mi> <mi>n</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>n</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>X</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>14</mn> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow>向量Y,λ都是已知的,利用最小二乘法,计算出向量X的估计值,过程如下:根据最小二乘法定义可以得到式(8): <mrow> <mi>F</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>,</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>,</mo> <msub> <mi>x</mi> <mn>14</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> <mo>-</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>-</mo> <msub> <mi>x</mi> <mn>14</mn> </msub> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>向量X的最小二乘估计值,应满足F取最小值时的解,对Q求偏导数,可以得到正规方程组如式(9): <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>nx</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <msub> <mrow> <mo>+</mo> <mi>x</mi> </mrow> <mn>14</mn> </msub> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>&Sigma;</mi> <msub> <mi>y</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>0</mn> </msub> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> <mo>+</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>+</mo> <msub> <mi>x</mi> <mn>14</mn> </msub> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> <msub> <mi>y</mi> <mi>i</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>0</mn> </msub> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>i</mi> </msub> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>4</mn> </msup> <mo>+</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>&CenterDot;</mo> <mo>+</mo> <msub> <mi>x</mi> <mn>14</mn> </msub> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>i</mi> </msub> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>=</mo> <mi>&Sigma;</mi> <msup> <msub> <mi>t</mi> <mi>i</mi> </msub> <mn>3</mn> </msup> <msub> <mi>l</mi> <mi>i</mi> </msub> <msub> <mi>y</mi> <mi>i</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>相应的矩阵形式为λTY=λTλX,λTλ列满秩,故向量X有唯一解,将其代入驾驶人瞳孔面积Q与驾驶人在隧道中视觉适应时间t和相应隧道环境照度l的二元四次函数式(4)中,即获得驾驶人视觉适应函数模型,如式(10)所示:Q=‑0.151t4‑35417×10‑10l4‑1.7093×1O‑5l3‑0.0146l2‑8.3614×10‑6tl3+0·0104tl2‑0·00l2t2l2+0·0151t3l    (10)步骤三:昼间长隧道入口段照明优化利用matlab软件做出驾驶人瞳孔面积Q、视觉适应时间t、相应环境照度l之间的三维曲面,驾驶人视觉适应函数模型中,瞳孔面积作因变量,分别对驾驶人在隧道中视觉适应时间t和相应隧道环境照度l求偏导数,建立式(11)方程组: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>Q</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>Q</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>l</mi> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>得到一系列数据组(t,l),即为驾驶人适合相应隧道环境照度变化视觉要求的控制参数;选取数据定为视觉适应时间tn对应的照度值ln;将驾驶人适和相应隧道环境照度变化视觉要求的照明控制参数,按照驾驶人视觉适应时间进行相应的回归拟合,其关系式为:l=4027t‑1.686    (12)依据式(12)建立函数模型求解驾驶人不同视觉适应时间下所需要的环境照度值,直至隧道入口,然后再根据 <mrow> <mi>d</mi> <mo>=</mo> <msubsup> <mo>&Integral;</mo> <mn>0</mn> <mi>t</mi> </msubsup> <mi>vtdt</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>vt</mi> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>式中,d为隧道纵深,t为视觉适应时间,v为车辆运行速度,确定隧道纵深d与驾驶人视觉适应时间t之间的关系;将式(12)代入式(13)得到公式(14),即可得到隧道纵深d、车速v与隧道优化照度l之间的关系:l=4027×(2d/v)‑0.843    (14)
地址 710064 陕西省西安市南二环中段