Ionospheric Correction and Analysis for Time-Series InSAR Based on Reformulating Range Split Spectrum Interferometry
-
摘要:
合成孔径雷达干涉技术（InSAR）受电离层延迟影响较大，特别是在低高纬地区，用于时序InSAR解算的干涉对数量大大减少，严重影响了低频合成孔径雷达（SAR）系统的监测精度. 针对该问题，引入重构距离向分谱法（RRSSI），将其扩展为时序InSAR电离层误差估计与校正方法. 为测试方法性能，获取2006年6月至2010年8月间覆盖阿拉斯加Anaktuvuk河区域的ALOS-1 PALSAR数据进行实验. 实验结果表明：本文方法能够有效估计时序InSAR结果中的电离层误差，且对小空间尺度电离层扰动具有较高的敏感性；经提出方法校正后，研究区域的年平均形变速率均值从校正前的1.46 cm/a降至0.49 cm/a，而标准差则从1.16 cm/a降至0.65 cm/a，精度显著提升；选取时序点的累积形变由校正前的大幅波动变得更加平稳，更符合形变规律.
-
关键词:
- 小基线集-InSAR /
- 电离层校正 /
- 距离向分谱法 /
- 时序
Abstract:Interferometric synthetic aperture radar (InSAR) is greatly affected by ionospheric delay, especially at low and high latitudes, which decreases the number of available interferometric pairs for time-series InSAR and the monitoring accuracy of low-frequency SAR systems. To this end, the reformulating range split spectrum interferometry (RRSSI) was introduced for the ionospheric error estimation and correction in time-series InSAR. The performance of the proposed method was tested by ALOS-1 PALSAR images covering the Anaktuvuk River area in Alaska from June 2006 to August 2010. The experimental results show that the proposed method can effectively estimate the ionospheric errors of time-series InSAR results and is sensitive to ionospheric disturbances at small spatial scales. After correction, the average annual deformation rate in the study area decreases from 1.46 to 0.49 cm/year, and the standard deviation from 1.16 to 0.65 cm/year, indicating an improvement in monitoring accuracy. Additionally, the cumulative deformation of the selected time-series points becomes more stable than the large fluctuation before correction, which is more consistent with the real deformation.
-
Key words:
- SBAS-InSAR /
- ionospheric correction /
- range split spectrum method /
- time series
-
中低速磁浮作为现有城市轨道交通的有益补充，具有安全低碳、静音舒适等优点，在我国具有广阔的应用前景[1]. 制动控制系统是磁浮列车自动驾驶系统（ATO）的重要子系统，在保障列车运行安全、停车精度和舒适性方面起到关键作用[2]. 为有效调节列车速度，中低速磁浮列车主要采用电-液混合制动实现正常到站停车. 由于电制动与液压制动转换过程状态强耦合、液压制动力强非线性等特点，传统的列车制动控制方法难以保障中低速磁浮列车的停车精度和舒适性. 因此，研究中低速磁浮列车电-液混合制动动态特性、混合制动优化控制方法，对提升磁浮列车停站精度和舒适性具有重要意义.
目前，针对传统列车制动过程的研究主要关注制动特性建模和控制算法方面. 精确的制动模型和良好的控制算法是列车平稳、精准停车的重要保障. 李中奇等[3]针对高速动车组动力分散特点建立多质点制动模型，基于径向基神经网络的比例-积分-微分（PID）控制策略与Smith预估器结合对给定速度跟踪控制，满足高速动车组制动过程精准停车要求. 周嘉俊等[4]分析列车实际减速度的影响因素，提出改进史密斯预估器对减速度进行预估控制，提高列车减速度的跟踪精度. 但该方法未识别造成列车减速度偏差原因，只针对跟踪误差对制动力进行调节的方式容易造成实际减速度的震荡，列车平稳性有待提升. 马天和等[5]采用参数估计算法对城轨列车制动过程受到的扰动进行估计补偿，通过在线修正目标制动力实现列车制动力闭环控制. 崔俊锋等[6]针对中低速磁浮列车制动力时滞问题，提出一种时滞补偿的广义预测控制器，有效降低列车制动过程时滞特性影响，提高控制精度. 传统基于物理参数建模方法难以准确表征磁浮列车动态特性，模型误差较大导致列车制动效果不佳. 为准确描述列车制动特性，机器学习算法被广泛用于列车数据驱动建模方面研究[7-8]. Yin等[9]针对高速列车制动系统延时问题，构建了一种滞后信息的长短期记忆网络（LSTM）模型，预测精度有很大提升. Liu等[10-11]结合回声状态神经网络建立高速列车状态预测模型，满足高速列车运行中安全、节能、舒适等多目标优化控制.
近年来，强化学习算法凭借其强大行动决策能力和解决不确定环境下控制问题的优势，在电力系统[12]、能源管理[13]、机器人[14]等领域取得了大量应用研究成果. 已有不少学者尝试将强化学习方法用于列车运行优化控制：张淼等[15]凭借优秀司机驾驶经验，设计一种基于策略梯度的强化学习算法，满足地铁变化场景下的能耗和舒适度等多目标需求； 高豪等[16]利用强化学习方法中的动态规划法优化地铁列车站间运行时间分配问题，有效降低列车牵引能耗；蒋灵明等[17]利用多智能体强化学习对大规模货运列车运行计划进行动态调整，有效降低货运列车延误影响； Shang等[18]提出一种带有参考系统的深度强化学习方法，有效提高智能体的学习效率； Liu等[19]提出一种基于深度Q网络（DQN）的重载列车智能控制方法，解决了列车在长陡坡情况下的循环空气制动控制问题.
本文针对中低速磁浮列车自动驾驶混合制动问题，提出一种混合制动特征自学习的深度强化学习（BFS-DQN）磁浮列车制动优化控制方法. 采用LSTM学习混合制动系统的动态特征，结合强化学习方法进行磁浮列车停站制动控制策略优化，并通过仿真实验验证其在中低速磁浮列车制动优化控制中的有效性.
1. 混合制动特征自学习方法
中低速磁浮列车混合制动特征因电液混合制动状态耦合、液压制动非线性而动态变化，目前广泛采用理论、静态制动特征模型，难以在仿真实验中准确反映强化学习的环境对agent （制动控制器）的奖励情况. 在此，采用LSTM学习混合制动的动态特征，以准确描述强化学习的奖励函数关系.
1.1 中低速磁浮列车混合制动原理
磁浮列车电-液混合制动原理如图1所示.
图 1 中低速磁浮列车电-液混合制动控制原理Figure 1. Principle of electro-hydraulic hybrid braking control of medium and low-speed maglev trainsATO发出一系列制动指令来操控磁浮列车减速停车. 磁浮列车制动过程遵循电制动优先、液压制动补偿的原则. 列车制动控制单元通过接收制动指令计算全列车所需制动力，按照电制动优先原则，向牵引控制单元发送电制动申请值，并根据反馈的电制动实际值计算所需补充的液压制动力. 然后，制动控制单元通过控制液压夹钳施加液压制动力，实现磁浮列车减速与停车控制.
1.2 混合制动力特征
中低速磁浮列车制动过程除了受到制动力控制外，还受到运行阻力的影响. 中低速磁悬浮列车的运行阻力比轮轨列车更为复杂，运行阻力包括空气阻力、涡流阻力以及附加阻力等.
空气阻力与列车运行速度及车辆编组相关，中低速磁浮列车空气阻力可表示为
Fa=(1.652+0.572N)v2， (1) 式中：$ N $为列车编组数，$ v $为列车速度.
集电器阻力Fc和电磁涡流阻力Fm为
{Fc=96N,Fm={3.354Mv,0<v⩽5.6m/s,(18.220+0.074v)M,v>5.6m/s, (2) 式中：$M$为列车质量.
附加阻力主要包括坡道阻力Fi和曲线阻力Fr，如式（3）.
{Fi=Mgsin(arctanβ),Fr=600RMg, (3) 式中：$\;\beta $为坡度值，$R$为曲线半径.
由于中低速磁浮列车电-液混合制动力实际值未知，本文将列车制动力特性用函数$ g(x) $表示. 根据中低速磁浮列车制动特性和运行阻力，建立基于牛顿方程的单质点列车制动模型，如式（4）.
{v=dsdt，a=dvdt,Ma=g(Fb,Fw)+ξ(t),Fw=Fa+Fc+Fm+Fr+Fi, (4) 式中：$s$和$a$为列车位置和加速度，${F_{\rm{b}}}$为列车制动力，${F_{\rm{w}}}$为列车运行阻力，$ \xi(t) $ 为未建模动态，t为时间.
1.3 混合制动动态特征自学习方法
中低速磁浮列车制动过程电制动、液压制动及总制动力变化特性如图2所示，图中：Bmax为最大理论制动力，d为制动距离，B为制动力. 设列车制动系统控制单元发出制动指令时列车制动距离为${d_0}$，由于电制动力响应存在短暂延时$\tau $，列车空走一段时间后距离为${d_1}$. 此时，电制动力开始减小，液压制动开始补充；制动距离为${d_2}$，电制动撤离，液压制动接管整车制动；当制动距离为${d_3}$时，液压制动力达到最大值，列车准备停车.
磁浮列车制动特性受到电制动响应延迟、液压制动离散及运行阻力变化等因素影响，列车实际制动力$B_{{\rm{s}}}$与理论制动力$ B_{\mathrm{n}} $存在较大差距，影响列车到站的停车精度. 同时，电制动和液压制动两者强非线性耦合特点导致磁浮列车电-液混合制动过渡不平稳，容易产生较大纵向冲击，对乘坐舒适性造成影响.
传统的机理模型难以准确描述中低速磁浮列车运行动态，无法保障列车制动效果. LSTM在处理各种非线性时序问题具有良好的效果，本文采用LSTM网络建立磁浮列车数据驱动的制动模型，利用列车历史制动数据进行制动特征自学习，以更准确地描述列车动态特性和预测列车运行状态. 基于LSTM的磁浮列车动态模型为
{xt+1=gLSTM(X(t),U(t),w)+ξ(t),X(t)={x1,x2,⋅⋅⋅,xT},U(t)={u0,u1,⋅⋅⋅,uT}, (5) 式中：${g_{{\rm{LSTM}}}}$(·)为基于磁浮列车动态模型网络；${x_{t + 1}} = \{ {d_{t + 1}},{v_{t + 1}}\} $，为列车下一时刻预测状态，即列车位置和速度；$U(t) $为列车前$T$时刻一系列的制动指令；$X(t) $为列车前$T$时刻的运行状态；$w = \{ {w_1},{w_2}, \cdot \cdot \cdot ,{w_n}\} $，为网络训练参数.
为了减小预测输出值与实际值之间的差距，采用均方根误差（RMSE）方法建立损失函数$l(\theta )$且使其最小化，如式（6）.
min (6) 式中：${x_t}$和${\hat x_t}$分别为列车时刻$t$的实际状态值和预测值.
本文建立基于LSTM网络的磁浮列车动态模型的目的是准确预测列车未来状态，如图3所示. 图中：Ct为时刻t时LSTM单元状态，ht为时刻t时LSTM的输出值. 制动指令${u_0}$和列车初始状态$\{ {d_0},{v_0}\} $作为已知输入，选取磁浮列车历史制动数据中一系列制动指令$\left\{ {{u_0},{u_1}, \cdot \cdot \cdot ,{u_T}} \right\}$对模型进行训练，模型可以输出一系列观测到的状态$\{ {d_1},{v_1}, \cdot \cdot \cdot ,{d_T},{v_T}\} $，最后将预测状态与磁浮列车现场采集的制动数据比较，验证模型准确性.
2. 磁浮列车BFS-DQN制动控制方法
深度强化学习算法融合了深度学习和强化学习各自优势，具有优秀的状态感知能力和强化的决策能力[20]. 本文主要关注提出一种BFS-DQN磁浮列车制动控制方法，通过学习磁浮列车制动特性以便准确获取运行状态和计算优化目标，提高端到端的强化学习控制效果.
本文提出的BFS-DQN控制方法如图4所示. 首先，通过列车历史运行数据训练LSTM网络，学习列车制动特性与动态环境的关系； 然后，将观测到的列车运行状态作为强化学习算法的输入，列车制动等级作为智能体的动作输出，同时设计舒适性和精准停车奖励函数，通过与环境交互的方式训练列车控制器学习最优控制策略；最后，输出合理的制动力以实现磁浮列车制动优化控制.
2.1 状态空间
根据式（4）、（5），将预测到的列车运行状态，即列车位置$d$、速度$v$和加速度$a$组成的三维状态空间作为强化学习智能体的观测状态，该状态空间S定义为
\begin{split} & S = \{ {\boldsymbol{s}} = {[d\;\;v\;\;a]^{\rm{T}}},\quad d \in [0,{d_{\mathrm{s}}}],v \in [0,{V_0}], \\ &\quad a \in [ - {A_{\max }},0]\}, \end{split} (7) 式中：${d_{\mathrm{s}}}$为列车目标停车位置，${V_0}$为列车电-液制动转换起始速度，${A_{\max }}$为列车最大减速度.
2.2 动作空间
磁浮列车制动过程中，ATO或司机根据运行状况操纵制动等级来调节列车制动力的大小，同样，本文选择列车制动等级作为智能体的动作，动作空间A定义为
{\boldsymbol{A}} = [{u_i}],\quad{u_i} \in [0,100{\text{%}} ], (8) 式中：${u_i}$为磁浮列车电-液混合制动等级，0%表示惰行，100%表示施加最大制动力.
2.3 奖励机制
BFS-DQN中，智能体通过每个时间步采取最佳行动来获得更多的奖励回报，可靠的奖励机制对于引导智能体学习最佳行动策略具有重要作用. 因此，本文根据磁浮列车制动优化控制目标，设计关于舒适性和精准停车的奖励函数，最后通过加权求和的方式得到总奖励，该奖励函数定义为
\left\{ \begin{array}{l} {r_{{\rm{total}}}} = {\lambda _1}{r_{{\rm{comfort}}}} + {\lambda _2}{r_{{\rm{stop}}}} + {r_{{\rm{fail}}}}, \\ {r_{{\rm{comfort}}}} = 1 - a_{{\rm{error}}}^2,\quad {a_{{\rm{error}}}} \leqslant 1.0\; {\rm{m/s}}^{3}, \\ {r_{{\rm{stop}}}} = 1 - d_{_{{\rm{error}}}}^2,\quad {d_{{\rm{error}}}} \leqslant 0.3\; {\rm{m}}, \\ {r_{{\rm{fail}}}}, \quad {a_{{\rm{error}}}} \gt 1.0\; {\rm{m/s}}^{3}\;{\mathrm{or}}\; {d_{{\rm{error}}}} \gt 0.3\;{\rm{m}}, \end{array} \right. (9) 式中：${r_{{\rm{total}}}}$为列车前后状态转移产生的即时总奖励；${\lambda _1}$,${\lambda _2}$为奖励权重；${r_{{\rm{confort}}}}$为乘客舒适性奖励函数；${a_{{\rm{error}}}}$为单位时间内减速度的变化激烈程度，用于描述乘客乘车体验[21]；${r_{{\rm{stop}}}}$为列车精准停车的奖励函数；${d_{{\rm{error}}}}$为列车实际制动距离与目标制动距离的差值，用于描述停车精准性；为节省时间成本，设置了惩罚函数${r_{{\rm{fail}}}}$对停车失败和减速度剧烈变化的情况给予训练失败惩罚.
2.4 学习策略
本文定义列车控制器作为BFS-DQN中的智能体，采用价值迭代方式学习最优的制动策略. 首先，定义一个价值函数${Q_\pi }(s,u)$用于描述智能体从开始到结束时的预期未来总奖励，列车制动控制器采用随机策略$\pi $进行动作学习时，从当前时间$t\;（t = 0,1, 2,\cdots）$观察环境状态${s_t}$（${s_t} \in S$），控制器再从制动级位的集合$A$中选择一个动作${u_t}$. 然后，环境根据列车当前状态信息和制动级位给予一定奖励${r_t}$且列车转移到新的状态${s_{t + 1}}$. 列车控制器通过评估每个动作的价值，并采取最大价值更新${Q_\pi }({s_t},{u_t})$，更新方式为
{Q_\pi }\left( {{s_t},{u_t}} \right) = {r_t} + \gamma \max {Q_\pi }\left( {{s_{t + 1}},{u_{t + 1}}} \right), (10) 式中：γ为目标值的折扣系数.
为了得到最大的奖励回报，假设状态转移概率$p$已知，根据贝尔曼最优方程和最优性原理求得最优动作策略${\pi ^ * }({s_t}|{u_t})$，智能体经过多次价值迭代后其价值函数${Q_\pi }({s_t},{u_t})$收敛于最优价值函数$ Q_\pi ^ * ({s_t},{u_t}) $，此时，最优价值函数$ Q_\pi ^ * ({s_t},{u_t}) $等于最优状态值函数$ V_\pi ^ * ({s_t}) $，策略更新如式（10）[22].
V_\pi ^ * ({s_t}) = Q_\pi ^ * ({s_t},{u_t}) = \mathop {\max }\limits_u \left({r_t} + \gamma \sum\limits_{s \in S} p({s_{t + 1}}|{s_t},{u_t}){V^ * }({s_t})\right), (11) 式中:$ {V^ * }({s_t}) $为时刻$t$的状态价值函数.
2.5 学习过程
2.5.1 数据生成
列车制动控制器通过观测当前时刻运行状态${s_t}$，主Q网络对当前状态根据动作选择策略输出对应列车制动等级的动作${u_t}$，列车状态预测模型根据动作和当前状态生成新的状态${s_{t + 1}}$；环境根据列车制动过程优化目标设计的奖励函数，计算状态${s_t}$转移到新状态${s_{t + 1}}$的奖励${r_t}$.
2.5.2 经验回放
将时刻$t$状态、动作和奖励数据$ \lt {s_t},{u_t},{r_t},{s_{t + 1}} \gt $存入经验池中，作为列车制动控制器的学习数据. 训练时，随机从经验池中抽取一定数量的样本$ \lt {s_j}, {u_j},{r_j},{s_{j + 1}} \gt $作为训练数据，将状态和动作$ \lt {s_j},{u_j} \gt $输入到主Q网络得到对应动作的$Q({s_j},{u_j};{\theta _j})$（$\theta_j $为当前神经网络权重系数），$ \lt {s_{j + 1}},{u_{j + 1}} \gt $输入到目标Q网络中得到目标值$\max \tilde Q({s_{j + 1}},{u_{j + 1}};{\theta _{j + 1}})$（$\theta_{j+1} $为目标网络的权重系数）.
2.5.3 网络更新
BFS-DQN中Q网络通过式（12）更新参数权重.
\left\{ \begin{gathered} L({\theta _j}) = {E}[({r_j} + \gamma \mathop {\max }\limits_u \tilde Q({s_{j + 1}},{u_{j + 1}};{\theta _{j + 1}}) - \\ \quad Q^2{({s_j},{u_j};{\theta _j})}], \\ {\theta _{j + 1}} = {\theta _j} + \alpha {\rm{loss}}(\theta )\nabla Q({s_j},{u_j};{\theta _j}), \\ \end{gathered} \right. (12) 式中： $\alpha $为网络学习率.
间隔一定训练次数复制主Q网络的权重参数到目标Q网络. 上述过程经过一定次数的循环训练，主Q网络收敛于目标Q网络，列车制动控制器学习到最优制动策略. BFS-DQN训练示意如图5.
图 5 强化学习制动优化控制算法流程Figure 5. Flowchart of optimization control algorithm for reinforcement learning braking3. 实验结果与分析
3.1 实验设置
为验证本文所提BFS-DQN方法有效性，以湖南省某条磁浮线路现场运营数据为例进行仿真验证，列车主要参数详见表1所示.
表 1 仿真列车参数Table 1. Simulation train parameters参数类别 参数特性 列车质量/t 75 线路最高限速/（km·h−1） 80 编组数量 3 最大常用制动力/kN 74.23 最大常用减速度/（m·s−2） 0.96 线路最大坡度/‰ 51.01 3.2 收敛性对比
本文采用传统DQN控制方法与BFS-DQN分别进行仿真，主要训练参数详见表2.
表 2 算法主要训练参数Table 2. Main training parameters for algorithm参数 BFS-DQN DQN LSTM 迭代次数/次 500 LSTM 学习率 0.001 LSTM 样本批量 50 单次训练最大步数/步 80 80 训练最大次数/次 20000 20000 Q 网络学习率 0.001 0.001 Q 网络更新频率 100 100 样本大小 32 32 经验池容量 2000 2000 折扣因子 0.96 0.96 贪婪率初始值 0.9 0.9 贪婪率最终值 0.1 0.1 同时，本文采用动态$ {\varepsilon {\text{-}} {\rm{greedy}}} $优化算法使智能体能够合理利用探索策略. 训练初期采用较大探索率，随着训练次数增加，探索率减小，智能体偏向选择奖励值最大的动作. 对列车控制器进行20000次训练，每250轮计算平均奖励值，通过观察奖励函数变化情况来判断算法是否收敛. BFS-DQN和DQN训练过程奖励变化情况详见图6所示.
从图6可以看出：训练初BFS-DQN和DQN的奖励值都呈现出明显上升趋势，说明智能体能够迅速找到合适的动作策略；训练次数达到12000次以上，2种优化算法平均奖励值缓慢上升并趋于稳定；BFS-DQN方法平均奖励高于传统DQN方法，制动特征自学习有助于强化学习智能体获取准确的状态信息和奖励，提高网络学习能力.
3.3 舒适性和停车精度
图7和图8分别为训练过程中平均加速度变化和平均停车误差变化情况. 统计最后5000次训练结果，加速度变化和停车误差情况如表3所示.
图 8 停车误差收敛情况（DQN和BFS-DQN分别取停车误差的绝对值和绝对值负值）Figure 8. Convergence of parking errors (DQN and BFS-DQN take the absolute value and negative value of absolute value of parking error, respectively)从图7和图8可以看出：初期加速度变化和停车误差都很大，但加速度变化很快收敛，由于停车奖励需要列车到站才能获得，智能体经过更多次的训练才寻得最优的停车策略. 由表3可得：BFS-DQN在舒适性和停车精度方面比传统DQN提高8.67%和33.33%. 训练结果表明：基于BFS-DQN磁浮列车制动优化控制方法取得了良好的控制效果，且比传统DQN方法在舒适性和停车精度方面更具有优越性.
表 3 算法训练结果Table 3. Training results for algorithm训练结果 BFS-DQN DQN 平均奖励值 33.5 27.8 平均状态转移次数/次 70 72 平均停车误差/m 0.10 0.15 平均加速度变化/（cm·s−3） 10.84 11.78 平均制动时间/s 14.0 14.4 3.4 实际运行数据验证实验
为验证BFS-DQN算法可行性，将训练好的列车控制器进行仿真测试，并记录50次测试的停车误差数据. 同时，本文分别采用均方根误差（RMSE）和标准差（SD）来评价不同制动算法的性能，结果如表4所示. 取其中一次数据的制动曲线如图9所示.
表 4 算法性能Table 4. Algorithm performance制动控制策略 RMSE SD BFS-DQN 0.099048 0.070652 DQN 0.142815 0.110446 传统 ATO 0.276103 0.140018 由表4可见：3种制动控制算法中，BFS-DQN制动控制算法的表现最好，停车误差小，且离散程度低，稳定性高. 由图9可知：采用传统ATO制动策略生成的制动曲线波动较大，而BFS-DQN和DQN制动策略生成的制动曲线较为平滑；停车精度方面，BFS-DQN、DQN和传统ATO 制动策略下的制动距离分别为24.26、24.19、24.45 m，停车误差分别为0.04、0.11、0.15 m，BFS-DQN和DQN方法都能够有效提高磁浮列车到站停车精度，且BFS-DQN方法更具优越性.
图10为3种制动策略下减速度变化情况，由图10分析可知：由于传统ATO制动策略采用电液混合制动-惰行-液压制动的方式进行制动停车，当磁浮列车制动距离为6.25、11.30、23.75 m时，由于电制动响应延迟和液压制动离散等特点，制动等级突变引起减速度剧烈变化，减速度最大改变量达到1.02 m/s3，较大的纵向冲击对乘客的安全性和舒适性造成影响；BFS-DQN方法在减速度控制方面更加合理，在列车停车的最后一刻的减速度最大改变量为0.6 m/s3，有效减少制动过程中产生的纵向冲击，舒适性上提高41.18%.
3.5 现场数据实验验证
为进一步检验制动控制器的性能和稳定性，将目标制动距离设为24.60、24.90、25.10 m，依次进行精准停车测试，50次测试结果如表5和图11所示.
表 5 停车误差分布情况Table 5. Distribution of parking errors次 停车误差/m BFS-DQN DQN ATO $ x \lt - 0.5 $ 0 0 3 $ - 0.5 \leqslant x \leqslant - 0.3 $ 0 2 3 $ - 0.3 \lt x \leqslant 0 $ 18 19 16 $0 \lt x \leqslant 0.3$ 32 29 20 $ 0.3 \lt x \leqslant 0.5 $ 0 0 6 $ x \gt 0.5 $ 0 0 2 图11可以看出：BFS-DQN和DQN方法的停车误差分别在±20 cm与±24 cm之间波动，表明BFS-DQN在不同目标制动距离干扰下稳定性要高于DQN. 由表3数据统计得：当前ATO制动策略下列车满足±30 cm精准停车要求39次，占比78.0%；DQN方法的停车精度满足±30 cm内48次，占比96.0%；本文提出的BFS-DQN算法满足停车精度±30 cm内50次，占比100%.
4. 结 论
1） 本文针对磁浮列车电-液混合制动问题提出了基于混合制动特征自学习的磁浮列车强化学习制动控制方法. 利用磁浮列车运行数据进行仿真实验，结果验证所提控制策略可以帮助磁浮列车平稳度过电液制动转换阶段，减少列车制动过程产生的纵向冲击，保障乘客舒适性的同时也提高了列车到站停车精度.
2） 对于不同目标制动距离的干扰，所提策略能够根据列车反馈的运行状态及时调整制动策略，控制列车输出合理的制动力，确保列车满足精准停车要求，具有一定的稳定性和自适应能力.
-
图 1 时序InSAR电离层延迟估计与校正
Figure 1. Delay estimation and correction for time-series InSAR ionospheric
图 3 电离层校正前后的重缠绕差分InSAR干涉图
Figure 3. Rewrapped differential InSAR interferograms before and after ionospheric correction
-
[1] WRIGHT P A, QUEGN S, WHEADON N S, et al. Faraday rotation effects on L-band spaceborne SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing,2003,41(12):2735-2744. [2] PI X Q. Ionospheric effects on spaceborne synthetic aperture radar and a new capability of imaging the ionosphere from space[J]. Space Weather,2015,13(11):737-741. [3] MEYER F, BAMLER R, JAKOWSKI N, et al. The potential of low-frequency SAR systems for mapping ionospheric TEC distributions[J]. IEEE Geoscience and Remote Sensing Letters,2006,3(4):560-564. [4] MAO W F, LIU G X, WANG X W, et al. An InSAR ionospheric correction method based on variance component estimation with integration of MAI and RSS measurements[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2021,14:1423-1433. doi: 10.1109/JSTARS.2020.3045267 [5] ZHANG B C, DING X L, ZHU W, et al. Mitigating ionospheric artifacts in coseismic interferogram based on offset field derived from ALOS-PALSAR data[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2016,9(7):3050-3059. doi: 10.1109/JSTARS.2016.2533441 [6] KAZUO O C. Recent trend and advance of synthetic aperture radar with selected topics[J]. Remote Sensing,2013,5(2):716-807. [7] ROTT H. Advances in interferometric synthetic aperture radar (InSAR) in earth system science[J]. Progress in Physical Geography: Earth and Environment,2009,33(6):769-791. doi: 10.1177/0309133309350263 [8] 刘国祥,陈强,罗小军,等. InSAR原理与应用[M]. 北京: 科学出版社, 2019. [9] AOKI Y, FURUYA M, DE ZAN F, et al. L-band synthetic aperture radar: current and future applications to earth sciences[J]. Earth, Planets and Space,2021,73(1):1-4. doi: 10.1186/s40623-020-01323-x [10] RAUCOULES D, DE MICHELE M. Assessing ionospheric influence on L-band SAR data: implications on coseismic displacement measurements of the 2008 Sichuan earthquake[J]. IEEE Geoscience and Remote Sensing Letters,2010,7(2):286-290. [11] JUNG H S, LEE D T, LU Z, et al. Ionospheric correction of SAR interferograms by multiple-aperture interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing,2013,51(5):3191-3199. doi: 10.1109/TGRS.2012.2218660 [12] WEGMULLER U, STROZZI T, WERNER C. Ionospheric path delay estimation using split-beam interferometry[C]//2012 IEEE International Geoscience and Remote Sensing Symposium. Munich: IEEE, 2012: 3631-3634. [13] FREEMAN A F. Calibration of linearly polarized polarimetric SAR data subject to faraday rotation[J]. IEEE Transactions on Geoscience and Remote Sensing,2004,42(8):1617-1624. [14] ZHU W, JUNG H S, CHEN J Y. Synthetic aperture radar interferometry (InSAR) ionospheric correction based on faraday rotation: two case studies[J]. Applied Sciences,2019,9(18):3871.1-3871.19. [15] GOMBA G, PARIZZI A, DE ZAN F, et al. Toward operational compensation of ionospheric effects in SAR interferograms: the split-spectrum method[J]. IEEE Transactions on Geoscience and Remote Sensing,2015,54(3):1446-1461. [16] FATTAHI H, SIMONS M, AGRAM P. InSAR time-series estimation of the ionospheric phase delay: an extension of the split range-spectrum technique[J]. IEEE Transactions on Geoscience and Remote Sensing,2017,55(10):5984-5996. [17] GOMBA G, DE ZAN F, PARIZZI A. Ionospheric phase screen and ionospheric azimuth shift estimation combining the split-spectrum and multi-squint methods[C]//Proceedings of EUSAR 2016: 11th European Conference on Synthetic Aperture Radar. Hamburg: VDE, 2016: 1-4. [18] LIANG C R, LIU Z, FIELDING E J, et al. InSAR time series analysis of L-band wide-swath SAR data acquired by ALOS-2[J]. IEEE Transactions on Geoscience and Remote Sensing,2018,56(8):4492-4506. doi: 10.1109/TGRS.2018.2821150 [19] LIANG C R, AGRAM P, SIMONS M, et al. Ionospheric correction of InSAR time series analysis of C-band sentinel-1 TOPS data[J]. IEEE Transactions on Geoscience and Remote Sensing,2019,57(9):6755-6773. [20] WEGMÜLLER U, WERNER C, FREY O, et al. Reformulating the split-spectrum method to facilitate the estimation and compensation of the ionospheric phase in SAR interferograms[J]. Procedia Computer Science,2018,138:318-325. doi: 10.1016/j.procs.2018.10.045 [21] BERARDINO P, FORNARO G, LANARI R, et al. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing,2002,40(11):2375-2383. [22] ZHANG Y, FATTAHI H, AMELUNG F. Small baseline InSAR time series analysis: Unwrapping error correction and noise reduction[J]. Computers & Geosciences,2019,133:104331.1-104331.29. [23] JONES B M, KOLDEN C A, JANDT R, et al. Fire behavior, weather, and burn severity of the 2007 anaktuvuk river tundra fire, north slope, Alaska[J]. Arctic, Antarctic, and Alpine Research,2009,41(3):309-316. doi: 10.1657/1938-4246-41.3.309 [24] LIU L, JAFAROV E E, SCHAEFER K M, et al. InSAR detects increase in surface subsidence caused by an Arctic tundra fire[J]. Geophysical Research Letters,2014,41(11):3906-3913. doi: 10.1002/2014GL060533 [25] LIAO H. Ionospheric correction of interferometric sar data with application to the cryospheric sciences[D]. Alaska Fairbanks: University of Alaska Fairbanks, 2018. [26] ZHANG B. Mitigation of ionospheric artifacts in InSAR data for estimating earthquake deformation[D]. Hong Kong: The Hong Kong Polytechnic University, 2020. -
下载:
下载:
点击查看大图
计量
- 文章访问数: 321
- HTML全文浏览量: 138
- PDF下载量: 48
- 被引次数: 0