用户名: 密码: 验证码:
基于组合测量的弹箭图像末修技术
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
为提高单兵武器对目标的打击精度和我方人员在战场上的生存概率,本文以单兵火箭弹为平台,研究了基于组合测量的弹箭图像末修技术。基于组合测量的弹箭图像末修技术以弹载CCD为目标探测器,利用地磁计和MEMS陀螺组合测量弹丸姿态,由微处理计算弹道偏差,并在弹道末段对弹丸进行脉冲修正。该技术不但可以用于单兵火箭弹,还可以用于其它野战火箭弹及各种炮兵弹药上。本文的主要研究内容有:
     1.以单兵火箭为研究背景,详细阐述了基于组合测量的图像末修弹药的系统组成和工作原理。通过对系统的三个关键子系统进行分析,重点研究了图像测量系统的目标测量原理、测量精度和姿态测量系统中的MEMS陀螺误差补偿方法。根据单兵火箭的运动特点,建立了适用于快速、准确计算的质点弹道方程和带修正的4自由度弹道方程。利用Simulink和6自由度弹道方程对某型单兵火箭弹进行了外弹道模型仿真。
     2.研究了利用CCD测量弹目相对位置的方法。根据弹丸运动方程和CCD成像特点,建立了弹丸运动时弹目相对运动方程和目标点在CCD上的成像轨迹方程;研究了利用CCD测量目标相对位置的估测方法,同时设计了利用单个CCD测量弹目相对位置的方法;分析了由于弹丸旋转对弹目测量所造成的误差;对各模型进行了Matlab仿真及转台实验验证。仿真及实验结果表明,利用文中模型可在一定条件下较准确测量弹目相对位置,但测量精度受姿态及基线距离测量精度的影响较大。
     3.研究了组合测量弹丸姿态的方法。针对单兵武器平台的弹丸运动特点,结合地磁矢量在弹体坐标系内的投影规律,利用椭圆拟合方法,建立了计算弹丸运动条件下的姿态角模型;通过引入中位数和修正系数,将中位法、UKF和椭圆理论相结合,提出了改进的椭圆参数拟合方法,提高了在奇异值干扰下的弹丸滚转角计算精度;设计了基于伺服稳定平台的弹丸姿态角测量方法;利用三轴地磁计和MEMS陀螺组成磁-陀螺系统,建立了组合测量模型,并进行了转台仿真实验。从对实验数据的处理结果可以得出:利用椭圆理论计算弹丸滚转角能得到较好的精度,而且实用性很强,但计算其余两姿态角的精度不高;磁-陀螺测量系统可互补磁强计和陀螺测量中的不足,测量弹丸姿态能得到较高的精度。
     4.以验证组合测量弹目相对位置模型为目的,设计了弹载测试系统和实验方案,并完成了用弹载测试系统组合测量弹目相对方位和弹丸姿态的原理性演示飞行实验。实验结果表明,本次实验较成功的测得了弹丸的姿态角和目标相对弹丸方位的变化规律,测量结果与文中仿真结果比较接近。
     5.研究了基于单兵火箭平台的修正策略。利用前文测量方法计算求取弹目相对偏差作为弹道修正量,将修正量分为距离偏差和角度偏差分别求取;根据某型脉冲发动机和单兵火箭运动特点,建立了弹道修正过程中的弹道变化模型;利用最优控制理论和弹道修正量,设计了以冲量为最优控制量的距离偏差、角度偏差的修正律;分析了修正过程中弹丸转速对修正的影响,提出了使弹丸保持稳定的方法,并计算出了本文修正方案的修正效率;最后设计了脉冲修正点火方案;并进行了仿真计算和转台实验进行验证。
In order to improve the firing accuracy to target and survival probabilities in battlefield, the technology of terminal trajectory correction is studied in this paper based on the combined measured method, applied in individual rocket. The technology makes use of CCD which is fixed in projectile as target detector and combines magnetometers and MEMS gyro to measure the projectile's attitudes, and then the trajectory error is calculated by microcomputer, and connect the trajectory by controlling small pulse motors at the end of trajectory at last. The technology not only can be used in individual rocket, but also in other artillery rockets and munitions.
     1. This paper expounds the system's composing and working principle of terminal trajectory correction munitions based on the combined measured technology and CCD. The principle detected target and accuracy of the image measured system and the error compensated method of MEMS gyro in attitude measured system are studied through analyzing to the three key sub-systems. According to the characteristics of projectile's motion, the projectile equations are built in this paper which can calculate the trajectory quickly and accurately. The trajectory of individual rocket is simulated by using of Simulink and 6-D equations.
     2. The method measured to target by CCD is studied in this paper. According to motion equations of projectile and the characteristics of imaging in CCD, the equations disturbed the relative motion between projectile and target is built and the method which is used to measure the relative position of target by CCD is studied in this paper. The measured error caused by the projectile rotation is analyzed. To test the above methods and equations, several experiments and simulations are designed. The results indicate that the relative positions of target can be measured accurately by the methods and equations.
     3. The method measured to projectile's attitudes is studied in this paper. According to the characteristics of individual rocket's motion, the attitudes measured model is built by using of B2AC and the projective geometry of geomagnetic vector in projectile coordinate. Through bringing several coefficients and UKF into B2AC, the improved B2AC is proposed which can calculated the attitude angle of projectile accurately when a part of the original data is not correct. A method is proposed which can measured the attitude based on the servo mechanism. The combine measured model is built using three-axes magnetometers and MEMS. To test the above methods and equations, several experiments and simulations are designed. The results indicate that the roll angle can be calculated accurately using improved B2AC, but the other two angles' accuracy is not very good while the combine measured model can compensate the shortage.
     4. To test the above models and methods, the experiment plan and detect system are designed. The test is finished to measure the attitude of projectile and the relative position between the projectile and target. The results show that the data is identical with the simulated data.
     5. The correct policy is studied in this paper. Using the trajectory error as correct parameter which can be calculated by the above models and methods, the trajectory equations which have been changed are built firstly, and then the correct policy is designed by use of optimized cybernetics. The influence of rotation to correcting is analyzed and the method to keep the projectile stabilized is proposed and the correct efficiency is calculated. Finally, the correct plan is designed and tested by simulation and experiments.
