CN112348793B - Method for automatically identifying and calculating antenna pattern pits - Google Patents
Method for automatically identifying and calculating antenna pattern pits Download PDFInfo
- Publication number
- CN112348793B CN112348793B CN202011222311.4A CN202011222311A CN112348793B CN 112348793 B CN112348793 B CN 112348793B CN 202011222311 A CN202011222311 A CN 202011222311A CN 112348793 B CN112348793 B CN 112348793B
- Authority
- CN
- China
- Prior art keywords
- antenna pattern
- pit
- data
- calculating
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000009499 grossing Methods 0.000 claims description 10
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method for automatically identifying and calculating antenna pattern pits. The invention provides a special algorithm which can eliminate the interference of small saw teeth of a pattern to pit identification caused by a test environment error, and accurately identify and calculate the pits; meanwhile, the calculation process is simple and clear, the data processing efficiency is high, and the software implementation is convenient. When the invention is used for engineering application, the pit processing period in the antenna manufacturing process is shortened by 96%, and the pit identification efficiency and quality are obviously improved.
Description
Technical Field
The invention relates to the technical field of antenna manufacturing, in particular to a method for automatically identifying and calculating antenna pattern pits.
Background
The antenna pattern is an important basis for measuring the performance of the antenna, various parameters of the antenna can be calculated and extracted from the antenna pattern, and the antenna pattern data processing service is widely applied to manufacturing activities of various antennas at present. At present, pits of an antenna pattern are mainly identified by means of naked eyes, so that the efficiency is low and the pits are easily influenced by the subjective view. The existing method related to pit identification cannot eliminate interference of small saw teeth of an antenna pattern to pit identification caused by test environment errors, so that the pit identification method cannot meet the pit identification and calculation requirements of the antenna pattern.
Disclosure of Invention
Aiming at the defects in the prior art, the method for automatically identifying and calculating the antenna pattern pits solves the problems that an algorithm for effectively identifying the antenna pattern pits is lacking at present, the manual identification efficiency is low, and the influence of the main view is large.
In order to achieve the aim of the invention, the invention adopts the following technical scheme: a method for automatically identifying and calculating antenna pattern pits, comprising the steps of:
s1, importing test data (A) of a two-dimensional antenna pattern into an xy coordinate axis n ，B n )(n＝1,2,3…m)，A n For angle data in test data, B n The m is the number of the test data;
s2, performing smoothing on the antenna pattern to obtain new amplitude data C n (n＝1,2,3…m)；
S3, passing new amplitude data C n And angle data A n Calculating slope D of amplitude n (n＝1,2,3…m-1)；
S4, slope D of amplitude n (n=1, 2,3 … m-1) traversing from left to right looking for the inflection point E of the antenna pattern i (i=1, 2,3 … j), j being the number of inflection points;
s5, taking out inflection point E i (i=1, 2,3 … j), finding the largest amplitude data and the smallest amplitude data, calculating the pit value by the largest amplitude data and the smallest amplitude data, judging whether the pit value meets the pit index requirement, and outputting the pit value if the pit value meets the pit index requirement.
Further: the specific steps of the smoothing process in the step S2 are as follows: according to the shape characteristics of the antenna pattern, natural spline interpolation of cubic spline interpolation is adopted, the value of the smooth point number is the test point number covered in a single small sawtooth, and the empirical value is 3.
Further: slope D of amplitude in step S3 n The calculation formula of (2) is as follows:
D n ＝(C n+1 -C n )/(A n+1 -A n )。
further: the identification condition of the inflection point in the step S4 is as follows:
D n ×D n+1 ≤0。
further: the pit index requirements in the step S5 are as follows: the empirical value of the pit index is less than 3dB.
Further: the calculation formula of the pit value in the step S5 is as follows:
B＝B max -B min
in the above, B is pit value, B max For maximum amplitude data, B min Is the smallest amplitude data.
The beneficial effects of the invention are as follows: the invention provides a special algorithm which can eliminate the interference of small saw teeth of a pattern to pit identification caused by a test environment error, and accurately identify and calculate the pits; meanwhile, the calculation process is simple and clear, the data processing efficiency is high, and the software implementation is convenient.
When the invention is used for engineering application, the pit processing period in the antenna manufacturing process is shortened by 96%, and the pit identification efficiency and quality are obviously improved.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a schematic diagram of eliminating small pattern serrations by smoothing in an embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating the location of inflection points of a smoothed antenna pattern in accordance with an embodiment of the present invention;
fig. 4 is a schematic diagram of pit calculation in an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and all the inventions which make use of the inventive concept are protected by the spirit and scope of the present invention as defined and defined in the appended claims to those skilled in the art.
As shown in fig. 1, a method for automatically identifying and calculating antenna pattern pits includes the steps of:
s1, importing test data (A) of a two-dimensional antenna pattern into an xy coordinate axis n ，B n )(n＝1,2,3…m)，A n For angle data in test data, B n The m is the number of the test data;
s2, performing smoothing on the antenna pattern to obtain new amplitude data C n (n=1, 2,3 … m); the smoothing is to filter small saw teeth caused by test environment errors and avoid interference with pit recognition. According to the shape characteristics of the antenna pattern, natural spline interpolation of cubic spline interpolation is adopted, the smooth points are related to angle stepping setting of a test environment and the precision of a test instrument, the smooth points are the test points covered in a single small saw tooth, and the empirical value is 3.
S3, passing new amplitude data C n And angle data A n Calculating slope D of amplitude n (n＝1,2,3…m-1)；
D n ＝(C n+1 -C n )/(A n+1 -A n )。
S4, slope D of amplitude n (n=1, 2,3 … m-1) traversing from left to right looking for the inflection point E of the antenna pattern i (i=1, 2,3 … j), j being the number of inflection points; the identification conditions of the inflection points are as follows: d (D) n ×D n+1 ≤0。
S5, taking out inflection point E i (i=1, 2,3 … j), the largest amplitude data and the smallest amplitude are found outThe data, the pit value is calculated by the maximum amplitude data and the minimum amplitude data, the maximum amplitude value is the main lobe peak value of the antenna pattern, the minimum amplitude value is the worst pit of the antenna pattern, and the calculation formula of the pit value is as follows:
B＝B max -B min
in the above, B is pit value, B max For maximum amplitude data, B min Is the smallest amplitude data.
And judging whether the pit value meets the pit index requirement, wherein the empirical value of the pit index is less than 3dB. And outputting the pit value if the pit value is satisfied.
In one embodiment of the invention, an angle array and an amplitude array of the antenna pattern are imported as per step S1.
And (2) smoothing the antenna pattern according to the step (S2), wherein the width of the small saw teeth introduced by the test environment does not exceed 4 degrees, the angle step of the test environment is 2 degrees, and the number of test points covered by a single small saw tooth is 3, so that the number of the smooth points takes a value of 3.
The SpInterp function of NI LabWindows CVI 8.5.5 software is called to perform three times of natural spline interpolation on the antenna pattern (other software can be used for realizing three times of natural spline interpolation), the obtained comparison pattern is shown in fig. 2, the effect of fitting the antenna pattern before and after smoothing is good as can be known from fig. 2, two small saw teeth caused by the testing environment error at-52 DEG are filtered through smoothing, and interference of the small saw teeth on pit identification is avoided.
Calculating the slope curve of the smoothed antenna pattern according to the step S3 to obtain a slope array D 1 、D 2 、D 3 …D m-1 。
Searching for inflection points of the antenna pattern according to the step S4, wherein the inflection point positioning condition is (D n ×D n+1 And 0), 5 inflection points can be identified, as shown in FIG. 3.
As shown in fig. 4, the maximum amplitude-inflection point 3 in the inflection points is positioned as the main lobe peak of the directional diagram according to the step S5, the minimum amplitude-inflection point 2 in the inflection points is positioned as the most serious pit of the directional diagram, the corresponding amplitude values-65.65 and-66.84 of the inflection points 3 and 2 are found, delta is 1.19dB, and 1.19dB meets the index condition of <3dB, so that the pit of the directional diagram of the antenna of the present frequency point is judged to be qualified.
Claims (6)
1. A method for automatically identifying and calculating antenna pattern pits, comprising the steps of:
s1, importing test data (A) of a two-dimensional antenna pattern into an xy coordinate axis n ，B n )(n＝1,2,3…m)，A n For angle data in test data, B n The m is the number of the test data;
s2, performing smoothing on the antenna pattern to obtain new amplitude data C n (n＝1,2,3…m)；
S3, passing new amplitude data C n And angle data A n Calculating slope D of amplitude n (n＝1,2,3…m-1)；
S4, slope D of amplitude n (n=1, 2,3 … m-1) traversing from left to right looking for the inflection point E of the antenna pattern i (i=1, 2,3 … j), j being the number of inflection points;
s5, taking out inflection point E i (i=1, 2,3 … j), finding the largest amplitude data and the smallest amplitude data, calculating the pit value by the largest amplitude data and the smallest amplitude data, judging whether the pit value meets the pit index requirement, and outputting the pit value if the pit value meets the pit index requirement.
2. The method for automatically identifying and calculating antenna pattern pits according to claim 1, wherein the specific steps of the smoothing process in step S2 are: according to the shape characteristics of the antenna pattern, natural spline interpolation of cubic spline interpolation is adopted, the value of the smooth point number is the test point number covered in a single small sawtooth, and the empirical value is 3.
3. The method for automatically identifying and calculating antenna pattern pits according to claim 1, wherein the slope D of the amplitude in step S3 n The calculation formula of (2) is as follows:
D n ＝(C n+1 -C n )/(A n+1 -A n )。
4. the method for automatically identifying and calculating antenna pattern pits according to claim 1, wherein the identification condition of the inflection point in the step S4 is:
D n ×D n+1 ≤0。
5. the method for automatically identifying and calculating antenna pattern pits according to claim 1, wherein the pit index requirement in step S5 is: the empirical value of the pit index is less than 3dB.
6. The method for automatically identifying and calculating antenna pattern pits according to claim 1, wherein the calculation formula of the pit value in step S5 is:
B＝B max -B min
in the above, B is pit value, B max For maximum amplitude data, B min Is the smallest amplitude data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011222311.4A CN112348793B (en) | 2020-11-05 | 2020-11-05 | Method for automatically identifying and calculating antenna pattern pits |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011222311.4A CN112348793B (en) | 2020-11-05 | 2020-11-05 | Method for automatically identifying and calculating antenna pattern pits |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112348793A CN112348793A (en) | 2021-02-09 |
| CN112348793B true CN112348793B (en) | 2023-04-25 |
Family
ID=74428298
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011222311.4A Active CN112348793B (en) | 2020-11-05 | 2020-11-05 | Method for automatically identifying and calculating antenna pattern pits |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112348793B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113300085A (en) * | 2021-05-25 | 2021-08-24 | 安徽大学 | 3G/4G/5G-oriented small broadband dual-polarized base station antenna unit |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103137093A (en) * | 2006-09-15 | 2013-06-05 | 株式会社半导体能源研究所 | Display device and method of driving the same |
| CA2821127A1 (en) * | 2012-07-16 | 2014-01-16 | Halliburton Energy Services, Inc. | Real-time casing detection using tilted and crossed antenna measurement |
| FR3031813A1 (en) * | 2015-01-16 | 2016-07-22 | Centre Nat D'etudes Spatiales (Cnes) | SYSTEM AND METHOD FOR MEASURING COMPONENTS OF THE RADIATION DIAGRAM OF AN ANTENNA |
| CN107657098A (en) * | 2017-09-15 | 2018-02-02 | 哈尔滨工程大学 | Perimeter antenna array Sparse methods based on quantum chicken group's mechanism of Evolution |
| CN108389181A (en) * | 2018-01-23 | 2018-08-10 | 南京大学 | It is a kind of based on the part defect of deep learning and normal direction figure detection and localization method |
| CN108957145A (en) * | 2018-06-04 | 2018-12-07 | 中国科学院大气物理研究所 | Local dodges the Return stroke of ground flash pulse checkout system and method for positioning net entirely |
| CN109001694A (en) * | 2018-05-31 | 2018-12-14 | 中国电子科技集团公司第二十九研究所 | A kind of dynamic self-adapting antenna scanning simulated behavior method and system |
| CN109670228A (en) * | 2018-12-12 | 2019-04-23 | 中国船舶重工集团公司第七二四研究所 | A kind of array antenna unit directional diagram Fast Interpolation calculation method |
| CN111798474A (en) * | 2019-04-05 | 2020-10-20 | 三星电子株式会社 | Image processing apparatus and image processing method thereof |
-
2020
- 2020-11-05 CN CN202011222311.4A patent/CN112348793B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103137093A (en) * | 2006-09-15 | 2013-06-05 | 株式会社半导体能源研究所 | Display device and method of driving the same |
| CA2821127A1 (en) * | 2012-07-16 | 2014-01-16 | Halliburton Energy Services, Inc. | Real-time casing detection using tilted and crossed antenna measurement |
| FR3031813A1 (en) * | 2015-01-16 | 2016-07-22 | Centre Nat D'etudes Spatiales (Cnes) | SYSTEM AND METHOD FOR MEASURING COMPONENTS OF THE RADIATION DIAGRAM OF AN ANTENNA |
| CN107657098A (en) * | 2017-09-15 | 2018-02-02 | 哈尔滨工程大学 | Perimeter antenna array Sparse methods based on quantum chicken group's mechanism of Evolution |
| CN108389181A (en) * | 2018-01-23 | 2018-08-10 | 南京大学 | It is a kind of based on the part defect of deep learning and normal direction figure detection and localization method |
| CN109001694A (en) * | 2018-05-31 | 2018-12-14 | 中国电子科技集团公司第二十九研究所 | A kind of dynamic self-adapting antenna scanning simulated behavior method and system |
| CN108957145A (en) * | 2018-06-04 | 2018-12-07 | 中国科学院大气物理研究所 | Local dodges the Return stroke of ground flash pulse checkout system and method for positioning net entirely |
| CN109670228A (en) * | 2018-12-12 | 2019-04-23 | 中国船舶重工集团公司第七二四研究所 | A kind of array antenna unit directional diagram Fast Interpolation calculation method |
| CN111798474A (en) * | 2019-04-05 | 2020-10-20 | 三星电子株式会社 | Image processing apparatus and image processing method thereof |
Non-Patent Citations (5)
| Title |
|---|
| Xu Liguo 等.Pulse doppler radar Seeker multi-target detection simulation research.《2017 First International Conference on Electronics Instrumentation & Information Systems》.2018,1-6. * |
| 任桐鑫 等.基于遗传算法和贪婪算法的新一代通讯网络设计研究.《数学建模及其应用》.2020,第9卷(第2期),41-52. * |
| 何大鹏 等.基于优化理论的网络部署与分配模型.《信息通信》.2020,(第5期),31-33. * |
| 张建民.基于卫星导航测量的零相差天线实现.《中国 优秀硕士学位论文全文数据库信息科技辑》.2018,I136-961. * |
| 陈凡.宽角度探测器天线的研究及其应用.《中国优秀硕士学位论文全文数据库信息科技辑》.2012,I136-28. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112348793A (en) | 2021-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112348793B (en) | Method for automatically identifying and calculating antenna pattern pits | |
| EP3896873A1 (en) | Field intensity prediction method and apparatus, and device and storage medium | |
| CN106096223B (en) | Five-hole probe data processing method | |
| CN109035363B (en) | Line-circle optimal fitting method for rapid iteration | |
| CN110797085B (en) | Method, system, equipment and storage medium for inquiring gene data | |
| CN107357959A (en) | A kind of high-precision blade profile characteristic parameter extraction method | |
| CN111307751A (en) | Spectrogram baseline correction method, system and detection method in tea near infrared spectrum analysis | |
| CN104155915A (en) | Optical element random processing path planning method | |
| CN110852652B (en) | Hydroelectric characteristic curve correction method and device | |
| CN109284062B (en) | Touch data processing method, device, terminal and medium | |
| WO2022127015A1 (en) | Rapid-converging ultrahigh-frequency microwave rainfall data discretization method | |
| CN113378433A (en) | Method, device, equipment and medium for determining forming technological parameters of composite aluminum plate | |
| CN111487440A (en) | Calibration method of five-hole probe | |
| CN117123922A (en) | Oval tube centering method, system, electronic equipment and storage medium | |
| CN109359360B (en) | Structural stress processing method based on local features | |
| CN108267743B (en) | Fast iterative underwater positioning method based on fitting interpolation | |
| KR102219745B1 (en) | Method and apparatus for processing biological sequence data | |
| CN110850812B (en) | Ion beam etching rate control method and device based on model | |
| CN108121955A (en) | A kind of smooth l based on adjusted mean approximation hyperbolic tangent function0Norm Method | |
| CN110095743B (en) | Distribution network terminal wave recording performance test module and waveform fitting method thereof | |
| CN107292220B (en) | Eyebrow position initialization method and electronic equipment | |
| CN112562012A (en) | Calibration plate, vehicle camera calibration method and related device | |
| CN118111358A (en) | Method for accurately detecting arc curvature radius of large-radius small-arc molded surface | |
| Ravi et al. | Compressor map prediction tool | |
| CN113673056B (en) | Method for determining cold test parameter limit value of engine |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |