Movatterモバイル変換


[0]ホーム

URL:


CN109862586A - A simulation system and method for realizing wireless network communication field strength coverage - Google Patents

A simulation system and method for realizing wireless network communication field strength coverage
Download PDF

Info

Publication number
CN109862586A
CN109862586ACN201910157744.7ACN201910157744ACN109862586ACN 109862586 ACN109862586 ACN 109862586ACN 201910157744 ACN201910157744 ACN 201910157744ACN 109862586 ACN109862586 ACN 109862586A
Authority
CN
China
Prior art keywords
field strength
module
signal
wireless network
network communication
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.)
Granted
Application number
CN201910157744.7A
Other languages
Chinese (zh)
Other versions
CN109862586B (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.)
Hunan City University
Original Assignee
Hunan City University
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 Hunan City UniversityfiledCriticalHunan City University
Priority to CN201910157744.7ApriorityCriticalpatent/CN109862586B/en
Publication of CN109862586ApublicationCriticalpatent/CN109862586A/en
Application grantedgrantedCritical
Publication of CN109862586BpublicationCriticalpatent/CN109862586B/en
Activelegal-statusCriticalCurrent
Anticipated expirationlegal-statusCritical

Links

Classifications

Landscapes

Abstract

Translated fromChinese

本发明属于无线网络通信技术领域,公开了一种实现无线网络通信场强覆盖仿真系统及方法,所述实现无线网络通信场强覆盖仿真系统包括:信号检测模块、通信场强测量模块、主控模块、场景模型构建模块、场强覆盖计算模块、数据优化模块、显示模块。本发明通过通信场强测量模块在无线电信号环境日益复杂的情况下,本发明具有更强的信号适应性;采用了信号带宽测量技术,可针对不同带宽的信号进行自动化场强测量;同时,通过场强覆盖计算模块将多个基站的经纬度位置点映射到平面区域上,构成平面基站点集,以平面几何的算法对所述基站点集进行处理,无须考虑基站的地形地貌因素,简单、快速地确定基站或小区的无线场强覆盖区域范围。

The invention belongs to the technical field of wireless network communication, and discloses a simulation system and method for realizing wireless network communication field strength coverage. The wireless network communication field strength coverage simulation system includes: a signal detection module, a communication field strength measurement module, a main Module, scene model building module, field strength coverage calculation module, data optimization module, display module. The present invention has stronger signal adaptability through the communication field strength measurement module under the condition that the radio signal environment is increasingly complex; the signal bandwidth measurement technology is adopted, and automatic field strength measurement can be performed for signals of different bandwidths; The field strength coverage calculation module maps the longitude and latitude position points of multiple base stations to a plane area to form a plane base station point set. The base station point set is processed by a plane geometry algorithm without considering the topographic and topographic factors of the base station, which is simple and fast. The wireless field strength coverage area range of the base station or cell is determined.

Description

A kind of realization wireless network communication field intensity emulation of coverage capability system and method
Technical field
The invention belongs to wireless network communication technique field more particularly to a kind of covering of realization wireless network communication field intensity are imitativeTrue system and method.
Background technique
Wireless network (wireless network) is the network realized using wireless communication technique.Wireless network both includedAllow user to establish global voice and data network that remote-wireless connects, also includes being optimized near radio connectionInfrared technology and radio-frequency technique, very similar with the purposes of cable network, maximum difference is the difference of transmission medium,Replace cable using radio technology, can be mutually backups with cable network.Wireless network, which is divided into, passes through public mobile networkThe wireless network (such as 4G, 3G or GPRS) and WLAN (WiFi) two ways of realization.GPRS surfing Internet with cell phone mode, is oneKind by mobile telephone network access Internet wireless networking mode, as long as therefore you where city opened GPRS onlineBusiness, you can be surfed the Internet in any one corner by laptop.However, prior wireless network communication field intensity measuresProcess is comparatively laborious, and the degree of automation is not high;Meanwhile wireless field density covering area range is determined, it is both needed to consider topography and geomorphologyDifference, calculating parameter is more, and method is complicated, low efficiency.
In conclusion problem of the existing technology is: prior wireless network communication field intensity measurement process is comparatively laborious, fromDynamicization degree it is not high and acquisition wireless communication signals noise it is big, inaccuracy;The field strength data degradation degree of acquisition is big, and availability is notIt is high.Meanwhile wireless field density covering area range is determined, it is both needed to consider the difference of topography and geomorphology, calculating parameter is more, and method is multipleIt is miscellaneous, low efficiency.
Summary of the invention
In view of the problems of the existing technology, the present invention provides a kind of realization wireless network communication field intensity emulation of coverage capability systemsSystem and method.
The invention is realized in this way a kind of realization wireless network communication field intensity emulation of coverage capability method, described to realize wirelesslyNetwork communication Field strength coverage emulation mode includes:
The first step detects wireless network using based on Optimum combining and the signal detector of the signal detection algorithm of characteristic valueSignal of communication;
Second step measures wireless network communication field intensity intensity data using field strength measurement device;
Third step utilizes the model of place of modeling software building wireless network communication field intensity distribution;
4th step calculates the wireless network communication field intensity area of coverage according to the field strength data of detection using calculation procedureDomain;
5th step optimizes processing to field strength data using the optimization program based on second order approximatioss;Second order Approximation NumberIt include: the Taylor expansion of the function of many variables according to optimization algorithm are as follows:
Second order expension is carried out to loss function:
The necessary condition that loss function obtains minimum value is:
Finally obtain the iterative formula of β:
6th step shows signal of communication, field strength, model of place, overlay area data using display.
Further, the signal detection algorithm based on Optimum combining and characteristic value specifically includes in the first step:
(1) collected wireless communication signal is done into discrete Fourier transform, and in frequency domain using maximization outputThe criterion of signal power determines a best initial weights, then merges multiple signals with best initial weights;
(2) real part of data after merging is taken, autocorrelation matrix R is constructedy(Ns);
(3) all characteristic values for seeking matrix Ry (Ns), then acquire detection statistic T according to corporal;
Tr (R in formulay(Ns)) be composite signal autocorrelation matrix mark, i.e. the sum of all characteristic values of matrix, by matrixThe ratio between mark and minimal eigenvalue are used as detection statistic;
(4) following formula is combined, the theoretical detection threshold γ under specific false-alarm probability is calculated;
F in formula1For the cumulative distribution function that Tracy-Widom the 1st is distributed, there is no expression formula form, the function is correspondingValue is acquired by tabling look-up;
(5) adjudicate: if T > γ, there are signal, affirmative judgements in channel;Otherwise signal, negative decision is not present.
The realization wireless network communication field intensity emulation of coverage capability method is realized another object of the present invention is to provide a kind ofRealization wireless network communication field intensity emulation of coverage capability system, the realization wireless network communication field intensity emulation of coverage capability system includes:
Signal detection module is connect with main control module, for detecting wireless communication signal by signal detector;
Communication field intensity measurement module, connect with main control module, for measuring wireless communication field by field strength measurement deviceStrong intensity data;
Main control module constructs module, Field strength coverage meter with signal detection module, communication field intensity measurement module, model of placeModule, data-optimized module, display module connection are calculated, is worked normally for controlling modules by single-chip microcontroller;
Model of place constructs module, connect with main control module, for constructing wireless network communication field intensity by modeling softwareThe model of place of distribution;
Field strength coverage computing module, connect with main control module, for the field strength number by calculation procedure according to detectionAccording to calculating wireless network communication field intensity overlay area;
Data-optimized module, connect with main control module, for optimizing processing to field strength data by optimizing program;
Display module is connect with main control module, for by display display signal of communication, field strength, model of place,Overlay area data.
Another object of the present invention is to provide a kind of application realization wireless network communication field intensity emulation of coverage capability methodsWireless communication platform.
Advantages of the present invention and good effect are as follows: the present invention is having existing field strength measurement by communication field intensity measurement moduleWhile technology carries out the ability of field strength measurement for radio frequency tone signal, it is also equipped with for widely used radio frequency tune in practiceSignal processed carries out the ability of field strength measurement, and in the case where radio signal environment is increasingly complicated, the present invention has stronger letterNumber adaptability;Using signal bandwidth measuring technique, this survey the present invention can with automatic measurement signal bandwidth for subsequent field strengthAmount provides the range of frequency bandwidths of power spectral integral, therefore the present invention can carry out automation field strength for the signal of different bandwidth and surveyAmount;Meanwhile the longitude and latitude position point of multiple base stations being mapped on plane domain by Field strength coverage computing module, constitute planeBase station point set is handled the base station point set with the algorithm of plane geometry, regardless of the topography and geomorphology factor of base station, letterWireless field density covering area range that is single, quickly determining base station or cell.
The present invention combines multi signal composition algorithm and characteristic value detection algorithm, constructs new detection statistic,Any prior information of undesired signal does not require multiple reception noises to meet the equal condition of power, while overcoming noisePower uncertainty influences, and belongs to total blindness's signal detection algorithm.The present invention carries out field strength data using second order approximate algorithm excellentChange, can effectively improve field strength data using degree.
Detailed description of the invention
Fig. 1 is realization wireless network communication field intensity emulation of coverage capability system structure diagram provided in an embodiment of the present invention;
In figure: 1, signal detection module;2, communication field intensity measurement module;3, main control module;4, model of place constructs module;5, Field strength coverage computing module;6, data-optimized module;7, display module.
Fig. 2 is realization wireless network communication field intensity emulation of coverage capability method flow diagram provided in an embodiment of the present invention.
Specific embodiment
In order to further understand the content, features and effects of the present invention, the following examples are hereby given, and cooperate attached drawingDetailed description are as follows.
Structure of the invention is explained in detail with reference to the accompanying drawing.
As shown in Figure 1, realization wireless network communication field intensity emulation of coverage capability system provided in an embodiment of the present invention includes: signalDetection module 1, communication field intensity measurement module 2, main control module 3, model of place building module 4, Field strength coverage computing module 5, numberAccording to optimization module 6, display module 7.
Signal detection module 1 is connect with main control module 3, for detecting wireless communication signal by signal detector;
Communication field intensity measurement module 2 is connect with main control module 3, for measuring wireless communication by field strength measurement deviceField strength data;
Main control module 3 is covered with signal detection module 1, communication field intensity measurement module 2, model of place building module 4, field strengthLid computing module 5, data-optimized module 6, display module 7 connect, and work normally for controlling modules by single-chip microcontroller;
Model of place constructs module 4, connect with main control module 3, for constructing wireless communication field by modeling softwareThe model of place being distributed by force;
Field strength coverage computing module 5 is connect with main control module 3, for the field strength by calculation procedure according to detectionData calculate wireless network communication field intensity overlay area;
Data-optimized module 6, connect with main control module 3, for optimizing processing to field strength data by optimizing program;
Display module 7 is connect with main control module 3, for showing signal of communication, field strength, scene mould by displayType, overlay area data.
2 measurement method of communication field intensity measurement module provided in an embodiment of the present invention is as follows:
(1) by rf modulated signal by being converted to 70MHz analog if signal after mixing and filtering;
(2) re-sampling obtains digital medium-frequency signal after 70MHz analog if signal being carried out variable gain amplification;
(3) zero-frequency digital baseband signal is obtained after digital medium-frequency signal being carried out Digital Down Convert;
(4) Fast Fourier Transform (FFT) is carried out to zero-frequency digital baseband signal and obtains its power spectrum, carried out further according to power spectrumBandwidth measurement obtains the bandwidth of zero-frequency digital baseband signal;
(5) power spectral integral in bandwidth is carried out to modulated signal, obtains rf modulated signal field strength.
5 calculation method of Field strength coverage computing module provided in an embodiment of the present invention is as follows:
(1) the full base station point set being mapped to the longitude and latitude of all base stations corresponding location point in plane;
(2) full base station point set is processed into subdivision triangle gridding;
(3) position in the center of circle of the circumscribed circle of the triangle in subdivision triangle gridding comprising same base point is determined, and willThe center of circle is sequentially connected the Thiessen polygon region that the region to be formed is determined as the base station;
(4) overlapping region of overlay area is initially justified in the Thiessen polygon region and the base station for determining base station, and will be overlappedThe range in region is determined as the wireless field density covering area range of base station;
(5) overlay area is cut using the azimuth of cell each in base station, determines the wireless of each cell in base stationField strength coverage region range.
In step (2), the step of full base station point set is processed into subdivision triangle gridding provided in an embodiment of the present inventionInclude:
1) three or three or more initial base station points are chosen and form current base station point set, determine the initial of initial base station pointConvex closure;It determines the triangle for meeting empty circle characteristic in initial convex closure, and triangle is added in triangle chained list;First primordium will be removedOther base sites other than website are sequentially inserted into current base station point set, per base sites are increased after being newly inserted into one, that is, are formed currentThe current convex closure of base station point set;
2) when after increase base sites inside the current convex closure, and at the inside for the triangle that convex closure includes, threeThe triangle is deleted in angular chained list, successively cuts open three in the triangle to purpose side when determining as purposeDivide processing;
3) when after increase base sites inside current convex closure, and when in the triangle edges that convex closure includes, in tripletThe triangle is deleted in table, successively by the triangle where rear increasing base sites in addition to while determine as purpose side, it is rightPurpose side carries out subdivision processing;
4) when rear increasing base sites are on current convex closure side, the triangle for increasing base sites after including is deleted in triangle chained listShape, by the triangle where rear increasing base sites in addition to while determines as purpose side, to the progress subdivision processing of purpose side;
5) when rear increasing base sites are at the outside of current convex closure, using any point on current convex closure as starting point, along convexPacket is found counterclockwise increases side of the base sites on the right side of it after making;When the quantity on side is one, this is determined while for purpose, by meshThe two-end-point on side be connected with the rear base sites that increase, form new current convex closure, and to the progress subdivision processing of purpose side;When number of edges amountWhen for two or more, the endpoint for the side chain that side is formed is chosen, endpoint is connected with rear increasing base sites, is formed new currentConvex closure, and successively determine while for purpose, subdivision processing successively is carried out to purpose side.
Purpose side subdivision processing step provided in an embodiment of the present invention specifically includes:
When purpose is in the current convex closure, by the purpose side and the rear triangle for increasing base sites and being formedIt is added in triangle chained list;
When the purpose side is inside the current convex closure, obtaining from triangle chained list includes the three of the purpose sideIt is angular, when the triangle is not in the purpose side and the rear circumscribed circle for increasing the triangle that base sites are formed, by the purposeSide and the rear triangle for increasing base sites formation are added in triangle chained list;Otherwise, when the triangle is in the purpose Bian YuhouWhen increasing in the circumscribed circle for the triangle that base sites are formed, cancels former purpose side, determine two except the purpose side in the triangleWhile for purpose, subdivision processing successively is continued to the purpose side, and the triangle is deleted from triangle chained list.
The determination step in the Thiessen polygon region of base sites provided by the invention includes:
Inquiry includes the triangle of same base point, association triangle is determined it as, when the base sites are described convexWhen the inside of packet, the center of circle of the circumscribed circle of the association triangle is sequentially connected the closed area to be formed and is determined as the base stationThe Thiessen polygon region of point;
When the base sites are on the side of the convex closure, by the center of circle successively phase of the circumscribed circle of the association triangleEven, side chain is formed;It determines the head end point of the side chain and the angle of tail point and base sites formation, determines that the angle of the angle is flatThe position of the intersection point of overlay area is initially justified in the extended line of separated time and base station, by the intersection point and the head end point and tail point phaseThe region even formed is determined as the Thiessen polygon region of the base sites.
As shown in Fig. 2, realization wireless network communication field intensity covering method provided in an embodiment of the present invention includes:
S101: logical using the signal detector detection wireless network of the signal detection algorithm based on Optimum combining and characteristic valueBelieve signal;
S102: wireless network communication field intensity intensity data is measured using field strength measurement device;
S103: the model of place of modeling software building wireless network communication field intensity distribution is utilized;
S104: wireless network communication field intensity overlay area is calculated according to the field strength data of detection using calculation procedure;
S105: processing is optimized to field strength data using the optimization program based on second order approximatioss;
S106: signal of communication, field strength, model of place, overlay area data are shown using display.
In step S101, the signal detection algorithm provided in an embodiment of the present invention based on Optimum combining and characteristic value is specifically wrappedIt includes:
(1) collected wireless communication signal is done into discrete Fourier transform, and in frequency domain using maximization outputThe criterion of signal power determines a best initial weights, then merges multiple signals with best initial weights;
(2) real part of data after merging is taken, autocorrelation matrix R is constructedy(Ns);
(3) all characteristic values for seeking matrix Ry (Ns), then acquire detection statistic T according to corporal;
Tr (R in formulay(Ns)) be composite signal autocorrelation matrix mark, i.e. the sum of all characteristic values of matrix, by matrixThe ratio between mark and minimal eigenvalue are used as detection statistic;
(4) following formula is combined, the theoretical detection threshold γ under specific false-alarm probability is calculated;
F in formula1For the cumulative distribution function that Tracy-Widom the 1st is distributed, there is no expression formula form, the function is correspondingValue is acquired by tabling look-up, and the 1st distribution values table of Tracy-Widom is as follows:
The practical false-alarm probability of 1 different theories value of table
t-3.90-3.18-2.78-1.91-1.27-0.590.450.98
F1(t)0.010.050.100.300.500.700.900.95
(5) adjudicate: if T > γ, there are signal, affirmative judgements in channel;Otherwise signal, negative decision is not present.
In step S105, second order provided in an embodiment of the present invention approaches data optimization algorithm and includes:
The numerical method that general loss function optimizes derives from Taylor expansion, the Taylor expansion of the function of many variablesAre as follows:
Second order expension is carried out to loss function:
The necessary condition that loss function obtains minimum value is:
Finally obtain the iterative formula of β:
When the invention works, believe firstly, detecting wireless communication using signal detector by signal detection module 1Number;Wireless network communication field intensity intensity data is measured using field strength measurement device by communication field intensity measurement module 2;Secondly, master controlModule 3 constructs module 4 by model of place and constructs the model of place that wireless network communication field intensity is distributed using modeling software;Pass throughField strength coverage computing module 5 calculates the wireless network communication field intensity area of coverage according to the field strength data of detection using calculation procedureDomain;Then, processing is optimized to field strength data using optimization program by data-optimized module 6;Finally, passing through display module7 show signal of communication, field strength, model of place, overlay area data using display.
The above is only the preferred embodiments of the present invention, and is not intended to limit the present invention in any form,Any simple modification made to the above embodiment according to the technical essence of the invention, equivalent variations and modification, belong toIn the range of technical solution of the present invention.

