CN103200608B - A kind of interference matrix compensation method and device based on measurement report - Google Patents

A kind of interference matrix compensation method and device based on measurement report Download PDF

Info

Publication number
CN103200608B
CN103200608B CN201210005734.XA CN201210005734A CN103200608B CN 103200608 B CN103200608 B CN 103200608B CN 201210005734 A CN201210005734 A CN 201210005734A CN 103200608 B CN103200608 B CN 103200608B
Authority
CN
China
Prior art keywords
interference
community
interfered cell
sector
compensation value
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.)
Active
Application number
CN201210005734.XA
Other languages
Chinese (zh)
Other versions
CN103200608A (en
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.)
China Mobile Group Zhejiang Co Ltd
Original Assignee
China Mobile Group Zhejiang 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 China Mobile Group Zhejiang Co Ltd filed Critical China Mobile Group Zhejiang Co Ltd
Priority to CN201210005734.XA priority Critical patent/CN103200608B/en
Publication of CN103200608A publication Critical patent/CN103200608A/en
Application granted granted Critical
Publication of CN103200608B publication Critical patent/CN103200608B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

A kind of interference matrix compensation method and device based on measurement report proposed, set up original interference matrix according to basic data and measurement report MR, in all communities, search the identical BCCH community pair of You Yuan community and interfered cell composition, for each identical BCCH community pair, determine interference compensation value according to position relation level between described former community and described interfered cell.

Description

