Planetary dynamic interactive simulation method and system

Abstract

The invention relates to the technical field of interactive simulation, in particular to a planet dynamic interactive simulation method and system, and the method comprises the following steps: collecting the position data and speed data of planets based on the data input at the initial stage of simulation starting, and measuring the orbit parameters and gravitational field of the planets by calculating the distance and the relative speed between the planets; and generating planet interaction feature data. Through real-time monitoring and dynamic data adjustment, the adaptability and accuracy of simulation are improved, the position and speed data of the planets are collected and analyzed in real time, and the simulation time step length is allowed to be adjusted in real time, so that any small change between the planets is reflected, and the adjustment ensures that when a key celestial body event, such as the planets, approaches, the planets, the planets can be accurately simulated. Simulation can reflect interaction and dynamic changes between planets in detail, waste of computing resources is reduced, simulation accuracy is improved, and the method is particularly more critical when processing complex celestial body interaction.

Description

Technical Field

[0001] The present invention relates to the technical field of interactive simulation, and in particular to a planetary dynamic interactive simulation method and system. Background Art

[0002] The field of interactive simulation technology mainly involves the technology of using computer systems to simulate real or fictional environments so that users can interact in real time in the environment. This technology is often used in training, education, research and entertainment, especially in the simulation of dynamic behavior of complex systems. Interactive simulation technology can provide a highly customizable and scalable environment in which physical, chemical or biological processes can be accurately simulated, and each operation of the user will instantly affect the results of the simulation.

[0003] Among them, the planetary dynamic interactive simulation method is a specific interactive simulation application, which is specially used to simulate and analyze the process of planetary motion and its dynamic interaction. This method can effectively simulate the gravitational effects between planets, orbital changes, and the interactions between planets and other celestial bodies (such as satellites and asteroids) in astronomical research and space mission planning. In addition, this method is also widely used in the field of education, helping students understand and explore the complex dynamics of the universe by providing an intuitive simulation interface.

[0004] Traditional methods cannot update data in a timely manner when simulating complex celestial dynamics, and it is difficult to accurately capture and reflect sudden events between planets, such as the rapid approach of asteroids to large planets. This slow response results in the simulation results being unable to reflect the real interaction between celestial bodies in real time, affecting the authenticity and scientificity of the simulation, and limiting the application effect in high-precision scientific research tasks and education. In addition, the lack of a mechanism to flexibly adjust the simulation accuracy limits the use in areas that require highly accurate control, such as deep space exploration mission planning and complex celestial behavior prediction. Summary of the invention

[0005] The purpose of the present invention is to solve the shortcomings of the prior art and to propose a planetary dynamic interactive simulation method and system.

[0006] In order to achieve the above object, the present invention adopts the following technical scheme: a planetary dynamic interactive simulation method, comprising the following steps:

[0007] S1: Based on the data input at the initial stage of simulation startup, collect the position data and velocity data of the planets, calculate the distance and relative velocity between the planets, measure the orbital parameters and gravitational field of the planets, and generate the planet interaction characteristic data;

[0008] S2: Based on the planetary interaction characteristic data, evaluating the computational load in the planetary interaction mechanical model, adjusting the simulation time step, optimizing the processing speed, and obtaining time step optimization data;

[0009] S3: Based on the time step optimization data, real-time monitoring and recording of planetary approach events and key celestial events, while detecting newly generated planetary position and velocity data during the simulation, and fine-tuning the time step to obtain refined and adjusted step data;

[0010] S4: Based on the refined and adjusted step size data, at the end of the simulation cycle, incremental data integration technology is used to merge new data generated in the simulation with existing data, and events in the simulation cycle are analyzed to generate comprehensive simulation cycle data.

[0011] The present invention is improved in that the steps for calculating the distance and relative speed between the planets are specifically as follows:

[0012] S111: extracting the position coordinates and velocity vector of each planet according to the data input at the initial stage of simulation startup, and obtaining a planet position and velocity information table;

[0013] S112: Based on the planetary position and velocity information table, use the formula:

[0014]

[0015] and

[0016]

[0017] Calculate the Euclidean distance d between the planets ij and relative velocity v rel，ij , where a1, a2 and a3 are weight coefficients for adjusting the accuracy of planetary position data, which are used to affect the accuracy of calculating the distance between planets; b1, b2 and b3 are weight coefficients for adjusting the accuracy of planetary velocity data, which are used to affect the accuracy of relative velocity calculation; x i ,y i and z i With x j ,y j and z j are the three-dimensional position coordinates of the two planets, used to represent the absolute positions of the planets in the simulation space, v x，i 、v y，i and v z，i With v x，j 、v y，j and v z，j It is the corresponding velocity vector component, which is used to calculate the relative motion speed between planets.

[0018] The present invention is improved in that the steps of measuring the planetary orbit parameters and the gravitational field are specifically as follows:

[0019] S121: Using the obtained distance and relative speed between the planets, use the formula:

[0020]

[0021] Calculate the gravitational force F ij , where c is the speed of light, v rel，ij is the relative velocity between the planets, d ij is the distance between the planets, m i and m j are the masses of the two planets, α is a parameter that adjusts the gravity calculation and incorporates the effects of motion, and G is the gravitational constant;

[0022] S122: Based on planetary position data and gravitational force F ij , using the formula:

[0023]

[0024] Calculate the orbital radius r of the planet orb and eccentricity e, where E is the total energy, L is the total angular momentum, m i and m j are the masses of the two planets, G is the gravitational constant;

[0025] S123: Combine the gravitational force between planets with the orbital parameters of the planets to analyze the stability and dynamic characteristics of the planetary system and obtain the planetary interaction characteristic data.

[0026] The present invention is improved in that the step of obtaining the time step optimization data is specifically as follows:

[0027] S211: Based on the planetary interaction characteristic data, the computational load of each simulation time step is calculated using the formula:

[0028]

[0029] Among them, d ij is the distance between the planets, v rel，ij is the relative velocity, v th is the speed threshold, α and β are adjustment coefficients, ∈ is a small adjustment coefficient of distance, which is used to adjust the effect of distance on load and obtain the calculated load L at each time step;

[0030] S212: According to the calculated load L and the set performance standard P, adjust the time step Δt using the formula:

[0031]

[0032] Where γ is a small constant added to avoid the denominator being zero, resulting in the adjusted time step Δt;

[0033] S213: Compare the adjusted time step Δt with the original time step, analyze the optimization effect, and obtain time step optimization data.

[0034] The present invention is improved in that the detection steps of the newly generated planet position and velocity data are specifically as follows:

[0035] S311: Based on the time step optimization data, the position and velocity data of the planet are collected in real time to obtain the current motion state of the planet;

[0036] S312: Analyze the continuously collected planet position and velocity data, and calculate the position change Δpos and velocity change Δvel in each simulation time step, using the formula:

[0037] Δpos=pos current -pos previous

[0038] and

[0039] Δvel=vel current -vel previous

[0040] Among them, pos current Andvel current Respectively represent the position and velocity of the current time step, pos previous Andvel previous Represents the position and velocity of the previous time step;

[0041] S313: Use the rate of change formula to calculate the trend of the planetary position and velocity and evaluate the stability and dynamic behavior of the simulation:

[0042]

[0043] The newly generated planetary position and velocity data are obtained, where α and β are weight factors used to control the impact of position and velocity changes, κ is an adjustment coefficient used to adjust the overall sensitivity of the trend, Δt is the current time step, r is the power of the time step, used to adjust the impact of time, and p and q are powers used to optimize the nonlinear impact of position and velocity changes.

[0044] The present invention is improved in that the step of obtaining the refined and adjusted step length data is specifically as follows:

[0045] S321: Analyze the newly generated planet position and velocity data with the current time step, and use the formula to adjust the time step to respond to the change of planet dynamics:

[0046]

[0047] Get the new time step Δtnew , where λ is the adjustment sensitivity parameter, controlling the sensitivity of trend changes to time step adjustment, σ is the index of trend impact strength, optimizing nonlinear impact, ν is the root number of time adjustment strength, used to refine the impact of time step adjustment, and Δt is the current time step;

[0048] S322: For the newly calculated time step Δt new Perform an evaluation and compare it with the current time step Δt to obtain the refined and adjusted step data.

[0049] The present invention is improved in that the step of merging the new data with the existing data is specifically as follows:

[0050] S411: Based on the refined and adjusted step size data, new data generated at the end of the simulation cycle is collected, and data existing before the start of the simulation cycle is extracted;

[0051] S412: Use the formula to integrate new and existing data:

[0052]

[0053] Among them, λ is the merging coefficient of new data, which is used to determine the weight of new data in the integration, α and β are adjustment parameters used to optimize the nonlinear effect when merging new data with existing data, and D new and D existing are the newly generated data and the existing data, respectively, combined The merged data set.

[0054] The present invention is improved in that the step of acquiring the comprehensive data of the simulation cycle is specifically as follows:

[0055] S421: Analyze the generated new data and existing data to identify key events and trends in the simulation cycle;

[0056] S422: Use data analysis techniques to classify and evaluate identified events, using the formula:

[0057]

[0058] Among them, θ is the analytical sensitivity parameter, which is used to adjust the degree of data analysis process, γ is the weighting factor, which is used to optimize the influence of raw data on the analysis results, k is the nonlinear response index, which optimizes the sensitivity of data features, κ is the normalization factor, and A results To analyze the results, D combined is the merged data set;

[0059] S423: According to the analysis result A results And the combined dataset D combined , using the formula:

[0060]

[0061] Get simulation cycle comprehensive data S cycle , where ξ and η are weighting parameters used to adjust the impact on the data, and ρ is the normalization coefficient used to balance the proportion of the data.

[0062] A planetary dynamic interactive simulation system, the system comprising:

[0063] The planetary data acquisition module collects the position and velocity data of the planets based on the data input at the initial simulation startup, calculates the distance and relative velocity between the planets, and measures the orbital parameters and gravitational field of the planets to obtain the planetary interaction characteristic data;

[0064] The load evaluation and optimization module evaluates the computational load of the planetary interaction mechanical model based on the planetary interaction characteristic data, adjusts the simulation time step according to the calculation result, reallocates the computing resources according to the adjusted time step, optimizes the processing speed, and obtains the time step optimization data;

[0065] The real-time monitoring and adjustment module uses the time step optimization data to monitor planetary approach events and key celestial body events, adjusts the time step according to the frequency of event occurrence, refines the time granularity of the simulation, and obtains refined and adjusted step data;

[0066] The data integration and analysis module applies incremental data integration technology based on the refined and adjusted step data, merges the newly generated data in the simulation with the existing data, analyzes the events in the simulation cycle, extracts key information, and obtains comprehensive data of the simulation cycle;

[0067] The system output and log module integrates the data of key celestial events and planetary dynamic interactions based on the comprehensive data of the simulation cycle, formulates simulation logs and charts, and forms a planetary dynamic analysis record.

[0068] Compared with the prior art, the advantages and positive effects of the present invention are:

[0069] In the present invention, the adaptability and accuracy of the simulation are improved through real-time monitoring and dynamic data adjustment. The position and velocity data of the planets are collected and analyzed in real time, allowing the simulation time step to be adjusted instantly to reflect any slight changes occurring between the planets. This adjustment ensures that during key celestial events such as planetary approaches, the simulation can reflect in detail the interactions and dynamic changes between the planets, which not only reduces the waste of computing resources but also improves the accuracy of the simulation, which is especially critical when dealing with complex celestial interactions. By merging new and old data and refining the adjusted step data, the consistency and reliability of the simulation results are further enhanced, providing more accurate data support. BRIEF DESCRIPTION OF THE DRAWINGS

[0070] Figure 1 A flow chart of a planetary dynamic interactive simulation method proposed by the present invention;

[0071] Figure 2 A schematic diagram of the calculation process of the distance and relative speed between planets in the planet dynamic interaction simulation method proposed by the present invention;

[0072] Figure 3 A schematic diagram of the measurement process of planetary orbit parameters and gravitational field in the planetary dynamic interactive simulation method proposed by the present invention;

[0073] Figure 4 A schematic diagram of the process of obtaining time step optimization data in the planetary dynamic interactive simulation method proposed by the present invention;

[0074] Figure 5 A schematic diagram of a detection process of newly generated planet position and velocity data in a planet dynamic interactive simulation method proposed by the present invention;

[0075] Figure 6 A schematic diagram of a process for obtaining step size data for refinement and adjustment in a planetary dynamic interactive simulation method proposed in the present invention;

[0076] Figure 7 A schematic diagram of the process of merging new data and existing data in the planetary dynamic interactive simulation method proposed by the present invention;

[0077] Figure 8 The present invention provides a schematic diagram of the process of obtaining comprehensive data of the simulation period in the planetary dynamic interactive simulation method. DETAILED DESCRIPTION

[0078] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0079] In the description of the present invention, it should be understood that the terms "length", "width", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside", etc., indicating positions or positional relationships, are based on the positions or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, in the description of the present invention, the meaning of "multiple" is two or more, unless otherwise clearly and specifically defined.

[0080] Example

[0081] See also Figure 1 The present invention provides a technical solution: a planetary dynamic interactive simulation method, comprising the following steps:

[0082] S1: Based on the data input at the initial stage of simulation startup, collect the position data and velocity data of the planets, calculate the distance and relative velocity between the planets, measure the orbital parameters and gravitational field of the planets, and generate the planet interaction characteristic data;

[0083] S2: Based on the planetary interaction characteristic data, evaluate the computational load in the planetary interaction mechanics model, adjust the simulation time step, optimize the processing speed, and obtain the time step optimization data;

[0084] S3: Based on the time step optimization data, the planetary approach events and key celestial events are monitored and recorded in real time. Meanwhile, the newly generated planetary position and velocity data are detected during the simulation, and the time step is fine-tuned to obtain the refined and adjusted step data.

[0085] S4: Based on the refined and adjusted step data, at the end of the simulation cycle, the incremental data integration technology is used to merge the new data generated in the simulation with the existing data, and the events within the simulation cycle are analyzed to generate comprehensive data for the simulation cycle.

[0086] The planetary interaction characteristic data include planetary trajectory data, velocity vector, and interaction force index. The time step optimization data include time adjustment factor, computational efficiency ratio, and resource consumption rate. The refined adjustment step data include adjustment accuracy information, step variation index, and response speed information. The simulation cycle comprehensive data include event frequency distribution data, simulation stability evaluation results, and data continuity evaluation results.

[0087] See also Figure 2 , the calculation steps of the distance and relative speed between planets are as follows:

[0088] S111: extracting the position coordinates and velocity vector of each planet according to the data input at the initial stage of simulation startup, and obtaining a planet position and velocity information table;

[0089] S112: Based on the planetary position and velocity information table, use the formula:

[0090]

[0091] and

[0092]

[0093] Calculate the Euclidean distance d between the planets ij and relative velocity v rel，ij, where a1, a2 and a3 are weight coefficients for adjusting the accuracy of planetary position data, which are used to affect the accuracy of calculating the distance between planets; b1, b2 and b3 are weight coefficients for adjusting the accuracy of planetary velocity data, which are used to affect the accuracy of relative velocity calculation; x i ,y i and z i With x j ,y j and z j are the three-dimensional position coordinates of the two planets, used to represent the absolute positions of the planets in the simulation space, v x，i 、v y，i and v z，i With v x，j 、v y，j and v z，j It is the corresponding velocity vector component, which is used to calculate the relative motion speed between planets.

[0094] Suppose there are two planets A and B, whose position coordinates and velocity vectors are as follows:

[0095] The position coordinates of Planet A: (x A ,y A 、z A ) = (1, 2, 3);

[0096] The position coordinates of planet B: (x B ,y B 、z B ) = (4, 6, 8);

[0097] The velocity vector of Planet A is: (v x，A 、v y，A 、v z，A ) = (2, 3, 1);

[0098] The velocity vector of planet B is: (v x，B 、v y，B 、v z，B )=(1,0,5).

[0099] Calculate the distance d between the planets AB

[0100] Calculate the difference in position coordinates between planets A and B:

[0101] Δx=x B -X A =4-1=3

[0102] Δy=y B -y A =6-2=4

[0103] Δz=z B -zA =8-3=5

[0104] Use the Euclidean distance formula to calculate the distance between the two planets, assuming that the weight coefficients a1, a2, and a3 are all set to 1 (no additional weighting):

[0105]

[0106] Calculate the relative speed v between the planets rel，AB

[0107] Calculate the velocity vector difference between planets A and B:

[0108] Δv x =v x，B -v x，A =1-2=-1

[0109] Δv y =v y，B -v y，A =0-3=-3

[0110] Δv z =v z，B -v z，A =5-1=4

[0111] The relative velocity between the two planets is calculated using the Euclidean formula for relative velocity, assuming that the weight coefficients b1, b2, and b3 are all set to 1 (no additional weighting):

[0112]

[0113] The distances and relative velocities between the planets are obtained from the position and velocity data input by the simulation.

[0114] See also Figure 3 , the measurement steps of planetary orbit parameters and gravitational field are as follows:

[0115] S121: Using the obtained distances and relative speeds between the planets, use the formula:

[0116]

[0117] Calculate the gravitational force F ij , where c is the speed of light, v rel，ij is the relative velocity between the planets, d ij is the distance between the planets, m i and m j are the masses of the two planets, α is a parameter that adjusts the gravity calculation and incorporates the effects of motion, and G is the gravitational constant;

[0118] S122: Based on planetary position data and gravitational force Fij , using the formula:

[0119]

[0120] Calculate the orbital radius r of the planet orb and eccentricity e, where E is the total energy, L is the total angular momentum, m i and m j are the masses of the two planets, G is the gravitational constant;

[0121] S123: Combine the gravitational force between planets with the orbital parameters of the planets to analyze the stability and dynamic characteristics of the planetary system and obtain the planetary interaction characteristic data.

[0122] Calculate the gravitational force F ij :

[0123] Assume that the masses of the two planets are:

[0124] m i =5.97×10 24 lg (earth-like mass) and m j =7.35×10 22 kg (similar to the mass of the moon);

[0125] The distance between the two planets d ij =384,400km (distance between the earth and the moon);

[0126] Gravitational constant G = 6.674 × 10 -11 m 3 kg -1 s -2 ;

[0127] Relative speed v rel，ij =1,022m / s;

[0128] The speed of light c = 3 × 10 8 m / s;

[0129] Adjust parameter α=0.001.

[0130] Use the formula to substitute specific values:

[0131]

[0132] F ij =1.982×10 20 (1+1.115×10 -8 )

[0133] F ij ≈1.982×10 20 N

[0134] Calculate the planetary orbit parameters:

[0135] Using the calculated gravitational force F ij , assuming that the total energy of the planetary system is E = -1.2×10 29 J, total angular momentum L = 2.9 × 10 34 Js.

[0136] Calculate the orbit radius r orb And the eccentricity e:

[0137]

[0138] Based on the calculated orbital parameters and gravitational forces, the stability and dynamic characteristics of the planetary system are analyzed and output as a comprehensive report.

[0139] See also Figure 4 , the specific steps for obtaining time step optimization data are:

[0140] S211: Based on the planetary interaction characteristic data, the computational load of each simulation time step is calculated using the formula:

[0141]

[0142] Among them, d ij is the distance between the planets, v rel，ij is the relative velocity, v th is the speed threshold, α and β are adjustment coefficients, ∈ is a small adjustment coefficient of distance, which is used to adjust the effect of distance on load and obtain the calculated load L at each time step;

[0143] S212: According to the calculated load L and the set performance standard P, adjust the time step Δt using the formula:

[0144]

[0145] Where γ is a small constant added to avoid the denominator being zero, resulting in the adjusted time step Δt;

[0146] S213: Compare the adjusted time step Δt with the original time step, analyze the optimization effect, and obtain time step optimization data.

[0147] Assumed parameters:

[0148] The distance between the planets d ij =10∧8m;

[0149] Relative speed v rel，ij =1000m / s;

[0150] Speed ​​threshold vth =500m / s;

[0151] Adjustment coefficient α = 0.5;

[0152] Speed ​​index β = 2;

[0153] Distance adjustment factor ∈ = 0.1;

[0154] Performance standard P = 0.01s;

[0155] Small constant γ = 0.001.

[0156] Substitute the specific values ​​and calculate the load L:

[0157]

[0158] Substitute specific values ​​and adjust the simulation time step Δt:

[0159]

[0160]

[0161] Δt≈0.31623s

[0162] Calculation load L: calculated using the output parameters of the simulation model. According to the distance and relative speed between the planets, the adjusted parameters α, β and ∈ are used to accurately control the sensitivity of the calculation load to the distance and speed;

[0163] Time step Δt: The time step is adjusted according to the calculated load L to ensure a balance between simulation efficiency and accuracy. The small constant γ ensures that the time step is not too large at extremely low loads.

[0164] See also Figure 5 , the detection steps of the newly generated planet position and velocity data are as follows:

[0165] S311: Based on the time step optimization data, the position and velocity data of the planet are collected in real time to obtain the current motion state of the planet;

[0166] S312: Analyze the continuously collected planet position and velocity data, and calculate the position change Δpos and velocity change Δvel in each simulation time step, using the formula:

[0167] Δpos=pos current -pos previous

[0168] and

[0169] Δvel=vel current -vel previous

[0170] Among them, pos current Andvel current Respectively represent the position and velocity of the current time step, pos previous Andvel previous Represents the position and velocity of the previous time step;

[0171] S313: Use the rate of change formula to calculate the trend of the planetary position and velocity and evaluate the stability and dynamic behavior of the simulation:

[0172]

[0173] The newly generated planetary position and velocity data are obtained, where α and β are weight factors used to control the impact of position and velocity changes, κ is an adjustment coefficient used to adjust the overall sensitivity of the trend, Δt is the current time step, r is the power of the time step, used to adjust the impact of time, and p and q are powers used to optimize the nonlinear impact of position and velocity changes.

[0174] Assumed parameters and data:

[0175] Assume that at a certain simulation time point, the position of planet A is pos current =(300,300,300) km;

[0176] The position pos of the previous time step previous =(290, 290, 290) km;

[0177] Current speed vel current = (2, 2, 2) km / s;

[0178] The velocity vel at the previous time step previous = (1.5, 1.5, 1.5) km / s;

[0179] Time step Δt = 10 seconds;

[0180] Weight factors α = 0.6, β = 0.4;

[0181] Power p = 1.5, q = 1.2;

[0182] Time power r = 0.8;

[0183] Adjustment coefficient κ=0.75.

[0184] Substitute the values ​​and calculate the position change (Δpos):

[0185] Δpos=pos current -pos previous= (300-290, 300-290, 300-290) = (10, 10, 10)km

[0186] Calculate the velocity change (Δvel):

[0187] Δvel=vel current -vel previous = (2-1.5, 2-1.5, 2-1.5) = (0.5, 0.5, 0.5)km / s

[0188] Substitute the values ​​and calculate the trend:

[0189]

[0190] trend≈2.28

[0191] See also Figure 6 , the specific steps for obtaining the refined and adjusted step length data are:

[0192] S321: Based on the newly generated planet position and velocity data, analyze the current time step and use the formula to adjust the time step to respond to changes in planetary dynamics:

[0193]

[0194] Get the new time step Δt new , where λ is the adjustment sensitivity parameter, controlling the sensitivity of trend changes to time step adjustment, σ is the index of trend impact strength, optimizing nonlinear impact, ν is the root number of time adjustment strength, used to refine the impact of time step adjustment, and Δt is the current time step;

[0195] S322: For the newly calculated time step Δt new Perform an evaluation and compare it with the current time step Δt to obtain the refined and adjusted step data.

[0196] Assume the following data:

[0197] Current time step (Δt): 0.1 seconds;

[0198] Dynamic trend (trend): 0.5;

[0199] Adjust the sensitivity parameter (λ): 0.05;

[0200] Index of trend influence strength (σ): 2;

[0201] Number of roots of time-adjusted intensity (ν): 3.

[0202] Substitute specific values:

[0203]

[0204] Δt new =0.1 0.99862

[0205] Δt new = 0.099862 seconds

[0206] The newly calculated time step Δt new It is 0.099862 seconds, which is slightly reduced compared to the original time step, indicating that it has been fine-tuned according to the dynamic trend of the planet. The new time step is compared with the original time step, and the specific value and proportion of the fine-tuning are recorded to provide data support for the accuracy and performance optimization of subsequent simulation cycles.

[0207] See also Figure 7 , the steps to merge new data with existing data are as follows:

[0208] S411: Based on the refined and adjusted step size data, new data generated at the end of the simulation cycle is collected, and data existing before the start of the simulation cycle is extracted;

[0209] S412: Use the formula to integrate new and existing data:

[0210]

[0211] Among them, λ is the merging coefficient of new data, which is used to determine the weight of new data in the integration, α and β are adjustment parameters used to optimize the nonlinear effect when merging new data with existing data, and D new and D existing are the newly generated data and the existing data, respectively, combined The merged data set.

[0212] Assume the following data:

[0213] New data (D new ): 50 (indicates the value of the newly added data);

[0214] Existing data (D existing ): 100 (indicates the value of existing data);

[0215] Merge coefficient (λ): 0.3 (determines the weight of new data in the merge);

[0216] Tuning parameters (α and β): 2, 1 (enhance the nonlinear effect of new data).

[0217] Substitute specific values ​​for calculation:

[0218]

[0219] Dcombined ≈273.33

[0220] D new New data is added to reflect the latest simulation results, D existing The existing data, including the accumulated results from previous simulation cycles, the merging coefficient λ, affects the relative importance of the new data in the integration process, and the nonlinear adjustment parameters α and β are used to strengthen the influence of new and old data in the data merging process. In this case, a higher α means that the new data is given a stronger nonlinear weighting, while a lower β means that the existing data is weighted in a linear manner.

[0221] Through this calculation process, it can be seen that the merged data D combined The value is about 273.33, which combines the influence of new data and existing data and increases the influence of new data by nonlinear adjustment parameters. This method ensures the flexibility and adaptability of data merging, so that the merged data can better reflect the latest development and historical accumulation of simulation.

[0222] See also Figure 8 ,The specific steps for obtaining comprehensive data of the simulation cycle are:

[0223] S421: Analyze the generated new data and existing data to identify key events and trends during the simulation cycle;

[0224] S422: Use data analysis techniques to classify and evaluate identified events, using the formula:

[0225]

[0226] Among them, θ is the analytical sensitivity parameter, which is used to adjust the degree of data analysis process, γ is the weighting factor, which is used to optimize the influence of raw data on the analysis results, k is the nonlinear response index, which optimizes the sensitivity of data features, κ is the normalization factor, and A results To analyze the results, D combined is the merged data set;

[0227] S423: According to the analysis result A results And the combined dataset D combined , using the formula:

[0228]

[0229] Get simulation cycle comprehensive data S cycle , where ξ and η are weighting parameters used to adjust the impact on the data, and ρ is the normalization coefficient used to balance the proportion of the data.

[0230] Assume the following data:

[0231] Merge dataset (D combined ): Assumed to be 120 (number of data points);

[0232] Analytical sensitivity parameter (θ): assumed to be 0.8 (used to adjust the sensitivity of the analysis);

[0233] Weighting factor (γ): assumed to be 1.5 (to enhance the influence of data features);

[0234] Nonlinear response exponent (k): assumed to be 2 (used to increase the nonlinear effect of the data);

[0235] Normalization factor (k): assumed to be 3 (for stabilizing the analysis results);

[0236] Weighting parameter (ξ): assumed to be 2 (to adjust the impact of the analysis results);

[0237] Weighting parameter (η): assumed to be 1 (adjusting for the effect of merging datasets);

[0238] Normalization coefficient (ρ): Assumed to be 5 (to balance the proportion of the comprehensive data).

[0239] Apply data analysis technology and substitute specific values ​​for calculation:

[0240]

[0241] Assume that the analyze function simply converts D combined Multiply by θ:

[0242]

[0243] Substitute the numerical values ​​for calculation and obtain the comprehensive data of the simulation cycle:

[0244]

[0245] Ensure the accuracy and suitability of the final output by processing and analyzing simulation data.