Claims (4)

CN201910157744.7A2019-03-022019-03-02Simulation system and method for realizing wireless network communication field intensity coverageActiveCN109862586B (en)

Priority Applications (1)

Application NumberPriority DateFiling DateTitle
CN201910157744.7ACN109862586B (en)2019-03-022019-03-02Simulation system and method for realizing wireless network communication field intensity coverage

Applications Claiming Priority (1)

Application NumberPriority DateFiling DateTitle
CN201910157744.7ACN109862586B (en)2019-03-022019-03-02Simulation system and method for realizing wireless network communication field intensity coverage

Publications (2)

Publication NumberPublication Date
CN109862586Atrue CN109862586A (en)2019-06-07
CN109862586B CN109862586B (en)2023-03-24

Family

ID=66899568

Family Applications (1)

Application NumberTitlePriority DateFiling Date
CN201910157744.7AActiveCN109862586B (en)2019-03-022019-03-02Simulation system and method for realizing wireless network communication field intensity coverage

Country Status (1)

CountryLink
CN (1)CN109862586B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN114585010A (en)*2022-03-042022-06-03广西润象信息网络工程有限公司Method for realizing wireless network communication field intensity coverage simulation

Citations (6)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
US5668562A (en)*1996-04-191997-09-16Lgc Wireless, Inc.Measurement-based method of optimizing the placement of antennas in a RF distribution system
CN101247613A (en)*2008-03-112008-08-20中兴通讯股份有限公司Wireless signal propagation model test approach and system used for honeycomb network
EP2227050A2 (en)*2009-03-062010-09-08Hitachi Ltd.Adjoining base station determining method and management device
CN103634810A (en)*2013-12-242014-03-12山东润谱通信工程有限公司Method of positioning problem regions covered with indoor wireless network
CN103841567A (en)*2012-11-262014-06-04中国移动通信集团广东有限公司Method for obtaining polygon area of base station and base station
CN106384023A (en)*2016-12-022017-02-08天津大学 Hybrid Field Strength Prediction Method Based on Main Path

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
US5668562A (en)*1996-04-191997-09-16Lgc Wireless, Inc.Measurement-based method of optimizing the placement of antennas in a RF distribution system
CN101247613A (en)*2008-03-112008-08-20中兴通讯股份有限公司Wireless signal propagation model test approach and system used for honeycomb network
EP2227050A2 (en)*2009-03-062010-09-08Hitachi Ltd.Adjoining base station determining method and management device
CN103841567A (en)*2012-11-262014-06-04中国移动通信集团广东有限公司Method for obtaining polygon area of base station and base station
CN103634810A (en)*2013-12-242014-03-12山东润谱通信工程有限公司Method of positioning problem regions covered with indoor wireless network
CN106384023A (en)*2016-12-022017-02-08天津大学 Hybrid Field Strength Prediction Method Based on Main Path

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘智勇: ""求解稳态系统优化与控制问题的递推二阶逼近法"", 《五邑大学学报》*

