建模题目：欢迎高手！大家进来讨论吧！

[KOBEBRYANT8]

车灯线光源的优化设计

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

[KOBEBRYANT8]

来吧，讨论！

[runnersun]

mark.

[xzzhq]

up

[joss]

算我一个数＞？！！

[longaway]

不懂

线光源的功率？怎么变化的？
业务内容太多

[liuyup]

up

[KOBEBRYANT8]

线光源的功率只是相对的！就是它的强度！用微分的方法可解！
第一道有点思路：
把旋转抛物曲面进行微分成每个小块，尽量均匀，减少误差！把每个点光源对小块曲面的反射再进行搜索、取舍！考虑物理上光的损失与距离有关后，再积分！
第二道题没固定的答案！自由发挥吧！
欢迎讨论第一道！！！！！！！！！！！

[btbuwangfei]

如果按照微分抛物面的方法，那也很不好搜索呀！你准备把抛物面分成多少份呢？
这也不好确定吧！这样想挺容易，，但做起来很难的

[ZeeGee]

呵呵，你所说的“证出是平面问题”具体指什么呢？

[KOBEBRYANT8]

是有的困难！但可以编程呀！

[KOBEBRYANT8]

我昨天晚9点就证明出是平面问题！
问题变的简单化了！
只要考虑在一个面上的光线取舍问题！编程选择、判断，结合光的强度，对B，C两点比较，找出与线光源的关系！OK！开始编程喽！

[KOBEBRYANT8]

还有10小时，大家一起努力吧！
先把求反射光的算法贴出！
#define L 4
#include "math.h"
main()
{int i=0,j;
double xm,ym,zm,zi;
float R[2700][3]={0};
double nx,ny,nz;
double I,N,t;
double Ix,Iy,Iz;


for(zm=0;zm<=5;zm=zm+1)
for(ym=0;ym<=5;ym=ym+1)
{ xm=(ym*ym+zm*zm)/60;
if (xm>21.6) continue;
for(zi=1;zi<=L;zi=zi+L/10)
{
I=sqrt((xm-15)*(xm-15)+(ym-0)*(ym-0)+(zm-zi)*(zm-zi));
Ix=(xm-15)/I; Iy=ym/I; Iz=(zm-zi)/I;
nx=-60; ny=2*ym; nz=2*zm;
N=sqrt(nx*nx+ny*ny+nz*nz);
nx=nx/N; ny=ny/N; nz=nz/N;
t=nx*Ix+ny*Iy+nz*Iz;
R[i][0]=Ix-2*t;
R[i][1]=Iy-2*t;
R[i][2]=Iz-2*t;
i++;
}

}
for(i=0;i<=500;i++)
for(j=0;j<=2;j++)
printf("\n%f\n",R[i][j]);


}

[KOBEBRYANT8]

欢迎提问！！

[zzitmf]

大哥,用VC能建模吗?