66 results about "Singular value decomposition method" patented technology

The invention discloses a water vapor chromatography modeling method and device based on GNSS-MET, and belongs to the technical field of atmospheric water vapor information inversion. According to the method, satellite signals and meteorological parameters are obtained through a GNSS-MET receiver, and the oblique path water vapor content (SWV) is calculated to serve as chromatography input; densely and uniformly distributed observation stations are selected to construct a chromatography area, and voxel grids are divided; a horizontal constraint is constructed by using Gaussian distance weighting, a vertical constraint is constructed in combination with a water vapor vertical index distribution characteristic, and ERA5 reanalysis data is introduced as a prior constraint, so that the limitation of a traditional sounding data constraint is solved; and resolving the tomographic equation through a singular value decomposition (SVD) method to obtain three-dimensional water vapor density distribution. According to the method, the high temporal-spatial resolution characteristic of ERA5 data is utilized, the precision and reliability of water vapor chromatography are remarkably improved, and the method is suitable for the fields of extreme weather prediction, meteorological monitoring and the like.

The invention provides an image restoration method based on adaptive weighted tensor completion, and relates to the technical field of image processing and application, and the method comprises the steps: obtaining to-be-restored image data, and carrying out the tensor of the to-be-restored image data, and obtaining input tensor data; constructing a tensor completion model based on an adaptive weighted tensor nuclear norm; wherein a weight matrix in the tensor completion model can be adaptively updated along with input tensor data; and based on an alternating direction multiplier method or an approximate singular value decomposition method based on tensor QR decomposition, solving the tensor completion model, and outputting restored tensor data to realize image restoration. According to the scheme, the image restoration quality can be improved.

The invention provides a Prony algorithm-based cable segment wave velocity acquisition method, device and system, and a medium, and the method comprises the steps: testing a to-be-tested cable, and obtaining a cable signal transfer function; on the basis of the cable signal transfer function, in combination with cable joint distribution, constructing a cable signal approximation function based on a Prony method; constructing a linear prediction model of the cable signal based on the cable signal approximation function; performing denoising processing on the linear prediction model based on a singular value decomposition method to obtain a denoised linear prediction model; solving the denoised linear prediction model to obtain an attenuation coefficient; and calculating the segmented wave velocity of the cable based on the attenuation coefficient. According to the method, the Prony estimation method is combined with the singular value decomposition noise reduction algorithm, high-precision extraction of attenuation constants and segmented wave velocity decoupling are achieved, and then the electrical distance positioning precision is improved.

The invention discloses a truss type load-bearing structure time domain dynamics topological optimization method and system based on an intrinsic orthogonal decomposition method, and relates to the technical field of truss type load-bearing structure topological optimization. The method is provided for solving the problems of high calculation power and the like caused by large-scale space-time response iteration of original transient problems and intensive calculation of dual problem sensitivity analysis in the prior art. According to the method, a local geometric control strategy for the cross interference problem of rod piece units is introduced, and a transient dynamics full-order topological optimization model of the truss structure is established through an equal geometric stiffness diffusion method. Based on an intrinsic orthogonal decomposition method, a transient dynamics reduced-order model of a truss structure is constructed, a reduced-order model base is updated in real time through an incremental singular value decomposition method, and then efficient solving of a transient dynamics equation and an adjoint equation in sensitivity analysis is achieved. And updating an iterative design variable by adopting a penalty function method, and seeking the optimal layout design of the truss structure. The invention provides a time-domain dynamics topological optimization method of the truss type load-bearing structure based on an intrinsic orthogonal decomposition method, and the method can be used for efficiently solving the dynamics layout optimization problem of the truss structure under the action of a transient load.

An electric vehicle charging station recommendation method and system, applied to the related field of charging station recommendation. The method comprises: acquiring user charging features and influence degree values of user charging evaluation features on user charging preferences, wherein the user charging evaluation features are constructed on the basis of the user's evaluation information on charging stations (S1); calculating the user's comprehensive scores for different charging stations on the basis of the influence degree values, and constructing a user preference matrix by using the user's comprehensive scores for different charging stations as elements (S2); and decomposing the user preference matrix by means of a singular value decomposition method, acquiring potential matching relationships between the user and charging station features, and recommending a charging station to the user on the basis of the potential matching relationships (S3). The method solves the problem in the prior art that charging stations cannot be personalizedly recommended for different users, thereby improving user experience.

