CN110196281A - A kind of detection method, system and the terminal device of signal wave transmission velocity of wave - Google Patents
A kind of detection method, system and the terminal device of signal wave transmission velocity of wave Download PDFInfo
- Publication number
- CN110196281A CN110196281A CN201910554692.7A CN201910554692A CN110196281A CN 110196281 A CN110196281 A CN 110196281A CN 201910554692 A CN201910554692 A CN 201910554692A CN 110196281 A CN110196281 A CN 110196281A
- Authority
- CN
- China
- Prior art keywords
- signal wave
- wave
- grid
- detection zone
- signal
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
Landscapes
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The present invention provides detection method, system and the terminal devices of a kind of signal wave transmission velocity of wave, and method includes: that the corresponding detection zone of carrier to be detected is divided into multiple grids；Based on signal wave producer in the position of position and signal wave receiver in the detection zone in the detection zone, propagation stroke of the signal wave in each grid is calculated；Propagation stroke based on signal wave in each grid determines transmission velocity of wave of the signal wave in the carrier to be detected.The present invention is by being divided into multiple grids for the corresponding detection zone of carrier to be detected, then the principle of image procossing is utilized, calculate propagation stroke of the signal wave in each grid, and then calculate transmission velocity of wave of the signal wave in the carrier to be detected, by calculating the propagation stroke in each grid, the stroke of signal wave transmission can be more precisely computed out, the error of signal wave velocity of wave detection is reduced.
Description
Technical field
The invention belongs to signal detection technique field more particularly to detection method, the systems of a kind of signal wave transmission velocity of wave
And terminal device.
Background technique
It is influenced by environmental factor and the long-term cyclic reverse loading of train, ballastless track of high-speed railway structure support layer concrete exists
Phenomena such as surface peeling, aggregate exposure, cement mortar are lost, causes concrete binding power to decline, and concrete strength reduces, and causes
Track irregularity jeopardizes traffic safety.Since supporting course major part region is blocked by track plates, inner concrete intensity cannot be straight
It connects and is detected using concrete bearing intension testings instrument such as reisilometers.
Concrete elastic wave velocity of wave is the important indicator for reflecting concrete strength, and concrete Elastic Wave velocity of wave detects at present
Error is larger, causes the estimation inaccuracy to concrete strength, cannot find the damage of supporting course in time, influence high-speed railway
Maintenance.
Summary of the invention
In view of this, the embodiment of the invention provides detection method, system and the terminals of a kind of signal wave transmission velocity of wave to set
It is standby, to solve the problems, such as that current signal wave velocity of wave detection error is larger.
The first aspect of the embodiment of the present invention provides a kind of detection method of signal wave transmission velocity of wave, comprising:
The corresponding detection zone of carrier to be detected is divided into multiple grids, the carrier to be detected is signal wave transmission
Carrier；
Based on signal wave producer in the detection zone position and signal wave receiver in the detection zone
Position, calculate propagation stroke of the signal wave in each grid, wherein the signal wave producer is located at the detection zone
In the first side, the signal wave receiver is located at second side in the detection zone, the first side of the detection zone and
Second side of the detection zone is the opposite two sides of the detection zone；
Propagation stroke based on signal wave in each grid determines biography of the signal wave in the carrier to be detected
Defeated velocity of wave.
The second aspect of the embodiment of the present invention provides a kind of detection system of signal wave transmission velocity of wave, comprising:
Division module, for the corresponding detection zone of carrier to be detected to be divided into multiple grids, the carrier to be detected
For the carrier of signal wave transmission；
First computing module, for based on position of the signal wave producer in the detection zone and signal wave receiver
Position in the detection zone calculates propagation stroke of the signal wave in each grid, wherein the signal wave producer
The first side in the detection zone, the signal wave receiver are located at second side in the detection zone, the inspection
Second side of the first side and the detection zone of surveying region is the opposite two sides of the detection zone；
Second computing module determines the signal wave in institute for the propagation stroke based on signal wave in each grid
State the transmission velocity of wave in carrier to be detected.
The third aspect of the embodiment of the present invention provides a kind of terminal device, including memory, processor and is stored in
In the memory and the computer program that can run on the processor, when the processor executes the computer program
The step of realizing the detection method of signal wave transmission velocity of wave as described above.
The fourth aspect of the embodiment of the present invention provides a kind of computer readable storage medium, the computer-readable storage
Media storage has computer program, and the computer program realizes signal wave transmission velocity of wave as described above when being executed by processor
The step of detection method.
Then the present invention utilizes image procossing by the way that the corresponding detection zone of carrier to be detected is divided into multiple grids
Principle calculates propagation stroke of the signal wave in each grid, and then calculates transmission of the signal wave in the carrier to be detected
Velocity of wave can more precisely compute out the stroke of signal wave transmission, reduce signal wave by calculating the propagation stroke in each grid
The error of velocity of wave detection.
Detailed description of the invention
It to describe the technical solutions in the embodiments of the present invention more clearly, below will be to embodiment or description of the prior art
Needed in attached drawing be briefly described, it should be apparent that, the accompanying drawings in the following description is only of the invention some
Embodiment for those of ordinary skill in the art without any creative labor, can also be according to these
Attached drawing obtains other attached drawings.
Fig. 1 is the structure for the grid that the corresponding detection zone of carrier to be detected that one embodiment of the present of invention provides divides
Schematic diagram；
Fig. 2 is the implementation process signal of the detection method for the signal wave transmission velocity of wave that one embodiment of the present of invention provides
Figure；
Fig. 3 is the structural schematic diagram of the detection system for the signal wave transmission velocity of wave that one embodiment of the present of invention provides；
Fig. 4 is the schematic diagram for the terminal device that one embodiment of the present of invention provides.
Specific embodiment
In being described below, for illustration and not for limitation, the tool of such as particular system structure, technology etc is proposed
Body details, to understand thoroughly the embodiment of the present invention.However, it will be clear to one skilled in the art that there is no these specific
The present invention also may be implemented in the other embodiments of details.In other situations, it omits to well-known system, device, electricity
The detailed description of road and method, in case unnecessary details interferes description of the invention.
Description and claims of this specification and term " includes " and other any deformations in above-mentioned attached drawing are
Refer to " including but not limited to ", it is intended that cover and non-exclusive include.Such as the process, method comprising a series of steps or units
Or system, product or equipment are not limited to listed step or unit, but optionally further comprising the step of not listing
Or unit, or optionally further comprising other step or units intrinsic for these process, methods, product or equipment.In addition, art
Language " first ", " second " and " third " etc. is for distinguishing different objects, not for description particular order.
In order to illustrate technical solutions according to the invention, the following is a description of specific embodiments.
As shown in Figure 1, ballastless track of high-speed railway supporting course two sides are respectively arranged the generation of signal wave (elastic wave) signal wave
The corresponding detection zone of carrier to be detected will be divided into (m × n) grid, wherein signal wave is in net by device and signal wave receiver
The coordinate of lattice is that the velocity of wave of the signal wave of (a, b) is vab；Conveying length when passing through grid (a, b) for signal wave.Assuming that signal
Wave producer is P, and signal wave receiver is Q, and the transmission time of signal wave is tpq:
Order matrix T is signal wave propagation time matrix T={ T1,T2,…,Tm}:
Wherein: Ti=[t11,t21,…,tm1] ', symbol " ' " indicate transposition.
It enablesL={ LpqIt is stroke matrix, V={ vabIt is rate matrices.In conjunction with above equation, T, V, L it
Between relationship are as follows:
L/V=T (5)
Since stroke matrix L matrix is that three-dimensional matrice ((m × m) × n × m) to simplify the calculation converts stroke matrix L
For two-dimensional matrix, size are as follows: L ((m × m) × (n × m))；Rate matrices V (n × m) be converted into one-dimensional matrix V ((n × m) ×
1)；Time matrix T (m × m) is converted into one-dimensional matrix T ((m × m) × 1).
Two-dimentional stroke matrix L at this time is that the solution of singular matrix speed is not unique, so the problem is ill-posed problem.
Embodiment 1:
Fig. 2 shows signal wave provided by one embodiment of the invention transmission velocity of wave detection method implementation flow chart,
For ease of description, only parts related to embodiments of the present invention are shown, and details are as follows:
As shown in Fig. 2, a kind of detection method of signal wave transmission velocity of wave provided by the embodiment of the present invention, comprising:
The corresponding detection zone of carrier to be detected is divided into multiple grids by S101, and the carrier to be detected is signal wave
The carrier of transmission.
S102, based on position of the signal wave producer in the detection zone and signal wave receiver in the detection zone
Position in domain calculates propagation stroke of the signal wave in each grid, wherein the signal wave producer is located at the detection
The first side in region, the signal wave receiver are located at second side in the detection zone, and the first of the detection zone
Second side of side and the detection zone is the opposite two sides of the detection zone.
S103, the propagation stroke based on signal wave in each grid determine the signal wave in the carrier to be detected
In transmission velocity of wave.
In the present embodiment, signal wave producer is at least one, and signal wave receiver is at least two.Optimal, it draws
The every row grid divided is equipped with a signal wave producer, and corresponding every row is equipped with a signal wave receiver.
In an embodiment of the present invention, after S101 further include:
S1101 obtains the abscissa resolution ratio of grid and the ordinate resolution ratio of grid.
In an embodiment of the present invention, S102 includes:
S201 obtains the coordinate of the signal wave producer and the coordinate of the signal wave receiver.
S202 sets all the points on the coordinate line of the coordinate of the signal wave producer and the signal wave receiver
Preset value, wherein the difference of the abscissa between two adjacent points is equal.
S203 obtains the number of preset value in each grid, and determines signal based on the number of preset value in each grid
Propagation stroke of the wave from signal wave producer into the communication process of signal wave receiver in each grid.
In an embodiment of the present invention, before S201 further include:
S2201 constructs initial pixel square according to the grid, the abscissa resolution ratio and the ordinate resolution ratio
Gust, the value of coordinate element is initial value in the initial pixel matrix.
In the present embodiment, abscissa resolution ratio is dx, and ordinate resolution ratio is dy.
In an embodiment of the present invention, after S201 further include:
S2202 calculates the signal using the coordinate of the signal wave producer and the coordinate of the signal wave receiver
The slope of straight line between wave producer and the signal wave receiver.
S2203, according to the coordinate of the signal wave producer, the abscissa resolution ratio, the ordinate resolution ratio and
The slope of the straight line obtains the seat of all the points on the straight line between the signal wave producer and the signal wave receiver
Mark.
In the present embodiment, coordinate of the signal wave producer on the grid are as follows:
Pi=(1, (mi-s)×dy)；
Wherein, PiFor the coordinate of the signal wave producer i；miFor the signal wave producer i grid m row；s
For position of the signal wave producer i on m row；Dy is the ordinate resolution ratio of the grid；
Coordinate of the signal wave receiver on the grid are as follows:
Qj=(n × dx, (nj-e)×dy)；
Wherein, QjFor the coordinate of the signal wave receiver j；njIt is arranged for the signal wave producer j the n-th of grid；e
For position of the signal wave producer j in line n；Dx is the abscissa resolution ratio of the grid；
The slope of straight line between the signal wave producer and the signal wave receiver:
Wherein, θ is the slope of the straight line；
The coordinate of all the points includes: on straight line between the signal wave producer and the signal wave receiver
The abscissa U of the pointkN × dx, the ordinate W of the point are arrived for 1kAre as follows:
Wk=(mi-s)×dy+tanθ×(Uk-1)；
Wherein, WkFor the ordinate of point k；UkFor the abscissa of point k.
In an embodiment of the present invention, S202 includes:
According to the coordinate of the point, coordinate member corresponding with the coordinate of the point is searched in the initial pixel matrix
Element, and the initial value of the coordinate element is replaced with preset value.
In the present embodiment, initial pixel matrix is m × dy row, the matrix of n × dx column, initial pixel matrix all elements
Initial value be 0 etc..
In an embodiment of the present invention, S203 includes:
Based on formulaCalculating signal wave is from signal wave producer into signal wave receiver communication process
Propagation stroke in each grid；
Wherein,When grid (a, b) passed through by the straight line between signal wave producer i and signal wave receiver j
Propagation stroke in the grid (a, b)；Number for the preset value for including in the grid (a, b)；θ is described
The slope of straight line between signal wave producer i and signal wave receiver j.
In an embodiment of the present invention, S103 includes:
S301, based on the detection zone after grid division, determine propagation stroke of the signal wave in each grid and
The loss function that transmission velocity of wave of the signal wave in the carrier to be detected is met, and the loss function is carried out just
Then change processing, the loss function after obtaining Regularization；
S302, obtains constraint condition, and the constraint condition is the biography for obtaining the signal wave in the carrier to be detected
The condition of the optimal solution of defeated velocity of wave；
S303, based on after regularization loss function and the constraint condition, determine the signal wave described to be detected
Transmission velocity of wave in carrier.
In an embodiment of the present invention, S301 includes:
The inverse of speed is indicated with slowness, rmn=1/vmn；vmnFor transmission speed of the signal wave in grid (m, n)；Slowly
Degree phasor is R={ r11,r12,…,,rmn}T；
It calculates, above formula can be simplified for simplicity are as follows:
th=tpq；
H=m2(m-p+1) × m+q, p ∈ [1, m], q ∈ [1, m], m are the line number of grid, and n is the columns of grid.
Wherein, tpqTransmission time of the signal wave in grid (p, q)；Wherein, m is the line number of the grid；N is the net
The columns of lattice；Respectively first signal wave producer divides to the straight line between first signal wave receiver
The length of transmission stroke not in grid (1,1), grid (1,2) ... in grid (m, n)；Respectively first
Signal wave producer to the straight line between second signal wave receiver respectively grid (1,1), grid (1,2) ... grid (m,
N) length of the transmission stroke in；Respectively m-th of signal wave producer is received to m-th of signal wave
The length of straight line between the device transmission stroke in grid (1,1), grid (1,2) ... in grid (m, n) respectively；t11For signal wave
From first signal wave producer to the time of first signal wave receiver；t12It is signal wave from first signal wave producer
To the time of second signal wave receiver；tmmFor signal wave from m-th of signal wave producer to m-th signal wave receiver
Time；R is slowness phasor；
R={ r11,r12,…,rmn}T；r11For in grid (1,1) at a slow speed；r12For in grid (1,2) at a slow speed；rmnFor
In grid (m, n) at a slow speed；thFor the transmission time of the signal wave on the h articles straight line.
fh(R)=th；h∈[1,m2]；
Wherein, D (R) is loss function；
By loss function D (R) regularization, obtain:
Wherein: F (R) is the loss function after regularization；λ is constant, generally 0.2.
By way of example it is assumed that m=3, n=2；Utilize h=m2(m-p+1) × m+q, p ∈ [1, m], q ∈ [1, m],And fh(R)=th；It can obtain:
Work as p=1, when q=1, h=1,
Work as p=1, when q=2, h=1,
.
.
.
Work as p=3, when q=3, h=9,
In an embodiment of the present invention, S302 includes:
It is R={ r for slowness phasor11,r12,…,,rmn}TDefine an initial value R0, the initial value of F (R) is F (R0)；
Constraint condition includes:
Rh+1=Rh-αF'(Rh)；
F(Rh+1)=F (Rh)+(Rh+1-Rh)F'(Rh)；
Wherein: F (Rh+1) it is objective function；F'(Rh) it is tangent slope；α is constant coefficient.
In an embodiment of the present invention, S303 includes:
S3301 constantly updates slowness phasor and loss function, until | F (Rh+1)-F(Rh) | < 0.01 terminates iteration.
S3302, finally according to slowness rmn=1/vmnFind out transmission velocity of wave v of the signal wave in the carrier to be detectedmn。
It should be understood that the size of the serial number of each step is not meant that the order of the execution order in above-described embodiment, each process
Execution sequence should be determined by its function and internal logic, the implementation process without coping with the embodiment of the present invention constitutes any limit
It is fixed.
Embodiment 2:
As shown in figure 3, the detection system 100 for the signal wave transmission velocity of wave that one embodiment of the present of invention provides, for holding
Method and step in embodiment corresponding to row Fig. 2 comprising:
Division module 110, for the corresponding detection zone of carrier to be detected to be divided into multiple grids, the load to be detected
Body is the carrier of signal wave transmission；
First computing module 120, for being connect based on position of the signal wave producer in the detection zone and signal wave
Position of the device in the detection zone is received, propagation stroke of the signal wave in each grid is calculated, wherein the signal wave hair
Raw device is located at the first side in the detection zone, and the signal wave receiver is located at second side in the detection zone, institute
Second side of the first side and the detection zone of stating detection zone is the opposite two sides of the detection zone；
Second computing module 130 determines that the signal wave exists for the propagation stroke based on signal wave in each grid
Transmission velocity of wave in the carrier to be detected.
In an embodiment of the present invention, the first computing module 120 includes:
Coordinate acquiring unit, for obtaining the coordinate of the signal wave producer and the coordinate of the signal wave receiver；
First computing unit, for by the coordinate line of the coordinate of the signal wave producer and the signal wave receiver
On all the points set preset value, wherein the difference of the abscissa between two adjacent points is equal；
Second computing unit, for obtaining the number of preset value in each grid, and based on preset value in each grid
Number determines propagation stroke of the signal wave from signal wave producer into the communication process of signal wave receiver in each grid.
In an embodiment of the present invention, the coordinate of signal wave producer are as follows: Pi=(1, (mi-s)×dy)；The signal wave
The coordinate of receiver are as follows: Qj=(n × dx, (nj-e)×dy)；
Wherein, PiFor the coordinate of the signal wave producer i；miFor the signal wave producer i grid m row；s
For position of the signal wave producer i on m row；Dy is the ordinate resolution ratio of the grid；QjFor the signal wave
The coordinate of receiver j；njIt is arranged for the signal wave producer j the n-th of grid；E is the signal wave producer j in line n
On position；Dx is the abscissa resolution ratio of the grid.
In an embodiment of the present invention, the second computing unit includes:
Based on formulaCalculating signal wave is from signal wave producer into signal wave receiver communication process
Propagation stroke in each grid；
Wherein,When grid (a, b) passed through by the straight line between signal wave producer i and signal wave receiver j
Propagation stroke in the grid (a, b)；Number for the preset value for including in the grid (a, b)；θ is described
The slope of straight line between signal wave producer i and signal wave receiver j.
In an embodiment of the present invention, the second computing module 130 includes:
Third computing unit, for determining the signal wave in each grid based on the detection zone after grid division
The loss function that is met of transmission velocity of wave in the carrier to be detected of propagation stroke and the signal wave, and to the damage
It loses function and carries out Regularization, the loss function after obtaining Regularization；
Condition acquiring unit, for obtaining constraint condition, the constraint condition is to obtain the signal wave described to be checked
Survey the condition of the optimal solution of the transmission velocity of wave in carrier；
4th computing unit, for based on after regularization loss function and the constraint condition, determine the signal wave
Transmission velocity of wave in the carrier to be detected.
In an embodiment of the present invention, third computing unit includes:
Wherein, fh(R)=th；h∈[1,m2]；D (R) is loss function；
R={ r11,r12,…,,rmn}T；
th=tpq；
Wherein, h=m2(m-p+1) × m+q, p ∈ [1, m], q ∈ [1, m], m are the line number of grid, and n is the column of grid
Number；tpqTransmission time of the signal wave in grid (p, q)；
Slowness rmn=1/vmn；vmnFor transmission speed of the signal wave in grid (m, n)；
Correspondingly, the loss function after the regularization are as follows:
Wherein: F (R) is the loss function after regularization；λ is constant.
In an embodiment of the present invention, constraint condition includes:
Rh+1=Rh-αF'(Rh)；
F(Rh+1)=F (Rh)+(Rh+1-Rh)F'(Rh)；
Wherein: F (Rh+1) it is objective function；F'(Rh) it is tangent slope；α is constant coefficient.
It is apparent to those skilled in the art that for convenience and simplicity of description, only with above-mentioned each function
The division progress of module can according to need and for example, in practical application by above-mentioned function distribution by different function moulds
Block is completed, i.e. the internal structure of the detection system of signal wave transmission velocity of wave is divided into different functional modules, with complete with
The all or part of function of upper description.Each functional module in embodiment can integrate in one processing unit, can also be with
It is that each unit physically exists alone, can also be integrated in one unit with two or more units, above-mentioned integrated mould
Block both can take the form of hardware realization, can also realize in the form of software functional units.In addition, each functional module
Specific name is also only for convenience of distinguishing each other, the protection scope being not intended to limit this application.Above-mentioned signal wave transmission wave
The specific work process of module in the detection system of speed, can be with reference to the corresponding process in preceding method embodiment 1, herein no longer
It repeats.
Embodiment 3:
Fig. 4 is the schematic diagram for the terminal device that one embodiment of the invention provides.As shown in figure 4, the terminal of the embodiment is set
Standby 4 include: processor 40, memory 41 and are stored in the meter that can be run in the memory 41 and on the processor 40
Calculation machine program 42.The processor 40 is realized in each embodiment as described in example 1 above when executing the computer program 42
The step of, such as step S101 to S103 shown in Fig. 2.Alternatively, reality when the processor 40 executes the computer program 42
The function of each module/unit in each system embodiment now as described in example 2 above, such as module 110 to 130 shown in Fig. 3
Function.
The terminal device 4 refers to the terminal with data-handling capacity, including but not limited to computer, work station, clothes
Business device, the smart phone more even haveing excellent performance, palm PC, tablet computer, personal digital assistant (PDA), intelligence electricity
Depending on (Smart TV) etc..Operating system is generally fitted on terminal device, including but not limited to: Windows operating system,
LINUX operating system, Android (Android) operating system, Symbian operating system, Windows mobile operating system, with
And iOS operating system etc..The specific example of terminal device 4 is enumerated in detail above, it will be appreciated by those of skill in the art that
Terminal device is not limited to above-mentioned enumerate example.
The terminal device may include, but be not limited only to, processor 40, memory 41.Those skilled in the art can manage
Solution, Fig. 4 is only the example of terminal device 4, does not constitute the restriction to terminal device 4, may include more or more than illustrating
Few component perhaps combines certain components or different components, such as the terminal device 4 can also include input and output
Equipment, network access equipment, bus etc..
Alleged processor 40 can be central processing unit (Central Processing Unit, CPU), can also be
Other general processors, digital signal processor (Digital Signal Processor, DSP), specific integrated circuit
(Application Specific Integrated Circuit, ASIC), ready-made programmable gate array (Field-
Programmable Gate Array, FPGA) either other programmable logic device, discrete gate or transistor logic,
Discrete hardware components etc..General processor can be microprocessor or the processor is also possible to any conventional processor
Deng.
The memory 41 can be the internal storage unit of the terminal device 4, such as the hard disk or interior of terminal device 4
It deposits.The memory 41 is also possible to the External memory equipment of the terminal device 4, such as be equipped on the terminal device 4
Plug-in type hard disk, intelligent memory card (Smart Media Card, SMC), secure digital (Secure Digital, SD) card dodge
Deposit card (Flash Card) etc..Further, the memory 41 can also both include the storage inside list of the terminal device 4
Member also includes External memory equipment.The memory 41 is for storing needed for the computer program and the terminal device 4
Other programs and data.The memory 41 can be also used for temporarily storing the data that has exported or will export.
Embodiment 4:
The embodiment of the invention also provides a kind of computer readable storage medium, computer-readable recording medium storage has meter
Calculation machine program is realized the step in each embodiment as described in example 1 above, such as is schemed when computer program is executed by processor
Step S101 shown in 2 to step S103.Alternatively, realizing when the computer program is executed by processor such as institute in embodiment 3
The function of each module/unit in each system embodiment stated, such as the function of module 110 to 130 shown in Fig. 3.
The computer program can be stored in a computer readable storage medium, and the computer program is by processor
When execution, it can be achieved that the step of above-mentioned each embodiment of the method.Wherein, the computer program includes computer program code,
The computer program code can be source code form, object identification code form, executable file or certain intermediate forms etc..Institute
State computer-readable medium may include: can carry the computer program code any entity or device, recording medium,
USB flash disk, mobile hard disk, magnetic disk, CD, computer storage, read-only memory (ROM, Read-Only Memory), arbitrary access
Memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium etc..
In the above-described embodiments, it all emphasizes particularly on different fields to the description of each embodiment, embodiment 1 to 4 can in any combination, group
The new embodiment formed after conjunction is also within the scope of protection of this application.There is no the portion for being described in detail or recording in some embodiment
Point, it may refer to the associated description of other embodiments.
Those of ordinary skill in the art may be aware that list described in conjunction with the examples disclosed in the embodiments of the present disclosure
Member and algorithm steps can be realized with the combination of electronic hardware or computer software and electronic hardware.These functions are actually
It is implemented in hardware or software, the specific application and design constraint depending on technical solution.Professional technician
Each specific application can be used different methods to achieve the described function, but this realization is it is not considered that exceed
The scope of the present invention.
In embodiment provided by the present invention, it should be understood that disclosed terminal device and method can pass through it
Its mode is realized.For example, system described above/terminal device embodiment is only schematical, for example, the module
Or the division of unit, only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple lists
Member or component can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point,
Shown or discussed mutual coupling or direct-coupling or communication connection can be through some interfaces, device or unit
INDIRECT COUPLING or communication connection, can be electrical property, mechanical or other forms.
Embodiment described above is merely illustrative of the technical solution of the present invention, rather than its limitations；Although referring to aforementioned reality
Applying example, invention is explained in detail, those skilled in the art should understand that: it still can be to aforementioned each
Technical solution documented by embodiment is modified or equivalent replacement of some of the technical features；And these are modified
Or replacement, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution should all
It is included within protection scope of the present invention.
Claims (10)
1. a kind of detection method of signal wave transmission velocity of wave characterized by comprising
The corresponding detection zone of carrier to be detected is divided into multiple grids, the carrier to be detected is the load of signal wave transmission
Body；
Based on signal wave producer in the position of position and signal wave receiver in the detection zone in the detection zone
It sets, calculates propagation stroke of the signal wave in each grid, wherein the signal wave producer is located in the detection zone
First side, the signal wave receiver are located at second side in the detection zone, the first side of the detection zone and described
Second side of detection zone is the opposite two sides of the detection zone；
Propagation stroke based on signal wave in each grid determines transmission wave of the signal wave in the carrier to be detected
Speed.
2. the detection method of signal wave transmission velocity of wave as described in claim 1, which is characterized in that described to be based on signal wave
Device calculates signal wave each in the position of position and signal wave receiver in the detection zone in the detection zone
Propagation stroke in grid includes:
Obtain the coordinate of the signal wave producer and the coordinate of the signal wave receiver；
All the points on the coordinate line of the coordinate of the signal wave producer and the signal wave receiver are set into preset value,
In, the difference of the abscissa between two adjacent points is equal；
The number of preset value in each grid is obtained, and determines signal wave from signal wave based on the number of preset value in each grid
Propagation stroke of the generator into the communication process of signal wave receiver in each grid.
3. the detection method of signal wave transmission velocity of wave as claimed in claim 2, which is characterized in that the signal wave producer
Coordinate are as follows: Pi=(1, (mi-s)×dy)；The coordinate of the signal wave receiver are as follows: Qj=(n × dx, (nj-e)×dy)；
Wherein, PiFor the coordinate of signal wave producer i；miFor signal wave producer i grid m row；S is signal wave
Position of the device i on m row；Dy is the ordinate resolution ratio of grid；QjFor the coordinate of signal wave receiver j；njFor signal wave hair
N-th column of the raw device j in grid；E is position of the signal wave producer j in line n；Dx is the abscissa resolution ratio of grid.
4. the detection method of signal wave transmission velocity of wave as claimed in claim 2, which is characterized in that described based in each grid
The number of preset value determine signal wave from signal wave producer into the communication process of signal wave receiver in each grid
Propagating stroke includes:
Based on formulaCalculating signal wave is from signal wave producer into signal wave receiver communication process every
Propagation stroke in a grid；
Wherein,In grid when grid (a, b) passed through by the straight line between signal wave producer i and signal wave receiver j
Propagation stroke in (a, b)；Number for the preset value for including in grid (a, b)；θ is signal wave producer i and signal
The slope of straight line between wave receiver j.
5. the detection method of signal wave transmission velocity of wave as described in claim 1, which is characterized in that the signal wave that is based on is every
Propagation stroke in a grid determines that transmission velocity of wave of the signal wave in the carrier to be detected includes:
Based on the detection zone after grid division, propagation stroke and the signal wave of the signal wave in each grid are determined
The loss function that transmission velocity of wave in the carrier to be detected is met, and Regularization is carried out to the loss function,
Loss function after obtaining Regularization；
Constraint condition is obtained, the constraint condition is to obtain transmission velocity of wave of the signal wave in the carrier to be detected most
The condition of excellent solution；
Based on after regularization loss function and the constraint condition, determine biography of the signal wave in the carrier to be detected
Defeated velocity of wave.
6. the detection method of signal wave transmission velocity of wave as claimed in claim 5, which is characterized in that the signal wave is in each net
The loss function that the transmission velocity of wave of propagation stroke and the signal wave in the carrier to be detected in lattice is met:
Wherein, fh(R)=th；h∈[1,m2]；D (R) is loss function；
R={ r11,r12,…,,rmn}T；
th=tpq；
Wherein, h=m2(m-p+1) × m+q, p ∈ [1, m], q ∈ [1, m], m are the line number of grid, and n is the columns of grid；tpq
Transmission time of the signal wave in grid (p, q)；
Slowness rmn=1/vmn；vmnFor transmission speed of the signal wave in grid (m, n)；
Correspondingly, the loss function after the regularization are as follows:
Wherein: F (R) is the loss function after regularization；λ is constant.
7. the detection method of signal wave as claimed in claim 6 transmission velocity of wave, which is characterized in that the constraint condition includes:
Rh+1=Rh-αF'(Rh)；
F(Rh+1)=F (Rh)+(Rh+1-Rh)F'(Rh)；
Wherein: F (Rh+1) it is objective function；F'(Rh) it is tangent slope；α is constant coefficient.
8. a kind of detection system of signal wave transmission velocity of wave characterized by comprising
Division module, for the corresponding detection zone of carrier to be detected to be divided into multiple grids, the carrier to be detected is letter
The carrier of number wave transmission；
First computing module, for based on position of the signal wave producer in the detection zone and signal wave receiver in institute
The position in detection zone is stated, calculates propagation stroke of the signal wave in each grid, wherein the signal wave producer is located at
The first side in the detection zone, the signal wave receiver are located at second side in the detection zone, the detection zone
First side in domain and second side of the detection zone are the opposite two sides of the detection zone；
Second computing module, for the propagation stroke based on signal wave in each grid, determine the signal wave it is described to
Detect the transmission velocity of wave in carrier.
9. a kind of terminal device, which is characterized in that in the memory and can be in institute including memory, processor and storage
The computer program run on processor is stated, the processor realizes such as claim 1 to 7 times when executing the computer program
The step of detection method of one signal wave transmission velocity of wave.
10. a kind of computer readable storage medium, which is characterized in that the computer-readable recording medium storage has computer journey
Sequence realizes the inspection of the signal wave transmission velocity of wave as described in any one of claim 1 to 7 when the computer program is executed by processor
The step of survey method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910554692.7A CN110196281A (en) | 2019-06-25 | 2019-06-25 | A kind of detection method, system and the terminal device of signal wave transmission velocity of wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910554692.7A CN110196281A (en) | 2019-06-25 | 2019-06-25 | A kind of detection method, system and the terminal device of signal wave transmission velocity of wave |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110196281A true CN110196281A (en) | 2019-09-03 |
Family
ID=67755074
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910554692.7A Pending CN110196281A (en) | 2019-06-25 | 2019-06-25 | A kind of detection method, system and the terminal device of signal wave transmission velocity of wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110196281A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114966855A (en) * | 2022-07-28 | 2022-08-30 | 煤炭科学研究总院有限公司 | Method, device, equipment and medium for detecting high stress area of coal rock mass |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1908651A (en) * | 2006-08-03 | 2007-02-07 | 长安大学 | Concrete ultrasound tomography algorithm |
| CN101581789A (en) * | 2009-06-23 | 2009-11-18 | 刘盛东 | Mine working face inter-lane seismic wave CT detection method |
| CN103258134A (en) * | 2013-05-14 | 2013-08-21 | 宁波大学 | Dimension reduction processing method of high-dimensional vibration signals |
| CN104749630A (en) * | 2015-04-01 | 2015-07-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for constructing microseism monitoring velocity model |
| CN106558086A (en) * | 2016-10-18 | 2017-04-05 | 江南大学 | A kind of velocity error correction method of forest stress wave fault imaging |
| CN106772597A (en) * | 2016-12-12 | 2017-05-31 | 中国矿业大学 | Take man-made explosion as the method for vibration signal inverting coal mine work area stress distribution |
| CN109003686A (en) * | 2018-05-25 | 2018-12-14 | 中广核工程有限公司 | Localization method, device, computer equipment and the storage medium of nuclear power station loosening part |
| CN109738525A (en) * | 2019-03-12 | 2019-05-10 | 石家庄铁道大学 | The high-speed railway bearing layer concrete velocity of wave detection device and estimation method of gridding |
-
2019
- 2019-06-25 CN CN201910554692.7A patent/CN110196281A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1908651A (en) * | 2006-08-03 | 2007-02-07 | 长安大学 | Concrete ultrasound tomography algorithm |
| CN101581789A (en) * | 2009-06-23 | 2009-11-18 | 刘盛东 | Mine working face inter-lane seismic wave CT detection method |
| CN103258134A (en) * | 2013-05-14 | 2013-08-21 | 宁波大学 | Dimension reduction processing method of high-dimensional vibration signals |
| CN104749630A (en) * | 2015-04-01 | 2015-07-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for constructing microseism monitoring velocity model |
| CN106558086A (en) * | 2016-10-18 | 2017-04-05 | 江南大学 | A kind of velocity error correction method of forest stress wave fault imaging |
| CN106772597A (en) * | 2016-12-12 | 2017-05-31 | 中国矿业大学 | Take man-made explosion as the method for vibration signal inverting coal mine work area stress distribution |
| CN109003686A (en) * | 2018-05-25 | 2018-12-14 | 中广核工程有限公司 | Localization method, device, computer equipment and the storage medium of nuclear power station loosening part |
| CN109738525A (en) * | 2019-03-12 | 2019-05-10 | 石家庄铁道大学 | The high-speed railway bearing layer concrete velocity of wave detection device and estimation method of gridding |
Non-Patent Citations (1)
| Title |
|---|
| 李勇 等: ""模型约束下的在线字典学习地震弱信号去噪方法"", 《地球物理学报》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114966855A (en) * | 2022-07-28 | 2022-08-30 | 煤炭科学研究总院有限公司 | Method, device, equipment and medium for detecting high stress area of coal rock mass |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108280685B (en) | Information acquisition method and device | |
| CN112598762B (en) | Three-dimensional lane line information generation method, device, electronic device, and medium | |
| Li et al. | A posteriori error estimates of recovery type for distributed convex optimal control problems | |
| CN111866775A (en) | Service arranging method and device | |
| Kobayashi et al. | Load balancing strategies for a parallel ray-tracing system based on constant subdivision | |
| CN110225453A (en) | Mobile terminal locating method, device, electronic equipment and storage medium | |
| Lai et al. | A parallel two-dimensional discontinuous galerkin method for shallow-water flows using high-resolution unstructured meshes | |
| Wang et al. | Testing a new adaptive grid‐based shallow flow model for different types of flood simulations | |
| CN110196281A (en) | A kind of detection method, system and the terminal device of signal wave transmission velocity of wave | |
| CN104573171B (en) | A kind of emulation mode of transmission line of electricity and system | |
| CN109699002A (en) | A kind of interior WiFi localization method, device and terminal device | |
| CN109815579A (en) | Determination method, computer readable storage medium and the terminal device of slope parameter | |
| CN114092642B (en) | Three-dimensional house type model generation method, device and equipment | |
| CN112816663A (en) | Method and device for monitoring soil water content of yellow river dam in flood control project | |
| CN108198173B (en) | A kind of online test method, device and the terminal device in distress in concrete region | |
| CN110991774A (en) | Electric quantity load prediction method and device | |
| US20100262804A1 (en) | Effective Memory Clustering to Minimize Page Fault and Optimize Memory Utilization | |
| Xiang et al. | Efficient shadowing determination at grazing incidence, for mesh‐based physical optics scattering analysis | |
| CN113030648B (en) | Power cable fault point position determining method and device and terminal equipment | |
| CN102339329A (en) | Method for dividing physical layout | |
| JP3415530B2 (en) | Electromagnetic field simulation method | |
| CN113127964A (en) | Multi-target equivalent static wind load calculation method and equipment for long-span roof structure | |
| CN113723712B (en) | Wind power prediction method, system, equipment and medium | |
| CN111224977A (en) | Process configuration method and system | |
| CN109612670A (en) | A kind of protective net monitoring method, system and terminal device |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190903 |