[0246] Planetary dynamic interactive simulation system, the system includes:

[0247] The planetary data acquisition module collects the position and velocity data of the planets based on the data input at the initial simulation startup, calculates the distance and relative velocity between the planets, and measures the orbital parameters and gravitational field of the planets to obtain the planetary interaction characteristic data;

[0248] The load assessment and optimization module assesses the computational load of the planetary interaction mechanics model based on the planetary interaction characteristic data, adjusts the simulation time step according to the calculation results, reallocates computing resources according to the adjusted time step, optimizes the processing speed, and obtains the time step optimization data;

[0249] The real-time monitoring and adjustment module uses time step optimization data to monitor planetary approach events and key celestial events, adjusts the time step according to the frequency of events, refines the time granularity of the simulation, and obtains refined and adjusted step data;

[0250] The data integration and analysis module applies incremental data integration technology based on the refined and adjusted step data to merge the newly generated data and existing data in the simulation, analyze the events in the simulation cycle, extract key information, and obtain comprehensive data of the simulation cycle;

[0251] The system output and log module integrates the data of key celestial events and planetary dynamic interactions based on the comprehensive data of the simulation period, formulates simulation logs and charts, and forms a record of planetary dynamic analysis.

[0252] The above are only preferred embodiments of the present invention and are not intended to limit the present invention in other forms. Any technician familiar with the profession may use the technical contents disclosed above to change or modify them into equivalent embodiments with equivalent changes and apply them to other fields. However, any simple modification, equivalent change and modification made to the above embodiments based on the technical essence of the present invention without departing from the technical solution of the present invention still falls within the protection scope of the technical solution of the present invention.

Claims

1. A planetary dynamic interactive simulation method, characterized in that: The following steps are involved: Based on the data input at the initial stage of simulation startup, the position and velocity data of the planets are collected, and the distance and relative velocity between the planets are calculated to measure the orbital parameters and gravitational field of the planets, thus generating the planetary interaction characteristic data; Based on the planetary interaction characteristic data, the computational load in the planetary interaction mechanical model is evaluated, the simulation time step is adjusted, and the processing speed is optimized to obtain time step optimization data; Based on the time step optimization data, the occurrence of planetary approach events and key celestial events is monitored and recorded in real time, while the newly generated planetary position and velocity data in the simulation process is detected, and the time step is fine-tuned to obtain refined and adjusted step data; Based on the refined and adjusted step size data, at the end of the simulation cycle, the incremental data integration technology is used to merge the new data generated in the simulation with the existing data, and the events in the simulation cycle are analyzed to generate comprehensive simulation cycle data.

2. The planetary dynamic interactive simulation method according to claim 1, characterized in that: The specific steps for calculating the distance and relative speed between the planets are: According to the data input at the initial stage of simulation startup, the position coordinates and velocity vector of each planet are extracted to obtain the planet position and velocity information table; Based on the planetary position and velocity information table, use the formula: and Calculate the Euclidean distance d between the planets ij and relative velocity v rel，ij , where a1, a2 and a3 are weight coefficients for adjusting the accuracy of planetary position data, which are used to affect the accuracy of calculating the distance between planets; b1, b2 and b3 are weight coefficients for adjusting the accuracy of planetary velocity data, which are used to affect the accuracy of relative velocity calculation; x i ,y i and z i With x j ,y j and z j are the three-dimensional position coordinates of the two planets, used to represent the absolute positions of the planets in the simulation space, v x，i 、v y，i and v z，i With v x，j 、v y，j and v z，j It is the corresponding velocity vector component, which is used to calculate the relative motion speed between planets.

3. The planetary dynamic interactive simulation method according to claim 2, characterized in that: The steps for measuring the planetary orbit parameters and gravitational field are specifically as follows: Using the obtained distances and relative speeds between the planets, use the formula: Calculate the gravitational force F ij , where c is the speed of light, v rel，ij is the relative velocity between the planets, d ij is the distance between the planets, m i and m j are the masses of the two planets, α is a parameter that adjusts the gravity calculation and incorporates the effects of motion, and G is the gravitational constant; Based on the planetary position data and the gravitational force F ij , using the formula: Calculate the orbital radius r of the planet orb and eccentricity e, where E is the total energy, L is the total angular momentum, m i and m j are the masses of the two planets, G is the gravitational constant; By combining the gravitational force between planets with the planetary orbit parameters, the stability and dynamic characteristics of the planetary system are analyzed to obtain the planetary interaction characteristic data.

4. The planetary dynamic interactive simulation method according to claim 1, characterized in that: The steps for obtaining the time step optimization data are specifically as follows: Based on the planetary interaction characteristic data, the computational load of each simulation time step is calculated using the formula: Among them, d ij is the distance between the planets, v rel，ij is the relative velocity, v th is the speed threshold, α and β are adjustment coefficients, ∈ is a small adjustment coefficient of distance, which is used to adjust the effect of distance on load and obtain the calculated load L at each time step; According to the calculated load L and the set performance standard P, the time step Δt is adjusted using the formula: Where γ is a small constant added to avoid the denominator being zero, resulting in the adjusted time step Δt; Compare the adjusted time step Δt with the original time step, analyze the optimization effect, and obtain the time step optimization data.

5. The planetary dynamic interactive simulation method according to claim 1, characterized in that: The detection steps of the newly generated planetary position and velocity data are specifically as follows: Based on the time step optimization data, the position and velocity data of the planet are collected in real time to obtain the current motion state of the planet; The continuously collected planet position and velocity data are analyzed to calculate the position change Δpos and velocity change Δvel in each simulation time step using the formula: Δpos=pos current -pos previous and Δvel=vel current -vel previous Among them, pos current Andvel current Respectively represent the position and velocity of the current time step, pos previous Andvel previous Represents the position and velocity of the previous time step; Use the rate of change formula to calculate the trend of the planetary position and velocity and evaluate the stability and dynamic behavior of the simulation: The newly generated planetary position and velocity data are obtained, where α and β are weight factors used to control the impact of position and velocity changes, K is the adjustment coefficient used to adjust the overall sensitivity of the trend, Δt is the current time step, r is the power of the time step, used to adjust the impact of time, and p and q are powers used to optimize the nonlinear impact of position and velocity changes.

6. The planetary dynamic interactive simulation method according to claim 5, characterized in that: The steps for obtaining the step length data for the refinement adjustment are specifically as follows: Based on the newly generated planet position and velocity data, and the current time step, the time step is adjusted using the formula to respond to changes in planetary dynamics: Get the new time step Δt new , where λ is the adjustment sensitivity parameter, controlling the sensitivity of trend changes to time step adjustment, σ is the index of trend impact strength, optimizing nonlinear impact, ν is the root number of time adjustment strength, used to refine the impact of time step adjustment, and Δt is the current time step; For the newly calculated time step Δt new Perform an evaluation and compare it with the current time step Δt to obtain the refined and adjusted step data.

7. The planetary dynamic interactive simulation method according to claim 1, characterized in that: The steps of merging the new data with the existing data are specifically as follows: Based on the refined and adjusted step size data, new data generated at the end of the simulation cycle is collected, and data existing before the start of the simulation cycle is extracted; Use the formula to combine new and existing data: Among them, λ is the merging coefficient of new data, which is used to determine the weight of new data in the integration, α and β are adjustment parameters used to optimize the nonlinear effect when merging new data with existing data, and D new and D existing are the newly generated data and the existing data, respectively, combined The merged data set.

8. The planetary dynamic interactive simulation method according to claim 7, characterized in that: The steps for obtaining the comprehensive data of the simulation cycle are specifically as follows: Analyze the new and existing data to identify key events and trends during the simulation cycle; Use data analysis techniques to classify and evaluate the identified events, using the formula: Among them, θ is the analysis sensitivity parameter, which is used to adjust the degree of data analysis process, γ is the weighting factor, which is used to optimize the influence of raw data on analysis results, k is the nonlinear response index, which optimizes the sensitivity of data features, K is the normalization factor, and A results To analyze the results, D combined is the merged data set; According to the analysis results A results And the combined dataset D combined , using the formula: Get simulation cycle comprehensive data S cycle , where ξ and η are weighting parameters used to adjust the impact on the data, and ρ is the normalization coefficient used to balance the proportion of the data.

9. Planetary dynamic interactive simulation system, characterized in that: According to any one of claims 1 to 8, the planetary dynamic interactive simulation method is implemented, and the system comprises: The planetary data acquisition module collects the position and velocity data of the planets based on the data input at the initial simulation startup, calculates the distance and relative velocity between the planets, and measures the orbital parameters and gravitational field of the planets to obtain the planetary interaction characteristic data; The load evaluation and optimization module evaluates the computational load of the planetary interaction mechanical model based on the planetary interaction characteristic data, adjusts the simulation time step according to the calculation result, reallocates the computing resources according to the adjusted time step, optimizes the processing speed, and obtains the time step optimization data; The real-time monitoring and adjustment module uses the time step optimization data to monitor planetary approach events and key celestial body events, adjusts the time step according to the frequency of event occurrence, refines the time granularity of the simulation, and obtains refined and adjusted step data; The data integration and analysis module applies incremental data integration technology based on the refined and adjusted step data, merges the newly generated data in the simulation with the existing data, analyzes the events in the simulation cycle, extracts key information, and obtains comprehensive data of the simulation cycle; The system output and log module integrates the data of key celestial events and planetary dynamic interactions based on the comprehensive data of the simulation cycle, formulates simulation logs and charts, and forms a planetary dynamic analysis record.

Images