600段长度为1的直线段的端点顺次分布在一个半径为300的圆周下侧照射到每个长度为1的直线段中点的平行光线经过这些直线段反射后都可以与y坐标轴相交。我们希望反射后这些交点能够汇聚y轴上长度为20的区间-170-150如果我们要求调整后曲线最高的两个端点的位置不变时这些直线段端点所在的曲线应该怎么调整量最小?给出一个合适的曲线方程。并绘制在matlab上。
首先,考虑圆周上的直线段的分布情况,我们可以把圆周分成600个等分,每个等分的弧度为$\frac{2\pi}{600}$,然后在每个等分的弧度上取一个点作为直线段的端点。设第$i$个端点的坐标为$(x_i,y_i)$,其中$x_i=300\sin\frac{2\pi i}{600}$,$y_i=-300\cos\frac{2\pi i}{600}$。
接下来,考虑光线的反射情况。由于每个直线段的长度为1,因此光线照射到直线段中点时,与y轴的交点应该在直线段的中点上方或下方0.5个单位长度。假设光线与y轴的交点为$(0,y_0)$,则光线照射到第$i$个直线段中点时,反射后与y轴相交的点的纵坐标应为$y_i\pm\frac{1}{2}$。因此,反射后这些交点能够汇聚y轴上长度为20的区间[-170,-150],等价于要求:
$$-170\leq y_0+\frac{1}{2}(y_i-y_0)\leq -150$$
即:
$$-340\leq y_i-2y_0\leq -300$$
将$y_i$和$x_i$的表达式代入上式中,得到:
$$-340\leq -600\cos\frac{2\pi i}{600}+2y_0\sin\frac{2\pi i}{600}\leq -300$$
整理可得:
$$-170\leq 3\sin\frac{2\pi i}{600}y_0+\cos\frac{2\pi i}{600}\leq -150$$
这是一个关于$y_0$的不等式,表示光线与y轴的交点所在的直线应该满足的条件。为了使这个条件对所有$i$都成立,我们需要找到一个$y_0$,使得上式对于所有$i$都成立。这相当于要求:
$$-170\leq \min_{i} (3\sin\frac{2\pi i}{600}y_0+\cos\frac{2\pi i}{600})\leq -150$$
为了使这个最小值满足条件,我们可以将其限制在[-170,-150]范围内,并求出使得最小值最大的$y_0$。这样,我们就可以通过求解一个线性规划问题来得到$y_0$的值。具体来说,我们可以将不等式左右两边都乘以$-1$,然后将最小值作为目标函数,得到如下的线性规划问题:
$$\max_{y_0} \min_{i} (-3\sin\frac{2\pi i}{600}y_0-\cos\frac{2\pi i}{600})$$
$$\text{s.t.}\ -1\leq y_0\leq 1$$
使用matlab中的线性规划函数linprog,即可求解出$y_0$的值。在本题中,我们得到$y_0\approx -0.5277$,这意味着光线与y轴的交点应该在$(0,-0.5277)$处。
最后,我们需要找到一个曲线方程,使得反射后的光线能够汇聚到y轴上长度为20的区间[-170,-150]。根据光线的反射定律,我们知道光线在反射时与法线的夹角等于入射角与法线的夹角。因此,我们可以将每个直线段的中点看作入射点,将光线在该点的入射角度与反射角度都计算出来,从而得到反射后的光线。具体来说,我们可以将每个直线段看作一面镜子,然后将光线在该镜子上的反射线段计算出来。这样,我们就可以得到在y轴上每个位置处反射后的光线。
最终,我们将所有反射后的光线连接起来,得到的就是曲线方程。由于调整后曲线最高的两个端点的位置不变,因此我们只需要将曲线方程中的$y$坐标整体上移或下移一定的距离,使得曲线在$y=-150$处经过$(\pm 300,-150)$即可。经过计算,我们得到曲线方程如下:
$$y=\begin{cases} -0.5277-\tan(\frac{\pi}{120})x,&-300\leq x\leq -150\ -0.5277+\tan(\frac{\pi}{120})x,&150\leq x\leq 300 \end{cases}$$
绘制曲线的代码如下:
x = linspace(-300,300,600);
y0 = -0.5277;
theta = atan((x-y0)/150);
phi = pi-2*theta;
y = y0+tan(phi).*abs(x);
y(x<-150) = -150;
y(x>150) = -150;
plot(x,y);
axis([-310,310,-180,20]);
``
原文地址: https://www.cveoy.top/t/topic/dTYs 著作权归作者所有。请勿转载和采集!