An adaptive temperature compensation method based on hemispherical resonator

Abstract

The present application relates to the technical field of temperature compensation, in particular to an adaptive temperature compensation method based on a hemispherical resonator, comprising the following steps: obtaining a node signal to construct a temperature sequence, extracting a strain distribution to form a chain path, generating a compensation channel and outputting a temperature compensation scheme. In the present application, a time sequence interference signal sequence facing the three-dimensional structure distribution of nodes is constructed, the strain dominant area is identified based on the difference in interference phase and reflectivity change, the chain path is constructed by fusing the polar coordinate position and the node connection relationship, the low-interference path segment in signal propagation is dynamically extracted, the continuous phase difference of the node response is obtained and sorted to generate a channel set, the multi-path driving instruction generation is realized by combining the mapping of interference characteristics, multi-point dynamic sensing and signal reconstruction can be realized under the three-dimensional thermal distribution field, the chain regulation and control is completed by guiding the compensation path to the thermal disturbance direction and the stress gradient, and the decoupling ability and response accuracy of temperature compensation in the non-steady thermal field are improved.

Description

Technical Field

[0001] This invention relates to the field of temperature compensation technology, and in particular to an adaptive temperature compensation method based on a hemispherical harmonic oscillator. Background Technology

[0002] The field of temperature compensation technology involves the dynamic control and precise correction of equipment performance deviations caused by temperature changes. This includes physical modeling of thermistor structures, construction of compensation curves, design of temperature sensing and response mechanisms, and error minimization in thermal environments. By establishing a temperature-performance relationship model, stable performance control of sensors, micromechanical structures, or inertial components under different thermal environments can be achieved. Traditional adaptive temperature compensation methods address the performance drift of resonant devices under temperature variations. These methods utilize real-time temperature values ​​obtained from surface thermistors and correct the resonant frequency shift using a preset linear or polynomial function. Typically, this involves fitting a temperature-frequency function relationship using static calibration data, and during operation, adjusting the resonant frequency based on the measured shell temperature value using a lookup table, and constructing a correction channel to adjust the resonant response.

[0003] Traditional adaptive temperature compensation mainly relies on surface thermistors to collect the shell temperature signal and construct a one-dimensional temperature-frequency mapping relationship based on static calibration data. In actual operation, the shell temperature and resonant frequency offset are compared by static function fitting. It lacks the ability to dynamically sense the temperature changes inside the resonator structure and cannot handle the dynamic changes of the three-dimensional temperature field inside the structure. It ignores the non-uniformity and directionality of temperature disturbance propagation in space, which makes it unable to accurately restore the thermal distribution state when there are multi-source coupled temperature disturbances. The compensation path is only constructed along a fixed direction and a single channel without considering the three-dimensional complex characteristics of the thermal distribution inside the structure. Furthermore, it does not model the thermal anisotropy of the material, which easily leads to difficulties in frequency offset decoupling in multi-modal response, thereby causing the amplification of higher-order thermal disturbance errors and affecting the stability and accuracy of temperature compensation. Summary of the Invention

[0004] To address the technical problems existing in the prior art, embodiments of the present invention provide an adaptive temperature compensation method based on a hemispherical harmonic oscillator, comprising the following steps:

[0005] To achieve the above objectives, the present invention adopts the following technical solution: an adaptive temperature compensation method based on a hemispherical harmonic oscillator, comprising the following steps:

[0006] S1: Obtain the three-dimensional coordinates and numbers of the sapphire fiber nodes in the spherical structure of the hemispherical harmonic oscillator, monitor the optical interference phase shift value and reflection intensity change signal within the period of each node, and combine them into a time-series signal sequence according to the acquisition order to obtain the temperature interference signal sequence structure.

[0007] S2: Based on the temperature interference signal sequence structure, read the interference phase shift value and reflection intensity change rate of adjacent nodes under the same period, perform slope ratio comparison, filter and judge the strain leading region, and record the corresponding node number and polar coordinate position to obtain the temperature strain dominant distribution map.

[0008] S3: Based on the node number and position in the dominant temperature strain distribution map, retrieve the corresponding communication path records within three periods, filter continuous path segments, construct the compensation signal forwarding signal response structure path according to the connection direction relationship between nodes, and generate a temperature compensation chain path template.

[0009] S4: Based on the signal response paths of each node in the temperature compensation chain path template, extract the interference phase data in the latest period, construct a continuous compensation channel, and generate an interference sequence temperature compensation channel set;

[0010] S5: Based on the set of temperature compensation channels of the interference sequence, map the corresponding nodes to the resonator structure to construct driving commands that suppress the resonance shift caused by temperature disturbances, and generate an adaptive temperature compensation response scheme.

[0011] As a further aspect of the present invention, the strain leading region specifically refers to the node pairs where the interference phase change is less than the reflection intensity increase.

[0012] As a further aspect of the present invention, the continuous path segment specifically refers to a path segment in which the difference in signal response fluctuation caused by temperature disturbance is less than a preset comparison amplitude.

[0013] As a further aspect of the present invention, during the construction of the continuous compensation channel, the phase difference between adjacent nodes is compared according to the order of nodes in the path, and the sorting order of the nodes is adjusted.

[0014] As a further aspect of the present invention, the temperature interference signal sequence structure includes phase perturbation characteristic indicators, reflectivity change characteristic indicators, and temporal node numbering relationships; the temperature strain dominant distribution map includes dominant node numbers, polar coordinate position labels, and strain dominant direction identifiers; the temperature compensation chain path template includes path segment sequence identifiers, connected direction chains, and signal forwarding control structures; the interference sequence temperature compensation channel set includes a phase difference sorting structure, a continuous node mapping sequence, and temperature perturbation compensation channel numbers; and the adaptive temperature compensation response scheme includes a control response position matching table, an interference feature mapping index, and driving command generation parameters.

[0015] As a further aspect of the present invention, the step of obtaining the temperature interference signal sequence structure is as follows:

[0016] S111: Obtain the number and corresponding three-dimensional coordinates of the sapphire fiber nodes in the spherical structure of the hemispherical harmonic oscillator, monitor the optical interference phase value and reflection intensity value of each node in the sampling period, and aggregate the data according to the node number and sampling order to generate an optical interference monitoring data sequence.

[0017] S112: Based on the optical interference phase value and reflection intensity value in the optical interference monitoring data sequence, the phase change of the same node in a continuous time period is differentially processed and arranged in chronological order to generate an interference phase change trend sequence.

[0018] S113: Based on the interference phase change trend sequence, determine whether the phase change trend and reflection intensity change of each node exceed the temperature disturbance detection threshold, extract all node time-series data segments that meet the conditions and recombine them to generate a temperature interference signal sequence structure.

[0019] As a further aspect of the present invention, the step of obtaining the dominant temperature strain distribution map is as follows:

[0020] S211: Based on the temperature interference signal sequence structure, read the interference phase shift value and reflection intensity change rate of all adjacent node pairs under the same sampling period, construct node pair data frames according to the node number order, extract the phase shift value and intensity change rate in the corresponding period for each pair of nodes, and aggregate and encode the value difference between node pairs to generate adjacent node slope comparison parameter set.

[0021] S212: Based on the adjacent node slope comparison parameter set, establish slope ratio values ​​for the interference phase shift difference and the reflection intensity change rate difference for each pair of nodes. If the phase shift slope is less than the intensity growth slope, mark the corresponding node pair as the strain response leading region, aggregate all marked node pair numbers and time segment information, and generate a strain-dominant node identifier set.

[0022] S213: Based on the node pair numbering information and their corresponding polar coordinate positions in the strain-dominant node identifier set, map all strain-dominant region identifier values ​​to the spherical coordinate graph structure in node order to construct a two-dimensional projection map and establish a temperature strain-dominant distribution map.

[0023] As a further aspect of the present invention, the step of obtaining the temperature compensation chain path template is as follows:

[0024] S311: Based on the node number and polar coordinate position recorded in the temperature strain dominant distribution map, retrieve the master controller communication path record data frames associated with each node in the most recent three sampling periods, extract all path segments of the corresponding node, and perform segment aggregation based on the timestamp field in the path to generate a periodic path segment sequence set.

[0025] S312: Extract the signal response intensity fluctuation data of each continuous path segment in the periodic path segment sequence set, calculate the maximum amplitude difference of three consecutive periods according to the node path order, and judge it with the preset comparison amplitude threshold. Select the path segments that meet the condition that the fluctuation difference is less than the preset comparison amplitude threshold, mark them as stable paths, organize and integrate them into a structured array according to the node number, and generate a set of stable path segments.

[0026] S313: Based on the node numbering order and the connection direction relationship between nodes in the stable path segment set, according to the downlink direction of the master controller in the network topology, perform directional consistency verification on each path segment, filter out path segments with inconsistent directions, and then splice them into a complete path sequence according to the communication direction to establish a temperature-compensated chain path template.

[0027] As a further aspect of the present invention, the step of obtaining the temperature compensation channel set of the interference sequence is as follows:

[0028] S411: Based on the path structure recorded in the temperature compensation chain path template, extract the interference phase value of each node in each path under the latest sampling period, retrieve the phase data of each node in all paths in the original order, and bind the node number and phase value to form a node phase pair list to establish a path node phase sequence group.

[0029] S412: Based on the phase values ​​corresponding to the adjacent numbers of each group of nodes in the path node phase sequence group, the interference phase values ​​between each two adjacent nodes are calculated by difference according to the node path order, and the absolute values ​​of the difference are sorted in ascending order. The node order in the original path is readjusted according to the sorting order to generate the minimum phase difference path sequence.

[0030] S413: Based on the adjusted node connection order in the minimum phase difference path sequence, connect the nodes sequentially according to their indexes to establish equally spaced path channels. Perform continuity judgment on the numbering and spatial distance between adjacent nodes, eliminate non-continuous segments, and establish an interference sequence temperature compensation channel set.

[0031] As a further aspect of the present invention, the step of obtaining the adaptive temperature compensation response scheme is as follows:

[0032] S511: Based on the node number sequence corresponding to each channel in the interference sequence temperature compensation channel set, retrieve the response region index of each channel termination node in the master controller mapping table, match the channel tail node with the master controller response table according to the path structure, and extract the corresponding control response coordinate field and pointing mark field to establish a channel control response mapping set.

[0033] S512: Based on the node index coordinates in the channel control response mapping set and the original harmonic oscillator structure coordinate index table, the node phase feature sequence in each control channel is projected onto the polar coordinate map of the harmonic oscillator spherical structure, and the spherical local interference feature region is reconstructed by interpolation according to the node number order to generate the node interference feature mapping structure.

[0034] S513: Based on the phase feature density and control response coordinates of each node in the node interference feature mapping structure, write all control response positions into the main controller input queue, and execute the main controller response instruction generation process in coordinate order to establish an adaptive temperature compensation response scheme.

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

[0036] In this invention, a time-series interference signal sequence oriented towards the three-dimensional structural distribution of nodes is constructed. Based on the difference between interference phase and reflectivity changes, strain-dominant regions are identified. A chain path is constructed by fusing polar coordinate positions and node connectivity. Low-interference path segments in signal propagation are dynamically extracted. The continuous phase difference of node responses is obtained and sorted to generate a channel set. Multi-path driving command generation is achieved by combining the mapping of interference features. This overcomes the limitations of a single correction channel and enables multi-point dynamic sensing and signal reconstruction in a three-dimensional thermal distribution field. At the same time, the compensation path is guided to complete chain-like control of thermal disturbance direction and stress gradient, effectively suppressing the nonlinear influence of thermoelastic coupling on the resonant frequency and improving the decoupling capability and response accuracy of temperature compensation in unsteady thermal fields. Attached Figure Description

[0037] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0038] Figure 1 This is a schematic diagram of the steps of the present invention;

[0039] Figure 2 This is a flowchart of the process for obtaining the temperature interference signal sequence structure of the present invention;

[0040] Figure 3 This is a flowchart of the process for obtaining the dominant temperature strain distribution map of the present invention;

[0041] Figure 4 This is a flowchart of the temperature compensation chain path template acquisition process of the present invention;

[0042] Figure 5 This is a flowchart of the process for obtaining the temperature compensation channel set of the interference sequence in this invention;

[0043] Figure 6 This is a flowchart of the adaptive temperature compensation response scheme of the present invention. Detailed Implementation

[0044] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0045] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0046] In the embodiments of this invention, the terms "image" and "picture" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, they convey the same meaning. Similarly, the terms "of," "corresponding (relevant)," and "corresponding" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, they convey the same meaning.

