标题 | 基于分治法的遍历搜索算法和粒子群算法的改进对太阳影子坐标的地理定位和事件确定的研究 |
范文 | 黄鸿基++钱圳冰++冯帆++周行洲 [摘要]研究首先分析了经度、纬度、杆子长度和日期这几个变量的关系,通过太阳影子长度计算公式,结合误差控制的相关理论,我们定义了最小误差函数,并建立了一个目标规划模型。对于多参数的情况,研究提出“基于分治法的遍历搜索算法”,用于求解出影子的地理位置和日期。后研究对在目标规划模型进行了合理的改进,建立了修正过的目标规划模型。通过对所收集的视频进行合理分析,以每三分钟为一个采样周期,得到了影长与时间的相关数据。同时,考虑到数据量大和算法复杂度高的因素,研究引入了粒子群(PSO)算法,并用遗传算法(GA)对其进行了改进,形成了GA-PSO算法。 [关键词]粒子群算法;遍历搜索;最小误差函数 [DOI]1013939/jcnkizgsc201718199 1引言 太阳影子的位置和时间的确定,在地理勘测和工程应用中有很高的应用价值。研究首先根据某固定直杆在水平地面上的太阳影子顶点静态坐标数据,建立数学模型确定直杆所处的地点和日期。后通过太阳影子的动态变化视频,判断该视频所发生的大致地点。从而推广这两种算法作为确定太阳影子发生地点和时间的方法。 2基于分治法的遍历搜索算法确定静态坐标下太阳影子发生地点和时间 21算法分析 研究根据相关数据给出的太阳影子顶点坐标数据,建立合理的数学模型来确定直杆所处的地点和日期。根据附件提供的坐标和其他数据,我们可以得到不同时间点所对应的太阳影子长度。很明显,这是一个目标规划问题。可以基于太阳影子长度公示,建立寻找地点和日期的规划模型。同时,考虑到研究涉及的未知量较多,直接求解很难得到最优的结果。因此,考虑纬度、经度和杆长等因素,研究采用先模块搜索后整体遍历的思想,建立一个新的“基于分治法的遍历算法”,寻找到准确的地理位置并确定它所对应的日期。 22算法模型的建立 步骤一:数据处理 定义四组变量α、β、h、N,分别代表经度、纬度、杆子高度和年份。由数学分析的相关理论和地理学的相关知识可知,一组连续的变量可以看成一组间隔无限小的离散型随机变量的线性组合。同时,在地球上,当两个经度之间相差1°时,它们之间的时间相差4分钟,因此,我们对数据进行了离散化处理,并且在误差范围内完全可以认为不会对模型的精确度造成影响。我们得到以下关系: 纬度:-90°≤α≤90° 高度:01≤001≤8,(m) 年份:1≤N≤365 步骤二:建立目标规划模型 由问题一,研究得到计算太阳影子的计算公式。若我们定义Lestimatei为与视频中所对应的第i个时刻代入经纬度、高度、日期和时间得到的太阳影子长度,L附件i为通过对附件中的数据进行处理得到的第i个时刻的太阳影子长度。当所有时刻这两个值的平方差最小时,这个地方将有最大的概率与视频中的地点相吻合,因为同一个地点之间由于有着相同的地理参数,它们在同一时刻的太阳影子长度必定完全吻合。得到以下模型: min=21i=1L附件i-Lestimatei)2 st-90°≤α≤90° (300-15t2)°≤β≤(15904-15t1)° 01≤001≤8,(m) 1≤N≤365 最后,对min函数数值求解的精度进行限制,当误差小于10-3时,停止遍历搜索,认为得到了最优的解决方案。 步骤三:建立基于分治法的遍历算法进行优化求解 (1)首先把四维向量空间划分为四个一维向量空间,分模块进行遍历搜索。 (2)接着,研究按照分治法的思想,分别对四个一维变量进行遍历搜索。同时进行全局搜索寻优。 (3)将第二步遍历寻优的结果与模型中研究要求的精度进行比较。若第二步中寻优的结果达到我们模型中所要求的精度要求时,遍历结束,否则返回第二步,进行递归的遍历求解。 23算法的求解 利用计算机模拟,我们得到已知数据中测量所在地数据如下表所示。 测量所在地 緯度(°N)经度(°E)杆高(m)日期 32248143223814/429 21573961331220 3基于GA-PSO算法对动态视频中太阳影子的大致地点的确定 31算法分析 研究在两种不同的情况下研究太阳影子的定位问题。首先,根据视频,研究可以得到在各个时间段所对应的太阳影子长度。接着,对于研究已有日期的视频部分,在目标规划模型的基础上进行改进,基于太阳影子长度公式,从而建立合理的规划模型。 32算法模型的建立 步骤一:数据处理 基于问题三建模的相关思路,我们定义两组变量α、β,分别代表经度和纬度。对数据进行了离散化处理,我们得到了以下关系: 纬度:-90°≤α≤90° 步骤二:建立目标规划模型 在问题三我们建立的目标规划模型的基础上,我们进行了适当的改进。若我们定义L′i为与视频中所对应的第i个时刻代入经纬度得到的太阳影子长度,Li为在视频中读出的第i个时刻的太阳影子长度。我们得到以下模型: min=22i=1(Li-L′i)2 st-90°≤α≤90° 675°≤β≤180° 同时,对于min的精度要求,我们定义,当误差小于10-3时,我们停止搜索,认为已经得到了最优解。 步骤三:用遗传算法优化粒子群算法(GA-PSO)以求得步骤二的最优解 在本算法中,为了得到最优解,我们设定了迭代次数为1000次。 Step1假定有一个S维目标搜索空间,其中第i个粒子表示为一个S维的向量 xi=xi1,x12,…,xiS,i=1,2,…,22 每一个粒子是一个潜在的解。将xi代入(3),我们可以算出它的适应值。第i个粒子飞翔的速度为S维向量,记为[AKV→]=Vi1,Vi2,…,ViS。在这里,我们设定每一个粒子存储了2个参数。同时,通过遗传算法的选择、交叉和遗传过程对所有变量进行初始化,记下第i个粒子迄今为止搜索到的最优位置为PiS[TX→]=PiS,…,PiS,整个粒子群搜索到的最优位置为PbestS[TX→]=PbestS,…,PbestS。 微粒i当前的最好位置可由下式确定: pi(t+1)= pi(t)→min(xi(t+1))≥min(xi(t)) Xi(t+1)→min(xi(t+1)) 根据Kennedy和Eberhart提出的相关理论,对粒子群进行以下操作: v=1min v(t+1)=v(t)+c1rand1(t)(piS(t)-xis(t))+c2rand2(pbestS(t)-xis(t)) xis(t+1)=xis(t)+v(t+1) 其中, i=[1,22],s=[1,S]; c1,c2分别为学习因子,rand(t)为产生随机数的函数,服从[0,1]区间的均匀分布。 根据上述方程组,分别对粒子群的速度和位置进行更新,若满足终止条件,则输出解,否则返回重新进行下一步的寻优。 最终,我们便可以得到视频拍摄地点的经度和纬度,从而进行比较精确的定位。 4算法推广 本模型所提出的“GA-PSO”算法对于解决大数据量的算法复杂度较高的问题有着较高的实用价值。该算法可以推广到人口相关数据的统计已经生产力的评估等多个领域。 参考文献: [1]司守奎,孙兆亮数学建模算法与应用北京国防工业出版社,2016(1) [2]卓金武,李必文,魏永生,等MATLAB在数学建模中的应用[M]北京北京航空航天大学出版社,2014(9) [基金项目]国家自然基金(项目编号:61322112)。 [作者简介]黄鸿基(1995—),男,汉族,江苏南京人,南京邮电大学本科生。研究方向:5G通信和智能算法;錢圳冰(1996—),男,汉族,江苏泰州人,南京邮电大学本科生。研究方向:计算机科学与技术;冯帆(1995—),男,汉族,江苏南京人,南京邮电大学本科生。研究方向:5G通信和智能算法;周行洲(1995—),男,汉族,江苏南京人,南京邮电大学本科生。研究方向:5G通信和智能算法。 |
随便看 |
|
科学优质学术资源、百科知识分享平台,免费提供知识科普、生活经验分享、中外学术论文、各类范文、学术文献、教学资料、学术期刊、会议、报纸、杂志、工具书等各类资源检索、在线阅读和软件app下载服务。