引文
[1]姚文进.防空弹药指令修正弹[D].南京理工大学博士论文,2006
    [2]王昆.单兵简易控制火箭修正特性研究[D].南京理工大学硕士论文,1999
    [3]孟新宇.指令修正弹药组合测量和自装定[D].南京理工大学博士论文,2008
    [4]贺彤蔚.弹载减旋机构设计[D].南京理工大学硕士论文,2008
    [5]闵俊峰.利用地磁场确定弹道修正弹的滚转姿态[D].南京理工大学,2006
    [6]王俊全.弹道修正弹舵机的结构设计和控制[D].南京理工大学硕士论文,2004
    [7]毛晓翔.弹道修正弹总体方案及关键技术研究.南京理工大学硕士论文,2005
    [8]王涛.脉冲式弹道修正弹的修正方法研究.南京理工大学硕士论文,2005
    [9]陈勇巍,程翔.弹道修正弹滚转角辨识系统模型与误差分析[J].探测与控制学报,2006,28(4):44-47,51
    [10]谭凤岗.弹道修正弹的概念研究[J].弹箭技术,1998,9(4):1-10
    [11]高敏,张强.弹道修正弹实际弹道探测技术综述[J].弹道学报,2003,15(1):87-92
    [12]赵金强,龙飞.弹道修正弹综述[J].制导与引信,2005,26(4):16-19
    [13]吴炎烜,范宁军.二维弹道修正引信总体方案和关键技术分析[J].战术导弹技术,2006(05)
    [14]洪元军,雷雅茹.国外弹道修正引信发展现状研究[J].探测与控制学报,2001,23(4):1-5
    [15]王雨时.我国引信技术发展刍议[J].国防技术基础,2005(12):31-34
    [16]陈科山,马宝华.一维弹道修正引信阻力器的研究现状分析及其设计原则探讨[J].探测与控制学报.2003,25(3).-24-29,2003,25(3):24-29
    [17]何光林,马宝华.一维弹道修正引信阻尼修正机构的改进设计与仿真[J].探测与控制学报,2005(1):1-4
    [18]徐劲祥,宋锦武.弹道修正弹末段脉冲推力控制研究[J].弹道学报,2005,17(2):19-23
    [19]陈科山,马宝华.迫弹一维弹道修正引信阻力器结构的空气阻力特性研究[J].北京理工大学学报,2004(6):477-480,491
    [20]王中原,王良明.修正弹道飞行稳定性分析[J].兵工学报,1998,19(4):298-300
    [21]王中原,艾东民.弹道修正中弹道诸元探测时间间隔的确定[J].弹道学报,2002,14(1):84-87
    [22]赵捍东,郭锡福.侧推矢量在提高火箭弹射击精度的技术研究[J].弹箭与制导学报,2005,25(2):72-74
    [23]桑保华,简易制导炸弹有控弹道的设计与研究[D],南京,南京理工大学博士论文,2006
    [24]孙磊,张河.具有弹道修正功能的引信设计[J].南京理工大学学报:自然科学版,2005(3):285-287
    [25]孙磊,张河.主备自动切换技术在弹道修正系统中的应用[J].战术导弹控制技术,2005(01):30-32
    [26]胡荣林,李兴国.遥控指令弹道修正弹外弹道最优预测模型[J].探测与控制学报,2006,28(4):39-43
    [27]吴金中.电视导引系统的视轴稳定跟踪技术[J].战术导弹技术,2002,(1):34-38
    [28]橱宜禾,周维真.成像跟踪技术导论.西安:西安电予科技大学出版杜,1991
    [29]陈明能,钟时俊.简述电视制导技术[J].空军装备,2003,(3):33-34.
    [30]刘松涛,周晓东.可见光电视/红外成像复合寻的制导技术[J].应用光学,2006,27(6):467-475
    [31]张中南,王富宾,李晓.发展中的红外成像制导技术[J].控制与制导,2006,(1)40-42
    [32]任宁.红外成像制导技术的发展[J].红外与激光工程,2007,36(增刊),99-102
    [33]赵强,刘隆和.红外成像制导及其目标背景特性分析[J].航天电子对抗,2006,1,27-29.
    [34]胡晓军,倪震明.毫米波/红外复合制导技术概述[J],制导与引信,2007,28(4),8-12.
    [35]胡体玲,李兴国.毫米波/红外复合导引头关键技术分析[J].激光与红外,2007,(2):1012103.
    [36]郭修煌.精确制导技术.北京:国防工业出版社,1999
    [37]宋爱群,黄元庆.CCD图像传感器的应用技术与发展趋势[J].电子测量与仪器学报,2007,增刊,301-304.
    [38]王庆有.CCD应用技术.天津大学出版社.2000
    [39]韩卫华,钱光第.CCD图像传感器新技术与发展方向[J].网络传播与技术,2007,(2):70-71.
    [40]张纯学.空地导弹红外导引头末制导技术[J].控制与制导,2004,(11):49-54.
    [41]钟坚.红外成像制导武器的装备与发展[J].船舶电子工程,2008,18(6):27-29.
    [42]高阳.成像制导中的图像预处理及目标识别技术研究[D].西安,西北工业大学博士论文,2006
    [43]刘岩.红外图像分析方法与目标识别研究[D].西北工业大学博士学位论文,2004
    [44]赵大政,张庆捷,李刚等.电视制导系统目标识别算法[J].弹道学报,20(1),95-97
    [45]祁载康.制导弹药技术.北京理工大学出版社,2002.
    [46]Michael S.L.Hollis.Preliminary Design of a Range Correction Module for an Artillery Shell.ARL-MR-298,1996.3
    [47]Michael S.L.Hollis.Design and Analysis of a Prototype Range Correction Device for a Mortar Projectile.ARL-MR-411,1998.8
    [48]Michael S.L.Hollis,Fred J.Brandon,Peter C.Muller.Design and Flight Test of a Prototype Range Control Module for an 81-mm Mortar.ARL-MR-463,1999.9
    [49]Michael S.L.Hollis,Fred J.Brandon.Design and Analysis of a Fuze-Configurable Range Correction Device for an Artillery Projectile.ARL-TR-2074,1999.12
    [50]A.Ziliani,C.Grignon,C.Trouillot,C.Jeannin.Diagnostic of the Behaviour of a Course Correction Ammunition During its Correction Phase.19th International Symposium of Ballistics.2001:479-487
    [51]王儒策,刘荣忠,苏玳,王晓鸣编.灵巧弹药的构造及作用.兵器工业出版社.2001.6
    [52]谭凤岗编译.采用鸭舵控制使先进炮弹射程增大和精度改进[J].弹箭技术.1998(1):15-22
    [53]曹鲜花等.Matlab在光纤陀螺随机漂移建模中的应用[M].航天控制,25(1):18-20,2007
    [54]马宝华.战争、技术与引信——关于引信及引信技术的发展[J].探测与控制学报.2001.23(1):1-6
    [55]《低伸弹道表》编制组.低伸弹道表.北京:国防工业出版社,1977
    [56]王小健.单兵反坦克火箭弹总体参量优化设计.南京理工大学硕士论文,1987
    [57]王连敏.低速低伸弹道解法[J].兵工学报,1981(1)
    [58]朱福亚.火箭弹构造与作用.北京:国防工业出版社,2005
    [59]徐明友.火箭外弹道.北京:北京航空航天大学出版社,2005
    [60]浦发.外弹道学.北京:国防工业出版社,1980
    [61]杨启仁.子母弹飞行动力学.北京:国防工业出版社,1999
    [62]许波,刘征.Matlab工程数学应用.清华大学出版社,2002,12.
    [63]郑阿奇.MATLAB实用教程.北京:电子工业出版社,2007
    [64]贾秋玲,袁冬.基于MATLAB 7.x/Simulink/Stateflow系统仿真、分析及设计.西安:西北工业大学出版社,2006.
    [65]周长省等.火箭弹设计理论.北京:北京理工大学出版社.2005
    [66]徐劲祥.靶场外弹道实验弹丸速度测定研究[J].战术导弹控制技术.2006,(2):59-61.
    [67]潘水泮.武器系统效力分析.北京:国防工业出版社,1985.
    [68]秦永元.惯性导航.科学出版社,2007,4
    [69]张学孚等.磁通门技术.北京:国防工业出版社,1995.
    [70]张淑琴.空间交会对接测量技术及工程应用.北京:中国宇航出版社,2006
    [71](苏)卡拉什尼科夫著.地磁及其运用.科学出版社.1956.6
    [72](苏)B.M扬托夫斯基.地磁学.地质出版社.1959
    [73]敏俊峰.利用地磁场确定弹道修正弹的滚转姿态
    [74]US/UK World Magnetic Model—Epoch 2005.2005
    [75]彭青玉.木星一上星边缘的椭圆拟合[J].云南天文台台刊,2003(4):43-48.
    [76]闫蓓.基于最小二乘法的椭圆拟合改进算法[J].北京航空航天大学学报,2008(3):295-297
    [77]Ahn SJ,Rauh W.Least-squares orthogonal distances fitting of circle,sphere,ellipse,hyperbola,and parabola[J].Pattern Recognition,2001,34(12):2283-2303
    [78]Gander W,GolubGH,StrebelR,Least squares fitting of circles and ellipses[J].BIT,1994,34:558-578
    [79]Spath H,Least squares fitting with spheres[J].Theory Application,1998,196(1):191-199
    [80]Spath H,Orthogonal distance fitting by circles and ellipses with given area.Computing Stat,1997,12:343-354
    [81]Yuen H K,Kitter IJ.Detecting partially occluded ellipses using the hough transform[J].Image Vision and Computing,1989,7(1):31-37
    [82]潘泉,杨峰.一类非线性滤波器——UKF综述[J].控制与决策,2005(5):481-489,494
    [83]Julier SJ,Uhlmann JK:Unscented filtering and nonlinear estimation[A],Proceedings of the IEEE,Institute of Electrical and Electronics Engineers Inc.,Piscataway,NJ 08855-1331,United States,2004,pp 401-422
    [84]Wan EA,Merwe Rvd:The unscented Kalman filter,Kalman Filtering and Neural Networks.Edited by Simon.Haykin,Wiley Publishing,2001
    [85]李博纳.概率论与数理统计.北京:清华大学出版社,2006
    [86]夏箐.椭圆拟合方法的比较[D].暨南大学硕士论文,2007
    [87]陈勇巍,程翔.弹道修正弹滚转角辨识系统模型与误差分析[J].探测与控制学报,2006,28(4):44-47,51
    [88]孟新宇等.弹道修正弹滚转角辨识和误差分析[J],计算机仿真,2008
    [89]方建芳.反直升机智能雷伺服机构设计[D],南京理工大学硕士论文,2008
    [90]曹红松,冯顺山.地磁陀螺组合弹药姿态探测技术研究[J].弹箭与制导学报,2006,26(3):142-145
    [91]K.马格努斯.陀螺理论与应用.北京:国防工业出版社,1983.
    [92]B.A.鲍德涅尔.飞行器的稳定与自动驾驶仪.北京:国防工业出版社,1965
    [93]崔敏.小型飞行体姿态测试研究.中北大学硕士论文,2005
    [94]Gebre-Egziabber D,Elkaim GH.Gyro-free quaternion-based attitude determination system suitable for implementation using low cost sensors[J].Record-IEEE PLANS,Position Location and Navigation Symposium,2000:185-192
    [95]李文亮.四元数矩阵.长沙,国防科技大学,2002年6月,
    [100]曹红松.地磁/GPS弹道探测技术研究[D].中北大学博士论文,2004
    [101]韩子鹏.外弹道学.北京:国防工业出版社,2000.
    [102]钱杏芳.导弹飞行力学.北京理工大学出版社,2003
    [103]刘松涛,周晓东.图像融合技术研究的最新进展[J].激光与红外,2006,36(8),627-631
    [104]Gabarda Salvador,Gabriel Gristo bal.Multfocus imagefusion through pseudo—Wigner distribution[J].OpticalEngineering,2005,44(4):047001(1-9).
    [105]王年,范益政,鲍文霞等.基于图割的图像匹配算法[J].电子学报,34(2),232-236
    [106]苗常青,汪勃,付梦印等.电视图像目标实时分割与识别算法[J].北京理工大学学报,2005,25(9),786-790
    [107]张鹏飞,朱鸿翔,王萧等.几种扫描搜索规律的对比研究[J].弹箭与制导学报,2005,25(2),45-47
    [108]张会凌,李战存.二次曲线.兰州,兰州大学出版社,2006
    [109]梁晋文,陈林才,何贡.误差理论与数据处理.北京,中国计量出版社,2001.
    [110]Set/Reset Pulse Circuits for Magnetic Sensors.Honeywell Application Note.2001
    [111]熊仲华,刘运生.地磁观测技术.北京,地震出版社,1997年10月,
    [112]倪永生.地磁学简明教程.北京,地震出版社,1990年8月,
    [113]徐文耀.地磁学.北京,地震出版社,2003年04月,
    [114]张红顶.法国研究磁场制导技术[J].现代军事,2005(6):11
    [115]冯顺山.地磁陀螺组合弹药姿态探测技术研究[J].弹箭与制导学报,2006,26(3):142-145
    [116]曹红松,陈国光.在利用地磁探测确定弹体滚转姿态时的使用域分析[J].弹箭与制导学报,2005,25(2):66-68
    [117]陈勇巍,程翔.弹道修正弹滚转角辨识系统模型与误差分析[J].探测与控制学报,2006,28(4):44-47,51
    [118]臧爱云,原魁.磁传感器在导航系统中的应用[J].高技术通讯.2004,14(10).-53-57,
    [119]王平,李铁寿.磁控小卫星周期时变PD控制设计方法[J].控制工程,1999(3):11-15,10
    [120]朱建丰,徐世杰.基于地磁场测量估计卫星姿态的UKF算法[J].宇航学报,2006,27(6)
    [121]李杨,涂歆滢.星体剩磁矩对姿控系统影响的仿真研究[J].系统仿真学报,2002,14(2)
    [122]田菁,吴美平.三轴磁强计姿态确定[J].国防科技大学学报,2001,23(5):17-21
    [123]黄琳,荆武兴.关于剩磁对卫星姿态确定与控制影响的研究[J].航空学报,2006,27(3)
    [124]Lee D-J,Park K-J.Sequential attitude estimation using particle filters[A].South Lake Tahoe,CA,United States,Univelt Inc.San Diego,CA92198,United States,2006,pp239-255
    [125]Crassidis JL,Markley FL.A Minimum model error approach for attitude estimation[J].Journal of Guidance,Control,and Dynamics,1997,20(6):1241-1247
    [126]de La Parra S,Angel J.Low cost navigation system for UAV's[J].Aerospace Science and Technology,2005,9(6):504-516
    [127]Jiang X-Y,Ma G-F.Gyroless satellite attitude estimation for square root unscented kalman filter[J].Nanjing Li Gong Daxue Xuebao/Journal of Nanjing University of Science and Technology,2005,29(4):399-402
    [128]De Ruiter AHJ,Damaren CJ.Extended Kalman filtering and nonlinear predictive filtering for spacecraft attitude determination[J].Canadian Aeronautics and Space Journal,2002,48(1):13-23
    [129]Comisel H,Forster M.Attitude estimation of a near-Earth satellite using magnetometer data[J].Acta Astronautica,1997,40(11):781-788
    [130]Basile G,Pierantoni M.:Attitude compensated electronic compass for aircraft navigation[A],1994,pp 126-131
    [131]Yuen H K,Kitter IJ.Detecting partially occluded ellipses using the hough transform[J].Image Vision and Computing,1989,7(1):31-37
    [132]杨忠根,任蕾.基于简化GEVD技术的椭圆拟合算法[J].模式识别与人工智能,2004,17(04):479-485
    [133]杨忠根,马彦.使用广义正交概念的K-RANSAC椭圆提取[J].自动化学报,2002,28(04):520-526
    [134]薛程,王士同.椭圆检测算法的比较[J].江南大学学报:自然科学版,2006,5(6):674-678
    [135]陈凯,刘青.一种随机化的椭圆拟合方法[J].计算机工程与科学,2005,27(06):48-49,80
    [136]Ahn SJ,Rauh W.Least-squares orthogonal distances fitting of circle,sphere,ellipse,hyperbola,and parabola[J].Pattern Recognition,2001,34(12):2283-2303
    [137]Rosin PL.Note on the least squares fitting of ellipses[J].Pattern Recognition Letters,1993,14(10):799
    [138]Rosin PL,West GAW.Nonparametric segmentation of curves into various representations[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,1995,17(12):1140-1153
    [139]陈基伟.椭圆直接拟合算法研究[J].工程勘察,2007(06):49-51
    [140]任蕾,杨忠根.基于子空间投影技术的椭圆拟合算法[J].上海海事大学学报:文理综合版,2006,27(03):90-94
    [141]刘书桂,李蓬.基于最小二乘原理的平面任意位置椭圆的评价[J].计量学报,2002(04)
    [142]刘延柱.高等动力学,高等教育出版社,2001
    [143]Fitzgibbon A,Pilu M.Direct least square fitting of ellipses[J].Pattern Analysis and Machine Intelligence,IEEE Transactions on,1999,21(5):476-480
    [144]Taubin G.Estimation of planar curves,surfaces,and nonplanar space curves defined by implicit equations with applications to edge and range image segmentation[J].Pattern Analysis and Machine Intelligence,IEEE Transactions on,1991,13(11):1115-1138
    [145]Bookstein FL:Fitting Conic Sections to Scattered Data[A],Computer Graphics and Image Processing,no.9,1979,pp 56-71
    [146]Caruso MJ:Applications of magnetic sensors for low cost compass systems[A],Position Location and Navigation Symposium,IEEE 2000,2000,pp 177-184
    [147]M.Moulin,J.C.Guodon.Process for compensating the magnetic disturbances in
    [148]Wan EA,Merwe Rvd:The unscented Kalman filter,Kalman Filtering and Neural Networks.Edited by Simon.Haykin,Wiley Publishing,2001
    [149]Changey S,Beauvois D.:A mixed extended-unscented filter for attitude estimation with magnetometer sensor[A].American Control Conference,2006,2006,pp 6 pp.
    [150]Yuamin W,Wu M.Unscented Kalman filtering for additive noise case:
    [151]Michael S.L.Hollis.Preliminary Design of a Range Correction Module for an Artillery Shell.ARL-MR-298,1996.3
    [152]Michael S.L.Hollis.Design and Analysis of a Prototype Range Correction Device for a Mortar Projectile.ARL-MR-411,1998.8.
    [153]Honeywell Magnetic Sensor Products user manual.Honeywell Application Note.2001

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700