1、- 1 -光照强度计算的优化模型摘 要文章把整个教室看成若干个点的集合，每一个点的光照强度为四个点光源照到该点的光照强度之和.用离散化的思想，以距离作为重要的桥梁连接目标点与光源和光照强度之间的关系： ( 是四点光源的光强).算PQJk2PQJ,是 目 标 点 的 光 照 强 度出任意一点的光照强度和距离地面 1 米处各点的光强，用 Maple8 画出位置与光照强度之间的曲面图,还作出密度图和等高线图进一步说明模型的合理性.对于第三问和第四问，重新考虑四个光源位置.根据在 1 米处的光照强度总和越大、光照越平均，学生的满意程度就越大的思想，用 MATLAB6 求出 1 米处的光照强度总和最大值
2、 M=0.1278 和光照最平均值 E=0.6 ，用最小二乘法找出一个最贴近 M 和 E310的值，其中目标函数关系式是 ，找出目标函数的最小值，确22)()( EJMLQw定四个光源的位置为（(80,47,25)、(80,103,25)、(40,103,25)、(40,44,25)）, 同理也能确定线光源的位置.关键词:光照强度;点光源;线光源;满意程度 - 2 -1 问题的提出现在已知一个教室长为 15 米，宽为 12 米，在距离地面高 2.5 米的位置均匀的安放 4 个光源（分别为 1、2、3、4） ，各个光源的光照强度均为一个单位，如下图要求：（1） 如何计算教室内任意一点的光照强度？
3、（光源对目标点的光照强度与该光源到目标点距离的平方成反比，与该光源的强度成正比）.（2） 画出距离地面 1 米处各个点的光照强度与位置（横纵坐标）之间的函数关系曲面图，试同时给出一个近似的函数关系式.（3） 假设离地面 1 米高正是学生桌面的高度，如何设计这四个点光源的位置，才能使学生对光照的平均满意度达到最高？（4） 若将题目中的点光源换成线光源，以上(2)、(3)问的结果又如何？（对于（1） 、 （2）问，假设横向（纵向）墙壁与光源、光源与光源、光源与墙壁之间的距离是相等的.）2 模型的假设21 光不会通过窗、门等外涉，也不考虑光在空气中的消耗，即光照强度和不变；2. 2 室内不受外界光源
4、影响；23 教室高度为 2.5 米；24 不考虑光的反射；25 线光源发光是均匀的.3 符号约定： 教室内任意一点Q： 四光源对点 的光照强度 JQ: 光源 P 的光照强度 (P=1,2,3,4) P： 光源 与点 之间的距离 (P=1,2,3,4)Q: 四个点光源位置 (P=1,2,3,4)Pw- 3 -: 1 米处各点的最平均光照强度E: 光源在 处时，1 米处各点的光源强度的总和wL： 1 米处各点光源强度总和的最大值M: 线光源的长度4 模型的建立由于光源对目标点的光照强度与该光源到目标点的距离的平方成反比，与该光源的强度成正比，故可作出任意一点光照强度函数：(1)412ppQJPkJ
5、其中 k 为正常数，建立右手直角坐标系，如图 1 所示，1.2.3.4 分别为四个光源的位置，由于横向（纵向）墙壁与光源、光源与光源、光源与墙壁之间的距离是相等的，可以得到此四光源在坐标系中的位置为: (单位:分米)1 (80, 50, 25)2 (80, 100, 25)3 (40, 100, 25)4 (40, 50, 25)由于各个光源的光照强度均为一个单位，且光都不外涉，不在空气中消耗，也不受外界光源影响，即所有目标点的光照强度之和等于四个点源的光强之和，故得： 41241pPQp JkJ 205 222222 )5()10()8()5()0()8(xyz zyxzyxk ( 4541
6、)4( 222222 zyxzyz25xy1501204 3 1 2图 1- 4 -由（2）式解得：k=0.01855612043故（1）式可化为： )3(QP04.185624ppQJJ5 模型的求解5.1 计算教室内任意一点的光照强度:把任意一点 Q 代入（3）式，可得教室内任意一点的光照强度：422232221 )5()0(x)-(43.861)5()10(x)-(44.86 JzyJzyJQ 5.2 距离地面 1 米处各个点的光照强度与位置（横纵坐标）之间的函数关系曲面图:根据函数关系式： 422322 21 5)50()4(3.8615)10()4(4.856 JyxJyxJQ 用
7、MAPLE8 软件作出相应的三维曲面图，如图 2 所示：用 MAPLE8 软件作出距离地面 1 米处各个点光照强度的密度图和等高线图，如图 3所示：（其中密度图能比较真实地反映四光源到离地面 1 米的水平面的亮暗程度，而等高线图也很好地反映在此平面上光照的分布.）图 2- 5 -（密度图） （等高线图）5.3 使学生对光照平均满意度最高的四个点光源的位置学生桌面所在平面的、越大,学生的满意度就越高, 要尽可能满足以上两个条件,以确定使学生满意度最大的光源位置.题目中没有关于优先权及权重的规定，可以设总光照强度的权重与光照平均的权重相同，要使在高为 1 米处光照越大就要使目标点尽可能地靠近光强总
8、和最大的位置,同理使光照越平均就要使目标点越靠近使光照最平均的位置,用最小二乘法得出它们的差值最小的位置即为所求. 22)()(min EJMLQws.t 80max- 9 -max=l;x11=x1;y11=y1;x12=x2;y12=y2;x13=x3;y13=y3;x14=x4;y14=y4;endp=l./(120.*150);Q=0;for x=0:0.1:12for y=0:0.1:15Q=Q+(k.*(1./(x1-x).2+(y1-y).2+(2.5-z).2)+1./(x2-x).2+(y2-y).2+(2.5-z).2)+1./(x3-x).2+(y3-y).2+(2.5-
9、z).2)+1./(x4-x).2+(y4-y).2+(2.5-z).2)-p).2.(1./2);endendif minQmin=Q;x21=x1;y21=y1;x22=x2;y22=y2;x23=x3;y23=y3;x24=x4;y24=y4;endendenddisp(最大值,x11=,num2str(x11), ,y11=,num2str(y11), ,x12=,num2str(x12), ,y12=,num2str(y12), ,x13=,num2str(x13), ,y13=,num2str(y13), ,x14=,num2str(x14), ,y14=,num2str(y14)
10、disp(最平均,x21=,num2str(x21), ,y21=,num2str(y21), ,x22=,num2str(x22), ,y22=,num2str(y22), ,x23=,num2str(x23), ,y23=,num2str(y23), ,x24=,num2str(x24), ,y24=,num2str(y24)附录二：clearclcmax=0;min=4;li=4;for i=0:0.1:3for j=0.1:0.1:4s=0;x1=8+i,y1=5-jx2=8+i,y2=10+jx3=4-i,y3=10+jx4=4-i,y4=5-j for x=0:0.1:12for
11、y=0:0.1:15for z=0:0.1:2.5if x1=x endend- 10 -endendk=4./s;l=0;z=1;e=0for x=0:0.1:12for y=0:0.1:15l=l+k.*(1./(x1-x).2+(y1-y).2+(2.5-z).2)+1/(x2-x).2+(y2-y).2+(2.5-z).2)+1./(x3-x).2+(y3-y).2+(2.5-z).2)+1./(x4-x).2+(y4-y).2+(2.5-z).2);r=k.*(1./(x1-x).2+(y1-y).2+(2.5-z).2)+1/(x2-x).2+(y2-y).2+(2.5-z).2)
12、+1./(x3-x).2+(y3-y).2+(2.5-z).2)+1./(x4-x).2+(y4-y).2+(2.5-z).2);e=e+(r-6*10(-32)2;endendS=(l-0.1278)2+eif Slili=Sx11=x1,y11=y1, x12=x2,y12=y2, x13=x3,y13=y3, x14=x4,y14=y4,en4en4en4disp(x11=,num2str(x11), ,y11=,num2str(y11), ,x12=,num2str(x12), ,y12=,num2str(y12), ,x13=,num2str(x13), ,y13=,num2str(y13), ,x14=,num2str(x14), ,y14=,num2str(y14)li