Cited By (1)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN114585010A (en)*2022-03-042022-06-03广西润象信息网络工程有限公司Method for realizing wireless network communication field intensity coverage simulation

Also Published As

Publication numberPublication date
CN109862586B (en)2023-03-24

Similar Documents

PublicationPublication DateTitle
Talvitie et al.Distance-based interpolation and extrapolation methods for RSS-based localization with indoor wireless signals
Abhayawardhana et al.Comparison of empirical propagation path loss models for fixed wireless access systems
Al-Ahmadi et al.On the approximation of the generalized-Κ distribution by a gamma distribution for modeling composite fading channels
US20200015047A1 (en)Wifi multi-band fingerprint-based indoor positioning
Ajose et al.Propagation measurements and modelling at 1800 MHz in Lagos Nigeria
CN103686999B (en)Indoor wireless positioning method based on WiFi signal
CN102480678B (en)Fingerprint positioning method and system
US20080085692A1 (en)Radio frequency coverage map generation in wireless networks
CN108989986B (en) Wi-Fi indoor positioning method based on iterative segmentation space method
CN102984745A (en)Combined estimation method for Wi-Fi AP (wireless fidelity access point) position and path loss model
CN108064075B (en)Method and apparatus for adjusting reference signal transmit power of one or more cells
CN109640262B (en)Positioning method, system, equipment and storage medium based on mixed fingerprints
Ozuomba et al.Characterisation of propagation loss for a 3G cellular network in a crowded market area using CCIR model
Racko et al.Wi-fi fingerprint radio map creation by using interpolation
CN111919221A (en)Method of processing image, computer-readable storage medium recording the method, and apparatus for processing image
CN105334493B (en)A kind of indoor orientation method based on WLAN
CN106211327A (en)A kind of method automatically generating location fingerprint data
CN108234686A (en)A kind of method and apparatus of indoor and outdoor judgement
CN109862586A (en) A simulation system and method for realizing wireless network communication field strength coverage
CN109462864A (en)A kind of 5G communication typical scene channel model adaptive matching method
CN112738838B (en)Method and device for determining propagation model
US20070093213A1 (en)Method and system for electromagnetic field evaluation
CN106707232A (en)WLAN (Wireless local area network) propagation model positioning method based on crowd sensing
Bhupuak et al.Minimizing path loss prediction error using k-means clustering and fuzzy logic
Gao et al.A New Approach for Wi‐Fi‐Based People Localization in a Long Narrow Space

Legal Events

DateCodeTitleDescription
PB01Publication
PB01Publication
SE01Entry into force of request for substantive examination
SE01Entry into force of request for substantive examination
GR01Patent grant
GR01Patent grant

[8]ページ先頭

©2009-2025 Movatter.jp