2016.11.7-双向时间比对仿真框架.docxVIP

双向时间同步算法的MATLAB仿真框架1.A和B的位置坐标A卫星随时间变化的位置坐标：B卫星随时间变化的位置坐标：2.B卫星向A发送时标信息，A卫星接收到之后，又将自己的信息返回给B的过程：在 t = 0s时刻，B向A发送带有时标信号的信息，并计算出A应该接收的时刻，计算过程是：for(i = 1; i = n; i++){if( ).....}说明A接收到信号的位置一定在和之间，采用位置插值法或最小二乘法拟合曲线，找出这个位置，即求出了相遇时刻。最小二乘的曲线拟合：为了使R最小，我们对系数a求偏导，得：...化成矩阵的形式：就可以求出系数.在 时刻 ，A接收到B发送过来的信号,同时记录下自己此时的钟面时刻，然后计算下行伪时延：。在 时刻，A向B发送自己的时标信号和上一步计算的下行伪时延，并用同样的方法算出B接收到信号的时刻。在 时刻，B接收A发送过来包含和的信号，并记录下此时刻的钟面时间，然后计算上行伪时延：。第一轮过程结束。然后重复上面的过程，就可以得到一系列伪时延数据：当B卫星接收到两组数据之后，就可以通过公式计算出与A的钟差：3.根据各个时刻列出事件调度表：时刻t(s)相应时刻应该处理的事件0sB发送A接收B的，记录，计算A向B发送和B接收和，记录，计算1B发送......B第二次接收到A的消息，B根据已知数据 计算出钟差......将各个事件时刻从小到大进行排列，即事件发生的先后顺序。例如，B第二次发送消息的时刻 1s，如果小于，就应该在其之前处理相应的事件。4.Matlab程序的框架：将事件调度表定义成一个元胞数组,每次从数组中读取时间数据。A和B的钟面时刻分别定义为读取当前计算机时间，和当前计算时间加上一个常数。While(1){t = t1（t0）;B发送信息：即记录钟面时间，B1 = current-clock;B计算时间：for(i = 1; i = n; i++){if( ).....}t = t2（）;A接收B的，记录（A的钟面时刻定义为current-clock + C）计算：= A1-B1;t = t3（）;A向B发送和:将A的两个数据定义为数组A = [ ,] ;t = t4();B接收和，记录，计算:=B2 - A2;再定义一个数组用来存放[ , ];当B接收到两组数据后，就可以计算出一个钟差：计算结果存储在一个数组中[ , , , ].If（当达到循环次数 || 事件列表里面没有事件 (t = null)）跳出循环 Break；}