CN103810358A - Technology for evaluating efficiency of constellation satellite communication system - Google Patents

Technology for evaluating efficiency of constellation satellite communication system Download PDF

Info

Publication number
CN103810358A
CN103810358A CN201210438777.7A CN201210438777A CN103810358A CN 103810358 A CN103810358 A CN 103810358A CN 201210438777 A CN201210438777 A CN 201210438777A CN 103810358 A CN103810358 A CN 103810358A
Authority
CN
China
Prior art keywords
satellite communication
communication system
constellation satellite
effectiveness
constellation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201210438777.7A
Other languages
Chinese (zh)
Inventor
耿亮
田卫平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
XI'AN YUANSHUO SCIENCE & TECHNOLOGY Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by XI'AN YUANSHUO SCIENCE & TECHNOLOGY Co Ltd filed Critical XI'AN YUANSHUO SCIENCE & TECHNOLOGY Co Ltd
Priority to CN201210438777.7A priority Critical patent/CN103810358A/en
Publication of CN103810358A publication Critical patent/CN103810358A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radio Relay Systems (AREA)

Abstract

In order to directly and visually evaluate efficiency of a constellation satellite communication system, and select a constellation system scheme which can best meet actual requirements and long-term development requirements to provide powerful decision support, the invention discloses a technology for evaluating the efficiency of the constellation satellite communication system. The technology includes the steps that firstly, an efficiency evaluation index system of the constellation satellite communication system is determined through a Delphi expert consultation method; secondly, weights of efficiency evaluation indexes of the constellation satellite communication system are determined through an analytic hierarchy process; eventually, the efficiency of the constellation satellite communication system is evaluated through an ADC efficiency evaluation model.

Description