The invention discloses a heterogeneous air quality data fusion method based on sparse matrix decomposition. The method comprises the following steps: S1, collecting air quality data; s2, preprocessing the collected air quality data, and constructing a preliminary fusion matrix; s3, marking missing data in the preliminary fusion matrix, and constructing a sparse matrix; s4, carrying out matrix decomposition by adopting a mode of combining an improved singular value decomposition method and a matrix completion method based on low-rank representation, filling missing data, and carrying out feature weighted optimization on a matrix decomposition result by utilizing an attention mechanism; and S5, error measurement and consistency check are carried out. According to the method, sparse matrix decomposition and low-rank completion technologies are adopted, the fusion process of the air quality data is optimized in combination with an attention mechanism, and high data precision, high consistency and improvement of prediction accuracy are achieved.

The invention relates to a singular value decomposition-based personalized federated learning method, which is used for solving the problem of uncertainty of a system caused by equipment isomerism and the problem that a client falls behind due to extra burden caused by a storage and transmission mechanism in the prior art. Firstly, a customer selection strategy based on resource awareness is determined, and a server side establishes a quantitative scoring system by comprehensively analyzing computing resources and communication conditions of clients. And screening a plurality of clients with optimal comprehensive performance from all online clients to participate in local training. During local training of the client, historical weight parameters are reserved by adopting self-knowledge distillation (SKDSVD) combined with a singular value decomposition method, residual errors after decomposition are personalized model parameters and are stored in the local client, and universal parameters are transmitted to the server for aggregation. The SKDSVD algorithm reduces the training time delay and the fall-behind rate of the local client while optimizing the global model of the server to obtain the optimal convergence precision.

The invention relates to a bolt looseness diagnosis method, system and device and a medium. The method comprises the following steps: acquiring an original sound signal for knocking a bolt; decomposing the original sound signal by using a CEEMDAN method to obtain a plurality of IMF components; calculating a plurality of preset types of quantitative indexes for each IMF component, performing dynamic weighting, and screening out a plurality of effective IMF components from the plurality of IMF components based on a weighting result of the IMF components as effective components; determining a dynamic function window and a frequency spectrum gravity center according to the plurality of effective components, decomposing the effective components in combination with a singular value decomposition method, and performing screening and reconstruction to obtain de-noised signals; performing feature extraction on the de-noised signal to obtain a diagnosis feature; and inputting the diagnosis characteristics into a bolt looseness diagnosis model to obtain a bolt looseness diagnosis result. According to the method, the noise reduction effect of the sound signals collected under the working condition can be effectively improved, the obvious bolt looseness diagnosis accuracy and superiority are achieved, and the bolt looseness diagnosis precision and efficiency can be improved.

The invention discloses a mechanical arm 3D printing normal calibration method based on singular value decomposition. The mechanical arm 3D printing normal calibration method comprises the steps that S1, a dial indicator is installed at the tail end of a mechanical arm; s2, respectively establishing coordinate systems by taking the center point of the flange at the tail end of the mechanical arm, the actual contact point of the dial indicator and the non-deformation tail end point of the dial indicator as original points; s3, taking the central point of the flange at the tail end of the mechanical arm as a movement track of the point in any set space curve, and sampling to obtain a space point set of sampling points; s4, fitting a plane in the space by using a singular value decomposition method, so that the sum of squares of distances from all sampling points in the space point set to the plane is minimum, and further obtaining a three-dimensional space plane equation and a normal vector of the fitted plane; and S5, coordinate system posture correction is conducted on the mechanical arm tail end flange center point coordinate system, so that the axis of the center of the mechanical arm tail end flange is parallel to the normal vector, and the parallelism between the printing track surface and the forming platform is improved. The forming quality of the first layer of 3D printing can be effectively improved.

The invention belongs to the technical field of health state assessment, and discloses a health state assessment method and device, electronic equipment and a storage medium, and the method comprises the steps: obtaining the physiological feature data of a user in a preset period, constructing the physiological feature data into tensor data containing an individual dimension, a feature dimension and a time dimension, complementing missing data in the tensor data by using a frequency domain transformation method and a singular value decomposition method to obtain complemented physiological feature data, inputting the complemented physiological feature data into a preset health state evaluation model, and calculating to obtain a health state evaluation result of the user in a preset period; by inputting the complemented physiological feature data obtained by complementing with the frequency domain transformation method and the singular value decomposition method into the preset health state assessment model, the health state assessment result of the user in the preset period is calculated, and the accuracy of health state assessment is improved.

