【数学建模】车灯线光源的优化设计之旋转抛物面反射问题求解

曲面的切平面与法线

(求梯度，对$x,y,z$分别求偏导)
1.曲面一般形式 : $F(x,y,z)=0$
算法 ： $gradF=(F_x,F_y,F_z)$
则有切平面 $F_x(x-x_0)+F_y(y-y_0)+F_z(z-z_0)=0$

法线$\frac{x-x_0}{F_x}=\frac{y-y_0}{F_y}=\frac{z-z_0}{F_z}$

2.曲面参数形式 $\{\begin{array}{l}x=X(u,v)\\ y=Y(u,v)\\ z=Z(u,v)\\ (u,v)\in D\end{array}$

---

空间曲线的切线和法平面

原理:切线就是割线无限逼近

$\Gamma =\{\begin{array}{l}x=x(t)\\ y=y(t)\\ z=z(t)\end{array}$

$M_0(x(t_0),y(t_0),z(t_0))$

$M(x(t_0+\Delta t),y(t_0+\Delta t),z(t_0+\Delta t))$

$\overrightarrow{M_{0}M}=(x(t_0+\Delta t)-x(t_0),y(t_0+\Delta t)-y(t_0),z(t_0+\Delta t)-z(t_0))$

$M_{0}M直线方程(点向式):\frac{x-x_0}{x(t_0+\Delta t)-x(t_0)}=\frac{y-y_0}{y(t_0+\Delta t)-y(t_0)}=\frac{z-z_0}{z(t_0+\Delta t)-z(t_0)}$

用$M$去逼近$M_0$，等价于 $\Delta t\to 0$

如把式子都除个$\Delta t$

式子变成$M_{0}M直线方程:\frac{x-x_0}{\frac{x(t_0+\Delta t)-x(t_0)}{\Delta t}}=\frac{y-y_0}{\frac{y(t_0+\Delta t)-y(t_0)}{\Delta t}}=\frac{z-z_0}{\frac{z(t_0+\Delta t)-z(t_0)}{\Delta t}}$

使$\Delta t\to 0$那么整个分母就变成了偏导

$M_{0}M直线方程L:\frac{x-x_0}{x^{\prime }(t_0)}=\frac{y-y_0}{y^{\prime }(t_0)}=\frac{z-z_0}{z^{\prime }(t_0)}$

空间平面的一般式表达式为：$Ax+By+Cz+D=0$,$D=0$时，该平面过原点
空间平面的点法式表达式为：$A(x-x_0)+B(y-y_0)+C(z-z_0)=0,(x_0,y_0,z_0)$是该平面上的点
包含了四个参数$A,B,C,D,(A,B,C)$是该平面的法向量。

所以法平面为$x^{\prime }(t_0)(x-x_0)=y^{\prime }(t_0)(y-y_0)=z^{\prime }(t_0)(z-z_0)$

---

空间曲面的切线和法平面

曲面方程 $\Sigma :F(x,y,z)=0$

点$M_0(x_0,y_0,z_0)$

任取经过$M_0$点的空间曲线$\Gamma =\{\begin{array}{l}x=x(t)\\ y=y(t)\\ z=z(t)\end{array}$

$M_0(x(t_0),y(t_0),z(t_0))$

一定会有$F(x(t),y(t),z(t))=0$
关于$t$求导后得到
$F_x\cdot x^{\prime }(t)+F_y\cdot y^{\prime }(t)+F_z\cdot z^{\prime }(t)=0$

带入$M_0(x_0,y_0,z_0)$
$F_x(x_0,y_0,z_0)\cdot x^{\prime }(t)+F_y(x_0,y_0,z_0)\cdot y^{\prime }(t)+F_z(x_0,y_0,z_0)\cdot z^{\prime }(t)=0$

用向量表示
$(F_x,F_y,F_z)\cdot (x^{\prime }(t_0),y^{\prime }(t_0),z^{\prime }(t_0))=0$

结论:
切平面
点$(x_0,y_0,z_0)$
法向量是$\vec{n}=(F_x(x_0,y_0,z_0),F_y(x_0,y_0,z_0),F_z(x_0,y_0,z_0))$
切平面方程为: $F_x(x_0,y_0,z_0)(x-x_0)+F_y(x_0,y_0,z_0)(y-y_0)+F_z(x_0,y_0,z_0)(z-z_0)$
法线为:
$\frac{x-x_0}{F_x}=\frac{y-y_0}{F_y}=\frac{z-z_0}{F_z}$

---

车灯线光源的优化设计

题目

A题 车灯线光源的优化设计

安装在汽车头部的车灯的形状为一旋转抛物面，车灯的对称轴水平地指向正前方, 其开口半径36毫米，深度21.6毫米。经过车灯的焦点，在与对称轴相垂直的水平方向，对称地放置一定长度的均匀分布的线光源。要求在某一设计规范标准下确定线光源的长度。
该设计规范在简化后可描述如下。在焦点F正前方25米处的A点放置一测试屏，屏与FA垂直，用以测试车灯的反射光。在屏上过A点引出一条与地面相平行的直线，在该直线A点的同侧取B点和C点，使AC=2AB=2.6米。要求C点的光强度不小于某一额定值（可取为1个单位），B点的光强度不小于该额定值的两倍（只须考虑一次反射）。

旋转抛物面求法

对于

$\{\begin{array}{l}y=x^2\\ z=0\end{array}$
绕$y$ 轴旋转
设旋转抛物面上一点为$A(x_1,y_1,z_1)$
设在$XOY$平面上的一点$B(x,y,z)$使得$y=y_1$

因为点$A$和点$B$距离$y$轴距离相等，通过两点间距离公式可得
$\sqrt{(x1-0{)}^2+(y_1-y_1{)}^2+(z_1-0{)}^2}=\sqrt{(x-0{)}^2+(y-y{)}^2+(z-0{)}^2}$

得到 $x_1^2+z_1^2=x^2+z^2$

因为$z=0$(点$B$在$XOY$平面上)
又因为$y=x^2$

所以$x^2+z^2=y$

如果$y=\frac{1}{2p}{x}^2$

那么$x^2+z^2=2py$

---

确定旋转抛物面方程

得到旋转抛物面方程为 $x^2+z^2=2py$
点$(0,0,0)$和点$(0,21.6,36)$在旋转抛物面上
带入得到$p=\frac{36^2}{2\ast 21.6}=30$

所以方程为$x^2+z^2=60y$

求反射到屏幕的点

假设线光源(焦点)$P(x_2,y_2,z_2)$射到抛物面上有某一反射点$Q(x_0,y_0,z_0)$
设反射光线上有一点$D(x_1,y_1,z_1)$
连接$PD$交经过反射点$Q$的切平面的法线于点$M(\frac{x_2+x_1}{2},\frac{y_2+y_1}{2},\frac{z_2+z_1}{2})$

将点$Q$带入分别对$x,y,z$求偏导，法线方程为

$\frac{x-x_0}{F_x}=\frac{y-y_0}{F_y}=\frac{z-z_0}{F_z}$

即法线为:
$\frac{x-x_0}{2x_0}=\frac{y-y_0}{-60}=\frac{z-z_0}{2z_0}$

点$M$在法线上
则有
$\frac{\frac{x_2+x_1}{2}-x_0}{2x_0}=\frac{\frac{y_2+y_1}{2}-y_0}{-60}=\frac{\frac{z_2+z_1}{2}-z_0}{2z_0}$

让公式等于$t$有
$\frac{\frac{x_2+x_1}{2}-x_0}{2x_0}=\frac{\frac{y_2+y_1}{2}-y_0}{-60}=\frac{\frac{z_2+z_1}{2}-z_0}{2z_0}=t$

解得
$\{\begin{array}{l}{x}_1=2x_0(2t+1)-x_2\\ y_1=-120t+2y_0-y_2\\ z_1=2z_0(2t+1)-z_2\end{array}$

根据$\overrightarrow{MQ}$与$\overrightarrow{PD}$垂直得到
$(\frac{x_2+x_1}{2}-x_0,\frac{y_2+y_1}{2}-y_0,\frac{z_2+z_1}{2}-z_0)\cdot (x_1-x_2,y_1-y_2,z_1-z_2)=0$

所以
$(\frac{x_2+2x_0(2t+1)-x_2}{2}-x_0,\frac{y_2+-120t+2y_0-y_2}{2}-y_0,\frac{z_2+2z_0(2t+1)-z_2}{2}-z_0)\cdot (2x_0(2t+1)-x_2-x_2,-120t+2y_0-y_2-y_2,2z_0(2t+1)-z_2-z_2)=0$

已知焦点的 $z_2=0,y_2=15$

但发现有$t^2$出现

其实$\overrightarrow{MQ}$就是法线向量，根据上面求法线为:
$\frac{x-x_0}{2x_0}=\frac{y-y_0}{-60}=\frac{z-z_0}{2z_0}$

得到法线向量为$\vec{n}=(2x_0,-60,2z_0)=(x_0,-30,z_0)$

那么关于垂直的等式就变为
$(x_0,-30,z_0)\cdot (x_1-x_2,y_1-y_2,z_1-z_2)=0$

即
$(x_0,-30,z_0)\cdot (2x_0(2t+1)-2x_2,-120t+2y_0-30,2z_0(2t+1))=0$

解得
$t=\frac{2x_0{x}_2-2x_0^2+60y_0-2z_0^2-900}{4x_0^2+4z_0^2+3600}$

所以解得点$D(x_1,y_1,z_1)$
$\{\begin{array}{l}{x}_1=2x_0(2t+1)-x_2\\ y_1=-120t+2y_0-y_2\\ z_1=2z_0(2t+1)-z_2\end{array}$

已知向量$\overrightarrow{QD}=(x_3,y_3,z_3)$，求向量$\overrightarrow{QD}$在$Y=25015$平面上的交点$G(x_4,y_4,z_4)$,其中
$\{\begin{array}{l}{x}_3=x_1-x_0\\ y_3=y_1-y_0\\ z_3=z_1-z_0\end{array}$

$QD:\frac{x-x_0}{x_3}=\frac{y-y_0}{y_3}=\frac{z-z_0}{z_3}$

带入点$G(x_4,25015,z_4)$即可求解