Constellation satellite communication system measures of effectiveness technology
Technical field
The present invention relates to constellation satellite communication system measures of effectiveness technology.
Background technology
In recent years, the development of terrestrial cellular wireless network is swift and violent, makes people's information interchange more convenient.But due to the impact of the factors such as geographical environment, China has the ground network in a lot of areas not yet to cover maybe and cannot cover, constellation satellite communication system can conveniently realize the seamless covering of earth surface, improve rapidly these regional communication conditions with its advantageous advantage, without the huge input of ground communication facilities.
Constellation satellite communication system is a complex structure, invests huge large scale system, and that its conceptual design comprises is interrelated, the many factors of restriction each other.From different angles, consider the impact of many factors on the practical usefulness of constellation satellite communication system, certainly will make scheme Design be full of diversity.Different design proposals makes to have systems technology feature separately, and effect is different, and How to choose can practical requirement, meets the constellation systems scheme of long term growth, avoids as far as possible great investment error, is the purpose of Effectiveness Evaluation.Under this background, the present invention proposes constellation satellite communication system measures of effectiveness technology.
ADC model is that American industry circle Weapons Systems Effectiveness advisory committee (WSEIAC) sets up, and obtains widespread use in a lot of fields.ADC model generalization taking into account system can be unavailable, can be unreliable while using and the three large key elements that can satisfy the demands, be a kind of comprehensively, objective synthetic effectiveness evaluation method.Therefore, the evaluation method of the present invention using ADC model as constellation satellite communication system measures of effectiveness.
Summary of the invention
The technical solution adopted in the present invention is: constellation satellite communication system measures of effectiveness technology, its step is as follows.
Step 1: set up constellation satellite communication system assessment indicator system.
By consulting pertinent literature, and sum up the influence factor that affects constellation satellite communication system, organize the authoritative expert of this respect, adopt Delphi expert consulting method to determine constellation satellite communication system assessment indicator system.
Step 2: the weight of determining constellation satellite communication system evaluation index.
Because the structure of constellation satellite communication system assessment indicator system is a multi-level structure, and analytical hierarchy process is processed the algorithms most in use of evaluation index weight in this class architecture just.Therefore, the present invention adopts analytical hierarchy process to determine the weight of constellation satellite communication system evaluation index.
Step 3: data acquisition.
In conjunction with constellation satellite communication system measures of effectiveness index, carry out data acquisition.
Step 4: adopt ADC model to carry out constellation satellite communication system measures of effectiveness.
ADC model formation is
E=ADC
Wherein, A is availability, the probability of expression system different conditions of living in the time starting to execute the task, and computing formula is as follows:
A=[a 11-a 1], wherein a 1=MTBF/ (MTBF+MTTR).
Here, MTBF is the average valuation of system interval time, and MTTR is mean repair time.
D is reliability matrix, represents known system residing state in the time starting working.If satellite mean failure rate probability is λ, spanning tree threshold value is T, and reliability matrix D can be expressed as:
D = 1 - λ ( S - T L ) λ ( S - T L ) 0 1
Here,
Figure BDA0000236316662
, L=S/N, N is system number.
C is ability matrix, represents that known system when residing state, completes measuring of assignment of mission in the process of executing the task, and computing formula is as follows:
C=[c 1c 2] t, c 1, c 2represent respectively the probability that under normal and malfunction, system is finished the work, generally task is nonserviceabled lower system cannot finish the work time, c 2=0.While assessment due to constellation satellite communication system, relate to index many, easily cause the calculation of complex of C.Therefore, in the present invention, index is carried out to linear weighted function, using the index attribute value of the reflection entire system ability obtaining as ability Matrix C.
Here, the index attribute value of entire system ability is:
c i = Σ j = 1 n w j x ij
Wherein, w jbe j (j=1,2 ..., the n) weight of individual constellation satellite communication system evaluation index, x ijbe i constellation satellite communication system about j (j=1,2 ..., the n) desired value of individual evaluation index.
The invention has the beneficial effects as follows: in conjunction with correlation theory research, proposed constellation satellite communication system assessment technology, through at practice test, shown that the method is reasonable in design, practical, for the assessment of constellation satellite communication system provides good method.Study constellation satellite communication system evaluation index system, analyzed, studied the method for constellation satellite communication system evaluation index weight, for constellation satellite communication planning provides theoretical reference.
Accompanying drawing explanation
Fig. 1 is process flow diagram of the present invention.
Embodiment
As described in Figure 1, constellation satellite communication system measures of effectiveness technology mainly comprises 6 steps: 1) tissue expert, adopts Delphi expert consulting method to determine constellation satellite communication system measures of effectiveness index system; 2) application level analytic approach is determined the weight of constellation satellite communication system measures of effectiveness index; 3) gather related data; 4) through ADC Effectiveness Evaluation Model, constellation satellite communication system usefulness is assessed.Idiographic flow is as follows:
Step 1: set up constellation satellite communication system assessment indicator system.
Organize the expert of constellation satellite communication system aspect, adopt Delphi expert consulting method to set up the assessment indicator system of constellation satellite communication system usefulness, in table 1.
Table 1 constellation satellite communication system index system
Figure BDA0000236316664
Step 2: the weight of determining constellation satellite communication system evaluation index.
Because the structure of constellation satellite communication system assessment indicator system is a multi-level structure, and analytical hierarchy process is processed the algorithm of evaluation index weight in this class architecture just.Therefore, the present invention adopts analytical hierarchy process to determine the weight of constellation satellite communication system evaluation index, and detailed process is as follows.
Step 2.1: expert seeks advice to multidigit, the data of obtaining are weighted on average, determine the judgment matrix of each layer, and result is as follows.
1) judgment matrix of rule layer
U = 1 3 4 6 1 / 3 1 1 / 3 3 1 / 4 3 1 3 1 / 6 1 / 3 1 / 3 1
2) judgment matrix of indicator layer, wherein covers, the judgment matrix of access, quality of service and system security is respectively matrix B, C, D, E.
U 1 = 1 3 5 1 / 3 1 1 / 2 1 / 5 2 1 U 2 = 1 2 5 1 / 3 1 / 2 1 3 1 / 5 1 / 5 1 / 3 1 1 / 7 3 5 7 1 U 3 = 1 5 1 / 3 3 1 / 5 1 1 / 7 1 / 2 3 7 1 7 1 / 3 2 1 / 7 1 U 4 = 1 1 / 3 3 2 3 1 4 3 1 / 3 1 / 4 1 1 / 2 1 / 2 1 / 3 2 1
Step 2.2: calculate above-mentioned judgment matrix U, U 1, U 2, U 3, U 4maximum characteristic root and corresponding maximal eigenvector, they are normalized respectively, and carry out consistency check, determine the weight of constellation satellite communication system evaluation index, it the results are shown in Table 2.
Table 2 constellation satellite communication system index weights allocation result
Figure BDA00002363166610
Step 3: data acquisition
Because constellation satellite communication system measures of effectiveness index system comprises quantitative target and qualitative index, therefore adopt diverse ways image data for different indexs, concrete evaluation index acquisition mode is in table 3.
Table 3 achievement data acquisition mode
Figure BDA00002363166611
Data after gathering are normalized.
Step 4: adopt ADC model to carry out constellation satellite communication system measures of effectiveness.
ADC model formation is
E=ADC (1)
Wherein, A is availability, the probability of expression system different conditions of living in the time starting to execute the task, and computing formula is as follows:
A=[a 1 1-a 1] (2)
Here a, 1=MTBF/ (MTBF+MTTR).
Wherein, MTBF is the average valuation of system interval time, and MTTR is mean repair time.
D is reliability matrix, represents known system residing state in the time starting working.If satellite mean failure rate probability is λ, spanning tree threshold value is T, and reliability matrix D can be expressed as:
D = 1 - λ ( S - T L ) λ ( S - T L ) 0 1 - - - ( 3 )
Here,
Figure BDA00002363166613
, L=S/N, N is system number.
C is ability matrix, represents that known system when residing state, completes measuring of assignment of mission in the process of executing the task, and computing formula is as follows:
C=[c 1c 2] t, c 1, c 2represent respectively the probability that under normal and malfunction, system is finished the work, generally task is nonserviceabled lower system cannot finish the work time, c 2=0.While assessment due to constellation satellite communication system, relate to index many, easily cause the calculation of complex of C.Therefore, in the present invention, index is carried out to linear weighted function, using the index attribute value of the reflection entire system ability obtaining as ability Matrix C.
Here, the index attribute value of entire system ability is:
c i = Σ j = 1 n w j x ij - - - ( 4 )
Wherein, w jbe j (j=1,2 ..., the n) weight of individual constellation satellite communication system evaluation index, x ijbe i constellation satellite communication system about j (j=1,2 ..., the n) desired value of individual evaluation index.
In order to verify the performance of algorithm of the present invention, IRIDIUM (Iridium) and GlobalStar (Globalstar) are assessed.Wherein, table 4 is the normalization evaluation index value of these two constellation satellite systems.
Table 4 evaluation index value
Index Iridium Globalstar
Single star coverage rate 1.0 0.83
Link mean down time 0.63 0.72
The minimum elevation angle that covers 0.85 0.8
Access capacity 0.81 0.76
Access rate 0.7 0.8
Average turn-on time 0.7 0.85
Access blocking rate 0.77 0.85
Transmission delay 0.7 0.75
Cutting off rate 0.75 0.85
Packet loss 0.8 0.82
The bit error rate 0.8 0.7
Maintainable 0.55 0.65
Anti-interference 0.6 0.75
Information encryption 0.75 0.65
Safety certification 0.8 0.7
On this basis, application ADC Effectiveness Evaluation Model is assessed Iridium system and Globalstar system effectiveness.As an example, the present invention provide Iridium system about in ADC Effectiveness Evaluation Model about parameter A, D, the analysis of C, specific as follows:
1) acquisition of parameter A
The mean time between failures (MTBF) of Iridium system is 10 years, and mean time between maintenance MTTR is 0.15, and application of formula (2) obtains the availability matrix of Iridium and is: [0.985,0.015].
2) acquisition of parameter D
By STK emulation, constantly reduce number of satellite, can obtain data shown in table 5:
Table 5 Iridium system coverage rate and satellite damage number relation table
Figure BDA00002363166615
Obviously, in the time that satellite damage reaches 12 stars, single star covers to reduce to 89.66%, approximately there is 10% possibility not have satellite to cover, system state now cannot meet the properties value measuring under normal operating conditions, thereby can, using this as thresholding, be about 9.6% according to rail satellite average year valuation rate.Computing system state transition probability, the result of calculation of Iridium system state transition probability is:
ρ = ( 0.096 ) 12 / 6 = 9.2 × 10 - 3
According to the reliable matrix D of the known Iridium system of formula (3) be
D = 0.9908 0.0092 0 1
3) acquisition of parameters C
In associative list 2, about the index evaluation value of Iridium and Globalstar in the weight of constellation satellite communication system measures of effectiveness index and table 4, should formula (4) obtain reflecting the index attribute value P of entire system ability, have:
P Iridium = 0.8011 , P Globalstar = 0.7756
The ability Matrix C of Iridium system is:
C = 0.8011 0.7756
Application of formula (1), determines that the measures of effectiveness value of Iridium system is: 0.8005.
Equally, according to the desired value about Globalstar system in table 4, application ADC Effectiveness Evaluation Model can determine that the measures of effectiveness value of Globalstar system is 0.769.
Therefore, can find out that through the present invention IRIDIUM comprehensive effectiveness will be higher than GlobalStar comprehensive effectiveness, this result meets generally acknowledges the fact.

