基于MQP模型电离图参数反演解析算法研究.docVIP

基于MQP模型电离图参数反演解析算法研究 　　摘要：针对现今普遍采用的偏导迭代电离图反演算法存在的求解依赖于初值选择的不足，给出了一种初值选择准则，保证了该算法的收敛性。同时，利用该准则，以及在对电离图虚高表达式进行合理近似的基础上,提出了一种新的电离图参数反演算法，该方法实现电离图参数反演的解析计算，较偏导迭代法计算量大大减少，更好满足工程实时性要求。 　　关键词:电离图；虚高；参数反演 　　中图分类号：TP381.1文献标识码：Adoi:10.3969/j.issn.1003-6970.2010.10.005 　　 　　An Analytic Algorithm on Inversion of Ionograms Parameters Based on MQP Model 　　 　　Hou Cheng-yu， Shen Yi-ying， Meng Wei-xiao 　　(School of Electronics information Engineering, Harbin Institute of Technology, Harbin 150001, China) 　　 　　Abstract:The algorithm of partial derivative iteration is often used to inverse the parameters of Ionograms. But the accuracy of this inversion results are depended on the selection of initial values. So a criterion for the selection of the initial value is proposed to solve this problem. Meanwhile, using this criterion, and the expressions of ionograms virtual high, based on a reasonable approximation, a new inversion algorithm on the ionograms’ parameters is presented. This method may calculate the ionogram’s parameters with analytical formula, which makes the calculation greatly reduced compared with the partial derivative iteration method. It will better meet engineering needs in real time. 　　Keywords: Ionogram; Virtual height; Parameters inversion 　　 　　1. 引 言 　　 　　 电离层作为地球大气的组成部分之一，对经其传播的电磁波将产生折射、吸收和时延等现象。无论是确保天波雷达系统正常工作[1]，还是反演地球大气参数（如利用GPS无线电掩星法）[2]，研究电离层的信息都是不可或缺。通常这些信息是利用电离图反演得到的[3]，因此，提高反演算法稳健性和自动化程度，降低计算量是实现上述工程应用所必需解决的问题。 　　 电离图是由电离层垂直测高仪扫频得到并利用测量回波和发射脉冲之间的延迟，得到电离层反射高度（称之为虚高）随扫描频率变化的曲线。根据该曲线，以及所建立的电离层模型，求解模型参数，进而反演得到电离层参数信息（主要是电离层各层的临界频率、峰值高度、半层厚度以及一些电离层不规则结构的信息）。如图1所示。 　　电离图反演的基础是电离层模型。在上个世纪九十年代，CCIR推荐采用修正B/D模型以及多准抛物线（MQP）模型。其中，MQP模型层数可调节，并采用准抛物线层与翻转抛物线层组合保证了层间电子浓度梯度的连续，较好表征了各层电子浓度分布特性，被人们广泛采用[4]。 　　基于MQP模型的电离图反演方法有切比雪夫多项式法[5]、偏导迭代法[6]等。后者相对前者简洁，参数含义明确而被普遍采用。对某一层参数反演时，偏导迭代法是通过输入至少两点经该层反射的电离图数据进行迭代而求解的。而在使用中发现，一旦数据选择不当，该算法将失效。所以，该方法要求数据录入者具备一定电离层知识，能选择合适的数据点，在算法失效时需人为干预，重新选择数据点，不便进行自动化处理。 　　本文将从MQP模型虚高表达式入手，在合理近似的基础上，推导出虚高与各电离层参数关系的简约表达式，进而寻求一种电