车灯线光源的优化设计.doc

车灯线光源的优化设计摘要本题主要以车灯线光源的优化设计为研究对象，在研究过程中建立了线性规化模型，给出算法并利用matlab软件求解。针对问题一，首先将线光源长度及其折射面离散化处理，将线光源的散射问题转换为光线间的折射问题，其次根据入射光线、法线和反射光线间的关系，推导出反射光线方程，从而可得反射光线与测试屏交点的表达式，然后规定测试屏上C、B两点周围10mm范围内均为其亮区，最后以满足B亮区光强度是A亮区光强度的两倍为约束条件，以线光源最短为目标，建立线性规划模型，给出算法并利用matlab求解，解得线光源长度最短为3.4mm，算法详见模型求解。针对问题二，首先将所求得的线光源长度及其折射面离散化处理，以一定步长分别逐步改变线光源长度及其折射点坐标，并利用问题一中反射光线与测试屏交点的表达式，得到每次改变后的交点坐标，记录下并画出其散点图，所得散点图即为车灯反射光亮区。对于问题三，车灯线光源设计规范的合理性主要体现在假设的科学性及量化标准的合理性上，应对其作出分析比较。关键词： 线性规划 离散化处理 1、问题重述问题背景：科学是第一生产力，设计经济适用的车灯需要融合光学、物理学等多方面的知识，基于节约能源和设计材料的原则，需要对车灯的线光源的长度进行优化设计。具体方案如下：明确车灯的构造考虑到安装在汽车头部的车灯的形状为一旋转抛物面，车灯的对称轴水平地指向正前方, 其开口半径36毫米，深度21.6毫米。经过车灯的焦点，在与对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。根据设计要求解决具体问题 该设计规范在简化后可描述如下。在焦点F正前方25米处的A点放置一测试屏，屏与FA垂直，用以测试车灯的反射光。在屏上过A点引出一条与地面相平行的直线，在该直线A点的同侧取B点和C点，使AC=2AB=2.6米。要求C点的光强度不小于某一额定值（可取为1个单位），B点的光强度不小于该额定值的两倍（只须考虑一次反射）。提出问题： 请解决下列问题：(1)在满足该设计规范的条件下，计算线光源长度，使线光源的功率最小。(2)对得到的线光源长度，在有标尺的坐标系中画出测试屏上反射光的亮区。(3)讨论该设计规范的合理性。2、问题分析题目中已给出了车灯的形状为抛物面，大小半径36毫米、深度21.6毫米，因此我们可以把车灯假设为三维空间坐标系中的几何图形进行分析，由以上给出的条件可以求出焦点在抛物面中的坐标位置，测试屏幕以及点A、B、C的坐标位置，据此可以对以下问题进行研究。 对于问题一，在焦点处放置一水平方向且与抛物面对称轴垂直的线光源，要求使得其满足C点的光强度不小于一个单位，B点的光强度不小于额定值的两倍。我们可以把光强度转化为反射光线数目，反射光线数目的多少就代表光强度的大小，反射到某一点的光线数目越多该点的光强度越大，B点的光强度为C点的两倍就相当于焦点处线光源照射到抛物面上一次反射到B点的光线数目是C点的两倍，在满足该条件的情况下，使线光源的功率最小，由发光功率公式得知线光源上的光分布均匀时，线光源的长度越短功率越小。对于问题二，在第一问的基础上我们知道了抛物面的表达式，算出了线光源的长度因此我们确定了线光源在坐标系中位置，根据空间解析几何的有关知识及光线的反射定律根据两点式得出反射光线的表达式，将题中给出的数据按原大小放在三维空间坐标坐标系中将线光源和抛物面离散化处理后将得到的反射光线在横坐标为25.015时，其余两坐标轴上的进行统计就可以得出亮区，将此数据做成平面图更加直观，具体图建下面题目求解中。对于问题三, 我们从测试屏位置的科学性、线光源位置的合理性、离散化处理的正确性、光强度近似量化的合理性四个方面对问题进行了讨论3、模型假设1、假设线光源可看做由离散的点组成的；2、假设不考虑光线在灯罩内发生二次折射的情况；3、假设光线经折射后到达B、C两处，光线直射的情况可忽略不计；4、假设光线在传播的过程中没有能量损失，即光线传播到光屏的距离不会影响到光能量的衰减；4、符号说明符号对符号的解释线光源的长度B位置的光线的条数C位置的光线的条数焦距的长度线光源分成等份离散化处理后每一份的长度抛物面与反射光线交点的坐标反射光线上由入射光线上一点得到你的对称点5、模型的建立与求解5.1构建优化模型计算线光源的长度5.1.1模型准备一 、基本问题的理解明确几个光学知识1、光线的传播路径在不经过反射的情况下，光线是沿直线传播的，其路径不会因外界环境发生改变。当光线经镜面反射后按一定角度折射射出。2、光线经反射后的位置关系光线的传播经反射屏后角度发生变化，入射光线与法线(垂直于反射屏的直线)的夹角与反射光线与法线的夹角相等。对线光源与灯罩的离散化处理结合光学知识，光的传播可看做是光线的传播；结合数学知识可将线近似看做很多点的集合，因此，我们将线光源离散化处理，为了便于研究，将线光源按一定标准分成很多段。同理可将灯罩看做很多点的集合，为了验证假设的合理性，可分别将线光源与灯罩分成不同数目的段研究。对B、C处光强度的理解由于在不同位置，经灯罩反射后到达测试屏上固定区域的光线条数是不同的，单位面积到达某一区域的光线数目的多少反映了光的强弱情况，结合实际情况，光线越密集光的强度越大，设一定区域内光线的条数为，则光强可近似表示为：由于要求B点的光强度不小于C点光强度的两倍，可将B、C附近面积为的相同区域作为研究对象，根据记录入射光线到反射屏上的点的个数，分别得到两区域中点的个数。将B、C区域看做直径为的圆形，则其表示区域的示意图如下：研究对象的坐标化处理为了定量地研究问题，需要在将线光源和灯罩离散化处理的基础上对其坐标化处理，考虑到安装在汽车头部的车灯的形状为一旋转抛物面，车灯的对称轴水平地指向正前方,将其放在空间直角坐标系中进行分析求解，在与车灯的对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。测试屏垂直于水平面。明确车灯内光线的传播路径 离散化的线光源发出的光一部分经过车灯内侧反射到达测试区域，一部分不经过反射直接到达测试区域，但是这样的光线相对于经过反射的光线数目是较少的，因此可忽略不计。二、目标的确定在满足设计规范的条件下，线光源的长度越大，在车灯正常工作时的功率越大，因此，线光源的长度与功率的大小是成正比的，所设计的线光源长度最小时其功率最小，综上所述，将功率最小这一目标转化为求解线光源长度的最小值，即：注：表示线光源离散化处理后得到的点光源个数；表示线光源分成等份离散化处理后每一份的长度四、约束条件的建立对B、C两点附近区域光强度的限制设计车灯线光源时要求B点的光强度不小于C点的光强度额定值的两倍，首先，对光强度合理近似量化，在单位面积S内光强度与照射在此区域内光线条数N有关，而B和C附近的区域面积可以看作相同的，因此，光强度和照射在B、C位置的光线的条数有关。根据光强度公式：设B、C点附近圆形区域的直径为，则其面积为，由于两者面积相同，可只考虑照射在B、C位置的光线的条数的倍数关系，设B位置的光线的条数为,C位置的光线的条数为，可将这一条件表示为： 对线光源长度的限制为了保证线光源位于车灯内部，线光源的长度不能大于旋转抛物面焦点处所对应的切面直径。结合车灯的开口半径和深度求得抛物线的曲线方程为：将其绕轴旋转得到旋转抛物面的方程为：设焦距的长度为，则切面半径的大小为，从而将线光源的长度限定为： 对离散化处理的点光源数目的约束线光源可以看作许多点光源的集合，由线光源分割得到的点光源的个数并不是越多越好，离散线光源的过程可近似看做单个点光源的移动，其移动经过的路径长度即为线光源的长度，将点光源的个数设为，仅需满足点光源个数的为整数的整数约束即可：五、反射光线方程的求解求解旋转曲面的法线方程建立坐标系求得抛物线的方程为，将抛物线沿轴旋转得到旋转抛物面方程为，设抛物面上一点，为法线上任意一点，根据旋转抛物面的方程得到法线的方向向量为，向量，根据直线的点向式方程得到曲面在已知点处的法线方程为： (1)为了得到反射光线的方程，需根据光线的反射定律，利用反射光线与入射光线对称的性质得到法线上点的坐标，示意图如下：【图示说明】上图中M为入射光线与反射面的交点，P点为入射光线上一点，过点P做法线的垂线，得到其与反射光线的交点Q，两点是关于法线对称的。根据抛物线的方程得到焦点坐标为,由于经过车灯的焦点，在与对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。设线光源的长度为，可求得线光源上点的坐标,设Q点的坐标为，根据对称性得到N的坐标。由公式(1)得到的参数方程将N点的坐标代入(1)式得到点向式方程，推得直线的参数方程为： (2)根据P、Q两点的连线与法线垂直得到： (3)将(2)式和抛物面的方程代入(3)式得： (4) 求得反射光线的方程已知抛物面与反射光线交点的坐标，结合反射光线上一点，由两点式方程得到反射光线的方程为： (5)将(4)式代入(5)式即可得到反射光线的方程。求解测试屏上点的坐标并记录根据测试屏的位置与距线光源的长度得到测试屏的表达式为：反射光线与测试屏的交点坐标表示为将其与反射光线的方程联立即可得到由线光源上的点经车灯反射后到达测试屏的点的坐标。分别记录折射点每次改变后与测试屏的交点是否在C、B两亮区范围内，从而记录下线光源每次改变长度后，折射在C、B两亮区内的光线条数。5.1.2模型的建立综上所述，建立以线光源长度最小为目标，B点的光强度不小于C点光强额定值的两倍、线光源长度不超过灯罩内径为为约束条件的最优化模型：目标函数：约束条件：【模型说明】目标函数满足条件的线光源长度为最小值约束条件(1)限定B点的光强度不小于C点光强度的两倍约束条件(2)保证线光源在车灯灯罩内约束条件(3)使线光源离散化处理的点光源的个数足够多约束条件(4)保证点光源的数目为整数5.1.3模型求解算法步骤：Step1. 首先将线光源离散化，以0.1mm为步长逐步增加线光源长度，其次将折射面离散化，分别以0.1mm，0.1度为步长，移动反射点的横纵坐标，从而将线光源的散射转化为点与点之间的反射问题。Step2. 将改变后的线光源长度及其反射点坐标代入所推导出的公式中，根据法线找出入射光线在反射光线上的对称点坐标，从而得到反射光线的直线方程，并求得反射光线与测试屏的交点即反射光线在测试屏上的亮点坐标。Step3. 分别以测试屏上C、B两点为圆心，以10mm为半径做圆，并规定反射在其内的光线均可表示C、B点亮度，并通过判断语句计算折射点每次改变后与测试屏的交点是否在C、B两亮区范围内，从而记录下线光源每次改变长度后，折射在C、B两亮区内的光线条数。Step4. 通过判断语句，从上述计算结果中找出B亮区光线条数是A亮区光线条数2倍以上的情况，并输出此时的线光源长度，从中找出最小值乘以2，即为所求结果。求解结果及分析：在求解过程中，采用“逐步逼近”法使线光源由原点沿z轴的正向和负向逐步变化，同时，求解出能够反射到达目标点的抛物面上的点，记录在该点反射时的值，并计算与值对应的光照度值，通过约束1和约束2的控制来选取满足题设的值，使功率最小。由编程实现，得到=3.4毫米。5.2绘制测试屏上反射光的亮区。5.2.1模型准备一、基本思路要绘制车灯反射光线的亮区，我们需要求解线光源上每个点的反射光线方程来确定这些光线所照射到的区域。，将题中给出的数据按原大小放在三维空间坐标坐标系中将线光源和抛物面离散化处理后将得到的反射光线在横坐标为25015毫米时，其余两坐标轴上的进行统计就可以得出亮区范围，用的作图功能得到亮区的示意图，具体图示见模型求解中。二、对求解用到公式的理解求解旋转曲面的法线方程设抛物面上一点，为法线上任意一点，根据旋转抛根据直线的点向式方程得到曲面在已知点处的法线方程为： (1)为了得到反射光线的方程，需根据光线的反射定律，利用反射光线与入射光线对称的性质得到法线上点的坐标，示意图如下：得到法线的参数方程由(1)式得到的点向式方程推得直线的参数方程为： (2)根据反射定律，反射光线与入射光线关于法线对称，因此入射光线上一点与对称点的连线同法线垂直，由这一关系得到： (3) 求得反射光线的方程已知抛物面与反射光线交点的坐标，结合反射光线上一点，由两点式方程得到反射光线的方程为： (4)将(3)式代入(4)式即可得到反射光线的方程。求解测试屏上点的坐标并记录根据测试屏的位置与距线光源的长度得到测试屏的表达式为：反射光线与测试屏的交点坐标表示为将其与反射光线的方程联立即可得到由线光源上的点经车灯反射后到达测试屏的点的坐标。分别记录折射点每次改变后与测试屏的交点是否在C、B两亮区范围内，从而记录下线光源每次改变长度后，折射在C、B两亮区内的光线条数。三、绘制亮区的基本方法结合求得的方程，将线光源的长度分解成n份，通过编程得到每一个点光源经过反射后能与反射屏相交的点的坐标，这些点的集合就是亮区的呈现形式，当点较密集时表示在这一区域内的光强度较高。5.2.2模型的建立本问题旨在求出在最小线光源下反射光线在测试屏上亮区，在将测试中各个因素坐标化处理的基础上，求得反射光线与测试屏的交点坐标表示为（单位为毫米），将其与反射光线联立可得到反射屏上的光点坐标：其中：【符号说明】表示该反射光线与抛物面的交点；反射光线上一点的坐标；为线光源的长度。【模型说明】根据空间解析几何的有关知识得出反射光线的表达式，将题给出数据按原大小放在三维空间坐标坐标系中将线光源和抛物面离散化处理后将得到的反射光线与测试屏的表达式联立得到测试屏上光点的坐标，将其记录下来并进行统计就可以得出亮区，具体图见下面模型求解中。5.2.3模型求解将线光源离散为点光源并求得每个点光源到达测试屏的光线数，由这些亮点构成的区域图即是测试屏上反射光的亮区。利用作图得到以下图形：结果分析：由上图可看出靠近车灯对称轴的区域内光强较强，越远离焦点所在平面光线越暗，这与实际情况中灯光照射的区域特点是一致的。5.3对设计规范合理性的分析一、测试屏位置的科学性题目中所给的设计规范是有一定物理意义的。首先，之所以把测试屏置于距离焦点25米处，是因为在25米的范围内，正常的车灯能够照亮足够大的区域使司机可以辨别行车路线、方向以及车辆前方的障碍物等。二、线光源位置的合理性设计车灯时将线光源置于焦点处，并且对称放置，我们要考察当线光源位置离开焦点时的情况，不妨将线光源沿光轴向两侧和上下移动，利用相同办法得到不同位置时测试屏上亮区的范围。可以发现，当光源位置沿主光轴方向偏离焦点时，中心亮度会急剧下降，这将导致司机在距离较远的正前方的人和物。因此，只有当光源对称地位于抛物面焦点处时，主亮度才会达到最大，才可以照亮一定范围内的所有区域。三、离散化处理的正确性 在解决问题的过程中，在对车灯的灯罩（即旋转抛物面）进行了离散化处理的同时将线光源进行离散化，线光源发出的光可以看做点光源的集合，通过离散化处理令线光源的长度从0开始累加处理即可得到线光源长度的最小值。离散点的集合形成面，因此对抛物面的离散化处理同样是正确的。四、光强度近似量化的合理性结合实际情况，在单位面积S内光强度与照射在此区域内光线条数n有关，因此我们可以将光强度量化为：光点个数的多少反映了光的强弱，因此这一近似量化是正确的。6、模型的评价与推广一、模型的优点将实际问题坐标化处理求解的得到反射光线的方程后运用了Matlab编程对模型进行了综合求解，简化了模型。模型具有较好的通用性，能够适应同类问题的各种变化，对不同大小的车灯都可求出其线光源功率最小时的长度，从而在满足条件的同时提高光线利用的效率。二、模型的缺点对线光源分割处理并累加的过程中，分割份数的不同可能会造成解出现一定的误差。对光强度的量化对光线的照射区域进行了等同化处理，这样可能在限定B、C处光强度倍数关系时出现一定的偏差。7附录问题一程序clc;clear;m=0;for x=0:1.2:21.6 for r=0:0.1:2*pi y=(60*x)(0.5)*cos(r); z=(60*x)(0.5)*sin(r); for d=-1.6:0.1:1.6 t=(30*x-y*d+450)/(30+2*x); y1=2*y-(y*t)/15-d; z1=2*z-(z*t)/15; x1=2*t+2*x-15; y2=(25015-x)*(y1-y)/(x1-x)+y; z2=(25015-x)*(z1-z)/(x1-x)+z; m=m+1; y3(m)=y2; z3(m)=z2; end endendplot(y3,z3,*)问题二程序clc;clear;m=0;n=0;l=0;k=0;l1=0;k1=0;for d=1.6:0.1:1.8 m1=0; n1=0; for x=0:0.1:21.6 for r=0:0.1:2*pi y=(60*x)(0.5)*sin(r); t=(30-y*d+450)/(30+2*x); y1=2*y-(y*t)/15-d; x1=2*t+2*x-15; f=(25015-x)*(y1-y)/(x1-x)+y; if(f2590&f1290&f=2) k(n)=d; k1(n)=n1; end end end endend2*l；k1经典婚庆主持词炮竹声声贺新婚,欢声笑语迎嘉宾.尊敬各位来宾,各位领导,各位亲朋好友,先生们,女士们,活泼可爱的小朋友们,大家好! 好歌好语好季节,好人好梦好姻缘.来宾们今天是公元*年*月*日(农历六月初八)是良辰吉日,在这大吉大利吉祥喜庆的日子里,我们怀着十二分的真诚的祝福相聚在*酒楼一楼婚宴大厅共同庆贺*先生与*小姐新婚典礼.(首先我们给予掌声的恭喜)大家都知道结婚是人生中的一件大事,而婚礼更是人生中最幸福神圣的时刻,尤其婚礼上浪漫温馨高雅别致的婚礼仪式以及亲朋好友的良好祝愿会给新人一生永远带来最美好的回忆.各位亲朋好友,我是本次婚礼庆典的主持人*.今天我十分荣幸地接受新郎新娘的重托,步入这神圣而庄重的婚礼殿堂为新郎*,新娘*的婚礼担任司仪之职.让我们在这里共同见证一对新人人生中最幸福神圣美好的一刻！真是:百鸟朝凤凤求凰,龙凤呈祥喜洋洋.让我们用掌声祝贺他们祝福新人凤凰展翅迎朝晖,恩爱鸳鸯比翼飞.携手同步知心人,共创宏图献真情.郎才女貌天作美,洞房花烛喜成双.在神圣的婚礼进行曲中一对新人手挽手,肩并肩缓缓步入婚礼大厅.脸上充满了无比幸福的笑容让我们用掌声与鲜花给予一对新人最诚挚的祝福.婚姻是人生大事,结婚典礼对青年男女来说是一生中最重要的时刻.你也笑,我也笑,亲朋好友齐来到.天也新,地也新,众星捧月迎新人.新郎新娘台上站,甜蜜感觉涌心间.风风雨雨牵手过,今天喜结美姻缘.亲朋好友齐相聚,欢欢喜喜来贺喜. * * * 天仙配,幸福的生活比蜜甜.在这个激动人心的美好时刻,作为婚庆司仪,首先请允许我代表新郎新娘以及新郎新娘的双方家长,对今天百忙当中来参加婚礼的各位来宾,各位亲朋好友的光临表示最诚挚的谢意和热烈的欢迎(谢谢大家)!欢迎你们!婚礼对每一个新婚的人而言,都是神圣,浪漫,唯美和经典的,随着神圣的婚礼进行曲奏响,英俊的新郎和美丽的新娘在掌声与祝福声中,缓缓的步上红地毯,那是万众瞩目的一瞬,那是梦寐以求的一瞬,那是凝结爱的万语千言的一瞬,那是最激动人心的一瞬,一同迷醉在尘世间最美妙的气氛里.爱情是古老而年轻的话题,也是不朽的人生主题.许多人已经拥有,更多人正