Claims (2)

1. constellation satellite communication system measures of effectiveness technology, its concrete steps are as follows: 1) application Delphi expert consulting method is set up constellation satellite communication system measures of effectiveness index system; 2) application level analytic approach is determined constellation satellite communication system measures of effectiveness index weights; 3) application ADC Effectiveness Evaluation Model assessment constellation satellite communication system usefulness.
2. constellation satellite communication system measures of effectiveness technology according to claim 1, its feature, in the feature in conjunction with constellation satellite communication system, has been set up ADC Effectiveness Evaluation Model, that is:
Figure 2012104387777100001DEST_PATH_IMAGE002
Wherein, 1) A is availability, the probability of expression system different conditions of living in the time starting to execute the task, and computing formula is as follows:
Figure 2012104387777100001DEST_PATH_IMAGE004
Here, MTBF is the average valuation of system interval time, and MTTR is mean repair time;
2) D is reliability matrix, represents known system residing state in the time starting working, and establishes satellite mean failure rate probability and is
Figure 2012104387777100001DEST_PATH_IMAGE006
, spanning tree threshold value is T, reliability matrix D can be expressed as:
Figure 2012104387777100001DEST_PATH_IMAGE008
Here,
Figure 2012104387777100001DEST_PATH_IMAGE010
, N is system number;
3) C is ability matrix, represents that known system when residing state, completes measuring of assignment of mission in the process of executing the task, and computing formula is as follows:
Figure 2012104387777100001DEST_PATH_IMAGE012
Wherein
Figure 2012104387777100001DEST_PATH_IMAGE014
represent respectively the probability that under normal and malfunction, system is finished the work, adopt the index attribute value of reflection entire system ability as ability Matrix C,
Figure 2012104387777100001DEST_PATH_IMAGE016
Here,
Figure 2012104387777100001DEST_PATH_IMAGE018
be the weight of j constellation satellite communication system evaluation index, be the desired value of i constellation satellite communication system about j evaluation index.
CN201210438777.7A 2012-11-06 2012-11-06 Technology for evaluating efficiency of constellation satellite communication system Pending CN103810358A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210438777.7A CN103810358A (en) 2012-11-06 2012-11-06 Technology for evaluating efficiency of constellation satellite communication system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210438777.7A CN103810358A (en) 2012-11-06 2012-11-06 Technology for evaluating efficiency of constellation satellite communication system

