51 results about "Simultaneous equations" patented technology

The invention discloses a subway tunnel fire behavior spreading prediction method, which comprises the steps of 1, fire source type digital modeling, 2, fire scene initialization, 3, dynamic fire source loading and ignition triggering, 4, fire behavior spreading multi-physics field coupling calculation, 5, ventilation system dynamic response, 6, visual risk thermodynamic diagram generation, and 7, fire behavior spreading prediction. Step 7, carrying out real-time early warning on the safety margin; the method has the beneficial effects that a total factor modeling capability is realized, multi-dimensional parameter input from fire source characteristics to a tunnel environment is supported, a scene range covering a typical subway fire disaster is expanded, solving is realized through a multi-physics field simultaneous equation, a fire spreading error is reduced, smoke diffusion simulation precision is improved, and a dynamic ventilation adjustment and safety early warning mechanism is realized. And more buffer time is provided for emergency response.

The invention provides a mountain slope transmission system catenary determining method and device, and the method comprises the steps: obtaining the arrangement parameters of a transmission system on a mountain slope, and the arrangement parameters at least comprise the horizontal distance between two fixed points of the transmission system, the vertical height difference, and the total length of a rope of the transmission system; based on the arrangement parameters, a simultaneous equation with the catenary opening coefficient as a core unknown quantity is established, and the simultaneous equation comprises the arc length constraint of the catenary and the height difference constraint of the two end points at the same time; solving the simultaneous equation to obtain a numerical solution of the catenary opening coefficient; and determining a catenary shape and a core engineering index of the transmission system on the mountain slope under the action of gravity according to the solved catenary opening coefficient. According to the method and the device, the engineering design precision is improved, and material redundancy and adjustment cost caused by shape uncertainty are saved.

The invention discloses a contact resistance-based method for synchronously measuring normal force and shear force of an interface, which comprises the following steps of: inspired by a yield criterion of material damage, establishing a relationship between normal stress and shear stress of the contact interface, thereby obtaining a relationship between a real contact area of the contact interface and the normal force and the shear force; obtaining the relationship between the contact resistance of the contact interface and the real contact area, and establishing the relationship between the contact resistance of the interface and the normal force and the shear force; according to the geometric deformation relation of the micro-convex bodies on the contact interface before and after the shearing force acts, the relation between the real contact area of the contact interface and the geometric deformation parameters of the micro-convex bodies is obtained, and the relation among the normal force of the contact interface, the shearing force and the geometric deformation parameters of the micro-convex bodies is obtained; and establishing a simultaneous equation set, measuring the contact resistance, and substituting the contact resistance into the equation set to calculate the normal force and the shear force. According to the invention, the normal force and shear force of the interface can be indirectly measured by measuring the contact resistance of the contact interface.

The invention discloses a power distribution network fault accurate positioning method based on negative-sequence component characteristics and two-end quantities, and the method comprises the steps: collecting fault voltage and current information, and extracting voltage and current fundamental frequency components; performing phase-mode transformation on the fundamental frequency phasors of the three-phase voltage and the three-phase current to obtain negative-sequence components of the voltage and the current; the amplitudes of the current negative-sequence components at the outlets of all the lines are extracted and compared, and the line with the maximum current negative-sequence component amplitude at the outlet is judged as a fault line; establishing head-end quantity and tail-end quantity expressions of the negative-sequence voltage amplitude of the fault point; and solving and judging a fault point distance by using a simultaneous equation set of fault point negative sequence voltage amplitude expressions established by the head-end quantity and the tail-end quantity. According to the method, only the negative-sequence parameters at the two ends of the line and the negative-sequence impedance parameters of the fault line are needed, the cost is low, the practicability is good, a new thought is provided in the aspect of fault distance measurement, and technical support can be provided for reducing the manual line patrol range on site and improving the fault maintenance efficiency.

The invention discloses a wind power plant current collection system topology design method and system based on an improved genetic algorithm, relates to the technical field of system topology design, and aims to solve the core problems of collaborative optimization of fan grouping and submarine cable connection and accurate matching of current-carrying capacity and cost under high-dimensional constraints. According to the method, a parameter calculation module is combined with a formula to determine the maximum feeder fan number and the number range, a grouping connection module combines polar coordinate grouping with a minimum spanning tree algorithm to construct an initial scheme, a cost modeling module calculates the cost through an exclusive function, and an iterative optimization module dynamically adjusts the grouping and fan number based on an improved genetic algorithm. And the constraint screening module eliminates violation schemes through a penalty formula, and finally generates optimal and suboptimal topological graphs containing core parameters. According to the method, the high-dimensional optimization problem is layered and disassembled, the constraint satisfaction degree and the economical efficiency are considered, the submarine cable investment cost is greatly reduced, the design efficiency and the engineering feasibility are improved, and the method is suitable for topology design of various offshore wind plant current collection systems.

