第 20卷 第 5期
2003年 03月
工 程 数 学 学 报
JOURNAL OF ENGINEERING MATHEMATICS
Vol .20 No.5
======================================================
Mar. 2003
文章编号:1005-3085(2003)05-0019-10
车灯线光源的优化设计
方案
气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载
薛 武, 杨铭和, 倪 冉
指导老师: 指导小组
(北京化工大学,北京 100029)
编者按:本文的主要特点是:在建立反射光强度的数学模型时,用反射点投影至测试屏目标点的变换的 Jacobi行列式来刻
画反射面元面积与该面元的反射光在测试屏照射到的微元面积之间的精确关系,并给出了相应的近似模型。
遗憾的是在实际离散化计算时,作者们做了过度的简化,未能真正应用上述模型,计算结果偏小。
摘 要:本问
题
快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题
是一个车灯线光源的优化设计问题。
首先,我们建立了一个连续的数学模型来描述这个优化问题,此模型研究了线光源上任意一点发出的光线
经过抛物面反射后到达光屏的情况,对于给定的考察点( B或 C)得到一个联系该点与光的发射点,反射点这 3
个点关系的方程组。由于光在传播过程中会有发散的现象,我们用 Jacobi行列式做了一个变换来描述这种散射。
根据 C点的光强度必须大于一个确定的值(文中设为 1),B点的光强度必须大于该值的 2倍的约束条件将这个
问题抽象成一个非线性规划问题。
由于解非线性规划问题是很复杂的过程,我们选择了将连续模型简化成一个离散模型。但是离散模型将光
离散成一条一条光线时,一般不能考虑光在传播过程中的散射问题。根据我们连续模型中用 Jacobi行列式算出
的结果知可以考虑光在传播过程中的散射问题,但是 Jacobi行列式是很难求出来的,为了解决这个困难我们将光
的散射用连续的方法做了一个简化。简化的方法是用向量投影的方法粗略作出了两个面积微元之间的关系,从
而得到了光线打在光屏上的散射效果与光线的起始单位方向向量和该光线经抛物面反射时的反射点的坐标的
关系。
运用以上的离散模型的算法,得到最优的线光源长度为 3.39mm。
关键词:线光源,非线性规划,反射,Jacobi行列式,向量投影。
分类号:AMS(2000)49K35 中图分类号:O224 文献标识码:A
1 问题的重述(略)
2 基本假设
1).不考虑光在空气中传播时的损耗。
2).点光源发光时,在各个方向上产生的光强度是一样的。
3 参数说明
P :旋转抛物面的焦参数
λ:线光源光能线密度,可以表示单位长度光源产生的光强度量
2R :线光源的长度
W :线光源半轴长发光的总强度量
I(B):B 点接收到的光强度(光强度量的面密度)
I(C):C 点接收到的光强度(光强度量的面密度)
d :旋转抛物面的深度
I(h):光在空间传播时在距离光源 h 处的光强度
Q :某点光源的发光强度量
n :每个发光点向空间发射的光线数
Δr :每个离散点光源的长度
Δm :每束光线可以产生的光强度量
E(P):抛物面上可以将光线反射到 P点的所有点的集合
O :光屏上每个小网格的面积
!v :测试屏的单位法向量(沿 z轴负方向)
4 问题的分析
本问题是一个面镜反射的几何光学问题。目的是要求在满足题目要求( I(B)" 2,
I(C)" 1)时对线光源的长度( 2R )和单位发光密度(λ)进行优化,使得线光源的发光总
强度量最小。
对于其中的考察点,例如 B ,首先要考虑光源上哪些点可以经过抛物面哪些反射点反
射到考察点 B 。然后根据他们之间的位置关系算出从线光源上某点射出的光线在到达考
察点 B 时,考察点 B 所接收到的光的强度。再对该考察点上所接收到的所有光的强度进行
积分,算出 I(B),将 I(B)" 2,作为问题的一个约束条件,同理算出关于 C 点的一个约束
条件。问题就变成了一个带约束条件的优化问题。
图 1 光反射路线
5 模型的构建
5 .1 反射点和发光点以及受光点之间的关系
考虑线光源上一点 H(0,r,p /2)发射一条光
线经过抛物面上一点 G(x1,y1,z1)射到光屏上
某点 P(x0,y0,z0)的情况。
如图 1,抛物面方程为
x2 + y2 = 2pz
G 点的单位内法向量为
{#lG}=
- x1,- y1,{ }p
x21 + y21 + p$ 2
将向量 #GP 和 #GH 单位化得
e#GP =
x0 - x1,y0 - y1,z0 - z{ }1
(x1 - x0)2 +(y1 - y0)2 +(z1 - z0)$ 2
02 工 程 数 学 学 报 第 20
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
卷
e→GH =
- x1,r - y1,
p
2 - z{ }1
x21 +( r - y1)2 +
p
2 - z( )1ヘ 2
若光线可以从 H出发经过抛物面上的G到达光屏上的P ,则向量 e→GP和 e→GH的和与G点
的单位内法向量同向,即向量 e→GP 和 e→GH 的和与 G 点的单位内法向量的向量积为 0,即
e→GP + e→( )GH × l→G = 0
令α =
→GP
→GH
,可以将原方程组化简为
x1(y0 +αr)- x0 y1 = 0
x1 z0 +α
p
2 -(1 +α)( z1 + p[ ]) + px0 = 0
y1 z0 +α
p
2 -(1 +α)( z1 + p[ ]) + p(y0 +αr)= 0
z1 ≤ d
x21 + y21 = 2pz
1
(*)
由于方程组(*)中前三个方程只有两个是独立的,则在前两个方程中可以任意去掉一
个,方程组就决定 H 、P 和 G 三点之间的位置关系。
5 .2 光线经光屏反射过程的数学描述
在线光源上取一点 H(0,r,p /2),过该点沿单位向量
u→1 = sinαcosβ,sinαsinβ,cos{ }α
做一条直线交抛物面于 G(x1,y1,z1),则
x1 - 0
sinαcosβ
=
y1 - r
sinαsinβ
=
z1 -
p
2
cosα =
t
又 G 在抛物面上,则
x21 + y21 = 2pz1
⇒ t =
pcosα - y1 sinαsinβ + p
2 cos2α - rpsinβsin(2α)- sin
2αcos2βr
2 + p2 sin2αsin2ヘ β
sin2α
则
图 2 向量加法示意
x1 = tsinαcosβ
y1 = r + tsinαsinβ
z1 =
p
2 + tcos
α
G 点的单位内法向量为
l→G =
- x1,- y1,{ }p
x21 + y21 + pヘ 2
则表示反射光线方向的单位向量 u→2 如图 2,
12第 5
======================================================
期 车灯线光源的优化设计方案
可知
- u→1 + u→2 = 2 l→G cosθ
又 cosθ =
- u→1· l→G
- u→1 l→G
= - u→1· l→G
则经抛物面反射的单位向量 u→2 = u→1 - 2 u→1· l→( )G l→G
5 .3 光在传播过程中的发散
光在空间传播时是在整个空间上均匀发散的。我们用光强度量来表示光源在整个空间
发射出去的光能,用光强度量的面密度来度量光发散到某个点在该点产生的光强度。当发
光强度量为 Q 的点光源发出的光传播到空间某点时,由于光波是球面波,则可知
I(h)= Q4πh2
那么当长度为 d r 的线光源(当成点光源看)从 H 发出的光传播到 G 时,G 点所接收到
的光强度为
d I(G)=
λd r
4π →GH 2
现在讨论光从 G 反射到 P 点时光强度的变化情况。在 G 周围取一个小的面积微元
dS = 1 + z21x1 + z
2
1yヘ 1 dx1 dy1
从方程组(*)中可以确定
x1 = x1(x0,y0)
y1 = y1(x0,y0{ )
则 x1,y1 对于 x0,y0 的 Jacobi行列式
J =
∂x1
∂x0
∂y1
∂x0
∂x1
∂y0
∂y1
∂y0
则可得
dx1 dy1 = J dx0 dy0
而 dx0 dy0 为光屏上的面积微元,记为 dS' ,即
dS' = dx0 dy0
则可以推出
dS
dS' = 1 + z
2
1x1
+ z21yヘ 1 J
G 点接收到的光强度量
dm(G)= - u→1· l→( )G d I(G)dS (取负号是由于u→1 和 l→G 夹角超过 90º)
又由于 dS 上的所有光线全反射到 dS' ,即它们两个面积微元上的光强度通量相等
dm(G)= dm(P)= d I(P)dS'
⇒d I(P)=
- u→1· l→( )G dS
dS' d I(G)
22 工 程 数 学 学 报 第 20
======================================================
卷
= - u→1· l→( )G 1 + z21x1 + z21yヘ 1 J d I(G)
=
- u→1· l→( )G 1 + z21x1 + z21yヘ 1 Jλd r
4π →GH 2
则 P 点的光强度
I(P)=∫
R
- R
(P)
- u→1· l→( )G 1 + z21x1 + z21yヘ 1 Jλd r
4π →GH 2
5 .4 非线性规划问题
通过以上的分析可以知道 B 点的光强度为
I(B)=∫
R
- R
(B)
- u→1· l→( )G 1 + z21x1 + z21yヘ 1 Jλd r
4π →GH 2
同理,C 点的光强度为
I(C)=∫
R
- R
(C)
- u→1· l→( )G 1 + z21x1 + z
2
1yヘ 1 Jλd r
4π →GH 2
根据题意,本问题的约束条件为
I(B)≥ 2
I(C)≥
{ 1
我们的目标是使线光源半轴发光的总强度量最小,即
minW =λR
那么本问题可以归结为一个求非线性规划最优解的问题
minW =λR
s . t
I(B)=∫
R
- R
(B)
- u→1· l→( )G 1 + z21x1 + z
2
1yヘ 1 Jλd r
4π →GH 2
≥ 2
I(C)=∫
R
- R
(C)
- u→1· l→( )G 1 + z21x1 + z
2
1yヘ 1 Jλd r
4π →GH 2
≥
1
6 模型的简化与求解
6 .1 光源发光的离散化
以上的模型最终是要求解一个非线性规划问题,解这类问题通常需要很大的计算量,而
且很容易陷入局部极小值。现在我们将模型中连续的线光源离散成许多点光源,而且将光
的均匀发散离散化,考虑点光源的光以 n 条光线的形式对外发散。将光屏也分割成许多小
网格,将经过抛物面反射到屏上各个网格内的光所具有的光强度量进行求和再除以该网格
的面积,定义为该网格内的光强。
若每个离散点光源的长度为Δr ,则该离散点发出的一条光线可以产生的光强度量
Δm =λΔ
r
n
32第 5
======================================================
期 车灯线光源的优化设计方案
6 .2 光传播过程中发散的一个近似处理
对于光在传播过程中的散射对射到光屏上后强度的影响模型中是用 Jacobi行列式体现
的。但是这个 Jacobi行列式的表达式是很复杂的,我们给出一种近似处理这种散射影响的
方法。这个近似的简化是在连续条件下考察面积微元得到的,但是这个结果对于离散的模
型同样适用。
如图 3,一束光线射到曲面上的一个单位法向量为 l→G 的面积微元 S1 后反射到一个单位
法向量为"v 的平面上形成一个面积微元 S2 。由于面积微元很小,我们近似的认为反射光线
沿一个方向射出。
图 3 面积微元投影
我们进行反向思考,由于 S2 是 S1 沿单位向量 u→2 投
影后在单位法向量为"v 的平面上的投影,那么 S2 和 S1
在 u→2 方向上的投影面积应该是相等的,即
S1 l→G·u→( )2 = S2 u→2·"v = - S2 u→2·"( )v (取负
号是由于u→2 和"v 夹角超过 90º) S1 面上接收到的光
强度量
Δm(S1)= - u→1· l→( )G ΔI(S1)S1
又由于 S2 和 S1 上通过的光强度通量相等,即
Δm(S1)= Δm(S2)= ΔI(S2)S2
#ΔI(S2)=
u→1· l→( )G S1
S2 Δ
I(S1)=
u→1· l→( )G u→2·"( )v
l→G·u→( )2
ΔI(S1)
光线从线光源上的 H 到达抛物面上的 G 点时的光强度根据[5 . 3光在传播过程中的发
散]可以得到,光从 G 传播至 P 时运用以上方法可以得到
ΔI(P)=
u→1· l→( )G u→2·"( )v
l→G·u→( )2 ΔI
(G)
=λ
u→1· l→( )G u→2·"( )v Δr
4π →GH 2 l→G·u→( )2
=λ
u→1· l→( )G u→1 - 2 u→1· l→( )G l→[ ]G ·"{ }v Δr
4π →GH 2 l→G· u→1 - 2 u→1· l→( )G l→[ ]{ }G
这时每条光线从 H 射到 P 时在 P 所在网格产生的光强度量
Δm(P)= ΔI(P)
4π →GH 2
n =
λ u→1· l→( )G u→1 - 2 u→1· l→( )G l→[ ]G ·"{ }v Δr
n l→G· u→1 - 2 u→1· l→( )G l→[ ]{ }G
那么 P 点的光强度为
I(P)=
$
G%E(P)
Δm(P)
O =
$
G%E(P)
λ u→1· l→( )G u→1 - 2 u→1· l→( )G l→[ ]G ·"{ }v Δr
n l→G· u→1 - 2 u→1· l→( )G l→[ ]{ }{ }G
O
6 .3 问题的求解
根据以上对模型的离散简化得到一个离散优化模型
42 工 程 数 学 学 报 第 20
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
卷
minW =λR
s . t . I(B)=
Σ
G∈E(B)
λ u→1· l→( )G u→1 - 2 u→1· l→( )G l→[ ]G ·${ }v Δr
n l→G· u→1 - 2 u→1· l→( )G l→[ ]{ }{ }G
O % 2
I(C)=
Σ
G∈E(C)
λ u→1· l→( )G u→1 - 2 u→1· l→( )G l→[ ]G ·${ }v Δr
n l→G· u→1 - 2 u→1· l→( )G l→[ ]{ }{ }G
O % 1
为了方便用计算机求解,我们对上面的优化模型进行转化,得到
min W = T·R
s . t . I(B)= T·NB % 2
I(C)= T·NC % 1
其中,T 为每条光线带的光强,NB 为打到以 B 为圆心,以 2cm为半径的圆内的光线条数,
NC 为打到以 C 为圆心,以 2cm为半径的圆内的光线条数,对于这个问题的求解我们采取分
步等步长搜索的方法。具体步骤如下:
步骤一 将线光源细化成许多个小的光源,每个小光源的长度Δr = 0.01mm。将整个
空间的立体角 4π分成 n = 40000份。在 B 和 C 点周围分别做一个半径为 2cm的圆。
步骤二 将每个小光源发射出的光离散为以 40000条光线的形式发射出去,每条可以
射到抛物面上的光线经抛物面反射后会在光屏上打下唯一的一个点。如果该点在以 B(或
C )为中心的圆内,就将打在这个区域的光线条数加 1。
步骤三 首先在 R ∈(0,30mm)内以为 0 .1mm步长对 R 进行搜索,得到的 W - R 的图
象如下图 4,对图中划在圆圈内部的 W - R 曲线进行放大,如图 5。
图 4 W-R关系图 图 5 W-R关系局部极大图
我们从图中可以看出最优的 R出现在 1 .5mm到 3mm区域内。现将进一步搜索的区域
定为 R ∈[1 .5mm,3mm]。
步骤四 在区间 R∈[1 .5mm,3mm]内以 0 .005mm为步长进行搜索,得到和关于线光源
半轴长 R 图形(如图 6)。
利用 NB 、NC 的值,对优化模型进行求解,得到在 R ∈[1 . 5mm,3mm]区间内 W - R 图
52第 5
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
期 车灯线光源的优化设计方案
线,如图 7。
图 6 线光源半径 R图形 图 7 W-R关系图
在图中圆圈(即最优化解)附近点 W - R - T的数据列表如下
R(mm) 1.670 1.675 1.680 1.685 1.690 1.695
T 0.00016145 0.00016072 0.00016 0.00015941 0.00015870 0.00015811
W(e - 004) 2.696215 2.69206 2.688 2.6860585 2.68203 2.6719645
R(mm) 1.700 1.705 1.710 1.715 1.720
T 0.00015777 0.00015723 0.00015679 0.00015632 0.00015586
W(e - 004) 2.68209 2.6807715 2.681109 2.680888 2.680792
从优化的结果数据,我们得出最优的线光源的半轴长为 R = 1. 695mm。此时这个线光
源的长度 2R = 3.39mm。
7 问题第二问的解答(略)
8
设计规范
民用建筑抗震设计规范配电网设计规范10kv变电所设计规范220kv变电站通用竖流式沉淀池设计
的合理性验证
8 .1 合理性的分析
本题中涉及到的灯是车的前照灯。前照灯是汽车在夜间行车或在能见度较低的条件
下,为驾驶员提供车道照明的重要设备,所以前照灯必须有足够的发光强度和正确照明区
域。不合格的前照灯会使驾驶员对前方道路情况辨认不清,或者在与来车交会时造成对方
驾驶员目眩等,从而导致事故的发生。
一个合理设计
规范
编程规范下载gsp规范下载钢格栅规范下载警徽规范下载建设厅规范下载
应该是在这个规范下做出的最合适的车灯具有以下照明效果:
1)车灯正前方为了照亮车道,应该亮度很大。
2)前照灯也要照亮邻道,使驾驶员能清楚把握邻道的情况,故车两侧附近也应该有较
高的亮度。
3)而更远一点的距离,如果太亮的话,就有可能使迎面而来车辆的驾驶员目炫,而使发
生事故的机率加大,故亮度要较低。
8 .2 设计规范的随机模拟验证
根据以上的合理性分析,我们需要用实验的方法检验这个规范的合理性。我们在以上
做的模型中要求点光源在空间各个方向上产生的光强度是一样的,但是考虑到光在实际传
62 工 程 数 学 学 报 第 20
======================================================
卷
播过程中,由于种种因素的影响,会使光源在发光时在各个方向上发光的强度是不一样的。
我们用计算机对线光源经过抛物面反射发光进行随机模拟实验。具体步骤如下:
步骤一 用 Matlab随机数函数 UNIFRND( )产生点光源向空间发射光线的随机单位方
向向量。
步骤二 用这些随机单位方向向量对上文已求出的最佳光源发光进行模拟,在离散成
许多小区域组成的光屏上统计每个区域内的光强度。在光屏上作出光强度分布图 8。
图 8 光强度分布图
上图中,光屏上某点的光强度用高度来体现,高度越大,亮度越大。从上图可以很明显的看
出来灯的正前方光的强度很大,而且两侧光的强度和正前方相比很小,但是和其他区域相比还是
比较强的,这就从实验的角度验证了这个设计规范的合理性。以下作出该图象的俯视图(图 9)。
图 9 上图的俯视图
72第 5
======================================================
期 车灯线光源的优化设计方案
从这个图中就可以更直接的看出光屏上各区域光强度的分布。中间部分的很亮可以为
夜间行车提供良好照明,一个车宽附近的区域也较亮,可以清楚看到邻道的情况,更远一点
的距离亮度迅速减弱,这样与来车交会时不易造成对方驾驶员目眩,从而有效地避免事故的
发生。
9 模型优缺点分析
对于本题可以建立连续和离散两种模型。前者可以将本题抽象成一个非线性规划问
题,但是求解的过程比较复杂,很容易陷入局部极小值的情况,而且对于问题第二问的求解
没什么帮助,但是优点是如果初值取的好可以很精确的求出本题的结果。后者的优点是对
于算法而言不是很复杂,求解结果很快,同时可以很方便的解决问题第二问,但是求解得到
的结果不是很精确,而且只能给出数值解,并且将光离散成光线以后就不能体现光在传播过
程中的散射。综合以上考虑,我们提出的这个模型既拥有离散模型的优点,又运用一个连续
模型的变换结果,考虑了光在传播过程中的散射,比较精确地刻画了光在传播过程中的散
射。
参考文献:
[1] 方得植编 .《解析几何》[M]. 北京:高等教育出版社,1986
[2] 张以谟编 .《应用光学》(修订本)[M].天津:机械工业出版社,1988
[3] 张凤林,孙学珠编 .《
工程
路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理
光学》[M]. 天津:天津大学出版社,1988
[4] 刘则毅编 .《科学计算技术与 Matlab》刘则毅[M].北京:科学出版社,2001
[5] 章学诚编 .《解析几何》[M]. 北京:北京大学出版社,1989
[6] 中国汽车检测维修网:http://www. auto - tester . com
An Optimal Design of Line Light Source
XUE Wu, YANG Ming-he, NI Ran
Advisor: Mathematical Modeling Tutor Group
(Beijing University of Chemical Technology,Beijing 100029)
Abstract:Firstly,we developed a continuous mathematic model to characterize this problem. It demonstrated how the beam reaches the
screen by the reflection of the paraboloid. It was a nonlinear programming problem. Because it was very difficult to be solved,we simpli-
fied it into a discrete model . Based on this model,we wrote a Matlab program and got the optimal length of the line light source was 3 .
39mm. Finally,we applied random simulation to verify that the design scheme was reasonable.
Keyword:line light source,Jacobi determinant,nonlinear programming,reflection,projection
82 工 程 数 学 学 报 第 20
======================================================
卷