Publications (1)

Publication Number Publication Date
CN103810358A true CN103810358A (en) 2014-05-21

Family

ID=50707122

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210438777.7A Pending CN103810358A (en) 2012-11-06 2012-11-06 Technology for evaluating efficiency of constellation satellite communication system

Country Status (1)

Country Link
CN (1) CN103810358A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104158609A (en) * 2014-08-13 2014-11-19 中国人民解放军理工大学 Satellite communication network system performance evaluation method
CN104732079A (en) * 2015-03-17 2015-06-24 中国空间技术研究院 Effectiveness determining method for remote sensing satellite system
CN106940446A (en) * 2017-03-15 2017-07-11 北京航空航天大学 A kind of satellite navigation system usability evaluation method and device
CN112149958A (en) * 2020-08-24 2020-12-29 成都天奥集团有限公司 Method for calculating multi-dimensional efficiency evaluation index of space-ground base station network resources
CN113052431A (en) * 2021-02-19 2021-06-29 中国人民解放军战略支援部队航天工程大学 Effectiveness evaluation method of aerospace test and launch system based on ADCP model
CN113189623A (en) * 2021-04-21 2021-07-30 中国电子科技集团公司第五十四研究所 Satellite navigation system signal quality grade evaluation method
CN117811961A (en) * 2024-03-01 2024-04-02 南京航空航天大学 Heterogeneous constellation observation efficiency evaluation method facing to offshore moving target

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101944160A (en) * 2010-08-31 2011-01-12 环境保护部华南环境科学研究所 Immediate offshore area ecological environment comprehensive evaluation method based on analytic hierarchy process and comprehensive evaluation method
US20120203369A1 (en) * 2011-02-04 2012-08-09 International Business Machines Corporation Manufacturing execution system (mes) including a wafer sampling engine (wse) for a semiconductor manufacturing process

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101944160A (en) * 2010-08-31 2011-01-12 环境保护部华南环境科学研究所 Immediate offshore area ecological environment comprehensive evaluation method based on analytic hierarchy process and comprehensive evaluation method
US20120203369A1 (en) * 2011-02-04 2012-08-09 International Business Machines Corporation Manufacturing execution system (mes) including a wafer sampling engine (wse) for a semiconductor manufacturing process

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨海涛: "星座卫星通信***效能评估研究", 《中国优秀硕士学位论文全文数据库 信息科技辑 》 *
王宇: "基于ADC模型的星座卫星通信***效能评估技术研究", 《中国优秀硕士学位论文全文数据库 经济与管理科学辑》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104158609A (en) * 2014-08-13 2014-11-19 中国人民解放军理工大学 Satellite communication network system performance evaluation method
CN104732079A (en) * 2015-03-17 2015-06-24 中国空间技术研究院 Effectiveness determining method for remote sensing satellite system
CN104732079B (en) * 2015-03-17 2017-06-27 中国空间技术研究院 A kind of remote sensing satellite System Effectiveness determines method
CN106940446A (en) * 2017-03-15 2017-07-11 北京航空航天大学 A kind of satellite navigation system usability evaluation method and device
CN106940446B (en) * 2017-03-15 2019-11-12 北京航空航天大学 A kind of satellite navigation system usability evaluation method and device
CN112149958A (en) * 2020-08-24 2020-12-29 成都天奥集团有限公司 Method for calculating multi-dimensional efficiency evaluation index of space-ground base station network resources
CN112149958B (en) * 2020-08-24 2024-02-13 成都天奥集团有限公司 Method for calculating multidimensional performance evaluation index of antenna base station network resources
CN113052431A (en) * 2021-02-19 2021-06-29 中国人民解放军战略支援部队航天工程大学 Effectiveness evaluation method of aerospace test and launch system based on ADCP model
CN113052431B (en) * 2021-02-19 2022-08-12 中国人民解放军战略支援部队航天工程大学 Effectiveness evaluation method of aerospace test and launch system based on ADCP model
CN113189623A (en) * 2021-04-21 2021-07-30 中国电子科技集团公司第五十四研究所 Satellite navigation system signal quality grade evaluation method
CN113189623B (en) * 2021-04-21 2022-08-12 中国电子科技集团公司第五十四研究所 Satellite navigation system signal quality grade evaluation method
CN117811961A (en) * 2024-03-01 2024-04-02 南京航空航天大学 Heterogeneous constellation observation efficiency evaluation method facing to offshore moving target