The application discloses a kind of wind farm power collection system topology design method and system based on improved genetic algorithm, it is related to system topology design technical field, to solve the core problem of wind turbine grouping and submarine cable connection collaborative optimization under high-dimensional constraint, accurate matching of load flow and cost.The method determines the maximum number of wind turbines and the number range of feeder by parameter calculation module simultaneous equations, grouping connection module constructs initial scheme by polar grouping combined with minimum spanning tree algorithm, cost modeling module calculates cost by exclusive function, iteration optimization module dynamically adjusts grouping and the number of wind turbines based on improved genetic algorithm, constraint screening module excludes illegal scheme by penalty formula, and finally generates optimal and suboptimal topology graph containing core parameters.The application divides high-dimensional optimization problem into layers, considers constraint satisfaction degree and economy, significantly reduces submarine cable investment cost, improves design efficiency and engineering feasibility, and is suitable for various offshore wind farm power collection system topology design.

This invention discloses a method, system, equipment, and medium for calculating pressure loss in perfluorohexanone (PFH) pipeline networks, relating to the field of computer-aided design technology. The method traverses the component tree of the pipeline network building information model using a three-dimensional geometry engine to form a set of pipeline network features; constructs a directed graph topology representing the flow direction of the pipeline network, calculates potential energy loss, converts local resistance into equivalent length, and constructs a weighted directional topology weight matrix; obtains hydraulic balance under single-phase flow conditions through iterative calculation of node simultaneous equations, and generates the corresponding static pressure distribution results of the pipeline network; recursively calculates the dryness and density changes of the gas-liquid mixture using the equation of state, introduces a two-phase flow friction multiplication factor for nonlinear compensation, outputs a three-dimensional dynamic compensation pressure field model, and performs rendering and digital judgment; this solution solves the problem of low accuracy in pressure loss calculation caused by the inability of existing technologies to accurately simulate the flash evaporation effect of PPHH, thus improving the reliability of pipeline network design.

The present application relates to the field of chemical process simulation, and particularly relates to a simultaneous equation simulation block solving method for process industry. The method comprises the following steps: extracting equations in an activated state and variables in an unfixed state from a chemical process model to form an original equation and variable set to be solved; converting the equation and variable set into a graph data structure reflecting system dependency according to a preset block strategy; using a strongly connected component identification algorithm to traverse the graph data structure to obtain a strongly connected component set, each strongly connected component in the strongly connected component set corresponding to a sub-equation block in the process equation; determining a solving order of the sub-equation blocks according to the dependency of the sub-equation blocks; solving each sub-equation block according to the solving order of the sub-equation blocks; and collecting the solving results of each sub-equation block to obtain a global solution of the complete chemical process model.

The application discloses a special-shaped aerodynamic contraction profile automatic design method based on a bicubic curve equation and belongs to the technical field of aviation aerodynamic wind tunnel test. The method solves the problems of many design parameters of the special-shaped aerodynamic profile, artificial calculation and the problem that existing mechanical software auxiliary design is difficult to complete due to a huge amount of calculation in the prior art. The method comprises the following steps: S1. calculating intersection coordinate values of an entrance profile control equation and calculating polar coordinate values corresponding to a feature point group; S2. dividing the aerodynamic contraction profile according to the feature point group, realizing equal phase and equal distance scattering of each zone in the polar coordinate, obtaining corresponding discrete coordinate point groups of entrance and exit sections and saving the discrete coordinate point groups; S3. using a simultaneous equation set, batch calculating coordinate values of each section space coordinate point group; S4. generating a curve according to the coordinate values of the space coordinate point group, generating a surface according to the curve and obtaining the required aerodynamic contraction profile. The application meets the actual engineering requirements and can be applied to the design of a wind tunnel special-shaped aerodynamic contraction profile.

