CN104766311B - Bin method for fast searching is blocked during a kind of sound wave irradiation random fluctuation interface - Google Patents
Bin method for fast searching is blocked during a kind of sound wave irradiation random fluctuation interface Download PDFInfo
- Publication number
- CN104766311B CN104766311B CN201510140918.0A CN201510140918A CN104766311B CN 104766311 B CN104766311 B CN 104766311B CN 201510140918 A CN201510140918 A CN 201510140918A CN 104766311 B CN104766311 B CN 104766311B
- Authority
- CN
- China
- Prior art keywords
- bin
- random fluctuation
- central point
- interface
- max
- 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
Abstract
The invention belongs to underwater target scattering field, and in particular to be a kind of scattering acoustic field at numerical computations random fluctuation interface when, search not by sound wave be irradiated to block bin sound wave irradiate random fluctuation interface when block bin method for fast searching.The present invention includes：Import N number of bin nodal information, bin central point information and the sound source position information at the random fluctuation interface for having carried out Triangular object model division；Ask for the maximum Z of z-component in all bin central point three-dimensional coordinatesmax, obtain z=ZmaxPlane etc..The region of search when present invention is due to judging whether bin is blocked reduces, and reduces amount of calculation, improves the search efficiency for blocking bin；Suitable for the occlusion issue at random fluctuation interface, and fluctuating interfacial area is bigger, and computational efficiency improves more obvious；The relevant calculating of point, line, surface in space geometry is only make use of, without the Processing Algorithm of complexity, is easy to program realization.
Description
Technical field
The invention belongs to underwater target scattering field, and in particular to be a kind of numerical computations random fluctuation interface dissipate
When penetrating sound field, the bin that the blocks when sound wave for blocking bin that search is not irradiated to by sound wave irradiates random fluctuation interface is quickly searched
Suo Fangfa.
Background technology
Using physical acoustic method estimate fluctuating interface or and its nearby target scattering acoustic field when, it is necessary to bin draw
Fluctuating INTERFACE MODEL after point block the search of bin, and then weeds out the bin being blocked.Traditional method is to carry out
Traversal search, judge whether mutually blocked between each bin, therefore amount of calculation is larger.Also used in electromagnetic scattering field
Some block the method for bin fast search, and such as (Tao Yang, Chen Yulin, aircraft engine complicated cavity block topological trace method
Algorithm research, radar science and technology, the 2nd phase of volume 2013,11), bin group technology (Ji Jinzu, Liu Zhanhe, based on bin
The electromagnetism of packet blocks algorithm and its optimization, BJ University of Aeronautics ＆ Astronautics's journal, the 4th phase of volume 2009,35) etc..Above method
The method for fast searching proposed both for regular geometric body target, the present invention is blocked for the bin at random fluctuation interface to be asked
Topic, bin method for fast searching is blocked when having invented a kind of sound wave irradiation random fluctuation interface.
The content of the invention
Solve that prior art search procedure is computationally intensive, calculates overlong time asks it is an object of the invention to provide a kind of
Bin method for fast searching is blocked during the sound wave irradiation random fluctuation interface of topic.
The object of the present invention is achieved like this：
(1) N number of bin nodal information, the bin central point letter at the random fluctuation interface for having carried out Triangular object model division are imported
Breath and sound source position information；
(2) the maximum Z of z-component in all bin central point three-dimensional coordinates is asked formax, obtain z=ZmaxPlane；
(3) for n-th of bin, connection sound source and its central point, line and z=Z are tried to achievemaxPlane intersection point (x', y',
Zmax), the subpoint found intersection in xoy planes is (x', y'), and n-th of bin central point is in the subpoint of xoy planes
(xn,yn)；
(4) if x'≤xn, make x'*=x'- Δs l, xn*=xn+Δl；If x'＞ xn, make x'*=x'+ Δs l, xn*=
xn-Δl；If y'≤yn, make y'*=y'- Δs l, yn*=yn+Δl；If y'＞ yn, make y'*=y'+ Δs l, yn*=yn-Δ
L, wherein Δ l are the bin length of side, and it is respectively (x'*, y'*), (x that summit is obtained in xoy planesn*,y'*)、(xn*,yn*)、
(x'*,yn* rectangular region)；
(5) bin central point subpoint in xoy planes is (x, y) in random fluctuation interface,All subpoints are taken to meet x
∈[x'*,xn] and y ∈ [y'*, y *n*] the bin corresponding to condition, obtain corresponding with the rectangular region tried to achieve in step (4)
M bin；
(6) for m-th of bin in M bin in step (5), plane where asking it and sound source in step (3) with
The intersection point P of n-th of bin central point line, tries to achieve βi
Wherein ri(i=1,2,3, r4=r1) for i-th of vertex position vector of m-th bin, rpFor intersection point P position
Vector, vmFor the normal vector of m-th of bin, βiFor the judgement parameter tried to achieve, if all of βi＞ 0, then n-th of bin quilt
M-th of bin blocks, otherwise is not blocked by m-th of bin；
(7) repeat step (6), the hiding relation of n-th of bin and whole M bins is judged, if be not blocked,
Retain n-th of bin, otherwise delete n-th of bin；
(8) repeat step (3)~step (7), all bins remained are incidence in N number of bin at fluctuating interface
The bin that sound wave illuminates, the bin of deletion are the bin being blocked.
The beneficial effects of the present invention are：
The present invention is directed to the bin occlusion issue at random fluctuation interface, there is provided one kind blocks bin method for fast searching, solution
The problem of existing method for fast searching of having determined can only be directed to regular geometric body target.Compared to traditional traversal search method, sheet
Invention has advantages below：
Region of search (bin quantity) when the 1st, due to judging whether bin is blocked reduces, and reduces amount of calculation, improves
Block the search efficiency of bin；
2nd, suitable for the occlusion issue at random fluctuation interface, and fluctuating interfacial area is bigger, and computational efficiency improves brighter
It is aobvious；
3rd, the relevant calculating of point, line, surface in space geometry is only make use of, without the Processing Algorithm of complexity, is easy to program reality
It is existing.
Brief description of the drawings
Fig. 1 is the random fluctuation interface for having carried out Triangular object model division.
Fig. 2 is the rectangular region in the xoy planes obtained by subpoint.
Fig. 3 is M bin corresponding with the rectangular region in xoy planes in random fluctuation interface.
Fig. 4 is to weed out the random fluctuation interface after the bin that is blocked.
The flow chart for blocking bin method for fast searching when Fig. 5 is sound wave irradiation fluctuating interface.
Embodiment
Below in conjunction with the accompanying drawings 1~5 and example the present invention is described in further detail.
The invention belongs to during the scattering acoustic field in underwater target scattering field, mainly numerical computations random fluctuation interface,
Search for the method for blocking bin not being irradiated to by sound wave.The present invention includes：The bin N for judging whether to be blocked for needs,
Seek the maximum Z of bin center point coordinate z-componentmax, obtain z=ZmaxPlane；Sound source and bin central point are connected, obtains line
With z=ZmaxThe intersection point of plane；Rectangular region is built in the subpoint of xoy planes according to intersection point and bin central point, and then obtained
All bin set corresponding with rectangular region on to random fluctuation interface；Bin N and the bin collection tried to achieve are judged successively
The hiding relation of bin in conjunction.Present invention decreases bin region of search, improves the search efficiency for blocking bin.
To reach above-mentioned purpose, the present invention uses following technical scheme：
1st, N number of bin nodal information and bin the central point letter at the random fluctuation interface for having carried out Triangular object model division are imported
Breath, and sound source position information.
2nd, the maximum Z of z-component in all bin central point three-dimensional coordinates is asked formax, obtain z=ZmaxPlane.
3rd, for n-th of bin, connection sound source and its central point, line and z=Z are tried to achievemaxPlane intersection point (x', y',
Zmax).The subpoint found intersection in xoy planes is (x', y'), and n-th of bin central point is in the subpoint of xoy planes
(xn,yn)。
If the 4th, x'≤xn, make x'*=x'- Δs l, xn*=xn+Δl；If x'＞ xn, make x'*=x'+ Δs l, xn*=
xn-Δl；If y'≤yn, make y'*=y'- Δs l, yn*=yn+Δl；If y'＞ yn, make y'*=y'+ Δs l, yn*=yn-Δ
l.Wherein Δ l is the bin length of side.It is respectively (x'*, y'*), (x that summit is obtained in xoy planesn*,y'*)、(xn*,yn*)、
(x'*,yn* rectangular region).
5th, bin central point subpoint in xoy planes is (x, y) in random fluctuation interface,All subpoints are taken to meet x
∈[x'*,xn] and y ∈ [y'*, y *n*] the bin corresponding to condition, that is, obtain corresponding with the rectangular region tried to achieve in step 4
M bin.
6th, for m-th of bin in M bin in step 5, plane where asking it and sound source in step 3 and n-th
The intersection point P of bin central point line.β is tried to achieve by following formulai
Wherein ri(i=1,2,3, r4=r1) for i-th of vertex position vector of m-th bin, rpFor intersection point P position
Vector, vmFor the normal vector of m-th of bin, βiFor the judgement parameter tried to achieve.If all of βi＞ 0, then n-th of bin quilt
M-th of bin blocks, otherwise is not blocked by m-th of bin.
7th, repeat step 6, judge the hiding relation of n-th of bin and whole M bins, if be not blocked, protect
N-th of bin is stayed, otherwise deletes n-th of bin.
8th, above step 3~7 is repeated, all bins remained are incident acoustic wave in N number of bin at fluctuating interface
The bin illuminated, the bin of deletion are the bin being blocked.
Embodiment
The first step, with reference to Fig. 1, random fluctuation interface x directions and y directions length are 3m, have divided N=15842 altogether
Bin, sound source coordinate are (0, -5,3), and unit m, the unit of all coordinates is m below.
Second step, the maximum for trying to achieve all bin central point z-components is z=0.4885m, and z=0.4885m is three-dimensional space
Between in a plane.
3rd step, for the n-th=1 bin, connection sound source coordinate (0, -5,3) and its center point coordinate (xn,yn,zn)=
(- 1.1097,0.8040, -0.0679), it is (x', y', Z to try to achieve line and the intersection point of z=0.4885 planesmax)=(-
0.9084,-0.2487,0.4885).Find intersection subpoint in x o y planes for (x', y')=(- 0.9084 ,-
0.2487), the n-th=1 bin central point is (x in the subpoint of x o y planesn,yn)=(- 1.1097,0.8040).
4th step, with reference to Fig. 2, according to following logical relation：If x'≤xn, make x'*=x'- Δs l, xn*=xn+Δl；
If x'＞ xn, make x'*=x'+ Δs l, xn*=xn-Δl；If y'≤yn, make y'*=y'- Δs l, yn*=yn+Δl；If
Y'＞ yn, make y'*=y'+ Δs l, yn*=yn- Δ l, and Δ l=0.0333 is taken, can obtain the summit in xoy planes is respectively
(x'*, y'*)=(- 0.8751, -0.2821), (xn*, y'*)=(- 1.1430, -0.2821), (xn*,yn*)=(- 1.1430,
0.8373)、(x'*,yn*)=(- 0.8751,0.8373) rectangle.
5th step, with reference to Fig. 3, the projection of all bin central points is taken to meet x ∈ [x'*, xn] and y ∈ [y'*, y *n*] institute
Corresponding bin, that is, obtain M=543 bin corresponding with the rectangular region tried to achieve in step 4.
6th step, for the m=1 bin in M bin in the 5th step, plane where asking it and sound in the 3rd step
Source and the intersection point P of the n-th=1 bin central point line.Judge parameter betaiTried to achieve by following formula
Wherein ri(i=1,2,3, r4=r1) for each vertex position vector of the m=1 bin,rpSweared for intersection point P position
Amount, vmFor the normal vector of the m=1 bin, βiFor the judgement parameter tried to achieve.If all of βi＞ 0, then the n-th=1 face
Member is blocked by the m=1 bin, otherwise is not blocked.In this example, rp=(- 0.9084, -0.2487,0.4885), r1=(-
0.9167,0.2833,0.2333), r2=(- 0.9167,0.3167,0.2804), r3=(- 0.9500,0.2833,0.2431),
r4=r1, vm=(0.1675, -0.8045,0.5698), tries to achieve β1=0.0052, β2=-0.0017, β3=-0.0016, not
All βi＞ 0, so the n-th=1 bin is not blocked by the m=1 bin.
7th step, the 6th step content is repeated, judge the n-th=1 bin and remaining m=2 ... M bin blocks pass
System, if be not blocked, retain the n-th=1 bin, otherwise delete the n-th=1 bin.In this example, the n-th=1 face
Member is blocked by the m=71 bin, then the n-th=1 bin is deleted in N number of bin at whole random fluctuation interface.
8th step, with reference to Fig. 4, repeat the step of the above the 3rd to the 7th, N number of bin that judges the n-th=2 ... whether by it
Its bin blocks, and all bins for remaining are the bin that incident acoustic wave illuminates, and the bin of deletion is the face being blocked
Member.
Fig. 5 is the flow chart that the present invention is " blocking bin method for fast searching when sound wave irradiates fluctuating interface ".
Claims (1)
1. block bin method for fast searching during a kind of sound wave irradiation random fluctuation interface, it is characterised in that：
(1) import the N number of bin nodal information at random fluctuation interface for having carried out Triangular object model division, bin central point information with
And sound source position information；
(2) the maximum Z of z-component in all bin central point three-dimensional coordinates is asked formax, obtain z=ZmaxPlane；
(3) for n-th of bin, connection sound source and its central point, line and z=Z are tried to achievemaxPlane intersection point (x', y',
Zmax), the subpoint found intersection in xoy planes is (x', y'), and n-th of bin central point is in the subpoint of xoy planes
(xn,yn)；
(4) if x'≤xn, make x'*=x'- Δs l, xn *=xn+Δl；If x'＞ xn, make x'*=x'+ Δs l, xn *=xn-Δl；
If y'≤yn, make y'*=y'- Δs l, yn *=yn+Δl；If y'＞ yn, make y'*=y'+ Δs l, yn *=yn- Δ l, wherein Δ
L is the bin length of side, and it is respectively (x' that summit is obtained in xoy planes*,y'*)、(xn *,y'*)、(xn *,yn *)、(x'*,yn *) length
Square region；
(5) bin central point subpoint in xoy planes is (x, y) in random fluctuation interface, takes all subpoints to meet x ∈
[x'*,xn *] and y ∈ [y'*,yn *] bin corresponding to condition, obtain M corresponding with the rectangular region tried to achieve in step (4)
Individual bin；
(6) for m-th of bin in M bin in step (5), plane where asking it and sound source in step (3) and n-th
The intersection point P of bin central point line, tries to achieve βi
Wherein riFor i-th of vertex position vector of m-th of bin, i=1,2,3,4, r4=r1, rpFor intersection point P position vector,
vmFor the normal vector of m-th of bin, βiFor the judgement parameter tried to achieve, if all of βi＞ 0, then n-th of bin is by m-th
Bin blocks, otherwise is not blocked by m-th of bin；
(7) repeat step (6), judge the hiding relation of n-th of bin and whole M bins, if be not blocked, retain
N-th of bin, on the contrary delete n-th of bin；
(8) repeat step (3)~step (7), all bins remained are incident acoustic wave in N number of bin at fluctuating interface
The bin illuminated, the bin of deletion are the bin being blocked.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510140918.0A CN104766311B (en) | 2015-03-27 | 2015-03-27 | Bin method for fast searching is blocked during a kind of sound wave irradiation random fluctuation interface |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510140918.0A CN104766311B (en) | 2015-03-27 | 2015-03-27 | Bin method for fast searching is blocked during a kind of sound wave irradiation random fluctuation interface |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104766311A CN104766311A (en) | 2015-07-08 |
| CN104766311B true CN104766311B (en) | 2017-12-19 |
Family
ID=53648120
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510140918.0A Active CN104766311B (en) | 2015-03-27 | 2015-03-27 | Bin method for fast searching is blocked during a kind of sound wave irradiation random fluctuation interface |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104766311B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5763837A (en) * | 1994-10-11 | 1998-06-09 | Davignon; Andre | Automatic biometric data system |
| CN103761416A (en) * | 2013-12-23 | 2014-04-30 | 哈尔滨工程大学 | Method for rapidly forecasting scattering sound field of underwater triangle corner reflector |
| CN103853914A (en) * | 2013-12-23 | 2014-06-11 | 哈尔滨工程大学 | Method for rapidly estimating scattering sound field of underwater circular corner reflector |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7034831B2 (en) * | 2001-02-05 | 2006-04-25 | National Instruments Corporation | System and method for generating a low discrepancy curve in a region |
-
2015
- 2015-03-27 CN CN201510140918.0A patent/CN104766311B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5763837A (en) * | 1994-10-11 | 1998-06-09 | Davignon; Andre | Automatic biometric data system |
| CN103761416A (en) * | 2013-12-23 | 2014-04-30 | 哈尔滨工程大学 | Method for rapidly forecasting scattering sound field of underwater triangle corner reflector |
| CN103853914A (en) * | 2013-12-23 | 2014-06-11 | 哈尔滨工程大学 | Method for rapidly estimating scattering sound field of underwater circular corner reflector |
Non-Patent Citations (2)
| Title |
|---|
| 基于时-空离散化的海底混响时间序列仿真方法;陈文剑 等;《声学技术》;20110430;第30卷(第2期);173-177 * |
| 基于面元分组的电磁遮挡算法及其优化;姬金祖 等;《北京航空航天大学学报》;20090430;第35卷(第4期);453-456 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104766311A (en) | 2015-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104318056B (en) | Soft tissue deformation and cutting simulation method based on position dynamics | |
| CN104794758B (en) | A kind of method of cutting out of 3-D view | |
| CN107633123B (en) | Method for simulating bleeding and accelerating treatment of smooth particle fluid dynamics | |
| CN105389850B (en) | A kind of observability generation method of extensive three-dimensional scenic | |
| Leeper et al. | Point clouds can be represented as implicit surfaces for constraint-based haptic rendering | |
| CN105261069B (en) | Soft tissue organs member spherical model based on GPU automatically generates the method with collision detection | |
| CN109859208A (en) | Scene cut and Target Modeling method based on concavity and convexity and RSD feature | |
| CN104134236A (en) | Boolean operation method of three-dimensional planar solids | |
| CN104766311B (en) | Bin method for fast searching is blocked during a kind of sound wave irradiation random fluctuation interface | |
| Qu | Research on parallel algorithm based on AABB bounding box coordinate chain method | |
| Chuang et al. | Drawing graphs with nonuniform nodes using potential fields | |
| CN103440684B (en) | A kind of method that humorous for ball lighting is applied to iso-surface patch | |
| Wang et al. | Design and application of bounding volume hierarchy collision detection algorithm based on virtual sphere | |
| Nemec et al. | Adjoint algorithm for CAD-based shape optimization using a Cartesian method | |
| Paiva et al. | Fluid-based hatching for tone mapping in line illustrations | |
| Li et al. | Roaming path generation algorithm and optimization based on Bezier curve | |
| Vyatkin et al. | Multi-level ray casting of function-based surfaces | |
| Audette et al. | A review of mesh generation for medical simulators | |
| Meng et al. | Anisotropic Cartesian Grid Generation Strategy for Arbitrarily Complex Geometry Based on a Fully Threaded Tree | |
| CN108182675A (en) | A kind of face element shadowing method during sound wave irradiation random fluctuation interface | |
| Lizier et al. | Comparing techniques for tetrahedral mesh generation | |
| Wang et al. | Improved simulation of virtual bone drilling surgery with a voxel-based model | |
| Pan et al. | Dissection of hybrid soft tissue models using position-based dynamics | |
| Fey | Towards Effective Graph Representation Learning | |
| Gan et al. | Liver cutting algorithm based on virtual plane |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |