11

两个二函数图像对应的函数分别为：

L 2=0.002T 2−0.0253T +1.2725 (3.2) L 3=0.0007T 2+0.0108T +3.5227 (3.3)

根据上述，附件2与附件3直杆影长最短时刻分别是北京时间12点36分、北京时间12点48分，而当地时间都是正午12点整，所以通过式（2.3）计算，可以基本确定附件2与附件3直杆所在地的地理经度分别为东经79度、东经108度。

3.2.2直杆的地理纬度φ∗、杆长H 和日期的确定

本问题的优化模型与问题2相比，增加了变量日期，所以，将优化模型中遗

传算法的初始化步骤中变量的设置变为——设置3个变量（直杆长度H 、地理纬度φ∗与日期N ），确定好变量上下界（0<H ≤10，0°≤φ∗≤90°，0<N ≤356）；随机生成40个个体作为初始种群X (0)

然后以相似的方法运算，可以得出两个附件的三个变量的一些最优值，如表

2与表3。

3.3问题3模型的结果分析

表2 附件2中的数据通过遗传算法所得的结果

日期 经度 纬度

杆长

误差

地点 4月3日 79 15.9328 7.3314 0.0656 印度

4月2月

79 9.5015 7.6344 0.0036 印度 5月12日

79

6.7742

8.915

0.0849

孟加拉湾

图 7 附件2二次函数拟合

图 8附件3目标函数随迭代次数的变化图像