A kind of interference matrix compensation method and device based on measurement report
Technical field
The present invention relates to a kind of interference matrix compensation method and device, or rather, relate to a kind of based on honeybeeInterference matrix compensation method and the device of the report of nest wireless network measurement.
Background technology
Building interference matrix based on measurement report MR (MeasureReport) is current automatic frequencyThe basis of optimizing, the accuracy of interference matrix plays vital effect to automatic frequency division. Adopt at present emptyPlan adjacent area is carried out comprehensive BCCH (BCCH) measurement and is made MR interference matrix substantially improve standardReally, but mobile phone reports MR measurement data to have a defect, at current service cell BCCH and itsIn the identical situation of its interfered cell BCCH, may interfered cell not reported to MR measurement report, byThis in MR measurement report, reporting measurement reports not of the community that in existing network, all BCCH are identical, because ofJust there is not mutual interference relationships, these identical BCCH community correspondences in the interference matrix of this foundation yetPoint be called blind spot.
The situation of blind spot on the impact of the dividing method based on MR is at present: if the BCCH before changing frequentlyPlan more rationally, may not exist interference to affect each other less; If instead before changing frequentlyBCCH plans unreasonablely, and so last interference matrix may lack interference relationships each other,Affect thus the effect of automatic frequency division. In the situation that dense city channeling is more and more higher, and existing networkUnder the strategy that BCCH and TCH (Traffic Channel) segmentation are used, this blind spot loss effect will be byFlaring is large, thereby last frequency division is brought to direct impact, simultaneously probably originally irrational after frequency divisionBCCH and TCH planning still exist, thereby cannot obtain the maximum improvement of performance. Therefore need blindPoint carries out necessary compensation, sets up complete interference matrix.
Summary of the invention
Based on foregoing, inventor notices and need to improve the foundation of interference matrix. In brief,Through necessary compensation and set up after complete interference matrix, this method and apparatus will improve a point yupin effect.
Utilize existing interference matrix can not solve the problems referred to above and the needs of this area completely. Correspondingly,Propose embodiment a kind of interference matrix compensation method and device based on measurement report is provided, can overcome and be permittedThe defect of many or all above-mentioned this areas of mentioning.
Correspondingly, the present invention proposes a kind of interference matrix compensation method based on measurement report, according to basic numberAccording to setting up original interference matrix with measurement report MR, in all communities, search You Yuan community and disturb littleThe identical BCCH community pair of district's composition, for each identical BCCH community pair, according to described former communityAnd between described interfered cell, position relation level is determined interference compensation value.
And, wherein determine interference compensation value according to position relation between described former community and described interfered cellFurther comprise: if described former community is omnidirectional sector, by nearest Dui Yuan community, interfered cell 3The mean value of individual area interference value is as interference compensation value; If described former community is directed sector, disturb littleDistrict is omnidirectional sector, and according to interfered cell, to other, the interference value with base station cell is determined described interference compensationValue; If the sector of the sector of described interfered cell and described former community is all directed sector, described in basisThe position relation of all sectors of the base station, place, He Yuan community, sector of interfered cell, determines level of interference, rootAccording to described level of interference, determine interference compensation value;
And, wherein saidly set up original interference matrix according to basic data and measurement report MR and comprise: pressDivide with frequency and adjacent frequency interference threshold according to C/I, and calculate and drop on different interval measuring object sampled points alwaysRatio in sampled point, obtains disturbing numerical value with frequency and/or adjacent frequency thus.
The present invention simultaneously proposes a kind of device that the interference matrix of measurement report is compensated, and comprises according to basePlinth data and measurement report MR set up the module of setting up of original interference matrix, in all communities, search byWhat the identical BCCH community of former community and interfered cell composition was right searches module, for each identical BCCHCommunity pair, determines interference compensation value really according to position relation level between described former community and described interfered cellCover half piece, and set up the compensating module of complete interference matrix according to described definite interference compensation value.
And wherein said determination module is further carried out following operation: if described former community is omnidirectional's fanDistrict, using the mean value of 3 nearest Dui Yuan community, interfered cell area interference values as interference compensation value;If described former community is directed sector, interfered cell is omnidirectional sector, same to other according to interfered cellThe interference value of base station cell is determined described interference compensation value; If the sector of described interfered cell and described former littleThe sector in district is all directed sector, owning according to the base station, place, He Yuan community, sector of described interfered cellThe position relation of sector, determines corresponding level of interference, according to described level of interference, determines interference compensation value;
And the wherein said module of setting up is divided with frequency and adjacent frequency interference threshold according to C/I, and calculates and drop onThe ratio of different interval measuring object sampled points in total sampled point, obtains with frequency and/or adjacent frequency dry thusDisturb numerical value.
In embodiment shown in the present, while utilizing described complete interference matrix to carry out frequency division, can obtainBe better than utilizing described original interference matrix to carry out the effect of frequency division.
Brief description of the drawings
Fig. 1 is the flow chart of the interference matrix compensation method based on measurement report;
Fig. 2 is the right flow chart in identical BCCH community of searching You Yuan community and interfered cell composition;
Fig. 3 is the interference matrix compensation arrangement block diagram based on measurement report;
Tu4Shi Yuan community and azimuth, interfered cell and position relation schematic diagram;
Fig. 5 calculates the dry of Dui Yuan community, interfered cell according to former community and the combination of interfered cell position relationDisturb the flow chart of offset.
Detailed description of the invention
Below in conjunction with accompanying drawing, embodiments of the present invention are described in detail.
Fig. 1 is the flow chart of the interference matrix compensation method based on measurement report. As shown in Figure 1, in stepS101, imports basic data and MR measurement report. Basic data can comprise for example base station latitude and longitude information,Azimuth information, angle of declination information, frequency configuration information, neighboring BS information, handover information etc., MR measuresReport comprises the statistical report of whole BCCH metrical informations of complete collection period, according to distinct deviceProducer defines different forms, metrical information wherein can comprise community, Measuring Time, measurement BCCH,Measure BCC, measure NCC, measuring object level Yu Zhe community level difference value is at different interval hits.Then, on the basic data importing and MR measurement report basis, set up original interference matrix. Specifically,Divide with frequency and adjacent frequency interference threshold according to C/I, then calculate and drop on different interval measuring object sampled pointsRatio in total sampled point, thus same, adjacent frequency interference numerical value obtained.
Subsequently, at step S102, in all communities, search the phase of You Yuan community and interfered cell compositionWith BCCH community pair. Search the identical BCCH community of You Yuan community and interfered cell composition right in detailProcess will be described later.
Then, at step S103, to searching each the identical BCCH community pair obtaining in step S102,Determine interference compensation value according to position relation between former community and interfered cell.
Then,, at step S104, set up complete interference matrix according to described definite interference compensation value.
Wherein, before setting up complete interference matrix according to described definite interference compensation value, can set up byCell basis information and original interference matrix form compensated information set jointly. Wherein cell basis information canComprise community CellID, community BCCH frequency, community longitude and latitude, azimuth, community, original interference squareBattle array comprises main plot, interfered cell, co-channel interference value, adjacent frequency interference value.
Fig. 2 is the right flow chart in identical BCCH community of searching You Yuan community and interfered cell composition. ?In Fig. 2, first in all communities, carry out community BCCH traversal, the i.e. BCCH of sequential search community.Subsequently, from the BCCH of the community of searching, find identical BCCH community to candidate. Community is then setTo between the distance factor thresholding row filter of going forward side by side, thereby find the community pair of identical BCCH.
Fig. 3 is the interference matrix compensation arrangement block diagram based on measurement report. In Fig. 3, interference matrix is mendedRepay device 300 and comprise the module of setting up of setting up original interference matrix according to basic data and measurement report MR301, right the looking in identical BCCH community of searching You Yuan community and interfered cell composition in all communitiesLook for module 302, for each identical BCCH community pair, according to described former community and described interfered cell itBetween position relation level determine the determination module 303 of interference compensation value, and according to described definite interference compensationValue is set up the compensating module 304 of complete interference matrix. In addition, wherein said determination module 303 is furtherCarry out following operation: if described former community is omnidirectional sector, by nearest Dui Yuan community, interfered cell 3The mean value of individual area interference value is as interference compensation value; If described former community is directed sector, disturb littleDistrict is omnidirectional sector, and according to interfered cell, to other, the interference value with base station cell is determined described interference compensationValue; If the sector of the sector of described interfered cell and described former community is all directed sector, described in basisThe position relation of all sectors of the base station, place, He Yuan community, sector of interfered cell, determines corresponding interference etc.Level, according to described level of interference, determines interference compensation value. The wherein said module 301 of setting up is drawn according to C/IDivide with frequency and adjacent frequency interference threshold, and calculate the measuring object sampled points that drop on different intervals in total sampled pointRatio, obtain thus with frequently and/or adjacent frequency disturb numerical value.
Tu4Shi Yuan community and azimuth, interfered cell and position relation schematic diagram. As shown in Figure 4, definition is dryDisturbing cell coordinate is: (x1,y1), the angle of interfered cell and standard coordinate axle horizontal direction is: β1(0≤β1<2π);Former cell coordinate is: (x2,y2), the angle of former community and standard coordinate axle horizontal direction is: β2(0≤β2<2π),The azimuth of definition relatively former community, interfered cell is α, by following formula computer azimuth angle α.
&alpha; = arctan y 2 - y 1 x 2 - x 1 ( x 1 < x 2 , y 1 < y 2 ) &pi; + arctan y 2 - y 1 x 2 - x 1 ( x 1 > x 2 ) 2 &pi; + arctan y 2 - y 1 x 2 - x 1 ( x 1 < x 2 , y 1 > y 2 ) &pi; / 2 ( y 2 > y 1 x 2 = x 1 ) 3 &pi; / 2 ( y 2 < y 1 , x 2 = x 1 ) 4
If (β2-α)(β1-α) >=0, He Yuan community in interfered cell is spatially relative, as in Fig. 4 14107Shown in 24124, be defined as " homonymy ";
If (β2-α)(β1-α) < 0, be defined as " heteropleural ";
If | &beta; 1 - &alpha; | < &pi; 2 | &pi; - | &beta; 2 - &alpha; | | < &pi; 2 , Be defined as " subtend ";
If | &beta; 1 - &alpha; | &GreaterEqual; &pi; 2 | &pi; - | &beta; 2 - &alpha; | | &GreaterEqual; &pi; 2 , Be defined as " dorsad ";
If | &beta; 1 - &alpha; | < &pi; 2 | &pi; - | &beta; 2 - &alpha; | | &GreaterEqual; &pi; 2 , Be defined as " greatly in the same way ";
If | &beta; 1 - &alpha; | &GreaterEqual; &pi; 2 | &pi; - | &beta; 2 - &alpha; | | < &pi; 2 , Be defined as " in the same way little ";
Draw following combination and level of interference according to above 6 kinds of position relations:
A Homonymy, relatively Very high
B Heteropleural, relatively High
C Homonymy, greatly in the same way In
D Homonymy is in the same way little In low
E Heteropleural, greatly in the same way Low
F Heteropleural is in the same way little Lower slightly
G Homonymy, dorsad Very low
H Heteropleural, dorsad Very low
Fig. 5 calculates the dry of Dui Yuan community, interfered cell according to former community and the combination of interfered cell position relationDisturb the flow chart of offset. In Fig. 5, start to compensate the identical of Dui Youyuan community and interfered cell compositionAfter the interference matrix of the right measurement report in BCCH community compensates, first at step S501, determine formerWhether community is directed sector. If the result of step S501 is "Yes", perform step S502, asFruit result is "No", performs step S506, also determines that former community is omnidirectional sector, and interfered cell isThe compensate for disturbances value of directed sector or omnidirectional sector. Specifically, if former community is omnidirectional sector, disturbCommunity is directed sector or omnidirectional sector, by the interference value of 3 nearest Dui Yuan community, interfered cell communitiesMean value be defined as interference compensation value. Then,, at step S502, determine whether interfered cell is directedSector. If the result of step S502 is "Yes", perform step S503, if result is "No",Perform step S507, also determine that former community is directed sector, interfered cell is the compensation of omnidirectional sectorInterference value. Specifically, if former community is directed sector in step S507, interfered cell is omnidirectionalSector, processes relatively according to position. That is to say, the in the situation that of relative in Yu Yuan community, interfered cell,Interference compensation value is defined as: max{I1,I2}(I1,I2Represent interfered cell interference with base station cell to otherValue, if do not have, is 0; ); In Yu Yuan community, interfered cell, dorsad in the situation that, interference compensation value isavg{I1,I2. Next,, at step S504, calculate other community, He Tong base station, Dui Yuan community, interfered cellLevel of interference. Specifically, if former community is directed sector, interfered cell is directed sector, rootAccording to the position relation of A-H noted earlier, judge the pass of all sectors, base station, place, He Yuan community, interfered cellSystem, draws corresponding level of interference. Then,, at step S505, carry out level of interference relatively and determine phaseThe interference compensation value of answering. Here the number of setting other community in same base station is 2, but is not limited to this. FixedJustice IntfOriFor the level of interference of Dui Yuan community, interfered cell, Intf1For other community 1, Dui Tong base station, interfered cellLevel of interference, Intf2For the level of interference of other community 2, Dui Tong base station, interfered cell. If satisfied condition:IntfOri>Intf1And IntfOri>Intf2, determine that interference compensation value is: sum{I1,I2}+0.001; If metCondition: IntfOri=Intf1And IntfOri>Intf2, or IntfOri>Intf1And IntfOri=Intf2, determine to disturb and mendThe value of repaying is: max{max{I1,I2}-0.001,0}; If satisfied condition: Intf2>IntfOri>Intf1OrIntf1>IntfOri>Intf2, interference compensation value is defined as to avg{I1,I2; If satisfied condition:IntfOri<Intf1And IntfOri<Intf2, interference compensation value is defined as:, byThis completes determining of level of interference under different situations. It is noted that above-mentioned 0.001 is exemplary,But those skilled in the art are known, this value is not limited to this. Then, all sectors are traveled through, calculateCorresponding " blind spot " area interference is related to offset, joins in whole interference matrix, finally forms wholeThe interference matrix that network is complete, thus, interference relationships is more rationally reliable concerning whole automatic frequency distributes.Although described by reference to the accompanying drawings embodiments of the present invention, those of ordinary skill in the art can be appendedIn the scope of claim, make various distortion or amendment.

Claims (8)

1. the interference matrix compensation method based on measurement report, is characterized in that, the method comprises:
Set up original interference matrix according to basic data and measurement report MR;
In all communities, search the identical BCCH community pair of You Yuan community and interfered cell composition;
For each identical BCCH community pair, determine interference compensation value according to position relation between described former community and described interfered cell;
Set up complete interference matrix according to described definite interference compensation value;
Wherein, while utilizing described complete interference matrix to carry out frequency division, can obtain and be better than utilizing described original interference matrix to carry out the effect of frequency division;
Wherein determine that according to position relation between described former community and described interfered cell interference compensation value further comprises:
If described former community is omnidirectional sector, using the mean value of 3 nearest Dui Yuan community, interfered cell area interference values as interference compensation value.
2. the interference matrix compensation method based on measurement report as claimed in claim 1, wherein determine that according to position relation between described former community and described interfered cell interference compensation value further comprises:
If described former community is directed sector, interfered cell is omnidirectional sector, and according to interfered cell, to other, the interference value with base station cell is determined described interference compensation value.
3. the interference matrix compensation method based on measurement report as claimed in claim 1, wherein determine that according to position relation between described former community and described interfered cell interference compensation value further comprises:
If the sector of the sector of described interfered cell and described former community is all directed sector, according to the position relation of all sectors of the base station, place, He Yuan community, sector of described interfered cell, determine level of interference, according to described level of interference, determine interference compensation value.
4. method in the interference matrix compensation method based on measurement report as described in as arbitrary in claim 1-3, wherein saidly set up original interference matrix according to basic data and measurement report MR and comprise: divide with frequently and adjacent frequency interference threshold according to C/I, and calculate and drop on the ratio of the interval measuring object sampled points of difference in total sampled point, obtain thus frequently same and/or adjacent frequency interference numerical value.
5. the device interference matrix of measurement report being compensated, comprising:
Set up module, set up original interference matrix according to basic data and measurement report MR;
Search module, in all communities, search the identical BCCH community pair of You Yuan community and interfered cell composition;
Determination module, for each identical BCCH community pair, determines interference compensation value according to position relation level between described former community and described interfered cell;
Compensating module, sets up complete interference matrix according to described definite interference compensation value;
Wherein, while utilizing described complete interference matrix to carry out frequency division, can obtain and be better than utilizing described original interference matrix to carry out the effect of frequency division;
Wherein said determination module is further carried out following operation:
If described former community is omnidirectional sector, using the mean value of 3 nearest Dui Yuan community, interfered cell area interference values as interference compensation value.
6. device as claimed in claim 5, wherein said determination module is further carried out following operation:
If described former community is directed sector, interfered cell is omnidirectional sector, and according to interfered cell, to other, the interference value with base station cell is determined described interference compensation value.
7. device as claimed in claim 5, wherein said determination module is further carried out following operation:
If the sector of the sector of described interfered cell and described former community is all directed sector, according to the position relation of all sectors of the base station, place, He Yuan community, sector of described interfered cell, determine corresponding level of interference, according to described level of interference, determine interference compensation value.
8. device as described in one of claim 5-7, the wherein said module of setting up is divided with frequency and adjacent frequency interference threshold according to C/I, and calculate and drop on the ratio of the interval measuring object sampled points of difference in total sampled point, obtain thus frequently same and/or adjacent frequency interference numerical value.
CN201210005734.XA 2012-01-09 2012-01-09 A kind of interference matrix compensation method and device based on measurement report Active CN103200608B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210005734.XA CN103200608B (en) 2012-01-09 2012-01-09 A kind of interference matrix compensation method and device based on measurement report

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210005734.XA CN103200608B (en) 2012-01-09 2012-01-09 A kind of interference matrix compensation method and device based on measurement report

Publications (2)

Publication Number Publication Date
CN103200608A CN103200608A (en) 2013-07-10
CN103200608B true CN103200608B (en) 2016-05-25

Family

ID=48722930

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210005734.XA Active CN103200608B (en) 2012-01-09 2012-01-09 A kind of interference matrix compensation method and device based on measurement report

Country Status (1)

Country Link
CN (1) CN103200608B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111278040A (en) * 2018-12-05 2020-06-12 中国移动通信集团四川有限公司 Interference source positioning method, device, equipment and computer storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004057900A1 (en) * 2002-12-20 2004-07-08 Nokia Corporation Method and system for allocating channels in a cellular communication network
CN101547449A (en) * 2009-05-04 2009-09-30 中国移动通信集团浙江有限公司 Frequency sweep and mobile phone measurement report-based method for automatic frequency optimization
CN101873623A (en) * 2010-07-09 2010-10-27 广州杰赛科技股份有限公司 Automatic frequency optimization method based on measurement report data
CN101917724A (en) * 2010-07-30 2010-12-15 深圳市科虹通信有限公司 Method and system for obtaining combined interference matrixes of broadcast control channels
CN102164379A (en) * 2010-02-23 2011-08-24 中兴通讯股份有限公司 Method and device for adjusting frequency

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004057900A1 (en) * 2002-12-20 2004-07-08 Nokia Corporation Method and system for allocating channels in a cellular communication network
CN101547449A (en) * 2009-05-04 2009-09-30 中国移动通信集团浙江有限公司 Frequency sweep and mobile phone measurement report-based method for automatic frequency optimization
CN102164379A (en) * 2010-02-23 2011-08-24 中兴通讯股份有限公司 Method and device for adjusting frequency
CN101873623A (en) * 2010-07-09 2010-10-27 广州杰赛科技股份有限公司 Automatic frequency optimization method based on measurement report data
CN101917724A (en) * 2010-07-30 2010-12-15 深圳市科虹通信有限公司 Method and system for obtaining combined interference matrixes of broadcast control channels

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GSM网络中同邻频干扰分析优化的方法研究与系统实现;柯文渊;《信息科技辑》;20111231;第19-43页 *
黄剑锋,侯优优.TD-SCDMA网优专家系统中的干扰矩阵应用.《电信工程技术与标准化》.2011, *

Also Published As

Publication number Publication date
CN103200608A (en) 2013-07-10

Similar Documents

Publication Publication Date Title
CN102149192B (en) Cellular network wireless positioning method based on cooperation of mobile stations
CN101426206B (en) Method and system for gridding communication network call distribution
WO2016165435A1 (en) Method and device for detecting overshoot coverage of radio network, and communication system
CN104144430B (en) A kind of Physical Cell Identifier optimization method and device
CN106973406B (en) Weak coverage area detection method and system
CN104469676A (en) Method and system for locating mobile terminal
CN101888640A (en) Positioning method of city mobile station
CN101790189A (en) Method for positioning mobile station
WO2022022486A1 (en) Processing method and processing apparatus for saving energy of base station
CN104038901A (en) Indoor positioning method for reducing fingerprint data acquisition workload
CN114363808B (en) Indoor positioning method based on RSSI ranging
CN103517310A (en) Method for automatically correcting positioning wireless base station dataset
CN101964981A (en) Method and device for pre-evaluating telephone traffic and configuring carrier frequency of new base station
CN101304595B (en) Method and apparatus for acquiring frequency plan
CN108281023A (en) A kind of method and system showing real-time road by mobile terminal
CN105866732A (en) Improved MK model and WKNN algorithm combined mixed indoor positioning method
CN106922017B (en) Positioning method and terminal
CN101201398B (en) Method for positioning indirect wave using auxiliary measuring angle of adjacent station ranging
CN106488555B (en) Indoor positioning method and device
CN103200608B (en) A kind of interference matrix compensation method and device based on measurement report
CN101674599A (en) Detection method and detection system of boundary roaming cell
Arai et al. Color radiomap interpolation for efficient fingerprint wifi-based indoor location estimation
CN109769297A (en) A kind of indoor orientation method merging multifrequency WIFI signal
CN105208582A (en) Method of signal optimization based on grid interference in LTE network and system thereof
CN110519685A (en) Indoor orientation method, device and medium based on WiFi

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant