CN115204317B - Sound velocity profile continuation method and system based on orthogonal empirical function decomposition - Google Patents
Sound velocity profile continuation method and system based on orthogonal empirical function decomposition Download PDFInfo
- Publication number
- CN115204317B CN115204317B CN202211118337.3A CN202211118337A CN115204317B CN 115204317 B CN115204317 B CN 115204317B CN 202211118337 A CN202211118337 A CN 202211118337A CN 115204317 B CN115204317 B CN 115204317B
- Authority
- CN
- China
- Prior art keywords
- sound velocity
- velocity profile
- depth
- sea
- full
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Computing Systems (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention provides a sound velocity profile continuation method and a sound velocity profile continuation system based on orthogonal empirical function decomposition, which relate to the technical field of ocean acoustic application and are used for solving the problem that a deep sea part after continuation by the existing method does not accord with the real sound velocity distribution rule of the deep sea in a region; and constructing an extended full-sea-depth sound velocity profile according to the characteristic vector coefficients and the sample data of the full-sea-depth historical sound velocity profile, and by the mode, matching the actually-measured sound velocity value part of the sound velocity profile to be extended with high precision, and simultaneously ensuring that the sound velocity distribution of the extended part conforms to the actual sound velocity distribution rule of the region.
Description
Technical Field
The invention belongs to the technical field of ocean acoustic application, and particularly relates to a sound velocity profile continuation method and system based on orthogonal empirical function decomposition.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art that is already known to a person of ordinary skill in the art.
With the development of the underwater sound velocity field construction technology, sound velocity profile extension can be widely applied to application systems which take sound waves as signal carriers, such as underwater communication, positioning, navigation and the like. Under the condition that shallow sea and deep sea acoustic velocity profile samples coexist, deep sea partial continuation under the condition of high-precision matching of shallow sea actual measurement acoustic velocity values is carried out on the shallow sea acoustic velocity profile samples, the number of reference samples of full-sea deep acoustic velocity profiles can be increased, and the real-time estimation precision of the regional acoustic velocity field is improved.
The existing underwater sound velocity profile continuation method comprises the following steps: liu Fuchen and the like provide a sound velocity profile continuation method based on polynomial fitting, fitting an actually measured sound velocity profile of a shallow sea part by using a polynomial, and performing continuation in a deep sea area; cheng Fang and the like provide a sound velocity profile continuation method based on sound velocity matching of a deep sea isothermal layer, a full-sea deep sound velocity profile sample and a sound velocity profile sample to be continued are matched at a common part of the deep sea isothermal layer, a sound velocity profile empirical orthogonal function reconstruction process is combined, a structural coefficient of the sound velocity profile to be continued in the deep sea isothermal layer is calculated, and the distribution condition of the sound velocity profile in the deep sea isothermal layer is constructed; qu Ke provides a method and a device for complementing sound velocity profiles based on historical data and machine learning, wherein the historical data is subjected to empirical orthogonal function analysis to obtain characteristic vectors and characteristic vector coefficients, neurons representing different reference classification information are obtained through unsupervised learning method training, the association degree between an actually measured sound velocity profile and each neuron is calculated according to the Euclidean distance, and the sound velocity profile is reconstructed by using the characteristic vector coefficients corresponding to the neurons with the highest association degree; cai Kuijie and the like provide a sound velocity profile continuation method based on low-rank matrix completion, a sound velocity profile to be continued and a historical full-sea deep sound velocity profile are utilized to construct a matrix to be supplemented, and a low-rank matrix completion algorithm is applied to supplement the matrix to be supplemented, so that sound velocity profile continuation is completed.
Although various sound velocity profile continuation methods exist at present, the following problems still exist: firstly, the matching precision of the actually measured part of the sound velocity profile to be extended is not high enough; and secondly, the actually measured part of the sound velocity profile is excessively matched, so that the extended deep sea part does not accord with the real sound velocity distribution rule of the regional deep sea.
Disclosure of Invention
In order to solve the problems, the invention provides a sound velocity profile continuation method and a sound velocity profile continuation system based on orthogonal empirical function decomposition, which are used for matching the actually measured sound velocity value part of the sound velocity profile to be continued with high precision, ensuring that the sound velocity distribution of the continuation part accords with the actual sound velocity distribution rule of the region, and providing more and more reliable reference samples for the sound velocity profile inversion.
In order to achieve the above object, the present invention mainly includes the following aspects:
in a first aspect, an embodiment of the present invention provides a sound velocity profile continuation method based on orthogonal empirical function decomposition, including:
acquiring sample data of a sound velocity profile to be extended, the maximum depth of an actually measured sound velocity value of the sound velocity profile to be extended and historical sound velocity profiles of the whole sea depth;
according to the maximum depth of the actually measured sound velocity value of the sound velocity profile to be extended, carrying out depth interception on the full-sea-depth historical sound velocity profile sample data to obtain historical sound velocity profile sample interception data;
performing empirical orthogonal function decomposition on the sample truncation data of the historical sound velocity profile to obtain a first feature vector of a preset order of the sample truncation data of the historical sound velocity profile, and performing matching field processing on the sound velocity profile to be extended by using the first feature vector to obtain a feature vector coefficient corresponding to the sound velocity profile to be extended in the first feature vector;
and constructing a prolonged full-sea deep sound velocity profile according to the characteristic vector coefficient and the sample data of the full-sea deep historical sound velocity profile.
In a possible implementation manner, the obtaining manner of the first feature vector includes:
determining an average sound velocity profile of the truncated data of the historical sound velocity profile sample according to the truncated data of the historical sound velocity profile sample, and constructing a first residual error matrix;
calculating an autocovariance matrix of the first residual error matrix, and solving an eigenvalue of the autocovariance matrix and a corresponding eigenvector;
and arranging the obtained eigenvalues from big to small, and selecting the eigenvector corresponding to the eigenvalue of the preset order arranged at the front as the first eigenvector.
In a possible implementation manner, the performing matching field processing on the sound velocity profile to be extended by using the first eigenvector to obtain an eigenvector coefficient of the sound velocity profile to be extended corresponding to the first eigenvector includes:
determining a residual vector between the sound velocity profile to be extended and an average sound velocity profile of sample truncation data of the historical sound velocity profile;
and determining a characteristic vector coefficient of the sound velocity profile to be extended corresponding to the first characteristic vector according to the residual vector and the first characteristic vector.
In one possible embodiment, the feature vector coefficients are products of transposes of the first feature vector and a residual vector.
In a possible implementation manner, the constructing a extended full-sea-depth acoustic velocity profile according to the eigenvector coefficients and sample data of the full-sea-depth historical acoustic velocity profile includes:
performing empirical orthogonal function decomposition on the sample data of the full-sea-depth historical sound velocity profile to obtain an average sound velocity profile of the full-sea-depth historical sound velocity profile and a second feature vector of a preset order;
and constructing the extended full-sea-depth sound velocity profile according to the average sound velocity profile, the characteristic vector coefficient and the second characteristic vector of the full-sea-depth historical sound velocity profile.
In a possible implementation manner, the constructing a extended full-sea-depth sound velocity profile according to an average sound velocity profile, a feature vector coefficient, and a second feature vector of the full-sea-depth historical sound velocity profile includes:
and calculating the product of the feature vector coefficient and the second feature vector, and determining the sum of the product and the average sound velocity profile of the whole-sea-depth historical sound velocity profile as the extended whole-sea-depth sound velocity profile.
In a possible implementation manner, the obtaining manner of the second feature vector includes:
determining an average sound velocity profile of the full-sea-depth historical sound velocity profile according to the sample data of the full-sea-depth historical sound velocity profile, and constructing a second residual matrix;
calculating an autocovariance matrix of the second residual error matrix, and solving an eigenvalue of the autocovariance matrix and a corresponding eigenvector;
and arranging the obtained eigenvalues from big to small, and selecting the eigenvector corresponding to the eigenvalue of the preset order arranged at the front as a second eigenvector.
In a possible implementation manner, preprocessing an acquired full-sea-depth historical sound velocity profile data set to obtain full-sea-depth historical sound velocity profile sample data, specifically including: and performing depth interception on the full-sea-depth historical sound velocity profile data set according to the maximum depth of the actually measured sound velocity value of the full-sea-depth sound velocity profile to obtain full-sea-depth historical sound velocity profile sample data.
In a second aspect, an embodiment of the present invention provides a sound velocity profile continuation system based on orthogonal empirical function decomposition, including:
the data acquisition module is used for acquiring a sound velocity profile to be extended, the maximum depth of the actually measured sound velocity value of the sound velocity profile to be extended and historical sound velocity profile sample data of the whole sea depth;
the depth interception module is used for carrying out depth interception on the sample data of the full-sea-depth historical sound velocity profile according to the maximum depth of the actually-measured sound velocity value of the sound velocity profile to be extended to obtain the intercepted data of the historical sound velocity profile sample;
the processing module is used for performing empirical orthogonal function decomposition on the sample truncation data of the historical sound velocity profile to obtain a first feature vector of a preset order of the sample truncation data of the historical sound velocity profile, and performing matching field processing on the sound velocity profile to be extended by using the first feature vector to obtain a feature vector coefficient of the sound velocity profile to be extended corresponding to the first feature vector;
and the construction module is used for constructing the extended full-sea-depth sound velocity profile according to the characteristic vector coefficient and the sample data of the full-sea-depth historical sound velocity profile.
In a possible implementation, the building module is specifically configured to: performing empirical orthogonal function decomposition on the sample data of the full-sea-depth historical sound velocity profile to obtain an average sound velocity profile of the full-sea-depth historical sound velocity profile and a second feature vector of a preset order;
and constructing the extended full-sea-depth sound velocity profile according to the average sound velocity profile, the characteristic vector coefficient and the second characteristic vector of the full-sea-depth historical sound velocity profile.
The above one or more technical solutions have the following beneficial effects:
according to the method, empirical orthogonal function decomposition is carried out on the truncated data of a historical sound velocity profile sample to obtain a first feature vector of a preset order of the truncated data of the historical sound velocity profile sample, and matching field processing is carried out on the sound velocity profile to be extended by utilizing the first feature vector to obtain a feature vector coefficient of the sound velocity profile to be extended corresponding to the first feature vector; and constructing a prolonged full-sea deep sound velocity profile according to the characteristic vector coefficient and the sample data of the full-sea deep historical sound velocity profile. By the method, the main characteristics of the sound velocity profile sample to be extended in the cut-off depth of the region can be fully captured, the characteristic vector coefficient and the characteristic vector can be used for high-precision matching reconstruction of the sound velocity profile sample to be extended in the cut-off depth region (shallow sea) with obvious sound velocity change, meanwhile, the data of the depth part of the historical full-sea sound velocity profile sample exceeding the sound velocity profile sample to be extended is used as reference constraint, the sound velocity distribution of the extended part is guaranteed to be in accordance with the actual sound velocity distribution rule of the region, and more reliable reference samples are provided for sound velocity profile inversion.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification, illustrate exemplary embodiments of the invention and together with the description serve to explain the invention and not to limit the invention.
Fig. 1 is a schematic flow chart of a sound velocity profile continuation method based on orthogonal empirical function decomposition according to an embodiment of the present invention;
FIG. 2 is an overall framework diagram of a sound velocity profile continuation method based on orthogonal empirical function decomposition according to an embodiment of the present invention;
fig. 3 is a schematic diagram illustrating a comparison between the sound velocity profile continuation provided by the first embodiment of the present invention.
Detailed Description
The invention is further described with reference to the following figures and examples.
It is to be understood that the following detailed description is exemplary and is intended to provide further explanation of the invention as claimed. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
Example one
Referring to fig. 1, the present embodiment provides a sound velocity profile continuation method based on orthogonal empirical function decomposition, which specifically includes the following steps:
s101: and acquiring sample data of the sound velocity profile to be extended, the maximum depth of the measured sound velocity value of the sound velocity profile to be extended and the historical sound velocity profile of the whole sea depth.
In the specific implementation, the maximum depth of the actually measured sound velocity value of the sound velocity profile to be extended is shallow seawater sound velocity profile data, the historical sound velocity profile sample data of the full-sea depth comprises complete seawater sound velocity profile data obtained through multiple measurements, andisample data S of sound velocity profile of whole sea i Is represented as follows:
wherein the content of the first and second substances,sin order to be the value of the speed of sound,dfor the purpose of the depth value,d D is the maximum depth of the measured sound velocity value of the full-sea-depth sound velocity profile,s D is the sound velocity value of the full-sea deep sound velocity profile,Iis the number of samples.
S102: and performing depth interception on the full-sea-depth historical sound velocity profile sample data according to the maximum depth of the actually measured sound velocity value of the sound velocity profile to be extended to obtain historical sound velocity profile sample interception data.
In specific implementation, the sound velocity profile sample data S of the whole sea depth history i Performing partial depth 0-d truncation:d d the maximum depth of the actually measured sound velocity value of the sound velocity profile to be extended is reserved for the part, the 0-d depth interval part of each historical sound velocity profile sample is reserved, and the truncated data of the historical sound velocity profile sample
Expressed as:
。
s103: and performing empirical orthogonal function decomposition on the sample truncation data of the historical sound velocity profile to obtain a first feature vector of a preset order of the sample truncation data of the historical sound velocity profile, and performing matching field processing on the sound velocity profile to be extended by using the first feature vector to obtain a feature vector coefficient corresponding to the sound velocity profile to be extended in the first feature vector.
In specific implementation, the historical sound velocity profile sample truncation data is analyzed by using an empirical orthogonal function analysis method
And performing empirical orthogonal function decomposition to obtain a first feature vector of a preset order of the truncated data of the historical sound velocity profile sample, and performing matching field processing on the sound velocity profile to be extended by using the first feature vector to obtain a feature vector coefficient corresponding to the sound velocity profile to be extended in the first feature vector.
The obtaining mode of the first feature vector comprises the following steps: determining an average sound velocity profile of the sample truncation data of the historical sound velocity profile according to the sample truncation data of the historical sound velocity profile, and constructing a first residual matrix; calculating an autocovariance matrix of the first residual error matrix, and solving an eigenvalue of the autocovariance matrix and a corresponding eigenvector; and arranging the obtained eigenvalues from big to small, and selecting the eigenvector corresponding to the eigenvalue of the preset order arranged at the front as the first eigenvector.
Specifically, the average sound velocity profile of the truncated data of the historical sound velocity profile sample is as follows:
constructing a first residual matrix:
the autocovariance matrix of the first residual matrix is:
wherein the content of the first and second substances,C S,d for a matrix of (d + 1) × (d + 1) order, calculateC S,d Eigenvalue matrix of
And eigenvector matrix V d And satisfies the following conditions:
wherein, the first and the second end of the pipe are connected with each other,
for the diagonal matrix of order N:
n is the order of the characteristic value, V d Is a (d + 1) x Nth order matrix, each column of which is a feature vector, V d Can be expressed as
Characteristic value
And feature vector
And (7) corresponding.
Arranging the obtained eigenvalues from big to small, and selecting eigenvectors corresponding to the top n eigenvalues
Truncating a first eigenvector of a preset order of the data for the historical sound velocity profile sample, where the order value is typically an integer between 3-5.
The matching field processing is performed on the sound velocity profile to be extended by using the first feature vector to obtain a feature vector coefficient of the sound velocity profile to be extended corresponding to the first feature vector, and the method comprises the following steps: determining a residual vector between the sound velocity profile to be extended and an average sound velocity profile of sample truncation data of the historical sound velocity profile; and determining a characteristic vector coefficient of the sound velocity profile to be extended corresponding to the first characteristic vector according to the residual vector and the first characteristic vector. Optionally, the feature vector coefficient is a product of a transpose of the first feature vector and the residual vector.
In particular, the sound velocity profile to be extended may be expressed as
The sound velocity value vector of the sound velocity profile to be extended is expressed as
Wherein, in the step (A),trepresenting the sound velocity profile to be extended as test data, and the residual vector between the test data and the average sound velocity profile of truncation data of 0-d depth historical sound velocity profile sampleS X,t Comprises the following steps:
performing a matching field processing process on the sound velocity profile to be extended by utilizing the first n-order feature vector of the depth of 0-d, calculating the projection of the sound velocity profile to be extended on the first n-order feature vector, and acquiring the feature vector coefficient of the projection of the first n-order feature vector of the sound velocity profile to be extended on the depth of 0-dCf d ：
Here, the first and second liquid crystal display panels are,
is the first feature vector
The transposing of (1).
S104: and constructing a prolonged full-sea deep sound velocity profile according to the characteristic vector coefficient and the sample data of the full-sea deep historical sound velocity profile.
Specifically, in S104, the constructing a extended full-sea-depth acoustic velocity profile according to the feature vector coefficient and the sample data of the full-sea-depth historical acoustic velocity profile includes: performing empirical orthogonal function decomposition on the sample data of the full-sea-depth historical sound velocity profile to obtain an average sound velocity profile of the full-sea-depth historical sound velocity profile and a second eigenvector of a preset order; and constructing the extended full-sea-depth sound velocity profile according to the average sound velocity profile, the characteristic vector coefficient and the second characteristic vector of the full-sea-depth historical sound velocity profile.
As shown in fig. 2, preprocessing the acquired full-sea-depth historical sound velocity profile data set to obtain full-sea-depth historical sound velocity profile sample data, specifically including: according to the measured maximum depth of the sound velocity value of the full-sea-depth sound velocity profiled D And performing depth interception on the full-sea-depth historical sound velocity profile data set to obtain full-sea-depth historical sound velocity profile sample data, namely reserving 0-D full-sea-depth sound velocity profile sample S in the part i 。
The second feature vector is obtained in a manner that: determining an average sound velocity profile of the full-sea-depth historical sound velocity profile according to the sample data of the full-sea-depth historical sound velocity profile, and constructing a second residual matrix; calculating an autocovariance matrix of the second residual error matrix, and solving an eigenvalue of the autocovariance matrix and a corresponding eigenvector; and arranging the obtained eigenvalues from large to small, and selecting the eigenvector corresponding to the eigenvalue of the preset order arranged at the front as a second eigenvector.
In a specific implementation, the average acoustic velocity profile of the 0-D full-sea-depth acoustic velocity profile sample is:
constructing a second residual matrix of the 0-D full-sea sound velocity profile sample:
the autocovariance matrix of the second residual matrix is:
wherein, C S,D For a matrix of order (D + 1) × (D + 1), C is calculated S,D Characteristic value of
And a feature vector V D And satisfies the following conditions:
wherein the characteristic value
For diagonal matrices:
where M is the order of the eigenvalues, V D Is a matrix of order (D + 1) xM, each column of which is a feature vector, V D Can be expressed as
Characteristic value
And feature vector
And (7) corresponding.
Arranging the diagonal elements of the eigenvalue matrix of the auto-covariance matrix corresponding to the second residual matrix in the 0-D full sea depth through empirical orthogonal function decomposition from large to small, and selecting the first m-order eigenvector
Here, the preset order m = n takes an integer of 3-5, and in the present embodiment, m = n =3 may be taken.
The constructing of the extended full-sea-depth sound velocity profile according to the average sound velocity profile, the characteristic vector coefficient and the second characteristic vector of the full-sea-depth historical sound velocity profile comprises the following steps: and calculating the product of the feature vector coefficient and the second feature vector, and determining the sum of the product and the average sound velocity profile of the whole-sea-depth historical sound velocity profile as the extended whole-sea-depth sound velocity profile.
In the concrete implementation, the characteristic vector coefficient of the 0-d depth front n-order characteristic vector when the sound velocity profiles to be extended are matched is utilizedCf d 0-D average sound velocity profile at full sea depth
And first m order eigenvectors
Combined construction of extended full-sea deep acoustic velocity profile
The method specifically comprises the following steps:
extended acoustic velocity profile
Expressed as:
Example two
The embodiment of the invention also provides a sound velocity profile continuation system based on orthogonal empirical function decomposition, which comprises:
the data acquisition module is used for acquiring a sound velocity profile to be extended, the maximum depth of the actually measured sound velocity value of the sound velocity profile to be extended and historical sound velocity profile sample data of the whole sea depth;
the depth interception module is used for carrying out depth interception on the sample data of the full-sea-depth historical sound velocity profile according to the maximum depth of the actually-measured sound velocity value of the sound velocity profile to be extended to obtain the intercepted data of the historical sound velocity profile sample;
the processing module is used for performing empirical orthogonal function decomposition on the sample truncation data of the historical sound velocity profile to obtain a first feature vector of a preset order of the sample truncation data of the historical sound velocity profile, and performing matching field processing on the sound velocity profile to be extended by using the first feature vector to obtain a feature vector coefficient of the sound velocity profile to be extended corresponding to the first feature vector;
and the construction module is used for constructing the extended full-sea-depth sound velocity profile according to the characteristic vector coefficients and the sample data of the full-sea-depth historical sound velocity profile.
As an optional implementation manner, the building module is specifically configured to: performing empirical orthogonal function decomposition on the sample data of the full-sea-depth historical sound velocity profile to obtain an average sound velocity profile of the full-sea-depth historical sound velocity profile and a second eigenvector of a preset truncation order;
and constructing the extended full-sea-depth sound velocity profile according to the average sound velocity profile, the characteristic vector coefficient and the second characteristic vector of the full-sea-depth historical sound velocity profile.
The sound velocity profile continuation system based on orthogonal empirical function decomposition provided in this embodiment is used to implement the sound velocity profile continuation method based on orthogonal empirical function decomposition, so that a specific implementation manner in the sound velocity profile continuation system based on orthogonal empirical function decomposition may be found in the foregoing example part of the sound velocity profile continuation method based on orthogonal empirical function decomposition, and details are not described here.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A sound velocity profile continuation method based on orthogonal empirical function decomposition is characterized by comprising the following steps:
acquiring sample data of a sound velocity profile to be extended, the maximum depth of an actually measured sound velocity value of the sound velocity profile to be extended and historical sound velocity profiles of the whole sea depth;
according to the maximum depth of the actually measured sound velocity value of the sound velocity profile to be extended, carrying out depth interception on the full-sea-depth historical sound velocity profile sample data to obtain historical sound velocity profile sample interception data;
performing empirical orthogonal function decomposition on the truncated data of the historical sound velocity profile sample to obtain a first feature vector of a preset order of the truncated data of the historical sound velocity profile sample, and performing matching field processing on the sound velocity profile to be extended by using the first feature vector to obtain a feature vector coefficient of the sound velocity profile to be extended corresponding to the first feature vector;
and constructing a prolonged full-sea deep sound velocity profile according to the characteristic vector coefficient and the sample data of the full-sea deep historical sound velocity profile.
2. The method according to claim 1, wherein the obtaining of the first eigenvector comprises:
determining an average sound velocity profile of the sample truncation data of the historical sound velocity profile according to the sample truncation data of the historical sound velocity profile, and constructing a first residual matrix;
calculating an autocovariance matrix of the first residual error matrix, and solving an eigenvalue of the autocovariance matrix and a corresponding eigenvector;
and arranging the obtained eigenvalues from big to small, and selecting the eigenvector corresponding to the eigenvalue of the preset order arranged at the front as the first eigenvector.
3. The method for extending a sound velocity profile based on orthogonal empirical function decomposition of claim 1, wherein the performing matching field processing on the sound velocity profile to be extended by using the first eigenvector to obtain eigenvector coefficients of the sound velocity profile to be extended corresponding to the first eigenvector comprises:
determining a residual vector between the sound velocity profile to be extended and an average sound velocity profile of sample truncation data of the historical sound velocity profile;
and determining a characteristic vector coefficient of the sound velocity profile to be extended corresponding to the first characteristic vector according to the residual vector and the first characteristic vector.
4. The method of claim 3, wherein the eigenvector coefficients are the product of a transpose of the first eigenvector and a residual vector.
5. The method for extending the sound velocity profile based on the orthogonal empirical function decomposition of claim 1, wherein the constructing the extended full-sea-depth sound velocity profile according to the eigenvector coefficients and the sample data of the full-sea-depth historical sound velocity profile includes:
performing empirical orthogonal function decomposition on the sample data of the full-sea-depth historical sound velocity profile to obtain an average sound velocity profile of the full-sea-depth historical sound velocity profile and a second feature vector of a preset order;
and constructing the extended full-sea-depth sound velocity profile according to the average sound velocity profile, the characteristic vector coefficient and the second characteristic vector of the full-sea-depth historical sound velocity profile.
6. The method for extending the sound velocity profile based on the orthogonal empirical function decomposition of claim 5, wherein the constructing the extended whole-sea-depth sound velocity profile according to the average sound velocity profile, the eigenvector coefficients and the second eigenvector of the whole-sea-depth historical sound velocity profile comprises:
and calculating the product of the feature vector coefficient and the second feature vector, and determining the sum of the product and the average sound velocity profile of the whole-sea-depth historical sound velocity profile as the extended whole-sea-depth sound velocity profile.
7. The method according to claim 5, wherein the second eigenvector is obtained by:
determining an average sound velocity profile of the full-sea-depth historical sound velocity profile according to the sample data of the full-sea-depth historical sound velocity profile, and constructing a second residual matrix;
calculating an autocovariance matrix of the second residual error matrix, and solving an eigenvalue of the autocovariance matrix and a corresponding eigenvector;
and arranging the obtained eigenvalues from big to small, and selecting the eigenvector corresponding to the eigenvalue of the preset order arranged at the front as a second eigenvector.
8. The method for extending a sound velocity profile based on orthogonal empirical function decomposition of claim 1, wherein the step of preprocessing the acquired full-sea-depth historical sound velocity profile data set to obtain sample data of the full-sea-depth historical sound velocity profile specifically comprises the steps of: and performing depth interception on the full-sea-depth historical sound velocity profile data set according to the maximum depth of the actually measured sound velocity value of the full-sea-depth sound velocity profile to obtain full-sea-depth historical sound velocity profile sample data.
9. An orthogonal empirical function decomposition-based sound velocity profile continuation system, comprising:
the data acquisition module is used for acquiring a sound velocity profile to be extended, the maximum depth of the actually measured sound velocity value of the sound velocity profile to be extended and historical sound velocity profile sample data of the whole sea depth;
the depth interception module is used for carrying out depth interception on the sample data of the full-sea-depth historical sound velocity profile according to the maximum depth of the actually-measured sound velocity value of the sound velocity profile to be extended to obtain the intercepted data of the historical sound velocity profile sample;
the processing module is used for performing empirical orthogonal function decomposition on the sample truncation data of the historical sound velocity profile to obtain a first feature vector of a preset order of the sample truncation data of the historical sound velocity profile, and performing matching field processing on the sound velocity profile to be extended by using the first feature vector to obtain a feature vector coefficient of the sound velocity profile to be extended corresponding to the first feature vector;
and the construction module is used for constructing the extended full-sea-depth sound velocity profile according to the characteristic vector coefficient and the sample data of the full-sea-depth historical sound velocity profile.
10. The orthogonal empirical function decomposition-based sound velocity profile continuation system of claim 9, wherein the construction module is specifically configured to: performing empirical orthogonal function decomposition on the sample data of the full-sea-depth historical sound velocity profile to obtain an average sound velocity profile of the full-sea-depth historical sound velocity profile and a second feature vector of a preset order;
and constructing an extended full-sea-depth sound velocity profile according to the average sound velocity profile, the characteristic vector coefficient and the second characteristic vector of the full-sea-depth historical sound velocity profile.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211118337.3A CN115204317B (en) | 2022-09-15 | 2022-09-15 | Sound velocity profile continuation method and system based on orthogonal empirical function decomposition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211118337.3A CN115204317B (en) | 2022-09-15 | 2022-09-15 | Sound velocity profile continuation method and system based on orthogonal empirical function decomposition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115204317A CN115204317A (en) | 2022-10-18 |
| CN115204317B true CN115204317B (en) | 2022-11-25 |
Family
ID=83573525
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211118337.3A Active CN115204317B (en) | 2022-09-15 | 2022-09-15 | Sound velocity profile continuation method and system based on orthogonal empirical function decomposition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115204317B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1062543A (en) * | 1996-08-26 | 1998-03-06 | Oki Electric Ind Co Ltd | Development function generation method and apparatus for ocean acoustic tomography |
| CN111307266A (en) * | 2020-02-21 | 2020-06-19 | 山东大学 | Sound velocity obtaining method and global ocean sound velocity field construction method based on same |
| CN112964231A (en) * | 2021-02-03 | 2021-06-15 | 广东海洋大学 | Method for obtaining depth of ocean mixing layer based on sound velocity disturbance modal matching |
| CN113051260A (en) * | 2021-04-27 | 2021-06-29 | 中国人民解放军国防科技大学 | High-resolution sound velocity profile data compression method based on empirical orthogonal function decomposition |
| CN113218493A (en) * | 2021-04-08 | 2021-08-06 | 中国人民解放军国防科技大学 | Sound velocity profile inversion method based on empirical orthogonal function method |
-
2022
- 2022-09-15 CN CN202211118337.3A patent/CN115204317B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1062543A (en) * | 1996-08-26 | 1998-03-06 | Oki Electric Ind Co Ltd | Development function generation method and apparatus for ocean acoustic tomography |
| CN111307266A (en) * | 2020-02-21 | 2020-06-19 | 山东大学 | Sound velocity obtaining method and global ocean sound velocity field construction method based on same |
| CN112964231A (en) * | 2021-02-03 | 2021-06-15 | 广东海洋大学 | Method for obtaining depth of ocean mixing layer based on sound velocity disturbance modal matching |
| CN113218493A (en) * | 2021-04-08 | 2021-08-06 | 中国人民解放军国防科技大学 | Sound velocity profile inversion method based on empirical orthogonal function method |
| CN113051260A (en) * | 2021-04-27 | 2021-06-29 | 中国人民解放军国防科技大学 | High-resolution sound velocity profile data compression method based on empirical orthogonal function decomposition |
Non-Patent Citations (2)
| Title |
|---|
| Precise Positioning of Underwater Static Objects without Sound Speed Profile;Fanlin Yang等;《Marine Geodesy》;20111231;全文 * |
| 水下声学定位关键技术问题研究;高茹茹;《万方数据》;20190118;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115204317A (en) | 2022-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101272168B (en) | Signal sources estimation method and its DOA estimation method | |
| CN109410917A (en) | Voice data classification method based on modified capsule network | |
| US8270632B2 (en) | Sound source localization system and method | |
| CN108828547A (en) | The high method of the low Elevation of metre wave radar based on deep neural network | |
| EP2232407B9 (en) | Method for separating mixed signals into a plurality of component signals | |
| CN113486574B (en) | Sound velocity profile completion method and device based on historical data and machine learning | |
| CN108375763A (en) | A kind of frequency dividing localization method applied to more sound source environment | |
| CN114022812B (en) | DeepSort water surface floater multi-target tracking method based on lightweight SSD | |
| CN110048416A (en) | The low-frequency oscillation modal identification method of S-G filtering and adaptive M P algorithm | |
| CN112285647B (en) | Signal azimuth high-resolution estimation method based on sparse representation and reconstruction | |
| CN114972339B (en) | Data enhancement system for bulldozer structural member production abnormity detection | |
| CN113109759A (en) | Underwater sound array signal direction-of-arrival estimation method based on wavelet transformation and convolutional neural network | |
| CN111785286A (en) | Home CNN classification and feature matching combined voiceprint recognition method | |
| CN111812581A (en) | Spherical array sound source direction of arrival estimation method based on atomic norm | |
| CN115204317B (en) | Sound velocity profile continuation method and system based on orthogonal empirical function decomposition | |
| CN106483193B (en) | A kind of wave based on High-order Cumulant reaches method for quick estimating | |
| CN110716203A (en) | Time-frequency analysis and tracking method of passive sonar target | |
| CN103837858A (en) | Far field direction of arrival estimation method applied to plane array and system thereof | |
| CN104407319A (en) | Method and system for finding direction of target source of array signal | |
| CN108107392B (en) | Multi-line TRL calibration method and terminal equipment | |
| CN112014791B (en) | Near-field source positioning method of array PCA-BP algorithm with array error | |
| CN105572629A (en) | Two-dimensional direction measuring method of low operation complexity and applicable to any array structure | |
| CN111175709B (en) | Error suppression-based large-range meteorological radar-oriented jigsaw method | |
| CN117169143A (en) | Full spectrum water quality detection method and system based on neural network | |
| CN112904270B (en) | Direction-of-arrival estimation method based on fitting model under array model error |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |