44results about How to "Short calculation time" patented technology

The invention discloses a charging pile dynamic power scheduling method, and relates to the technical field of charging pile power scheduling, and the method specifically comprises the following steps: obtaining node state information of each remaining normal node in real time, and generating a remaining bearing capability index and a load balancing index of each normal node; and constructing a power distribution evaluation model for the generated residual bearing capacity index and the load balancing index of each normal node, generating a power distribution evaluation coefficient of each normal node, analyzing the power distribution evaluation coefficient after generation, and evaluating the power bearing capacity of each normal node. Dividing each normal node into a high-bearing node, a middle-bearing node and a low-bearing node according to an evaluation result; and respectively adopting a corresponding power scheduling strategy according to a division result of each normal node. According to the invention, the problem that the power resource cannot be quickly and efficiently allocated after the fault node exits is solved, and efficient and stable dynamic balanced power scheduling and system operation are realized.

The invention discloses a remote sensing image target detection method based on an accelerated projection recurrent neural network. The method comprises the following steps: acquiring an original remote sensing image and calculating a remote sensing image autocorrelation matrix; constructing a linear constraint optimization mathematical model of filtering output based on the remote sensing image autocorrelation matrix, and obtaining a linear equation mathematical model through conversion; and solving the linear equation mathematical model by adopting an accelerated projection recurrent neural network algorithm, carrying out vector decomposition on a solving result to obtain a filtering coefficient vector, obtaining a remote sensing image after filtering output based on the filtering coefficient vector, and completing remote sensing image target detection. According to the remote sensing image target detection method based on the accelerated projection recurrent neural network, the convergence speed is higher, the calculation time is shorter, the image threshold classification capability is higher and the target detection precision is higher when the optimization problem is solved.

The invention discloses a high-robustness demosaicing bilinear interpolation method for polarization hyperspectral imaging, relates to the field of polarization spectral imaging, and solves the problems that in the prior art, when an optical system is located on a platform, vibration occurs, frames of collected images are lost, and spectrum aliasing exists in the optical system. The method comprises the following steps: collecting a group of polarized hyperspectral images, and carrying out high-robustness demosaicing reconstruction on the collected group of polarized hyperspectral images; the image I is recorded as the image of the obtained Kth spectral band, and the mosaic image II is an empty matrix and is used as an interpolation space of the image I; performing assignment on the mosaic image II in the horizontal direction; filling by using boundary conditions to complete interpolation in the horizontal direction; extracting polarization channel intensity images in four different directions in the mosaic image II, and performing vertical direction interpolation on the polarization channel intensity images; and polarization demosaicing reconstruction of the Kth spectral band image is completed. The method is suitable for the field of suppressing image artifacts generated in reconstruction.

The invention relates to an MIMO system modulation type identification method and system under communication interference, and belongs to the technical field of signal identification, and the method comprises the steps: carrying out the filtering, denoising and whitening preprocessing of a received MIMO signal; a five-dimensional fourth-order mutual cumulant tensor is constructed by using the preprocessed signal, and a separation matrix is estimated through parallel factorization, so that suppressing interference reconstruction separation is realized; constructing a three-dimensional tensor again for the separated and suppressed signal, and performing dimensionality reduction and feature extraction on the tensor by using a generalized tensor discriminant analysis algorithm; and performing Tucker decomposition on the tensor after dimension reduction to establish a kernel tensor and a factor matrix, and identifying the similarity of the kernel tensors of different modulation types according to the Euclidean distance criterion to complete MIMO signal modulation type identification. Simulation results show that the method has good recognition accuracy under the condition of low signal to interference ratio, and prior information such as channel coefficients and noise distribution is not needed.

The invention relates to the technical field of tire dynamics, in particular to a finite element simulation method, device, medium and program for lateral deviation and longitudinal slip composite slip characteristics of a tire. The invention provides a tire lateral deviation and longitudinal slip composite slip characteristic finite element simulation method which comprises the steps of finite element pretreatment, two-dimensional inflation and three-dimensional loading modeling analysis, steady-state simulation analysis, dynamic lateral deviation and longitudinal slip composite slip characteristic simulation modeling analysis and simulation result extraction, the method is high in simulation precision, and meanwhile the simulation efficiency is considered.