[0047] In this embodiment of the invention, sometimes a subscript such as W1 may be written in a non-subscript form such as W1. When the difference is not emphasized, the meaning they express is the same.

[0048] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0049] Please see Figure 1 This invention provides an adaptive temperature compensation method based on a hemispherical harmonic oscillator, comprising the following steps:

[0050] S1: Obtain the three-dimensional coordinates and numbers of the sapphire fiber nodes within the spherical structure of the hemispherical harmonic oscillator, monitor the optical interference phase shift value and reflection intensity change signal caused by temperature disturbance within each node period, and combine them into a time-series signal sequence according to the acquisition order to obtain the temperature interference signal sequence structure.

[0051] S2: Based on the temperature interference signal sequence structure, read the interference phase shift value and reflection intensity change rate of adjacent nodes under the same period, perform slope ratio comparison operation, filter out node pairs with interference phase change less than reflection intensity increase, judge them as strain advance region under temperature disturbance, and record the corresponding node number and polar coordinate position to obtain temperature strain dominant distribution map;

[0052] S3: Based on the node number and position in the temperature strain dominant distribution map, retrieve the corresponding master controller communication path records within three cycles, filter the continuous path segments where the difference in signal response fluctuation caused by temperature disturbance is less than the preset comparison amplitude, and construct the signal response structure path for compensating signal forwarding according to the connection direction relationship between nodes, and generate a temperature compensation chain path template.

[0053] S4: Based on the signal response path of each node in the temperature compensation chain path template, extract the interference phase data in the latest period, compare the phase difference between adjacent nodes according to the node order in the path, adjust the sorting order of the nodes, and construct a continuous compensation channel under temperature disturbance to generate an interference sequence temperature compensation channel set.

[0054] S5: Based on the temperature compensation channel set of the interference sequence, determine the control response position pointed to by each channel, map the corresponding node interference feature sequence to the resonator structure, trigger the master controller to construct the driving command to suppress the resonance offset caused by temperature disturbance, and generate an adaptive temperature compensation response scheme.

[0055] The temperature interference signal sequence structure includes phase perturbation characteristic indicators, reflectivity change characteristic indicators, and temporal node numbering relationships. The dominant temperature strain distribution map includes dominant node numbers, polar coordinate position labels, and strain dominant direction identifiers. The temperature compensation chain path template includes path segment sequence identifiers, connected direction chains, and signal forwarding control structures. The temperature compensation channel set of the interference sequence includes a phase difference sorting structure, a continuous node mapping sequence, and temperature perturbation compensation channel numbers. The adaptive temperature compensation response scheme includes a control response position matching table, an interference feature mapping index, and driving command generation parameters.

[0056] Please see Figure 2 The specific steps of S1 are as follows:

[0057] S111: Obtain the number and corresponding three-dimensional coordinates of the sapphire fiber nodes in the spherical structure of the hemispherical harmonic oscillator, monitor the optical interference phase value and reflection intensity value of each node in the sampling period, and aggregate the data according to the node number and sampling order to generate an optical interference monitoring data sequence.

[0058] To obtain the node numbers and corresponding three-dimensional coordinates of the sapphire fiber nodes within the spherical structure of a hemispherical harmonic oscillator, a point-by-point scan of the oscillator surface is required using a spatial positioning device. Each node is represented in spherical coordinates, which are converted to Cartesian coordinates using x = rsinθcosΦ, y = rsinθsinΦ, and z = rcosθ. Here, r = 30 mm, and θ and Φ represent the polar angle and azimuth angle, respectively. Let the angle values ​​of a given node be θ = 60° and Φ = 45°.