The invention provides a spectral data cube reconstruction method, device and equipment based on compressed physical prior, a storage medium and a program product, and relates to the technical field of spectral reconstruction, and the method comprises the steps: obtaining a two-dimensional measurement image formed after a to-be-observed target is coded and modulated by a spectral imaging chip; obtaining an original transmission spectrum matrix corresponding to each pixel in the spectral imaging chip; performing principal component analysis on the original transmission spectrum matrix in a spectrum dimension, and extracting a principal component feature vector through a singular value decomposition method; projecting the original transmission spectrum matrix into an orthogonal subspace formed by the principal component feature vectors to obtain a compressed physical prior vector; and inputting the two-dimensional measurement image and the compressed physical prior vector into a deep expansion neural network, and outputting a high-dimensional spectral data cube of the to-be-observed target. The reconstruction efficiency of the spectral data cube can be remarkably improved.

The invention provides an optical guiding method for recycling a deep-sea unmanned underwater vehicle, and belongs to the technical field of underwater vehicle recycling. Four blue-green LED guiding light sources are non-uniformly and annularly distributed and installed on the front end face of a recycling cage, the three-dimensional coordinates of the guiding light sources are calibrated in advance, and guiding light source images are collected in real time through a binocular camera; an adaptive exposure control technology is adopted to adjust image quality, statistical information calculation is performed on an identified light source contour, an effective guide light source target is determined through convex constraint configuration judgment and ellipse fitting screening, four guide light sources are constructed as a traveling salesman problem, and a nearest neighbor heuristic algorithm is adopted to complete number matching. A rotation matrix and a translation vector are solved through a singular value decomposition method according to the pre-measured coordinates and the real-time calculated coordinates, and collaborative optimization is carried out at the same time, so that the technical problems of insufficient optical guidance positioning precision and poor real-time performance in the recovery process of the unmanned underwater vehicle in the deep sea environment are solved.

This invention discloses a high-precision FMCW laser ranging method based on the ESPRIT algorithm, belonging to the field of laser ranging technology. First, the laser beam is split into two beams by an optical beam splitter. One beam serves as the local oscillator beam, and the other as the signal beam. The signal beam hits the target object and returns, carrying the distance information of the target. The local oscillator beam and the signal beam are coupled in an optical coupler to obtain a difference frequency signal, which is then acquired by a photodetector and a digital acquisition card. A Hankel matrix is ​​constructed from the difference frequency signal. The Hankel matrix is ​​decomposed using singular value decomposition (SVD) to extract the signal subspace. Based on the decomposition of the signal subspace, an invertible matrix rotation operator is obtained and optimized to obtain the signal frequency, thereby calculating the distance to the target object. This invention is not only simple to implement and easy to operate, but also highly practical and suitable for widespread use.

The application belongs to the technical field of three-dimensional point cloud registration, and specifically provides a point cloud registration method based on second-order spatial compatibility measurement, comprising the following steps: extracting corner points of source point cloud and target point cloud as key points by using Harris-3D algorithm; describing and preliminarily matching the extracted key points by using FPFH algorithm to obtain matching pairs; calculating a second-order spatial compatibility matrix between the matching pairs, and removing false matching pairs to obtain correct matching pairs; calculating a coarse registration result of two point clouds by using singular value decomposition method in combination with the correct matching pairs to obtain an initial pose; and calculating a fine registration result of point clouds according to the initial pose of two point clouds in combination with K-dimensional tree accelerated ICP algorithm to finally complete the point cloud registration work. The application effectively suppresses the influence of outliers on registration accuracy, and improves the accuracy and speed of point cloud registration.

The present application relates to aero-engine test bench structure health monitoring and performance evaluation technology, provide a kind of test bench performance attenuation evaluation method, system and equipment based on defect detection, the method comprises: based on at least two types of non-destructive testing method to construct the finite element numerical model of test bench;By singular value decomposition method to finite element numerical model order reduction and solution obtain order reduction model;The predicted structure displacement of test bench is obtained in order reduction model, the actual structure displacement of test bench under test load is collected, and the structure performance attenuation evaluation of test bench is carried out by predicted structure displacement and actual structure displacement.This application method can carry out online structure health identification and performance early warning to test bench by the degradation model constructed, with good accuracy, scalability and engineering practical value, can be widely applied in the state evaluation and fault prediction of aero-engine test bench.