Similar Documents

Publication Publication Date Title
CN103810358A (en) Technology for evaluating efficiency of constellation satellite communication system
Solangi et al. Assessing the solar PV power project site selection in Pakistan: based on AHP-fuzzy VIKOR approach
Santos et al. Regional frequency analysis of droughts in Portugal
Fan et al. Reshaping the sustainable geographical pattern: A major function zoning model and its applications in China
AU2006324395B2 (en) Location determination of power system disturbances based on frequency responses of the system
CN108280553A (en) Regional Torrent Risk Zonation based on GIS- Artificial neural network ensembles and prediction technique
Du et al. Optimal sensor placement and measurement of wind for water quality studies in urban reservoirs
Boland Spatial-temporal forecasting of solar radiation
Guo et al. Spatial and temporal variation of extreme precipitation indices in the Yangtze River basin, China
CN104376413A (en) Power grid planning scheme evaluation system based on analytic hierarchy process and data envelopment analysis
CN104517037A (en) Remote sensing estimation method of biological bearing capacity
CN103179582B (en) A kind of room subsystem site selecting method and system
CN109919514B (en) Transformer area line loss abnormity detection method and device based on line loss marker post threshold value
CN105222787A (en) Based on the location fingerprint base construction method of matrix fill-in
CN104732079A (en) Effectiveness determining method for remote sensing satellite system
CN106788815A (en) A kind of short wave communication reliability estimation method based on many system detection datas
Barbosa et al. Applying fuzzy logic to comparative distribution modelling: a case study with two sympatric amphibians
De Mello et al. Interpolation methods for improving the RUSLE R-factor mapping in Brazil
CN105930645A (en) Communication station equipment maintenance support capability assessment method based on principal component analysis
Masseran et al. On spatial analysis of wind energy potential in Malaysia
CN106202491A (en) A kind of based on seasonal effect in time series data digging method and system
CN111144637A (en) Regional power grid geological disaster forecasting model construction method based on machine learning
CN114723283A (en) Ecological bearing capacity remote sensing evaluation method and device for urban group
CN115330086A (en) Photovoltaic land automatic site selection method under multi-condition constraint of spatial big data
CN104050254B (en) The method that feature database in 3D rooms is built using house data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB02 Change of applicant information

Address after: 710077, block 13, building A, Jiayu building, No. 58, Kam Yip Road, Xi'an hi tech Zone, Shaanxi, China

Applicant after: XI'AN YUANSHUO SCIENCE & TECHNOLOGY CO., LTD.

Address before: 710077 Shaanxi city of Xi'an province high tech Zone Jinye Road No. 69 business R & D Park C District No. 1 gazelle Valley E room 503

Applicant before: XI'AN YUANSHUO SCIENCE & TECHNOLOGY CO., LTD.

CB03 Change of inventor or designer information

Inventor after: Yue Rong

Inventor before: Geng Liang

Inventor before: Tian Weiping

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: GENG LIANG TIAN WEIPING TO: YUE RONG

C41 Transfer of patent application or patent right or utility model
CB03 Change of inventor or designer information

Inventor after: Deng Zhenghong

Inventor before: Yue Rong

COR Change of bibliographic data
TA01 Transfer of patent application right

Effective date of registration: 20160128

Address after: Beilin District Shaanxi province Xi'an City friendship road 710072 No. 127 box 205

Applicant after: Northwestern Polytechnical University

Address before: 710077, block 13, building A, Jiayu building, No. 58, Kam Yip Road, Xi'an hi tech Zone, Shaanxi, China

Applicant before: XI'AN YUANSHUO SCIENCE & TECHNOLOGY CO., LTD.

WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140521