[0059] Therefore, x = 30 × sin60° × cos45° ≈ 18.37 mm, and so on, the three-dimensional coordinates of the node can be obtained. During the sampling process, a 1550 nm wavelength laser is injected into the sapphire fiber through the light source coupling module at a period of 0.5 seconds. The intensity of the reflected light signal and the interference phase value are recorded. The reflection intensity value (a measure of the relative intensity of the reflected light signal, obtained by dividing the node's reflected light power value by a normalization factor) is recorded as a ratio between 0 and 1. The phase value is in radians, ranging from -π to +π. After the signal value of each node is collected, it is aggregated by the node number and timestamp to construct a structured data frame. For example, if the phase values ​​of node 1 are 0.85 rad and 1.15 rad at t = 0.5 s and 1.0 s respectively, the phase change is 0.30 rad. This value and the reflection intensity are simultaneously bound to the sequence array. After sampling, the nodes are sorted from 1 to n according to their numbers, and then a complete sequence is constructed in chronological order.

[0060] For example, the node sequence is: P1=[0.85,1.15,0.92], R1=[0.70,0.74,0.69], which represents the phase and reflectance sequence of the first node at three time points (P1 represents the sequence of interference phase values ​​collected by the first monitoring node at multiple time points, and R1 represents the sequence of reflection intensity values ​​collected by the first monitoring node at the same time point). All nodes are processed sequentially, and the aggregated data forms a two-dimensional time-node structure data matrix, which is then recorded as the input to the main system. The number of samples can be referenced from the typical data given in Table 1.

[0061] Table 1. Optical signal sampling table of monitoring nodes

[0062] Node number Sampling time (s) Phase value (rad) Reflection intensity (ratio) 1 0.5 0.85 0.70 1 1.0 1.15 0.74 2 0.5 -0.20 0.62 2 1.0 0.10 0.65 3 0.5 2.50 0.78 3 1.0 2.60 0.80

[0063] As shown in Table 1, the phase value of node 1 changes from 0.85 rad to 1.15 rad, a change of 0.30 rad, corresponding to a change in reflection intensity of 0.04 (according to Table 1, 0.74 - 0.70 = 0.04). Data sequences for all nodes are extracted and constructed in this manner for subsequent trend analysis, ultimately generating an optical interferometry monitoring data sequence.

[0064] S112: Based on the optical interference phase value and reflection intensity value in the optical interference monitoring data sequence, the phase change of the same node in a continuous time period is differentially processed and arranged in chronological order to generate an interference phase change trend sequence.

[0065] The phase and reflection intensity values ​​are read from the optical interferometry monitoring data sequence. The phase change trend of each node at continuous sampling times is calculated using a differential processing method, i.e., subtracting the phase values ​​between two adjacent times. For example, if the phase of node 2 at 0.5s and 1.0s is -0.20rad and 0.10rad respectively, the change is 0.10rad - (-0.20rad) = 0.30rad. This difference is the phase shift of the node during that time period. This operation is performed on all nodes, and the differences are arranged in chronological order to form a vector sequence. For example, if the phase of node 2 at three sampling times (0.5s, 1.0s, 1.5s) is -0.20rad, 0.10rad, and 0.25rad respectively, the differential sequence is [0.30, 0.15]. The reflection intensity is similarly differentially processed to record its change trend and serve as an auxiliary analysis parameter. If the sampling time interval is 0.5s, the average rate of change can be calculated by dividing the difference value by the time interval. That is, the first rate of change for node 2 is 0.30rad / 0.5s = 0.60rad / s. The trend vectors of all nodes are organized into a two-dimensional matrix using node numbers. Each row of the matrix represents the phase change trend sequence of a node, and the columns correspond to time periods. For example, in the node sequence matrix:

[0066]

[0067] This represents the phase change trend of three nodes over two time periods. This matrix structure can be used for subsequent temperature perturbation judgment and trend continuity analysis, ultimately generating an interferometric phase change trend sequence.

[0068] S113: Based on the interference phase change trend sequence, determine whether the phase change trend and reflection intensity change of each node exceed the temperature disturbance detection threshold, extract all node time-series data segments that meet the conditions and recombine them to generate a temperature interference signal sequence structure.

[0069] Based on the interference phase change trend sequence, it is determined whether the phase change direction of a node remains consistent within a continuous time period. Combined with the change in reflection intensity, it is checked whether it exceeds the temperature disturbance detection threshold. The disturbance threshold is set at 0.25 rad. If two consecutive phase differences of a node are 0.30 rad and 0.32 rad, respectively, and both are positive, it indicates that the trend direction is consistent and continuously exceeds the set threshold. Therefore, it is determined that the node is affected by temperature disturbance within that time period. To quantify continuity, consistent direction is defined as the difference having the same sign, and continuous values ​​are defined as the absolute value of the difference being greater than 0.25 rad. Regarding reflection intensity, a fluctuation threshold of 0.03 is set. If the change exceeds this value, it is considered an intensity anomaly. Taking node 2 as an example, assuming its phase difference is [0.30, 0.32] and reflectivity difference is [0.04, 0.05], both exceeding the threshold and increasing in the same direction, this node is marked as a disturbance node. Nodes meeting this condition and their continuous time period data are selected from all nodes. Their numbers and corresponding time periods are extracted, and time-series aggregation is performed. If the time period during which node 2 satisfies the perturbation condition is from 1.0s to 2.0s, then its number and the sequence segment are output as a set of structures. By combining all segments that satisfy the condition in this way, a complete data structure is formed, and finally, a temperature interference signal sequence structure is generated.

[0070] Please see Figure 3 The specific steps of S2 are as follows:

[0071] S211: Based on the temperature interference signal sequence structure, read the interference phase shift value and reflection intensity change rate of all adjacent node pairs under the same sampling period, construct node pair data frames according to the node number order, extract the phase shift value and intensity change rate in the corresponding period for each pair of nodes, and aggregate and encode the value difference between the node pairs to generate the adjacent node slope comparison parameter set.

[0072] Based on the temperature interference signal sequence structure, the interference phase shift value and reflection intensity change rate under the same period are extracted from each pair of adjacent nodes. First, the order of the adjacent node numbers needs to be retrieved. Let the 5th and 6th nodes be considered a pair of adjacent nodes. Then, their phase shift value and reflectivity change rate are read at the period t1 = 1.0s. If the phase shift of the 5th node is 0.22 rad and the intensity change rate is 0.05, and the phase shift of the 6th node is 0.28 rad and the intensity change rate is 0.08, then the corresponding difference pair for this node pair is a phase shift difference of 0.06 rad and an intensity change rate difference of 0.03. This difference pair is used as the basis for comparison. When performing this extraction operation, all consecutively numbered node pairs need to be traversed. In each pair, the phase shift value and reflectivity change rate data of each node within the current period are read. The data storage structure adopts a triplet format, recording the node number, phase shift value, and reflectivity change rate respectively. Multiple structured triplet data are then combined into a node pair data frame. During the combination process, the period timestamp is used as the primary key to match the data of adjacent nodes, thereby constructing a node pair data matrix within the period, and organizing and storing it row by row according to the node number order. Based on the above, the differences corresponding to each node pair are extracted and aggregation is performed. The aggregation method is to sequentially concatenate the difference groups within the period to form a two-dimensional structure, where each row represents the difference data of a node pair. To ensure the uniformity of this structure, the node number order must be strictly incremental, the time period must be synchronized, and the sampling frequency must be consistent. For example, if the data of nodes 6 and 7 in period t1 are 0.35 rad, 0.09 and 0.31 rad, 0.06 respectively, then the difference is 0.04 rad, 0.03, which is taken as the data item of node pair (6, 7) in period t1. After aggregating the corresponding differences of all node pairs in all periods, a standardized data set is formed, as shown in Table 2:

[0073] Table 2. Data on the periodic difference between adjacent nodes

[0074] Node pairs Periodicity (s) Phase offset difference (rad) Difference in intensity （5，6） 1.0 0.06 0.03 （6，7） 1.0 0.04 0.03 （7，8） 1.0 0.09 0.05

[0075] As shown in Table 2, the phase offset difference of node pair (5, 6) is 0.06 rad, and the intensity change difference is 0.03 rad. This set of differences will be used as the parameter input for subsequent slope ratio comparison, and finally generate the slope comparison parameter set of adjacent nodes.

[0076] S212: Based on the slope comparison parameter set of adjacent nodes, establish slope ratio values ​​for the difference in interference phase shift and the difference in reflection intensity change rate for each pair of nodes. If the phase shift slope is less than the intensity growth slope, mark the corresponding node pair as the strain response leading region, aggregate all marked node pair numbers and time segment information, and generate a strain-dominant node identifier set.

[0077] The slope comparison parameter set of adjacent nodes was extracted. The ratio of the phase shift difference and the reflection intensity change difference in each pair of nodes was calculated and compared with a set reference slope ratio threshold. The threshold was set to 1.0. When the phase shift difference of node pair (5, 6) was 0.06 rad and the intensity change difference was 0.03 rad, the slope ratio was 0.06 / 0.03 = 2.0, which is greater than the threshold of 1.0, thus not meeting the screening condition. However, for node pair (6, 7), the phase shift difference was 0.04 rad and the intensity change difference was 0.06 rad, so the slope ratio was 0.04 / 0.06 ≈ 0.667, meeting the condition. This threshold was set based on experimental correspondence. When the reflectivity change rate was significantly higher than the phase shift rate and the slope ratio was less than 1.0, it was considered that temperature perturbation was dominant. This process requires calculating the slope ratio of all node pairs one by one, and setting the judgment expression as follows: if the slope ratio is less than the reference slope ratio threshold, then it is marked as a node pair in the strain leading region. All node pairs that meet this logical judgment are recorded in a structured manner with their corresponding numbers and sampling period information, organized into a two-dimensional array according to the period. Each row of the array is a node pair and its corresponding period. In the example, node pair (6, 7) meets the judgment condition at a period of 1.0s, and it is written into the result set as an effective strain-dominant segment. The structure record is (6, 7, 1.0). After aggregating all node pairs that meet the condition, the strain-dominant node identifier set is obtained.

[0078] S213: Based on the node pair numbering information and their corresponding polar coordinate positions in the strain-dominant node identifier set, map all strain-dominant region identifier values ​​to the spherical coordinate map structure in node order to construct a two-dimensional projection map and establish a temperature strain-dominant distribution map.

[0079] Based on the strain-dominant node identifier set, the numbering information of each node pair is extracted, and the polar coordinates of each node are retrieved from the system node spatial positioning database. Let the polar coordinates of nodes 6 and 7 be respectively... The corresponding spherical region is a local area near the equator of the sphere. After recording its number and corresponding polar coordinates, it is matched and mapped onto the spherical distribution map. The mapping method is as follows: read the polar angle values ​​of the nodes to the two endpoints in numerical order, construct local patch regions using trigonometric interpolation, and represent the perturbation level of the region in the spherical distribution map with color intensity. The perturbation level can be linearly mapped according to the inverse of the slope ratio, for example, a ratio of 0.667 corresponds to a perturbation level of 1.5. Then, all mapped regions are projected onto the two-dimensional unfolded map of the spherical map to construct an image matrix. Each element of the image matrix represents the perturbation level at a spherical coordinate point. The matrix is ​​assembled into a spherical image frame structure according to the point matrix sequence, and the output is a high-resolution map file for visualization. The polar coordinate information and perturbation level values ​​of all strain-dominant regions are embedded as pixel metadata into the image structure, and finally, a temperature-strain dominant distribution map is established.

[0080] Please see Figure 4 The specific steps of S3 are as follows:

[0081] S311: Based on the node number and polar coordinate position recorded in the temperature strain dominant distribution map, retrieve the master controller communication path record data frames associated with each node in the most recent three sampling periods, extract all path segments of the corresponding node, and perform segment aggregation based on the timestamp field in the path to generate a periodic path segment sequence set.

[0082] Based on the node numbers and polar coordinate positions recorded in the temperature strain dominant distribution map, the master controller communication path records for each node within three consecutive sampling periods are retrieved sequentially. The master controller communication path records are constructed using timestamps and communication hop counts as basic fields. Each path record is expressed in the form of a structure containing node number, start time, path sequence, and signal response value. When performing path retrieval, the set of node numbers extracted from the map is traversed first, and data frames with the same number and corresponding sampling period range are called from the path database. Assuming the node number is 7 and the retrieval period is from t0=1.0s to t2=3.0s, the corresponding path segments are 7, 1.0s, [7→9→11], 0.72; 7, 2.0s, [7→9→11], 0.74; and 7, 3.0s, [7→9→11], 0.76. The acquisition method judges the consistency of the starting point number and path by using the path matching field. At the same time, the path segment records of different periods are aggregated and sorted in ascending order by the time field to construct a continuous communication path structure. To ensure the continuity of path segments, it is necessary to determine whether the endpoint of each path segment is consistent with the starting point of the next one during the aggregation process. If they are consistent, they are marked as periodically continuous segments; otherwise, they are discarded. Path segments with inconsistent path node structures or sampling time jumps must be excluded to avoid interfering with chain recognition. Successfully aggregated path segments will be organized into a three-dimensional array, with dimensions of node number, period number, and path structure, ultimately establishing a structured set of path segment sequences. Table 3 shows an example of a three-period path segment for a typical node.

[0083] Table 3. Main Controller Three-Cycle Path Segment Table

[0084] Node number Sampling period (s) Path structure Signal response value 7 1.0 7→9→11 0.72 7 2.0 7→9→11 0.74 7 3.0 7→9→11 0.76

[0085] As shown in Table 3, node 7 has a consistent path and continuous time in three cycles, and the signal response value structure is stable, which satisfies the condition of continuity of periodic path segments. This path sequence will be used as the input for subsequent response fluctuation judgment, and finally a set of periodic path segment sequences will be generated.

[0086] S312: Extract the signal response intensity fluctuation data of each continuous path segment in the periodic path segment sequence set, calculate the maximum amplitude difference of three consecutive periods according to the node path order, and compare it with the preset comparison amplitude threshold. Select the path segments that meet the condition that the fluctuation difference is less than the preset comparison amplitude threshold, mark them as stable paths, organize and integrate them into a structured array according to the node number, and generate a set of stable path segments.

[0087] The signal response value sequences of each path segment in the periodic path segment sequence set are extracted. The change in response intensity of each path over three consecutive periods is compared. The signal response value sequences of the path segments are set to 0.72, 0.74, and 0.76. The difference between the maximum and minimum values ​​is calculated as ∆ = 0.76 - 0.72 = 0.04. This difference is compared with a preset comparison threshold of 0.05. Based on experimental statistical analysis, if the difference in response value fluctuation is less than 0.05, the path is considered stable over the consecutive periods. The threshold discrimination logic expression is used in the judgment process: if ∆ < Y, where Y = 0.05, the path is marked as stable. All path segments that meet this condition are numbered, recorded, and labeled. The label structure is node number, start time, path structure, and stability indicator. For example, the path segment at node 7 changes to 0.04 over three periods, satisfying the stability condition. It is recorded as 7, 1.0, 7→9→11, 1, indicating that the path segment started at 1.0s and satisfied the stability condition after three periods. All records meeting the criteria are grouped and organized into a structured array by node number. The array structure supports subsequent path concatenation. Record precision is retained to two decimal places. The fluctuation calculation process is based on the absolute difference value, and minute jitters with an average error of less than 0.005 are removed. The stability indicator is a Boolean code, 1 for stable and 0 for unstable, which facilitates subsequent conditional filtering. Finally, a set of stable path segments is generated.

[0088] S313: Based on the node numbering order and the connection direction relationship between nodes in the stable path segment set, according to the downlink direction of the master controller in the network topology, the direction consistency of each path segment is checked. After filtering out path segments with inconsistent directions, they are spliced ​​into a complete path sequence according to the communication direction to establish a temperature-compensated chain path template.

[0089] Based on a set of stable path segments, the starting and ending node numbers of each path segment are extracted. The connectivity direction relationships between all nodes in that path segment are retrieved from the network topology table of the master controller (i.e., the central control unit, which undertakes key tasks such as coordination, control, and data aggregation; typically an embedded control device with communication interfaces, data processing capabilities, and network scheduling logic, or an industrial computing control unit). The direction relationship is based on the downstream direction; that is, if the master controller number is less than the subordinate node number, it is a downstream path. For example, in the path segment 7→9→11, the node numbers monotonically increase, conforming to the downstream rule. During connectivity direction verification, each pair of adjacent nodes in the path is traversed, and a number comparison operation is performed on each pair of nodes. If the current node number is less than the next node number, it is a valid downstream connection; otherwise, the path segment is not included in the chained path group. The judgment expression is: For all If the path is valid, it is a path with consistent direction, and the identifier is set to 1; otherwise, it is set to 0. This is recorded as the path structure and direction verification identifier. For stable path segments with consistent direction, they are spliced ​​according to the logical order from the master controller to the terminal node. The condition for continuous splicing is that the terminal node of the path segment has the same number as the starting node of the next segment. If the subsequent segment of path segment 7→9→11 is 11→13→15, then it is spliced ​​as 7→9→11→13→15, forming a complete communication path structure. The spliced ​​chain structure is output in sequence form. Each chain path structure contains a sequence of path node numbers, path length, and direction consistency identifier. After splicing, this path structure is embedded into the signal compensation strategy template as a key parameter field, and finally, a temperature compensation chain path template is established.

[0090] Please see Figure 5 The specific steps of S4 are as follows:

[0091] S411: Based on the path structure recorded in the temperature-compensated chain path template, extract the interference phase value of each node in each path under the latest sampling period, retrieve the phase data of each node in all paths in the original order, and bind the node number and phase value to form a node phase pair list to establish a path node phase sequence group.

[0092] Based on the path structure recorded in the temperature-compensated chain path template, interference phase data of all nodes in the latest sampling period are extracted from the phase data acquisition module according to the node numbering order in each chain path. The sampling period is uniformly set to t=3.0s based on the system timestamp. Each node corresponds to a unique phase value, and the sampling data precision is set to 0.01rad. The extraction method is to directly match the period time after locating the node by index, and construct an array of binding relationships between node numbers and their phase values. Suppose a path is 7→9→11→13, and the corresponding node phase values ​​are 1.12rad, 1.08rad, 1.05rad, and 1.02rad respectively. After binding, the node pair sequence (7, 1.12), (9, 1.08), (11, 1.05), (13, 1.02) is constructed. This binding operation is performed independently for each path, constructing a mapping structure between the path and the node phase values. After all paths have been processed, the aforementioned node phase pairs are organized into a two-dimensional matrix. Each row represents a path, and the columns represent the path node numbers and their corresponding phase values. This structure is used for subsequent phase difference comparisons and node order reordering operations. For ease of access, a hash index structure is built for each group of path node phase data according to the path number, supporting path-level access operations, ultimately resulting in a group of path node phase sequences.

[0093] S412: Based on the phase values ​​corresponding to the adjacent numbers of each group of nodes in the path node phase sequence group, calculate the difference between the interference phase values ​​between each pair of adjacent nodes according to the node path order, sort the absolute values ​​of the difference in ascending order, and readjust the order of the nodes in the original path according to the sorting order to generate the minimum phase difference path sequence.

[0094] Based on the path node information in the path node phase sequence group, the interference phase values ​​between adjacent node pairs are extracted in node order. A phase difference calculation is performed by subtracting the phase value of the subsequent node from the current node's phase value. For example, if the path is 7, 9, 11, 13, with corresponding phase values ​​of 1.12 rad, 1.08 rad, 1.05 rad, and 1.02 rad, the phase differences are ∆1 = 1.12 - 1.08 = 0.04, ∆2 = 1.08 - 1.05 = 0.03, and ∆3 = 1.05 - 1.02 = 0.03, respectively. These three phase differences are combined into a difference sequence of 0.04 rad, 0.03 rad, and 0.03 rad. This sequence is then sorted in ascending order of absolute value to obtain the node connection priority corresponding to the new sequence index. The node order is then reorganized according to the priority, connecting node pairs with the smallest phase differences first to ensure a smooth phase transition in the path under temperature perturbations. The reordering rule is as follows: each time, select the unconnected node pair with the smallest phase difference and connect them to the new path sequence until all nodes form a closed link. Based on the example path, since the smallest phase difference is 0.03, the corresponding node pairs are 9, 11 and 11, 13. These are first combined and then connected to the remaining nodes to form a new sequence 9→11→13→7. This new sequence replaces the original path node order, and the new sorted structure is saved as the path reordering result, ultimately generating the path sequence with the smallest phase difference.

[0095] S413: Based on the adjusted node connection order in the minimum phase difference path sequence, connect the nodes sequentially according to their indexes to establish equally spaced path channels. Perform continuity judgment on the numbering and spatial distance between adjacent nodes, eliminate non-continuous segments, and establish an interference sequence temperature compensation channel set.

[0096] Based on the arrangement order of new nodes corresponding to each path in the minimum phase difference path sequence, the spatial coordinate information of each node is extracted sequentially according to its node number. The spatial distance value between each pair of adjacent nodes is obtained from the node topology table. The distance calculation is completed according to the Euclidean distance formula for three-dimensional coordinates between nodes. If a certain path is 9→11→13→7, the corresponding coordinates are... The distance between each pair of nodes is calculated, and a continuity threshold of 3 mm is set. If the spatial distance between a pair of nodes exceeds this threshold, it is considered a discontinuous connection segment and is discarded. In the example, if the distance between nodes 11 and 13 is 2.4 mm and the distance between 13 and 7 is 3.2 mm, the last segment is discarded, and 9→11→13 is retained. Subsequently, all retained segments are numbered and labeled, and a new channel connection structure is constructed. The channel structure is a sequential node sequence with node coordinates and phase value information, and is uniformly archived according to path number, forming a set of paths with continuity and phase stability, thus establishing an interference sequence temperature compensation channel set.

[0097] Please see Figure 6 The specific steps of S5 are as follows:

[0098] S511: Based on the node number sequence corresponding to each channel in the temperature compensation channel set of the interference sequence, retrieve the response region index of each channel termination node in the master controller mapping table, match the channel tail node with the master controller response table according to the path structure, and extract the corresponding control response coordinate field and pointing mark field to establish a channel control response mapping set.

[0099] Based on the path structure of each channel in the temperature compensation channel set of the interference sequence, the termination node number is extracted, and the control response record of that node in the master controller mapping table is retrieved. The mapping table uses the node number as the primary key index field and the control response coordinates and device ID as data fields. The retrieval method is to match the channel termination node number with the primary key field of the response table to obtain the control response position corresponding to each channel. For example, in the channel path 9→11→13→15, the termination node is 15. If the control coordinates corresponding to node 15 in the response mapping table are x=18.0, y=22.0, z=14.5, then this control response position is the control action point of that channel. After traversing all channels, the path number and control response coordinates are bound to form a binary dataset with the structure (path ID...). i Control coordinates i Simultaneously, the type number of the control device to which the response point belongs is recorded for easy retrieval during subsequent command generation. After the dataset is formed, it is uniformly sorted according to the path number, invalid channels without corresponding control response points are removed, and missing fields are marked. All structured response point data will be organized into a list structure, and an index will be built to accelerate the retrieval of response coordinate information by path number, ultimately generating a channel control response mapping set.

[0100] S512: Based on the positional correspondence between the node index coordinates in the channel control response mapping set and the original harmonic oscillator structure coordinate index table, the node phase feature sequence in each control channel is projected onto the polar coordinate map of the harmonic oscillator spherical structure, and the spherical local interference feature region is reconstructed by interpolation according to the node number order, generating the node interference feature mapping structure.

[0101] The process reads the path number and corresponding response coordinates of each control channel in the channel control response mapping set. Then, it sequentially refers to the three-dimensional polar coordinate index table of the resonator structure to complete the spatial mapping between the control channel nodes and the resonator structure nodes. The polar coordinate index table uses the node number and its spherical coordinates (θ, Φ) as its main fields. After retrieving all node numbers in each control path, its phase feature values ​​and corresponding polar coordinates are extracted sequentially. The node phase features are then distributed in a spherical map using a polar coordinate interpolation function to form a continuous surface. For example, if the polar coordinates of the nodes in paths 9→11→13→15 are (60°, 45°), (65°, 50°), (70°, 55°), and (75°, 60°), with corresponding phase values ​​of 1.02 rad, 1.08 rad, 1.12 rad, and 1.18 rad, then the interpolation generates a spherical region mapping image and generates a matrix structure M. i,j This is used to represent the phase feature density corresponding to the polar coordinate position. To unify the map construction method, the interpolation operation adopts two-dimensional polar coordinate sequence regularization processing, constructs a node interpolation array with a step size of 5°, removes null value regions, and finally constructs a corresponding spherical region feature map segment for each path. After all map segments are stitched together, a feature space mapping record between nodes is established. The above structure is archived according to channel number, and a structure format record field is constructed, including: path number, node list, phase sequence, mapped polar coordinate sequence, and interpolation matrix, finally generating a node interference feature mapping structure.

[0102] S513: Based on the phase feature density and control response coordinates of each node in the node interference feature mapping structure, write all control response positions into the main controller input queue, and execute the main controller response instruction generation process in the order of coordinates to establish an adaptive temperature compensation response scheme.

[0103] Based on the interpolation matrix structure corresponding to each control channel in the node interference feature mapping structure, the phase density value of each node at its polar coordinate position is extracted and the difference is calculated with the interference disturbance detection threshold. Let the threshold be 0.03 rad. If the phase density value of a node in the mapping is 1.12 rad, the difference is 0.07 rad, which is greater than the threshold judgment criterion of 0.03, and thus marked as a trigger region. After judging all nodes in each channel one by one, the control response coordinates corresponding to all nodes that meet the threshold triggering condition are written into the master controller response input buffer. The writing structure is (path ID...). i Control coordinates iThe controller reads the buffer data and records the timestamp and channel number to form a response sequence. After reading the buffer data, the controller parses the control coordinate information into device ID and command parameters. It then sets a compensation command intensity factor based on the density difference. The factor calculation rule is α = (density value - threshold) × λ, where λ = 10. In the example above, the factor is α = 0.07 × 10 = 0.7. This factor is then bound to the control command parameter field. The command structure consists of device ID, action coordinates, and compensation factor. All commands are aggregated into a command queue according to the channel number. The set of structures is stored sequentially in the controller's command buffer, ultimately establishing an adaptive temperature compensation response scheme.

[0104] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. An adaptive temperature compensation method based on hemispherical resonator, characterized in that, Comprise the following steps: S1: Obtain the three-dimensional coordinates and numbers of sapphire optical fiber nodes in the hemispherical harmonic oscillator spherical structure, monitor the optical interference phase shift value and reflection intensity change signal in each node cycle, combine into a time sequence signal sequence according to the collection sequence, and obtain a temperature interference signal sequence structure; S2: Based on the temperature interference signal sequence structure, read the interference phase shift value and reflection intensity change rate of adjacent nodes under the same cycle, perform slope ratio comparison, screen and judge the strain leading area, record the corresponding node number and polar coordinate position, and obtain a temperature strain dominant distribution map; S3: According to the node number and position in the temperature strain dominant distribution map, search the corresponding three-cycle communication path record, screen the continuous path segment, construct the compensation signal forwarding signal response structure path according to the connection direction relationship between nodes, and generate a temperature compensation chain path template; S4: Based on the signal response path of each node in the temperature compensation chain path template, extract the interference phase data in the corresponding latest cycle, construct a continuous compensation channel, and generate an interference sequence temperature compensation channel set; S5: According to the interference sequence temperature compensation channel set, map the corresponding node to the harmonic oscillator structure to construct a driving instruction for inhibiting the resonance shift caused by temperature disturbance, and generate an adaptive temperature compensation response scheme.

2. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that: The strain leading area specifically refers to the node pair whose interference phase change is less than the reflection intensity growth.

3. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that: The continuous path segment specifically refers to the path segment whose signal response fluctuation difference caused by temperature disturbance is less than the preset comparison amplitude.

4. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that: In the process of constructing the continuous compensation channel, the phase difference between adjacent nodes is compared according to the node order in the path, and the order of the nodes is adjusted.

5. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that: The temperature interference signal sequence structure includes phase disturbance characteristic index, reflectivity change characteristic index, and time sequence node number relationship, the temperature strain dominant distribution map includes dominant node number, polar coordinate position label, and strain dominant direction identifier, the temperature compensation chain path template includes path segment sequence identifier, connection direction chain, and signal forwarding control structure, the interference sequence temperature compensation channel set includes phase difference sorting structure, continuous node mapping sequence, and temperature disturbance compensation channel number, and the adaptive temperature compensation response scheme includes control response position matching table, interference feature mapping index, and driving instruction generation parameter.

6. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that, The acquisition step of the temperature interference signal sequence structure is: S111: Obtain the number and corresponding three-dimensional coordinates of sapphire optical fiber nodes in the hemispherical harmonic oscillator spherical structure, monitor the optical interference phase value and reflection intensity value of each node in the sampling cycle, aggregate the data according to the node number and sampling sequence, and generate an optical interference monitoring data sequence; S112: Based on the optical interference phase value and reflection intensity value in the optical interference monitoring data sequence, difference process the phase change amount of the same node in the continuous time period, and arrange in time sequence, generate an interference phase change trend sequence; S113: According to the interference phase change trend sequence, it is judged whether the phase change trend of each node and the reflection intensity change exceed the temperature disturbance detection threshold, all node time sequence data segments meeting the condition are extracted and recombined, and a temperature interference signal sequence structure is generated.

7. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that, The temperature strain dominant distribution map acquisition step is: S211: Based on the temperature interference signal sequence structure, the interference phase offset values and the reflection intensity change rates of all adjacent node pairs in the same sampling period are read, a node pair data frame is constructed in the node number order, the phase offset values and the intensity change rates in the corresponding period are extracted for each node pair, the value difference between the node pairs is aggregated and encoded, and a neighboring node slope ratio comparison parameter set is generated; S212: According to the neighboring node slope ratio comparison parameter set, the interference phase offset difference value and the reflection intensity change rate difference value of each node pair are respectively established as a slope ratio, if the phase offset slope is less than the intensity growth slope, the corresponding node pair is marked as a strain response leading area, and all marked node pair numbers and time segment information are aggregated to generate a strain dominant node identification set; S213: According to the node pair number information in the strain dominant node identification set and the corresponding polar coordinate position, all strain dominant area identification values are mapped into a spherical coordinate graph structure in the node order, a two-dimensional projection map is constructed, and a temperature strain dominant distribution map is established.

8. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that, The temperature compensation chain path template acquisition step is: S311: According to the node number and the polar coordinate position recorded in the temperature strain dominant distribution map, the communication path record data frame associated with each node in the last three sampling periods is retrieved, all path segments of the corresponding node are extracted, and the path segments are segmented and aggregated according to the timestamp field to generate a cycle path segment sequence set; S312: The signal response intensity fluctuation data of each continuous path segment in the cycle path segment sequence set is extracted, the maximum amplitude difference of three continuous cycles is calculated in the node path order, and a judgment is made with a preset comparison amplitude threshold, path segments meeting the fluctuation difference less than the preset comparison amplitude threshold are screened and marked as stable paths, and are organized and integrated into a structured array according to the node number to generate a stable path segment set; S313: Based on the node number order and the node interconnection direction relationship in the stable path segment set, the direction consistency of each path segment is verified according to the downward direction of the master controller in the network topology, after excluding the path segments with inconsistent directions, the complete path sequence is spliced in the communication direction in turn, and a temperature compensation chain path template is established.

9. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that, The interference sequence temperature compensation channel set acquisition step is: S411: Based on the path structure recorded in the temperature compensation chain path template, the interference phase values of each node in each path in the latest sampling period are extracted, the phase data of the nodes in all paths are retrieved one by one in the original order, and the node number and the phase value are bound as a node phase pair list to establish a path node phase sequence group; S412: Based on the adjacent number corresponding phase value of each group of nodes in the path node phase sequence group, the interference phase value between each two adjacent nodes is calculated by difference in node path order, and the absolute value of the difference is sorted in ascending order, and the node order in the original path is adjusted according to the sorting order to generate a minimum phase difference path sequence; S413: According to the adjusted node connection order in the minimum phase difference path sequence, the equal interval path channel is established by connecting the nodes in sequence according to the node index, and the continuity of the number and spatial distance between adjacent nodes is judged, and the non-continuous paragraphs are removed to establish the interference sequence temperature compensation channel set.

10. The self-adapting temperature compensation method based on hemispherical resonator according to claim 1, characterized in that, The adaptive temperature compensation response scheme acquisition step is: S511: According to the node number sequence corresponding to each channel in the interference sequence temperature compensation channel set, the response area index of each channel termination node in the host controller mapping table is searched, the tail node of the channel is numbered and matched with the host response table according to the path structure, and the corresponding control response coordinate field and pointing mark field are extracted to establish the channel control response mapping set; S512: Based on the node index coordinates in the channel control response mapping set and the original resonant substructure coordinate index table, the position corresponding to each control channel node phase feature sequence is projected into the polar coordinate map of the resonant sub-spherical structure, and the spherical local interference feature area is reconstructed by interpolation according to the node number order to generate the node interference feature mapping structure; S513: According to the phase feature density and control response coordinates of each node corresponding position in the node interference feature mapping structure, all control response positions are written into the host controller input queue, and the host controller response instruction generation process is executed according to the coordinate order to establish the adaptive temperature compensation response scheme.

Images