The application discloses a low-cost rotor displacement recognition method for an eight-tooth one-pair-pole bearingless permanent-magnet sheet motor, four programmable Hall sensors are sequentially installed on radial stator slots, and the Hall sensors are orthogonally distributed on the radial stator slots in pairs; the application proposes a four-phase Hall solving algorithm based on an extreme decoupling coefficient; firstly, an output signal expression of each Hall element is theoretically derived; secondly, other coefficients are separated through offline rotor edge sticking, and a variation extreme value of a single coefficient is obtained to measure and calculate a Hall signal coefficient, so that the constraint of the coefficient in the Hall signal is solved, and a demodulation process of multiplying a sine and a cosine to obtain the same coefficient proportional displacement information in the original method is cancelled; finally, rotor displacements x and y are solved through simultaneous equations and data summation and difference equations; the application is suitable for an eight-tooth one-pair-pole bearingless two-phase motor, and only four Hall sensors are needed to recognize the displacement, so that the cost is greatly reduced.

This invention discloses a method for generating two-dimensional closed curves using NUBRS interpolation, belonging to the field of computer-aided geometric design. The method first obtains parameter values ​​and a knot array using the interpolation points passed by the closed path curve and the order of the closed path curve. The coordinates of all control points of the closed path curve are then solved using a system of simultaneous equations. Finally, the parameter values ​​are constrained according to the NUBRS generation formula to generate a closed NUBRS curve. Based on NUBRS curves, the method accurately represents various complex free-form closed path curves.

The application discloses a kind of measurement and calculation method for long-distance positioning electric power conductor and tower connecting point, it is characterized by comprising the following steps: step a, using the walking robot mounted on the conductor, combined with satellite positioning technology, sampling conductor space coordinates;Step b, according to the sampling coordinates, fitting the straight line equation under the overhead projection of conductor and the catenary equation under the orthographic projection of conductor;Step c, in the front of tower, find the auxiliary sampling point P that forms vertical line with connecting point P in orthographic projection plane s , and obtain the latitude and longitude coordinates of auxiliary sampling point P s , and fitting PP s Straight line overhead projection equation;Step d, by simultaneous equations and solving, obtain the space coordinates of connecting point P.The application does not need to climb to tower and conductor connecting place, can obtain the geographical space coordinates of tower cross arm conductor connecting point, ensure the safety of electric power operation, reduce work difficulty, reduce workload.

The invention relates to a simplified calculation method and device for settlement of a rigid pile composite foundation and a medium. The calculation method is simplified, on the basis of the assumption theory of an equivalent pier body and bilinear frictional resistance, an equivalent pier body model assumes that relative slippage details of pile soil are ignored, but the core of overall bearing and deformation is grasped, static balance analysis can be conducted on the whole pier body according to the theory, and calculation is greatly simplified; on the basis of the hypothesis theory of bilinear frictional resistance, the frictional resistance which is traditionally considered to be in triangular or nonlinear distribution is simplified to be in a distribution mode that the upper portion and the lower portion of the equal sinking face are both constant values, calculation of the axial force of the pile body is made to be linear, and therefore an equilibrium equation can be directly established; complex simultaneous equation solving in a traditional method is avoided, all parameters have clear physical significance and determination methods, and the method is suitable for different types of rigid pile composite foundations and can adapt to different engineering conditions through parameter adjustment.

The application discloses a cable cross-connection grounding system connection state live-line testing device and method. A joint or terminal or grounding box copper bar containing an independent lead is selected as a position for coupling injection and current induction testing, a live-line testing device is installed, loop impedance of any phase branch plus parallel impedance of other two phase branches under two different frequencies is tested and obtained, simultaneous equations are established by using Ohm's law and characteristics that resistance and inductance are irrelevant to frequency, impedance and resistance of each branch are calculated, and whether the grounding system has a connection defect is judged according to the calculated impedance and resistance, so that the electrical connection state of the cable cross-connection grounding system can be live-line detected, and the method is simple, convenient and efficient.

The application relates to a ship propulsion shafting alignment method and system, which comprises the following steps: establishing a first propulsion shafting alignment calculation model of a rear bearing of a stern shaft, defining a dimensionless proportion; obtaining an average specific pressure parameter based on the first propulsion shafting alignment calculation model; obtaining an equivalent support point axial distance proportion of the rear bearing of the stern shaft based on the average specific pressure of each support point and the dimensionless proportion; obtaining a second propulsion shafting alignment calculation model of single-point support of the rear bearing of the stern shaft based on the first propulsion shafting alignment calculation model and the equivalent support point axial distance proportion of the rear bearing of the stern shaft; fitting a distribution relationship formula of the wear of the rear bearing of the stern shaft along the axial direction based on the average specific pressure of each support point, the dimensionless proportion and the second propulsion shafting alignment calculation model, and calculating the equivalent wear of the rear bearing of the stern shaft; and adjusting the alignment state of the shafting to meet the specification requirements based on the equivalent wear of the rear bearing of the stern shaft and the wear of the front bearing of the stern shaft calculated according to the simultaneous equations of the front and rear bearings.

The invention belongs to the technical field of wind power generation, and provides a method, a system and equipment for estimating the effective wind speed of a fan with an unknown Cp curve and a medium. Deducing a pneumatic torque expression based on an extended Kalman filtering algorithm; and establishing simultaneous equations of the relation function of the pneumatic torque and the wind speed and the pneumatic torque expression, and solving a wind speed estimated value by using a Newton iteration method. According to the method, a known Cp curve is not needed, calculation is fast, robustness is high, high precision can be kept even under the complex turbulence working condition, only conventional operation data are needed, and the fan control precision and the power generation efficiency are greatly improved.

The present application relates to a kind of mountainous highway intelligent car anthropomorphism speed and trajectory planning method and device, wherein the method comprises the following steps: obtaining vehicle operating data, road geometry data and environmental data;Data preprocessing;With road geometry data and environmental data as explanatory variable, based on simultaneous equations model, establish theoretical vehicle operating speed and trajectory model;Based on the result of vehicle operating speed and trajectory model, establish anthropomorphism motion planning basic rule;Anthropomorphism motion planning model is constructed;Anthropomorphism motion planning model is corrected, and the behavior that due to driver manually drives vehicle and leads to deviation out of lane is avoided;The output of corrected anthropomorphism motion planning model is used as a kind of soft constraint, for intelligent car motion control module reference, realize anthropomorphism speed and trajectory planning.Compared with prior art, the complexity of driver driving behavior in mountainous highway is considered, more in line with the motion expectation of intelligent car driver.

The invention relates to a tunnel surrounding rock crustal stress rapid determination and verification method based on secant modulus, and belongs to the technical field of geotechnical engineering testing. The method comprises the following steps: firstly, sampling from a drill hole, sampling a directional thick rock core, and drilling cylindrical rock samples in all directions in the thick rock core; then preloading the rock sample for specified time; the method comprises the following steps: performing uniaxial compression on a sample at a higher level, synchronously acquiring axial and radial strain data by using a strain gauge in the compression process, and drawing a stress-strain curve; obtaining a strain value of a corresponding stress increment by using a secant modulus; calculating a strain difference by using a strain difference function to obtain a characteristic SMM curve; determining a normal stress through a strain difference function, establishing a simultaneous equation to solve a stress tensor, and obtaining the size and direction of a principal stress through eigenvalue analysis; the crustal stress measured by the secant modulus method is substituted into the calculated displacement, and degradation correction is carried out by considering the influence of construction disturbance; and comparing the actually measured displacement with a theoretical calculation value for verification.

The application discloses a bearing mechanical property calculation method considering rotor deformation, first, according to rotor load and rotor structure data, a rotor force-deformation coordination equation is established, and the external load borne by the bearing is calculated; a bearing load deformation equation considering bearing centrifugal force and gyroscopic force is established; then, the rotor force-deformation coordination equation and the bearing load deformation equation are solved by using a step-by-step approximation method, and the mechanical property of the bearing is obtained. The combined equation is solved step by step, the solving efficiency of the combined equation is improved, the mechanical property calculation precision of the bearing is improved, and a new method is provided for structure design and performance evaluation of a bearing-rotor system. The bearing and the rotor are regarded as a system, the nonlinear coupling relationship between the additional torque generated by the rotor deformation and the internal load distribution of the bearing is considered, and a more accurate method is provided for mechanical property analysis of the bearing.

This invention discloses a method for simultaneously measuring interfacial normal force and shear force based on contact resistance, comprising the following steps: Inspired by the yield criterion for material failure, establishing the relationship between normal stress and shear stress at the contact interface, thereby obtaining the relationship between the actual contact area of ​​the contact interface and the normal force and shear force, obtaining the relationship between the contact resistance of the contact interface and the actual contact area, and establishing the relationship between the interface contact resistance and the normal force and shear force; Based on the geometric deformation relationship of the micro-protrusions on the contact interface before and after the application of shear force, obtaining the relationship between the actual contact area of ​​the contact interface and the geometric deformation parameters of the micro-protrusions, obtaining the relationship between the normal force and shear force of the contact interface and the geometric deformation parameters of the micro-protrusions; Solving a system of equations, and calculating the normal force and shear force by measuring the contact resistance and substituting it into the system of equations. This invention enables the indirect measurement of the magnitude of the interface normal force and shear force by measuring the contact resistance of the contact interface.

The application relates to the field of backscattering technology and discloses a channel decoupling and equalization method and system for WiFi backscattering, which comprises the following steps: constructing an equation group based on a reference signal and a backscattering signal to solve a forward channel and a backscattering channel, wherein the equation group comprises 56 equations constructed by using a VHT-LTF field and 64 equations constructed by using a reference symbol; obtaining a matrix M through frequency domain to time domain conversion, time domain modulation and time domain to frequency domain conversion based on the reference symbol and a known tag sequence; solving a particular solution of the forward channel and a particular solution of the backscattering channel by simultaneously solving the equation group and combining the matrix M; and using a double-end channel equalization method to decode tag data through obtained first equalization signals and second equalization signals. The application designs a channel evaluation and equalization method, so that the modulation granularity of the tag can be further reduced, and the system transmission rate is improved.

A code conversion method includes a conversion step of converting a first calculation model applicable to an Ising model or a QUBO in which options are represented by a first representation method into a second calculation model applicable to an Ising model or a QUBO in which options are represented by a second representation method, in which the conversion step includes solving a simultaneous equation in which each component of an energy calculated on the basis of the first calculation model matches each component of an energy calculated on the basis of the second calculation model, and obtaining a submatrix of the Hamiltonian used in computation of the second calculation model.

The invention discloses an olefin conversion operation condition optimization method based on an OCU device mechanism model, and relates to the field of chemical engineering and industrial informatization. According to the method, the whole OCU device is decomposed into a plurality of sub-process models, and a whole-process model is built by using process simulation software. The model is calculated by adopting a simultaneous equation method, and a physical property database based on an SRK state equation is integrated. Model parameters are calibrated through real-time operation data of the device, a convergence environment and the model parameters are adjusted according to heat and material balance, and an accurate OCU device optimization model is established. And the optimization platform performs calculation by combining the optimization model, and provides optimal operation setting reference for operators, so that the improvement of the operation performance of the device is realized.

The embodiment of the invention provides a coupling method and device for a flexible bending wing and a distributed electric propulsion system, computer equipment and a medium, and relates to the technical field of coupling control. The method comprises the following steps that aircraft basic parameters are obtained, and the flexible bending wing and a ducted air inlet are obtained through sensor arrays arranged on the flexible bending wing and the ducted air inlet; collecting flight state data in real time; based on the fan characteristic model, calculating an equivalent jet flow momentum coefficient of the whole wing in the current state; a coupling total lift force model is constructed, the target lift force coefficient and the jet flow momentum coefficient are input into the coupling total lift force model, and a trailing edge camber instruction meeting the lift force requirement is obtained through simultaneous solving; performing dynamic response matching and symmetric consistency constraint processing on the trailing edge camber instruction and the rotating speed instruction of the distributed ducted fan; and performing real-time closed-loop correction on the trailing edge camber and the rotating speed of the distributed ducted fan based on the actual lift coefficient. According to the scheme, through the coupled cooperative control method, the propelling stability under the complex working condition is improved.

The application discloses a method for correcting test results of tower grounding resistance by using a clamp-on method, and comprises the following steps: judging whether the grounding resistance of a tower meets the use condition of the clamp-on method; if yes, the clamp-on method is used; if not, a three-pole method is used to test all towers one by one to obtain the real grounding resistance of all towers; the clamp-on method is used to test the grounding resistance of all towers meeting the condition to obtain the grounding resistance of all towers tested by using the clamp-on method; the three-pole method is used to test any tower to obtain the grounding resistance of the tower tested by using the three-pole method; a new set of simultaneous equations is obtained according to the relationship between the grounding resistance of all towers tested by using the clamp-on method and the grounding resistance of any tower tested by using the three-pole method; and the grounding resistance of all towers is obtained by solving the new set of simultaneous equations. The grounding resistance of all towers is obtained by solving the set of simultaneous equations, and the test efficiency of the grounding resistance of the tower is remarkably improved.

This invention discloses a testing device and method for the nanosecond pulse response delay of a metal oxide surge arrester. The device mainly includes a nanosecond pulse generator, a high-voltage power divider, an absorption resistor, the surge arrester under test, and an oscilloscope. The overall structure is simple, the platform size is flexible, and the measurement location is unrestricted. The testing method can accurately obtain the surge arrester's response delay with a single pulse current, avoiding errors caused by multiple pulses generated by the nanosecond pulse generator. Furthermore, this method obtains the wave velocity and attenuation coefficient on the transmission line, and then solves simultaneous equations based on the current waveform measured at the measurement point to derive the nanosecond pulse waveform and the current waveform flowing through the surge arrester. The surge arrester's response delay is quantified using specific formulas, resulting in high-precision calculations and reliable test results.

The invention relates to the technical field of power system simulation, in particular to an electromagnetic transient rapid representation method and system for energy storage equipment in power simulation. According to the method, topology and boundary information of energy storage equipment is acquired, and an internal high-frequency switching process is equivalent to a controlled source, so that topology simplification is realized; mapping a high-dimensional state space to a low-dimensional subspace to obtain a dimensionality reduction state description; and forming a simultaneous equation by associating a port variable with a power grid node variable, solving the simultaneous equation, obtaining a port time domain response, and performing dynamic correction based on a state error. According to the method, the calculation complexity of electromagnetic transient simulation of the energy storage equipment is remarkably reduced on the premise of ensuring the simulation precision.

The invention relates to the technical field of power electronic device testing, in particular to an interelectrode capacitance measuring system for an IGBT (Insulated Gate Bipolar Transistor) and a calibration and calculation method thereof, the measuring system is integrated with a measuring-calibration integrated clamp, a standard load is arranged in the clamp, and impedance error calibration can be performed on a whole measuring channel through a mode switching switch before measurement. Based on the method, a set of mandatory three-step sequential measurement process is established: measurement is completed under three specific physical topologies in sequence, error compensation is performed by using a calibration result, and finally, three independent interelectrode capacitors of the IGBT are accurately solved through simultaneous equations. Therefore, the problems that a related high-precision scheme (a VNA method) is high in cost and difficult to operate, a simple scheme (a traditional LCR method) is low in precision, irregular in process and the like are solved, IGBT interelectrode capacitance measurement close to laboratory precision can be achieved on the premise that engineering practicability and economical efficiency are guaranteed, and high precision, high repeatability and high safety are achieved.

The invention relates to a decision intention estimation method based on different strategy data adaptive inverse reinforcement learning, belongs to the technical field of artificial intelligence and machine learning, solves the problems of too high data sampling requirements and too large calculation burden in the prior art, and comprises the following steps: S1, obtaining expert demonstration data, and obtaining an optimal decision matrix in the expert demonstration data; s2, initializing an algorithm, giving an initial moment and an initial data set, establishing a simultaneous Lyapunov equation set about potential variables, and obtaining parameter estimation and related variable initial values; s3, constructing an adaptive online inverse reinforcement learning algorithm to obtain real-time parameter estimation and related variables; and S4, repeating S3, iterating to obtain parameter estimation and related variables at each online moment, taking the parameter estimation and related variables as an intelligent agent decision-making intention estimation result, and outputting a performance index function parameter matrix flow for describing the intelligent agent decision-making intention in real time.

The present application relates to remote sensing inversion technical field, disclose a kind of coverage inversion method fusing multi-angle model and space-time algorithm, comprising the following steps: extracting near infrared and red band value, calculate normalized difference vegetation index and integrate into multi-angle vegetation index data, using ESTARFM model to obtain basic prediction result by downsizing processing, extract mixed pixel reflectivity variation ratio, combined with moving window method to generate predicted high-resolution image, eliminate porosity by simultaneous equation of pixel two-part formula, component summation formula and Lambert Beer formula to obtain parameter simultaneous equation formula, derive double-angle solution formula by substituting specific angle, set sliding window with target pixel as center, the predicted high-resolution image is substituted into the double-angle solution formula to construct overdetermined equation group, obtain optimal pure vegetation NDVI and optimal pure bare soil NDVI, calculate output final coverage.The present application realizes parameter dynamic extraction, eliminates the interference of vegetation three-dimensional structure, improves inversion accuracy.