The invention relates to an aero-engine test bed structure health monitoring and performance evaluation technology, and provides a test bed performance attenuation evaluation method, system and equipment based on defect detection, and the method comprises the steps: building a finite element numerical model of a test bed based on at least two types of non-destructive detection methods; performing order reduction and solution on the finite element numerical model through a singular value decomposition method to obtain an order reduction model; and acquiring the predicted structural displacement of the test bench in the reduced-order model, acquiring the actual structural displacement of the test bench under the test load, and evaluating the structural performance attenuation of the test bench through the predicted structural displacement and the actual structural displacement. According to the method, on-line structure health identification and performance early warning can be carried out on the test bed through the constructed degradation model in combination with real-time monitoring data, and the method has good accuracy, expandability and engineering practical value and can be widely applied to state evaluation and fault prediction of the aero-engine test bed.

The invention discloses a multi-scale space target relative pose estimation method and system using parallel double encoders, and belongs to the technical field of space target pose estimation. The invention aims to solve the problems caused by distance change, illumination difference and insufficient multi-scale feature expression of a space target in a dynamic environment. According to the technical key points, a parallel double-encoder architecture is adopted, and local and global features of a space target are extracted; constructing a feature fusion module to fuse features extracted by the CNN and the Transform; multi-scale features are fused through FPN, and angular points are predicted; training and optimizing the network; and solving the pose by using a random consistency perspective n-point method. A plurality of non-coplanar control points are selected, other points are expressed as a linear combination of the control points, a linear equation set is constructed, and a rotation matrix and a translation vector between a target and a camera coordinate system are solved by using a singular value decomposition method. The deep learning network algorithm provides an efficient and accurate solution for space target pose estimation, the algorithm can deeply mine rich information of image data by adopting an encoder-decoder structure, and finally high-precision calculation of the six-degree-of-freedom pose of the target is achieved.

The utility model provides a kind of power station parameter matrix analysis method, belongs to power station fault diagnosis technical field, solve how to design a kind of power station parameter matrix analysis method, solve the numerous power station parameters, data effective information density is low, lead to the problem that fault cause analysis is difficult, and the massive data stream generated occupies a large amount of storage space, by power station parameter data matrixization representation, constructs parameter data matrix, by matrix operation, the rank of matrix is as quantization index, effectively distinguishes the information density of parameter, the correlation between different parameters can be analyzed by the multi-parameter matrix constructed, not limited by parameter number and time sequence length, and by matrix similarity operation, judge the similarity of unit state, provide quick judgment basis for similar fault in later period;Dimensionality reduction is carried out to operation data matrix using matrix singular value decomposition method, retain the part of large singular value proportion, can extract fault feature, realize compressed storage, save storage space.

The present invention provides a bearing degradation trend prediction method based on multi-resolution feature extraction and Bi-LSTM network. It aims to solve the problem that the bearing vibration signal contains a large amount of random noise in complex working conditions and environments, which affects the bearing degradation trend prediction accuracy. It fully considers the impact of historical vibration state on current vibration. The main steps are as follows: applying the full life cycle vibration signal of the accelerated bearing degradation test, performing time domain analysis on it to extract multi-angle eigenvalues ​​to form a feature set, screening sensitive features according to the trend and correlation of the features, using the multi-resolution singular value decomposition method to decompose and reconstruct the sensitive features, inputting the normalized sensitive features into the constructed Bi-LSTM network for training, and predicting the bearing degradation trend. Its purpose is to provide a prediction method that fully considers the impact of the historical vibration state of the bearing on the current moment, eliminates random interference during operation, retains important fault information in the bearing degradation process, and improves the accuracy of bearing degradation trend prediction.

The application relates to the technical field of radars, and provides a multi-target parameter dimension reduction estimation method and device for a frequency control array radar, which comprises the following steps: determining a receiving signal matrix of the frequency control array radar based on a transmitting steering vector, a receiving steering vector of each target and a transmitting signal of the frequency control array radar at each snapshot moment; solving the receiving signal matrix by using a high-order singular value decomposition method to obtain a noise subspace matrix of the receiving signal matrix; constructing a multi-target parameter dimension reduction estimation model based on the noise subspace matrix and the transmitting steering vector and the receiving steering vector of each target; solving the multi-target parameter dimension reduction estimation model by using a polynomial root-finding method to estimate an angle estimation value of each target; and estimating a distance estimation value of each target corresponding to the angle estimation value of each target based on the angle estimation value of each target. The application not only reduces the calculation complexity under the condition of few snapshots, but also improves the precision of parameter estimation.

This application discloses an adaptive singular value decomposition method for underwater acoustic signal denoising and its application. The method includes: Step 1) acquiring the noisy underwater acoustic signal to be analyzed; Step 2) constructing an m×n dimensional Hankel matrix based on the underwater acoustic signal, performing singular value decomposition on the Hankel matrix to obtain a singular value diagonal matrix; Step 3) obtaining the curvature of each singular value in the singular value diagonal matrix, determining a threshold based on the curvature, and acquiring the singular values ​​within the threshold as valid singular values; Step 4) retaining the valid singular values ​​in the singular value diagonal matrix, setting the other singular values ​​to zero, and updating the singular value diagonal matrix; Step 5) reconstructing the denoised Hankel matrix using the updated singular value diagonal matrix, selecting all elements of the first row of the denoised Hankel matrix and the m-1 elements from the nth column of the second row to the nth column of the mth row for restoration, to obtain the denoised underwater acoustic signal.

The invention relates to the technical field of traffic flow analysis, in particular to a road network traffic flow analysis method based on robust principal component tracking. The method comprises the steps of collecting traffic information of a target area road network; performing compensation processing and protocol processing on the traffic information to obtain an original traffic matrix of the traffic information; performing singular value decomposition on the original traffic matrix to obtain an original noise matrix; obtaining a random matrix constraint parameter of the relaxation principal component tracking method according to the original noise matrix; decomposing the original traffic matrix by using a relaxation principal component tracking method to obtain a low-rank matrix of traffic information; and decomposing the low-rank matrix by adopting a singular value decomposition method to obtain a time-varying feature component and a spatial feature component, and carrying out traffic flow analysis after clustering the spatial feature component. According to the method, sudden change traffic flow components contained in traffic information and noise components caused by various reasons can be effectively eliminated, and the accuracy of traffic flow analysis is improved.

The invention belongs to the technical field of rainfall forecasting, and particularly provides a Yangtze River upstream flood season inter-monthly rainfall forecasting method based on SVD, and the method comprises the steps: data preprocessing: processing historical data in an inter-adult increment form according to an original time sequence; determining a critical sea temperature period influencing the upstream rainfall of the Yangtze River; performing principal component decomposition and reconstruction on the Yangtze river upstream rainfall data corresponding to the to-be-measured month and the screened-out sea temperature data of the key sea temperature period to obtain an average pitch percentage and a global sea temperature average pitch corresponding to the reconstructed Yangtze river upstream to-be-measured month; establishing a rainfall prediction model based on a singular value decomposition method; constructing a system error correction model based on a singular value decomposition method; and predicting in real time. According to the method, the current Yangtze river upstream rainfall objective forecasting level is improved by establishing the Yangtze river upstream flood season monthly rainfall machine learning forecasting model.

The application discloses a method for correcting time base jitter of a high-speed sampling oscilloscope, and the method comprises the following steps: constructing a function model of an oscilloscope waveform, processing an equation group of data by using a matrix, introducing a correction vector and a matrix, expressing a constraint matrix and solving; wherein the step of expressing the constraint matrix and solving comprises the following steps: calculating jitter errors generated by the oscilloscope by using an orthogonal distance regression algorithm; the calculation and solving step is to adopt different methods according to the number of solutions of the matrix, when the matrix is a single solution, a singular value decomposition method is adopted; when the matrix is a multiple solution, a gradient descent method is adopted. The technical scheme disclosed by the application corrects the time base jitter problem of the high-speed high-bandwidth sampling oscilloscope, and makes the measurement result more accurate.

A canopy morphology 3D feature analysis method based on lightweight drones and point cloud singular value decomposition belongs to the field of plant morphology analysis. The method is as follows: a drone-mounted RGB camera is combined with the SfM‑MVS algorithm to perform high-precision three-dimensional reconstruction of the soybean canopy structure. At the same time, a label-based method introduces the singular value decomposition (SVD) method for the point cloud coordinate matrix, and 9 canopy morphology 3D features at multiple scales are established as a new interpretable feature at the crop plot level. The results show that the global canopy 3D features are significantly correlated with soybean yield, lodging, canopy height, and vegetation index (p <= 0.05). 3D feature extraction based on local range can extract more detailed structural information of the canopy surface, which helps to characterize the differences in more growth stages of soybeans.

This application relates to a method for evaluating the stability of roadway roof based on the response characteristics of drilling parameters. The method includes: selecting any point on the roof within the roadway and continuously collecting drilling parameters during the drilling process using sensors; using singular value decomposition to denoise the drilling parameters to obtain a denoised drilling dataset, analyzing the energy consumption per unit volume of rock during drilling, and obtaining the uniaxial compressive strength of the rock based on the energy consumption per unit volume of rock; analyzing the thickness, number of fractures, and width of intact rock strata in the roadway roof based on drilling energy consumption or borehole images at different drilling depths; analyzing the cumulative percentage of intact rock strata in the roadway roof based on the thickness of intact rock strata; and determining the stability evaluation result of the roadway roof based on the uniaxial compressive strength of the rock, the thickness of intact rock strata, the number of fractures, the width of fractures, and the cumulative percentage of intact rock strata. This allows for a simple and rapid effective rating and evaluation of roof stability.

The invention discloses a bamboo strip character denoising method based on a fractional order diffusion equation inverse problem POD algorithm, and belongs to the field of image processing and inverse problem calculation, and the method comprises the steps: carrying out the preprocessing of a bamboo strip image, and obtaining a terminal observation image; based on the set diffusion time and the fractional order, a finite difference method is adopted to solve a positive problem of a time fractional order diffusion equation, and a snapshot data set is generated; based on the snapshot data set, a singular value decomposition method is adopted to extract a main mode, and a reduced-order model is established; based on the terminal observation image and the reduced-order model, constructing a regularization least square optimization problem and adopting a gradient iteration algorithm for solving, and performing inversion to obtain a clear image at an initial moment; and processing the clear image at the initial moment by adopting a threshold segmentation method, and outputting a finally recovered high-definition bamboo strip image. The invention provides an efficient and stable treatment means with clear physical significance for digital protection of the fragile cultural relics such as the bamboo strips.

The application provides a bearing wear degree diagnosis method, device, medium and product, the method comprises the following steps: decomposing an acceleration vibration signal based on a variational mode decomposition method to obtain K sub-signals of different vibration modes, the acceleration vibration signal is an acceleration vibration signal of a target bearing; extracting feature information of each sub-signal of the different vibration modes based on a singular value decomposition method to obtain a feature component corresponding to each sub-signal of the different vibration modes; and sending the feature component to a gated recurrent unit to diagnose the wear degree of the target bearing through the gated recurrent unit. Therefore, the application improves the precision of bearing wear degree diagnosis.

The invention discloses a particle number state prediction method and device, equipment and a medium. The method comprises the steps of performing data prediction on current state information based on a target system model corresponding to a target particle counter and target noise parameter information to determine a basic state prediction value and a basic measurement prediction value corresponding to a target sample; performing covariance calculation on the basic state prediction value based on a preset singular value decomposition method to determine a prediction covariance error, and performing residual calculation on the basic measurement prediction value based on a preset residual rule to determine a measurement prediction error; and performing gain calculation on the prediction covariance error based on a preset adjustment fusion strategy to determine a target gain, and performing state updating on the basic state prediction value based on the target gain and the measurement prediction error to obtain a target state prediction value corresponding to the target sample. According to the technical scheme, the influence caused by multiplicative noise can be effectively reduced, and the accuracy of a particle number state result is improved.

The application provides a method and device for determining the deformation of an open-pit mine dump based on time-series InSAR technology, which aims to effectively fuse continuous multi-period observations, eliminate the interference of surface elevation on SAR interference phase, and improve the determination accuracy of the deformation of the open-pit mine dump. The method comprises the following steps: obtaining SAR satellite images and a digital surface model of a target area in a preset period; constructing a plurality of multi-master image differential interference queues based on the SAR satellite images; calculating differential interference phases based on the plurality of multi-master image differential interference queues and the corresponding digital surface model, and obtaining unwrapped phases after atmospheric filtering processing by combining a singular value decomposition method; and calculating the deformation of the target area based on a sequential adjustment method, wherein after constructing an error estimation matrix containing elevation errors and unwrapped phases, the sequential adjustment method is used to obtain a velocity matrix corresponding to the time-series phase by using the error estimation matrix, and the deformation of the target area is calculated based on the velocity matrix.