The invention discloses an SRAF rule construction method and device and a computer readable storage medium. The construction method comprises the steps that a test pattern is designed; establishing a photoetching model; and according to the test pattern and the photoetching model, constructing an optimal SRAF rule through a simulated annealing algorithm. And according to the test pattern and the photoetching model, constructing an optimal SRAF rule through a simulated annealing algorithm. The simulated annealing algorithm can obtain a better effect in the solution of the multi-extremum problem, and the SRAF rule is constructed through the simulated annealing algorithm, so that a better result can be obtained in the multi-extremum solution problem of the shape of the auxiliary graph and the distance between the auxiliary graph and the main graph, and the obtained SRAF rule has a better improvement effect on the process window of the main graph. Moreover, the method is higher in calculation efficiency and shorter in calculation time.

The invention relates to the field of cloud computing and deep reinforcement learning, and discloses a large-scale scientific computing workflow scheduling method based on deep reinforcement learning. According to the method, the graph attention network is introduced to extract workflow structure features, and a deep reinforcement learning sampling strategy based on a genetic algorithm is designed to generate an initial sample trajectory for a deep reinforcement learning agent, so that the training process is accelerated. Experiments show that a faster training speed can be obtained in a simulated cloud environment, the workflow is executed on the cloud computing platform according to the generated large-scale scientific computing workflow scheduling strategy, and better completion time, energy consumption and reliability can be reflected.

The invention belongs to the technical field of wind measurement radars, and relates to a weak Doppler signal extraction method based on random modulation narrow pulse, which comprises the following steps: generating a random signal by adopting a random sequence generator, generating a narrow pulse signal by adopting a sine wave generator, modulating the random signal on the narrow pulse signal to obtain a modulation signal, and transmitting the modulation signal; receiving an echo signal after the modulation signal is transmitted, carrying out frequency mixing and amplification filtering, and carrying out integral accumulation after random signal autocorrelation test to obtain a sampling sequence; doppler frequency shift is extracted according to the sampling sequence, wind speed calculation is carried out after preprocessing, decimation filtering and FFT filtering are carried out in sequence, and the wind speed and the wind direction are obtained. According to the invention, weak Doppler frequency shift signals can be extracted in a short near-space range.

The application is a kind of automation testing and application fluency testing method based on control matching, relating to the field of Android end automation testing.The application method comprises: setting a test script containing behavior simulation and fluency determination in advance, reading the corresponding test script after the Android end receives a test command; acquiring interface control structure positioning target control, executing simulation operation; starting timing after each simulation operation, and acquiring all control elements of the current interface when reaching the set time threshold, judging whether the interface is loaded or not in combination with the interface structure information recorded in the fluency determination script; obtaining application fluency determination result according to the interface loading condition at different time thresholds and returning to the server end.The application does not need to acquire ROOT permission, does not need to connect PC end for testing, is safer, is more convenient to use, reduces the adaptation cost, and is higher in testing accuracy and testing speed.

The application discloses a creep life prediction method for a nickel-based single crystal alloy considering a primary orientation deviation angle, and comprises the following steps: S1, acquiring a primary orientation deviation angle of a sample; S2, performing a creep test on the sample to obtain creep test data; S3, establishing a creep constitutive damage equation, constructing a correlation formula of a primary orientation deviation angle influence factor and the primary orientation deviation angle, and introducing the primary orientation deviation angle influence factor into a material degradation term in the creep constitutive damage equation for correction; S4, calibrating undetermined parameters in the corrected creep constitutive damage equation by using the creep test data; and S5, establishing a finite element model of a nickel-based single crystal high-temperature alloy sample based on the calibrated creep constitutive damage equation, and calculating the creep life of the sample by using the finite element method. The application changes the influence of the primary orientation deviation on the creep life from "systematic theory coupling" to "engineering precise quantification", and is more suitable for actual engineering application requirements.

The invention provides a method and system for calculating the opening number of eyelets, computer equipment and a medium. The method comprises the steps that a seam pressure calculation equation is obtained through coupling; carrying out two adjacent discharge capacities and two seam pressures, and carrying out coupling to obtain a calculation equation of the opening number of the eyelets; according to the obtained hole opening number calculation equation, the number N1 of opened holes in the pump stopping stage before the temporary plugging material is added, the number N2 of remaining opened holes in the pump starting stage after the temporary plugging material is added and the number N3 of newly opened holes in the pump stopping stage after the temporary plugging material is added are calculated in the fracturing construction process of each section; according to the hole number N1 and the hole number N2, the hole number N4 plugged by the temporary plugging material is obtained; and according to the hole number N3 and the hole number N4, the total hole opening number N5 after construction is completed is obtained. The system comprises a data collection module, a data processing module and a digital display module which are connected in sequence. The method has the beneficial effect that the method can be popularized and applied on a large scale.

The present invention discloses a method for eliminating the influence of BLE direction finding carrier frequency deviation to obtain phase delay, so that the phase delay only related to the direction of the incident signal source can be obtained with extremely low computing resources without estimating the carrier frequency deviation, thereby improving the accuracy of the BLE direction finding system in the presence of carrier frequency deviation and improving the stability of the BLE direction finding system under various operating temperature conditions.

The invention discloses a vehicle cluster cooperative positioning method based on an MDS-IDBO algorithm. The method comprises the following steps: step 1, constructing a multi-node cooperative positioning model; step 2, calculating relative coordinates by using an MDS-MAP algorithm; step 3, anchor node positioning correction is carried out by using an IDBO algorithm; 4, recalculating the distance square matrix and the relative coordinates of the nodes according to the step 1 and the step 2 by using the anchor nodes positioned and corrected in the step 3; 5, the distance square matrix and the node relative coordinates are utilized, a positioning error function serves as a fitness function to be substituted into an IDBO algorithm, an optimal coordinate transformation matrix is obtained, and absolute coordinates are calculated. The method is suitable for the technical field of communication, node positioning correction is carried out based on the MDS-MAP algorithm, errors caused by BDS / GPS positioning can be corrected to a great extent, and the accuracy of node positioning correction is improved. The provided IDBO algorithm is used for optimizing the MDS-MAP algorithm, the algorithm complexity is low, the calculation time is short, and the positioning precision is high.

The invention belongs to the field of transformers, and particularly relates to an LPST (linear phase shift transformer) air-gap magnetic field modeling method and an inductance parameter acquisition method thereof, and the method comprises the steps: analyzing the relation between the primary and secondary sides of the LPST and an air gap, and building an air-gap magnetic field analysis model considering the saturation of a ferromagnetic material; analyzing the relationship between the LPST cogging and the air gap, and building an air gap magnetic field correction model considering the cogging effect; and analyzing the relationship between the LPST side end and the air gap, and building an air gap magnetic field correction model considering the on-off of the iron core. The established LPST air-gap magnetic field model can consider the influence of ferromagnetic material saturation, cogging effect and iron core breaking at the same time, and the calculation precision of each inductance parameter is high, the calculation dimension is small, and the calculation time is short based on the LPST air-gap magnetic field model.

The invention relates to the field of power systems and automation thereof, in particular to a wind power cluster equivalent aggregation method, system and equipment suitable for subsynchronous oscillation analysis, and the method comprises the steps: S1, obtaining oscillation-considered fan grouping indexes suitable for actual engineering; s2, sampling and constructing sample data by a grid-connected side according to a link sequence, generating a clustering feature tree through a BIRCH algorithm according to a grouping index, and realizing grouping of wind turbine generators; and S3, aggregating different units in the same wind turbine generator cluster according to a topological structure and an algorithm result, aggregating the units in the same cluster according to a power consistency and capacity consistency principle to obtain an aggregation model, and performing equivalent aggregation on the wind turbine generator cluster according to a subsynchronous oscillation working condition. Compared with the prior art, the method has the advantages that the grouping indexes with high fitness are set and determined, and the BIRCH algorithm is applied to group the wind turbine generator sets, so that the technical problems of high grouping difficulty, low algorithm operation efficiency and inadaptability to subsynchronous oscillation operation state analysis of a specific wind power system in the prior art are solved.

The invention provides a capacitance current calculation method under distributed compensation. The method comprises the following steps: S1, fitting a zero-sequence transient voltage of a bus into a sine function of which a rectangular coordinate is a time function; s2, setting a line terminal zero sequence voltage initial value; s3, calculating bus zero-sequence voltage and line outgoing line zero-sequence current according to the line distribution parameter circuit model; s4, comparing the calculated value of the bus voltage with an actual value, and calculating an error; s5, judging whether the error meets the precision requirement or not, and if yes, outputting line outgoing line zero-sequence current; and if not, correcting the initial value of the zero-sequence voltage of the line terminal and repeating the steps S3-S5. According to the method, the calculation error can meet the requirement of capacitance current compensation, the calculation time of the capacitance current is short, and the requirement of online adjustment of the arc suppression coil can be met.

The application discloses a kind of self-adapting polishing system and method of single type deflection reach stacker girder, comprising: line sweep laser sensor, flexible tooling, controller, polishing robot end;Polishing method is: collecting girder laser scanning data and analysis, obtain girder pose data and adjust girder to standard pose according to girder pose data;Scan girder surface under standard pose, respectively obtain the deflection deformation of girder surface, the appearance of weld surface;And according to the appearance of weld surface, obtain weld width, height data and automatically adjust polishing processing parameters;Polishing trajectory is generated based on adjusted polishing processing parameters for controlling polishing robot end;The application is by integrating workpiece pose adjustment algorithm, polishing trajectory guiding algorithm and laser line sweep sensing technology, to solve adaptive polishing problem, improve artificial polishing effect, to realize the automation of adaptive polishing process with deflection deformation girder guiding and improve the environmental problem caused by polishing.

This application provides a method and apparatus for cutting quantum circuits. The method includes dividing the quantum circuit to obtain two bits of data; generating cutting position data based on the quantum circuit and the two bits; cutting the quantum circuit based on the cutting position data to obtain circuit cutting information, which includes two initial quantum sub-circuits; updating the wire cutting position and gate cutting position to obtain new cutting position data if either initial quantum sub-circuit does not meet the bit number constraint; determining the circuit cutting cost for generating new circuit cutting information based on the new cutting position data if the new initial quantum sub-circuit meets the bit number constraint; and determining the target circuit cutting information from multiple new circuit cutting information or from existing circuit cutting information based on the circuit cutting cost.

This invention belongs to the technical field of robotic arm motion planning and discloses a robotic arm motion planning method and system based on an error-driven adaptive gradient model. The method includes: S1 establishing the joint velocity at two adjacent moments with respect to the adaptive coefficient λ based on a Lyapunov function. n The gradient relationship is obtained; S2. The joint velocity at time n+1 is calculated using the gradient relationship, and the error between the calculated joint velocity at time n+1 and the adaptive coefficient is updated using the central finite difference method; S3. n = n+1; return to step S2 until the joint velocities at all times are obtained; S4. The joint position is calculated using the relationship between the joint velocity and position; the velocity of the end effector is calculated using the robotic arm motion planning model and the joint velocities at each time, thereby realizing the motion planning of the robotic arm. This invention solves the problem of low computational efficiency in robotic arm motion planning.

The invention discloses an intelligent calculation method and system for muddy water pressure of the whole line of an excavation face of a shield tunnel, and the method comprises the steps: processing a CAD drawing of the shield tunnel, reading the information of the drawing, and dividing the whole line of the shield tunnel; basic property parameters of the shield tunnel stratum are extracted, the drawing information and the basic property parameters are matched, and a shield tunnel whole-line stratum parameter database is established; a JavaScript programming language is adopted, and a shield tunnel excavation face muddy water pressure calculation theory library is established; a shield tunnel whole-line stratum parameter database is called, the basic property parameters are assigned to a shield tunnel excavation face muddy water pressure calculation theory library for calculation, a shield tunnel excavation face whole-line muddy water pressure theoretical value is obtained, and the theoretical value and the calculation ring number are matched and output in a table mode. The muddy water pressure calculation time of the excavation face of the shield tunnel can be shortened, the muddy water pressure calculation efficiency is improved, and the method has the advantages of being wide in application, high in reliability and convenient to apply.

This application belongs to the field of radar reconnaissance technology and provides a single-source radiation determination method based on temporal information entropy. It first calculates the Time Difference of Arrival (TDOA) of all adjacent pulses, then sets a clustering threshold and a first-class cluster center, and uses a clustering algorithm to obtain multiple clusters based on the TDOA. Next, it reduces the number of clusters and removes invalid elements by merging and removing harmonics. The distribution probability of each TDOA is calculated based on the number of indices in the merged clusters after harmonic removal. The temporal information entropy is calculated based on the distribution probability of each TDOA. A set information entropy discrimination threshold is set; if the temporal information entropy does not exceed the threshold, the pre-sorted signal is considered to be composed of a single radiation source; otherwise, it is composed of multiple radiation sources. This method can accurately determine the composition of radiation source signals in the radar pre-sorting clustering results, avoiding further PRI main sorting processing of single-source signals and reducing redundant sorting calculations.

The present invention discloses a method for initial hydraulic regulation of a heating metering system based on resistance identification, comprising the following steps: S1, obtaining a "resistance coefficient-valve opening" characteristic curve of each balancing valve; S2, forming a search range S of the resistance coefficient of each pipe section at the thermal inlet; n ; S3, adjust the balance valve of each thermal inlet pipe section to the fully open state, adjust the water pump to the design head and operate it in a constant pressure difference mode, and obtain the flow data of each thermal inlet pipe section as the measured flow; S4, use the particle swarm resistance optimization identification algorithm to search for the resistance coefficient in the range S n The thermal inlet resistance coefficient is generated and brought into the hydraulic calculation model to obtain the calculated flow rate of each thermal inlet pipe section; S5 substitutes the measured flow rate and the calculated flow rate into the objective function of the particle swarm resistance optimization identification algorithm to obtain the actual resistance coefficient S t ; S6. Substitute the calculated flow rate and the design flow rate into the objective function of the particle swarm resistance optimization identification algorithm to obtain the resistance coefficient S of the thermal inlet pipe section under the design flow rate c ; S7, S t and S c The difference is the change in the resistance of the balancing valve. The characteristic curve of the "resistance coefficient-valve opening" of the balancing valve is referred to to adjust each balancing valve to complete the initial hydraulic adjustment. Compared with the prior art, the beneficial effect of the present invention is that the flow is distributed completely according to the resistance coefficient, which is more accurate, fast and convenient.

A robust demosaic bilinear interpolation method for polarization hyperspectral imaging is presented, addressing issues such as optical system platform vibration, frame loss in acquired images, and spectral aliasing in existing technologies. The method includes acquiring a set of polarization hyperspectral images and performing robust demosaic reconstruction on them. Image I is the acquired image of the Kth spectral band, and mosaic image II is an empty matrix used as the interpolation space for image I. Horizontal values ​​are assigned to mosaic image II, and boundary conditions are applied for filling to complete the horizontal interpolation. The intensity images of the four polarization channels in different directions in mosaic image II are extracted, and vertical interpolation is performed on these polarization channel intensity images. This completes the polarization demosaic reconstruction of the Kth spectral band image. This method is applicable to suppressing image artifacts generated during reconstruction.

The invention discloses a method for optimizing and adjusting the demagnetization current of a distributed demagnetization system. The method comprises the following steps: searching for the flight speed and position of an individual, and calculating the fitness function value of each individual at the current position; for each individual, comparing the fitness function value of the current position with the fitness function value of the current individual historical optimal position, and determining the current individual historical optimal position; for all individuals, comparing the fitness function value of the current individual historical optimal position of each individual with the fitness function value of the current global historical optimal position, and determining the current global historical optimal position; judging whether a termination condition is met according to the current global historical optimal position, and terminating the adjustment of the demagnetization current if the demagnetization current corresponding to the current global historical optimal position reaches a preset value; the invention optimizes the calculation algorithm, simplifies the calculation procedure, and shortens the calculation time, so as to minimize the peak error of the residual magnetic field after the ship is demagnetized.

The present invention proposes a vibration signal recognition method and system based on improved multi-scale permutation entropy. The method includes: collecting vibration signals under different working states, performing multi-scale analysis based on sliding windows and sample quantiles on the collected vibration signals under different working states, and obtaining a sequence after multi-scale analysis; obtaining the permutation entropy of the sequence after multi-scale analysis, obtaining a multi-scale permutation entropy value, and using the value of the multi-scale permutation entropy as a feature vector of the vibration signal under different fault states; inputting the feature vector into a classifier based on machine learning, thereby realizing the recognition of vibration signals under different fault states and obtaining a fault state recognition result; and verifying the rationality of the recognition method according to the fault state recognition result. The present invention solves the defect that the conventional multi-scale permutation entropy fault diagnosis algorithm has insufficient ability in fault feature extraction, and has a short calculation time, thereby improving the fault diagnosis efficiency.

The invention discloses a satellite attitude super-agile sweep neural network minimum parameter learning sliding mode control method, which belongs to the technical field of microsatellite attitude control, and comprises the following steps: S1, in a satellite control process, establishing an actual output model of a neural network; s2, introducing a minimum parameter learning method to replace weight adjustment of a radial basis function neural network, and converting calculation of a neural network weight matrix W into solution of a single parameter # imgabs0 #; according to the method, the system state is converged to zero within finite time by using fast Terminal sliding mode control, and a switching item in a sliding mode is approached through a radial basis function neural network, so that the aim of weakening buffeting is fulfilled; meanwhile, a minimum parameter learning method is designed to replace weight adjustment of the radial basis function neural network, a large amount of operation of a weight matrix is avoided, the efficiency of the controller is effectively improved, and the method has the advantages of being good in convergence, imaging and control effect.

The present invention discloses a method for quickly calculating the cavitation fraction distribution of gas-liquid two-phase flow in a rod bundle sub-channel, comprising the steps of: dividing the flow channel cross section into a grid, discretizing the flow channel cross section into a plurality of grid regions, and calculating the cavitation fraction distribution of the gas-liquid two-phase flow in a rod bundle sub-channel according to the average axial flow velocity of the flow channel. <j f >, preset cavitation fraction distribution α1 and total viscosity μ T , use the numerical solution method to obtain the axial velocity distribution u of the liquid phase: Step 2: Calculate the turbulent kinetic energy distribution TKE according to the axial velocity distribution u of the liquid phase and the preset cavitation fraction distribution α1; Step 3: Calculate the different radial interface forces of the bubble according to the liquid phase velocity distribution u and the turbulent kinetic energy distribution TKE; Step 4: Determine the cavitation fraction distribution α2 according to the mutual balance of different interface forces; Step 5: Determine whether the difference between α2 and α1 exceeds the threshold ε. If not, output α2 and complete the calculation; if exceeded, correct α1: and return to step 1 to recalculate until the iteration converges.

The invention discloses a shallow sleeve type buried pipe soil heat conductivity distribution measuring method, and belongs to the technical field of ground source heat pumps. The method comprises the following steps: carrying out a distributed thermal response test on the shallow sleeve type buried pipe, and layering the shallow sleeve type buried pipe in the vertical direction; sequentially carrying out linear fitting on the measured value of the average temperature of the annular fluid of each layer and the natural logarithm of the time to obtain the slope of a linear relation; constructing a one-dimensional heat exchange model for each layer in sequence, wherein the one-dimensional heat exchange models are respectively used for simulating the average temperature of annular fluid of each layer and calculating the slope of a linear fitting relational expression of natural logarithms about time; for any layer, taking the absolute deviation between the slope based on the actual measurement value and the slope based on the simulation data of the one-dimensional heat exchange model of the layer as an objective function of the layer, and further calculating the soil heat conductivity corresponding to the minimum objective function, namely the measured soil heat conductivity of the layer; through the steps, the thermal conductivity of each layer of soil is measured in sequence. The method is high in precision, simple and suitable for engineering application.

The embodiment of the invention discloses a gradient information storage method and device, electronic equipment and a storage medium. The method comprises the following steps: acquiring first capacity information corresponding to a target video memory space for storing gradient information of a neural network model; obtaining video memory size information corresponding to each group of gradient information of the neural network model, wherein each group of gradient information corresponds to a unique computing node of the neural network model; and according to the first capacity information and video memory size information corresponding to each group of gradient information, storing to-be-reserved gradient information in the neural network model in a target video memory space, and deleting gradient information except the to-be-reserved gradient information, the video memory size occupied by the to-be-reserved gradient information is smaller than or equal to the storage capacity of the target video memory space. According to the scheme, gradient information generated in the neural network model training process can be adaptively reserved according to the storage capacity of the video memory space.

The invention belongs to the field of heat conduction equation solving, and particularly relates to a novel calculation method for solving internal temperature distribution of an annular region, which has the characteristic of rapid calculation. The main thought of the numerical calculation method is as follows: 1, reducing the order of a differential equation by using two new unknown numbers; 2, deriving the relationship of the temperatures of four adjacent nodes through integration; 3, deriving the relationship of the temperatures of any node and a boundary; and 4, solving the temperature of each node by using a boundary condition. The method can perform global solution, has the characteristics of small calculation amount and small occupied computer memory while ensuring the calculation precision, and can quickly obtain the internal temperature distribution condition of the annular region.