WO2023077760A1 - 一种面向工业无线与tsn融合的跨网时间同步方法 - Google Patents

一种面向工业无线与tsn融合的跨网时间同步方法 Download PDF

Info

Publication number
WO2023077760A1
WO2023077760A1 PCT/CN2022/091457 CN2022091457W WO2023077760A1 WO 2023077760 A1 WO2023077760 A1 WO 2023077760A1 CN 2022091457 W CN2022091457 W CN 2022091457W WO 2023077760 A1 WO2023077760 A1 WO 2023077760A1
Authority
WO
WIPO (PCT)
Prior art keywords
node
time
synchronization
clock
sending
Prior art date
Application number
PCT/CN2022/091457
Other languages
English (en)
French (fr)
Inventor
魏旻
李小云
牛爽
付清云
邢益铭
王平
Original Assignee
重庆邮电大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 重庆邮电大学 filed Critical 重庆邮电大学
Publication of WO2023077760A1 publication Critical patent/WO2023077760A1/zh

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J3/00Time-division multiplex systems
    • H04J3/02Details
    • H04J3/06Synchronising arrangements
    • H04J3/0635Clock or time synchronisation in a network
    • H04J3/0638Clock or time synchronisation among nodes; Internode synchronisation
    • H04J3/0644External master-clock
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J3/00Time-division multiplex systems
    • H04J3/02Details
    • H04J3/06Synchronising arrangements
    • H04J3/0635Clock or time synchronisation in a network
    • H04J3/0638Clock or time synchronisation among nodes; Internode synchronisation
    • H04J3/0658Clock or time synchronisation among packet nodes
    • H04J3/0661Clock or time synchronisation among packet nodes using timestamps
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Definitions

  • the invention belongs to the field of industrial Internet and relates to a cross-network time synchronization method oriented to the integration of industrial wireless and TSN.
  • Intranets such as industrial control networks and management networks in traditional manufacturing plants can be considered as closed networks that cannot be accessed by the outside world.
  • Industrial scenarios such as inconvenient wiring for old factory upgrades and long distances to production sites such as oil and gas pipelines must use wireless communication technology; more industrial interconnection applications such as machine-to-machine interaction and collaboration, factory logistics, remote maintenance of large-scale manufacturing equipment, and AR/ VR, etc. require wireless communication to meet its high-efficiency requirements and various performance indicators.
  • research on converged networks including wired and wireless technologies has a certain foundation and has begun to be widely deployed in automation systems, smart cities and home area networks.
  • wired networks have higher reliability and performance due to topology stability and smaller packet loss rate; on the other hand, in terms of spatial range mobility, energy consumption, data throughput, timeliness, reliability and other characteristics Network connectivity requirements vary, so it is necessary to combine wired and wireless network technologies to provide the required quality of service for factory networks.
  • TSN Industrial control services need to be carried by a network with definite delay and reliable guarantee.
  • TSN has attracted widespread attention because of its definite delay forwarding capability and compatibility with Ethernet protocols.
  • wired TSN cannot meet the extensive deployment requirements of massive sensors and new equipment such as AGV in smart factories.
  • the coordinated transmission of industrial wireless networks and TSN has become a research hotspot in the field of industry and automation.
  • Clock synchronization refers to the process of limiting the deviation between the clock of the device in the system and a specific time information source within a certain range. It is a particularly important technology in a distributed system. Its purpose is to ensure that any two nodes finite maximum offset between to establish a global notion of time with a predefined precision. Especially in industrial control networks, many basic operations in the network, such as object tracking, deterministic scheduling, and resource management, are inseparable from clock synchronization, which means that all nodes in the network need to share a unified time scale.
  • the purpose of the present invention is to effectively synchronize time across networks, taking heterogeneous networks as objects, and combining statistical algorithms through synchronization methods such as IEEE802.1AS, beacon frame synchronization, and time stamp-free synchronization to study heterogeneous networks.
  • the time modeling and time error compensation model between synchronization nodes in the network improve the accuracy of cross-network time synchronization.
  • the present invention provides the following technical solutions:
  • a cross-network time synchronization method for industrial wireless and TSN integration comprising the following steps:
  • TSN module of the border gateway performs clock synchronization with the TSN switch of TSN network 2 in the slave clock state;
  • the TSN module is used as the main clock of the industrial wireless module, and the industrial wireless module synchronizes the clock with the TSN module through the serial port, so as to forward the time information from the TSN network 1 to the industrial wireless network;
  • the routing device in the industrial wireless network performs clock synchronization with the slave clock status and the industrial wireless module of the border gateway through beacon frame synchronization. At the same time, the routing device acts as the master clock of the node device to synchronize the clock of the node device, thus ensuring Clock synchronization for industrial wireless networks;
  • the conversion node on the terminal side performs clock synchronization between the slave clock state and the node devices of the industrial wireless network using a time-stamp-free synchronization method, does not rely on a special synchronization protocol, and only follows the existing network data flow to complete the synchronization parameter estimation, so as to This forwards the time information from the industrial wireless network to the TSN network 1. Subsequently, the TSN switch of TSN network 1 completes clock synchronization in slave clock state and transition node at the terminal side.
  • step S1 the following steps are specifically included:
  • step S1 the delay measurement mechanism is used to measure the link delay between the master and slave clocks, and the corresponding timestamp is obtained by exchanging the delay request message Delay_rep and the delay response message Delay_resp between the slave clock and the master clock,
  • Delay_resp(t2 ) indicates that the time t2 is carried in the delay response message sent by the master clock to the slave clock at time t3;
  • Follow_up_delay_resp(t3) indicates that the time t3 is carried in the delay response message sent by the master clock to the slave clock; in the Follow_up_delay_resp message
  • the slave clock receives timestamps t1, t2, t3, t4;
  • the link delay delay_resp between the slave clock and the master clock is calculated as shown in expression (1):
  • the processing time is the local time, which is the duration of the time interval between the master clock receiving the time synchronization event message and the master clock sending the next subsequent time synchronization event message, and the processing time is recorded as time_process, which is t3 here;
  • the border gateway includes an industrial wireless module and a TSN module, uses the SLIP protocol to encapsulate IP data packets, and ensures the transparency of data transmission through the SLIP serial port, making the data seamlessly connected between the two networks.
  • step S2 the synchronization pulse signal output by the TSN module clock is 1PPS pulse, and the industrial wireless module adjusts its own local time by detecting the synchronization pulse signal output by the TSN module clock; including the following steps:
  • the TSN module sends its own initial clock T 0 to the industrial wireless module through the serial port, and the industrial wireless module sets its own local initial time local_time0 according to this time;
  • S22 Start the interrupt program and start to detect the 1PPS pulse signal output by the TSN module. Every time the industrial wireless module detects a 1PPS signal, record the time when the 1PPS signal is received as local_time1;
  • the corrected time local_time satisfies the following expression:
  • S represents the length of the serial port message
  • R represents the serial port baud rate
  • N represents the number of pulses received by the industrial wireless module.
  • step S3 the industrial wireless module is marked as A, the routing device is marked as P 1 , and the node device is marked as P 2 .
  • node P 1 calculates the clock frequency offset through the timestamp information to adjust;
  • node P 1 After node P 1 adjusts, it sends a time message to node P 2. The moment when node P 1 sends the time message to node P 2 to the time when node P 2 receives the message is an interaction between node P 1 and node P 2 . Synchronization period, this process between node P 1 and node P 2 repeats N times to generate N synchronization periods; after each two nodes perform one-way interaction of time information (that is, a synchronization period), node P 2 passes the time stamp The information calculates the clock frequency offset for adjustment.
  • the adjusted new clock offset of node P 1 is obtained from formula (7)-(5) for:
  • node A sends its current time information Sent to node P 1 , the receiving time recorded by node P 1 Expressed as:
  • Equation (10) Equation (10)
  • P 1 only needs to obtain the parameters in formula (19), then it can be estimated by selecting the appropriate value can be adjusted
  • the frequency of node P 1 and node A is gradually approached so as to complete the time synchronization of P 1 and A.
  • step S3 the first layer node P 1 and the second layer node P 2 .
  • node P 1 receives the time information from reference node A, it adjusts its own local clock according to the time stamp calculation and then serves as the reference node of node P 2 .
  • node P 2 takes node P 1 as a reference node, and node P 1 sends time information to node P 2 ;
  • node P 1 sends its current time information sent to node P 2 , and the time when node P 2 receives the time information sent by node P 1 is recorded as Expressed as:
  • P2 is the correction time of the first synchronization period between node P1 and node P2 . Then node P2 is at Timestamp difference Adjust the local clock, the adjusted local time is used said, and has
  • the reference node P 1 needs to use the difference between the sending and receiving timestamps of node A and node P 1 in the first synchronization cycle in each cycle to correct their own local clock, so for nodes P 2 and P 1 , the difference between their sending and receiving timestamps in the second synchronization cycle Expressed as:
  • node P2 is the correction moment of P2 in the second period of node P1 and node P2 . Then node P2 is at Timestamp difference Adjust the local clock, the adjusted time is used said, and has
  • node P 2 is at The adjusted time is The time corresponding to this time to the time of this correction During this period, the true clock skew between node P 2 compared to P 1 for:
  • the adjusted clock offset of node P 2 in the second synchronization period between node P 1 and node P 2 is obtained from formulas (26)-(25) for:
  • the reference node P 1 needs to use the difference between the sending and receiving timestamps of node A and node P 1 in the second synchronization cycle in each cycle To correct their own local clock, so for nodes P 2 and P 1 , the difference between the sending and receiving timestamps of nodes P 1 and node P 2 in the third synchronization cycle Expressed as:
  • the reference node P1 For the i-th synchronization period of node P1 and node P2 , the reference node P1 needs to use the difference between the sending and receiving timestamps of node A and node P1 in the i-th synchronization period in each period to correct their own local clock, so for nodes P 2 and P 1 , the difference between the sending and receiving timestamps in the i-th synchronization cycle Expressed as:
  • P 2 only needs to obtain the parameters in formula (33), then it can be estimated by selecting the appropriate value can be adjusted Make the frequency of node P2 and node P1 approach gradually so as to complete the time synchronization of P2 and P1 .
  • step S4 the node device of the industrial wireless network is used as the master clock, denoted as M, and the terminal side converts the node as the slave clock, denoted as S, node S periodically sends ordinary data packets to node M, and by controlling the receiver
  • the response time interval ⁇ i to the sender implicitly conveys synchronization information
  • the sending node S is at Send an ordinary data packet to the receiving node M at any time, and the receiving node M is at The data packet is received at time, and after waiting for the response time interval ⁇ i , at Return a response message to the sending node S at all times, and the sending node S is in The response message is received at all times, and then the local time of sending and receiving is recorded and stored in the local cache.
  • the node can use the sending and receiving time of the local cache to estimate the clock frequency offset.
  • f S and f M be the crystal oscillator frequencies of node S and node M respectively, and define the clock frequency offset of node S relative to node M
  • the time when node M sends the response message can be obtained according to the RSA algorithm, namely The time when node S receives the response message sent by node M is recorded as but:
  • node S is at Send an ordinary data packet to node M at all times, and node M is in The data packet is received at any time, and the value of the response time interval ⁇ 2 of the second synchronization cycle is obtained according to the RSA algorithm.
  • the time for node M to send the response message but:
  • d 2 is the fixed time delay that the data packet produces during the second synchronous period data packet sending and receiving process
  • Y 2 is the downlink random time delay in the second synchronous period
  • node S is at Send an ordinary data packet to node M at all times, and node M is in Receive the data packet at any time, and get the value of the response time interval ⁇ 3 of the third synchronization cycle according to the RSA algorithm.
  • the time for node M to send the response message but:
  • Y 3 represents the random time delay of the downlink in the third synchronization period
  • node S is at Send an ordinary data packet to node M at all times, and node M is in Receive the data packet at any time, according to the RSA algorithm, the value of the response time interval ⁇ i of the i-th cycle can be obtained.
  • the time for node M to send the response message Node S is at Receive the response message returned by node M at all times. Available,
  • X i represents the random time delay generated during the sending and receiving process of the i-th cycle data packet
  • Y i represents the random time delay of the downlink in the i-th cycle.
  • the sending node S obtains a set of local timestamp pairs
  • the receiving node M obtains a set of local timestamp pairs Combined with the response time interval ⁇ i , use the maximum likelihood estimation method to estimate the frequency offset and fixed delay between nodes;
  • Equation (43) is expressed in the following matrix form:
  • Equation (44) expresses the maximum likelihood function of parameters ( ⁇ , d 1 , ⁇ 2 ) as:
  • Equation (47) and Equation (48) be 0, the maximum likelihood function of ⁇ (MS) and d 1 can be obtained by simultaneous solution:
  • the TSN switch of TSN network 1 completes synchronization with the conversion node on the terminal side in a slave clock state, and the synchronization method is consistent with step S1.
  • the method extends the time synchronization information from the two-level network to the three-level network in the wireless network, and the time adjustment of the first-level node after each synchronization period must consider the correction amount of the upper-level node, And the influence of fixed delay is added to the traditional correction amount. At the same time, it is not necessary to jointly estimate the frequency offset and clock offset. Only the frequency offset can be estimated separately to achieve synchronization with the multi-level nodes of the network. When optimizing the correction amount while reducing the complexity of the algorithm.
  • the present invention transmits the synchronization information from the master clock of the fusion network to the slave clock through the industrial wireless network through multi-stage forwarding of the time synchronization information, and completes the cross-network time synchronization of the industrial wireless network and the TSN network.
  • Fig. 1 is a cross-network network architecture diagram for industrial wireless and TSN integration according to the present invention
  • Fig. 2 is a schematic diagram of a cross-network time synchronization method
  • Fig. 3 is a flow chart of TSN network synchronization
  • FIG. 4 is a diagram of a peer-to-peer delay measurement mechanism
  • FIG. 5 is a schematic structural diagram of a border gateway
  • Fig. 6 is a flowchart of interaction of corrected time synchronization messages
  • Figure 7 is a timestamp-free synchronization protocol based on the RSA algorithm.
  • this solution innovatively proposes a cross-network network architecture oriented to the integration of industrial wireless and TSN.
  • This architecture transforms both ends of the industrial wireless network into TSNs. It consists of a border gateway and a terminal-side conversion node, where the border gateway includes an industrial wireless module and a TSN module.
  • TSN network 1 and TSN network 2 use the same protocol, industrial wireless network and the two TSN networks use different protocols.
  • the TSN network includes TSN terminals and TSN switches; industrial wireless networks can be wireless communication networks such as WIA-PA, wireless sensor networks, WirelessHART, and 5G; border gateways and terminal-side conversion nodes are used to convert protocols between industrial wireless networks and TSN networks , which acts as a bridge between the two networks.
  • industrial wireless networks can be wireless communication networks such as WIA-PA, wireless sensor networks, WirelessHART, and 5G; border gateways and terminal-side conversion nodes are used to convert protocols between industrial wireless networks and TSN networks , which acts as a bridge between the two networks.
  • Asymmetric delay in mixed wireless and wired networks is one of the main reasons for the inaccuracy of the synchronization process.
  • the synchronization accuracy is affected by the asymmetric delay, and the degree of influence depends on the accuracy of the calculated path delay.
  • the schematic diagram of the cross-network time synchronization method of this solution is shown in FIG. 2 .
  • the entire converged network uses the TSN switch of TSN network 2 as the master clock, and the TSN switch of TSN network 1 as the slave clock.
  • the time information is forwarded from the TSN switch of TSN network 2 to the In the TSN switch of TSN network 1, the master clock and the slave clock span the industrial wireless network.
  • Both TSN network and industrial wireless network use existing relatively mature clock synchronization protocols such as IEEE 802.1AS and beacon frame synchronization for synchronization, without introducing additional protocol overhead.
  • the TSN module of the border gateway performs clock synchronization with the TSN switch of TSN network 2 in the slave clock state; then inside the border gateway, the TSN module acts as the master clock of the industrial wireless module, and the industrial wireless module performs clock synchronization with the TSN module through the serial port , thus forwarding the time information from the TSN network 2 to the industrial wireless network.
  • the routing device in the industrial wireless network performs clock synchronization with the slave clock state and the industrial wireless module of the border gateway through beacon frame synchronization.
  • the routing device acts as the master clock of the node device, thereby ensuring the clock synchronization of the industrial wireless network.
  • the conversion node on the terminal side performs clock synchronization between the slave clock state and the node devices of the industrial wireless network using a time-stamp-free synchronization method, does not rely on a special synchronization protocol, and only follows the existing network data flow to complete the synchronization parameter estimation. This forwards the time information from the industrial wireless network to the TSN network 1. Subsequently, the TSN switch of TSN network 1 completes clock synchronization in slave clock state and transition node at the terminal side.
  • the synchronization information is transmitted from the master clock of the fusion network to the slave clock through the industrial wireless network step by step, and the cross-network time synchronization between the industrial wireless network and the TSN network is completed.
  • the TSN switch of TSN network 2 is used as the master clock, and the TSN module inside the border gateway is used as the slave clock.
  • the link delay and processing time of the master-slave clock are mainly considered.
  • the link delay is calculated according to IEEE 802.1AS.
  • the synchronization flow chart is shown in Figure 3 .
  • the slave clock and the master clock exchange the delay request message Delay_rep and the delay response message Delay_resp to obtain the corresponding timestamp.
  • Delay_resp(t2) indicates that the master clock replies to the slave clock at time t3. The time t2 is carried in the text.
  • Follow_up_delay_resp(t3) indicates that the time t3 is carried in the follow-up delay response message sent by the master clock to the slave clock.
  • the slave clock receives time stamps such as t1, t2, t3, t4, etc.
  • the link delay delay_resp between the slave clock and the master clock is calculated as the expression (1 ) as shown:
  • the processing time is the local time (here the duration of the time interval between the master clock receiving a time synchronization event message and the master clock sending the next subsequent time synchronization event message).
  • the processing time is recorded as time_process, which is t3 here.
  • the structure of the border gateway is shown in Figure 5, which includes an industrial wireless module and a TSN module.
  • the SLIP protocol is used in this solution to encapsulate IP data packets, and the SLIP serial port is used to ensure The transparency of data transmission makes data seamlessly connected between the two networks.
  • the TSN module synchronizes to the main clock of TSN network 2 in the state of the slave clock, and then forwards the time synchronization information to the industrial wireless module through serial communication in the state of the main clock, thus forwarding the time information from the TSN network 2 to the industrial wireless network .
  • the synchronization pulse signal output by the TSN module clock is 1PPS pulse, and the industrial wireless module adjusts its own local time by detecting the synchronization pulse signal output by the TSN module clock.
  • the TSN module For the clock synchronization between the two modules, the TSN module first needs to send its own initial clock T 0 to the industrial wireless module through the serial port, and the industrial wireless module sets its own local initial time local_time0 according to this time, and then starts the interrupt program to start detecting the TSN module For the output 1PPS pulse signal, every time a 1PPS signal is detected, record the time when the 1PPS signal is received as local_time1, and then perform a correction on its own local clock. Assuming that the industrial wireless module is corrected at local_time2, the corrected time local_time satisfies the following expression:
  • S represents the length of the serial port message
  • R represents the baud rate of the serial port
  • N represents the number of pulses received by the industrial wireless module.
  • Industrial wireless networks can run protocols such as WIA-PA, ISA100.11a, 5G or WirelessHART.
  • the clock is often adjusted only by calculating the frequency offset from the difference between the time stamps of the sending end and the receiving end, without considering the error effect caused by the fixed delay; or considering the fixed time delay. delay, but only in single-hop networks.
  • this paper estimates and compensates the frequency offset in a multi-hop network based on the synchronization mechanism of one-way message propagation, and adds a fixed delay to the correction amount.
  • the synchronization process is as follows:
  • the industrial wireless module has been synchronized to the TSN network 2 in the slave clock state according to the above method, and then the industrial wireless module forwards the time synchronization information to the routing device and the node device in the master clock state.
  • the entire industrial wireless network first uses the industrial wireless module in the border gateway as a reference node for clock synchronization.
  • the industrial wireless module transmits the synchronization information to the routing device, and the routing device adjusts its local clock according to this information; then the routing device acts as the clock of the node device
  • the reference node forwards the synchronization information to the node device, and the node device adjusts its local clock with this time synchronization information.
  • the industrial wireless module is marked as A
  • the routing device is marked as P 1
  • the node device is marked as P 2 .
  • node P 1 calculates the clock frequency offset through the timestamp information to adjust;
  • node P 1 After node P 1 adjusts, it sends a time message to node P 2. The moment when node P 1 sends the time message to node P 2 to the time when node P 2 receives the message is an interaction between node P 1 and node P 2 . Synchronization period, this process between node P 1 and node P 2 repeats N times to generate N synchronization periods; after each two nodes perform one-way interaction of time information (that is, a synchronization period), node P 2 passes the time stamp The information calculates the clock frequency offset for adjustment.
  • node A sends its current time information Sent to node P 1 , received time recorded by P 1 Can be expressed as:
  • the difference between the sending and receiving time stamps between the reference node A and the node P1 to be synchronized Can be expressed as:
  • node A sends its current time information Send to node P 1 .
  • node P 1 records the reception time Can be expressed as:
  • the difference between the sending and receiving time stamps of the reference node A and node P1 can be expressed as:
  • the difference between the sending and receiving timestamps of reference node A and node P1 can be obtained for:
  • Table 1 shows the main parameters in the process of node clock synchronization, including reference time, initial offset, correction time and correction amount, etc.
  • P 1 only needs to obtain the parameters in formula (19), then it can be estimated by selecting the appropriate value can be adjusted
  • the frequency of node P 1 and node A is gradually approached so as to complete the time synchronization of P 1 and A.
  • node P 1 After node P 1 receives the time information from reference node A, it adjusts its own local clock according to the time stamp calculation and then serves as the reference node of node P 2 . Similarly, node P2 uses node P1 as a reference node, and node P1 sends time information to node P2 .
  • node P 1 sends its current time information sent to node P 2 , and the time when node P 2 receives the time information sent by node P 1 is recorded as Can be expressed as:
  • the difference between the sending and receiving time stamps of node P 1 and node P 2 in the first synchronization period of node P 1 and node P 2 can be obtained for:
  • node P2 is at Timestamp difference Adjust the local clock, the adjusted local time is used said, and has
  • the reference node P 1 needs to use the difference between the sending and receiving timestamps of node A and node P 1 in the first synchronization cycle in each cycle To correct their own local clock, so for nodes P 2 and P 1 , the difference between the sending and receiving timestamps of nodes P 1 and node P 2 in the second synchronization cycle Can be expressed as:
  • node P2 is at Timestamp difference Adjust the local clock, the adjusted time is used said, and has
  • node P 2 is at The adjusted time is The time corresponding to this time to the time of this correction During this period, the true clock skew between node P 2 and node P 1 for:
  • the adjusted clock offset of node P 2 in the second synchronization period between node P 1 and node P 2 is obtained from formulas (26)-(25) for:
  • the reference node P 1 needs to use the difference between the sending and receiving timestamps of node A and node P 1 in the second synchronization cycle in each cycle To correct their own local clock, so for nodes P 2 and P 1 , the difference between the sending and receiving timestamps of nodes P 1 and node P 2 in the third synchronization cycle Can be expressed as:
  • the receiving time when node P 2 receives the synchronization message for the third synchronization period of node P 1 and node P 2 ; is the random time delay generated during the sending and receiving process of the time information of the third synchronization cycle between the node P1 and the node P2 ;
  • the reference node P 1 needs to use the difference between the sending and receiving timestamps of node A and node P 1 in the i-th synchronization cycle in each cycle to correct their own local clocks, so for nodes P 2 and P 1 , the difference between the sending and receiving timestamps of nodes P 1 and node P 2 in the i-th synchronization cycle Can be expressed as:
  • node P 1 in Represented as the sending time of node P 1 and node P 2 's i-th synchronization cycle node P 1 sending a synchronization message to node P 2 ;
  • P 2 only needs to obtain the parameters in formula (33), then it can be estimated by selecting the appropriate value can be adjusted Make the frequency of node P2 and node P1 approach gradually so as to complete the time synchronization of P2 and P1 .
  • This method extends the time synchronization information from the two-level network to the three-level network.
  • the time adjustment of the first-level node after each synchronization cycle must consider the correction amount of the upper-level node, and a fixed delay is added to the traditional correction amount.
  • the conversion node on the terminal side is an important node spanning the industrial wireless network and the TSN network 1, and acts as a bridge connecting the two networks.
  • most synchronization methods used hardware clock modules to obtain time stamps, which not only brought high hardware complexity, but also caused additional time synchronization because different protocol networks use different time synchronization mechanisms protocol overhead.
  • the node device and the terminal-side conversion node are synchronized here.
  • This paper uses a time stamp-free synchronization method based on the RSA algorithm, which does not rely on a special synchronization protocol and is not limited to different time synchronization protocols used by different networks.
  • the timestamp synchronization method relies on following the response message, and only follows the transmission of the existing network data stream to complete the synchronization parameter estimation, and does not need to carry timestamp information during the data stream interaction process, thereby implicitly obtaining the synchronization required for cross-network time synchronization information.
  • the node device of the industrial wireless network is used as the master clock, denoted as M, and the terminal-side conversion node is used as the slave clock, denoted as S.
  • Node S periodically sends ordinary data packets to node M.
  • the receiver Synchronization information is conveyed implicitly for the sender's response time interval ⁇ i .
  • the sending node S is at Send an ordinary data packet to the receiving node M at any time, and the receiving node M is at The data packet is received at time, and after waiting for the response time interval ⁇ i , at Return a response message to the sending node S at all times, and the sending node S is in The response message is received at all times, and then the local time of sending and receiving is recorded and stored in the local cache.
  • the node can use the sending and receiving time of the local cache to estimate the clock frequency offset.
  • the receiving node M determines the response time interval ⁇ i according to the RSA algorithm.
  • the RSA algorithm is an asymmetric cryptographic algorithm. Given a number, it can generate Different random numbers, this article only draws on the characteristics of the RSA algorithm that can generate random numbers. The specific mechanism is as follows:
  • f S and f M be the crystal oscillator frequencies of node S and node M respectively, and define the clock frequency offset of node S relative to node M
  • the time when node M sends a response message to node S can be obtained according to the above RSA algorithm, namely The time when node S receives the response message sent by node M is recorded as but:
  • node S is at Send an ordinary data packet to node M at all times, and node M is in Receive the data packet at any time, according to the RSA algorithm, the value of the response time interval ⁇ 2 of the second synchronization cycle can be obtained, and the time for node M to send the response message at this time but:
  • d 2 is the fixed time delay generated by the data packet in the process of transmitting and receiving the data packet in the second synchronous cycle
  • Y 2 is the random time delay of the downlink in the second synchronous cycle.
  • node S is at Send an ordinary data packet to node M at all times, and node M is in Receive the data packet at any time, according to the RSA algorithm, the value of the response time interval ⁇ 3 of the third synchronization cycle can be obtained.
  • the time for node M to send the response message but:
  • Y 3 represents the downlink random time delay in the third synchronization period.
  • node S is at Send an ordinary data packet to node M at all times, and node M is in Receive the data packet at any time, according to the RSA algorithm, the value of the response time interval ⁇ i of the i-th synchronization cycle can be obtained, and the time for node M to send the response message at this time Node S is at Receive the response message returned by node M at all times. Available,
  • X i represents the random time delay generated during the sending and receiving process of the i-th synchronization cycle data packet
  • Y i represents the random downlink time delay of the i-th synchronization cycle.
  • the sending node S can obtain a set of local timestamp pairs
  • the receiving node M can obtain a set of local timestamp pairs Combined with the response time interval ⁇ i , the frequency offset and fixed delay between nodes can be estimated by using the maximum likelihood estimation method.
  • Equation (43) can be expressed in the following matrix form:
  • Equation (47) and Equation (48) be 0, the maximum likelihood function of ⁇ (MS) and d 1 can be obtained by simultaneous solution:
  • the TSN switch of TSN network 1 completes synchronization with the conversion node on the terminal side in the slave clock state.
  • the synchronization method is consistent with the time synchronization process of the TSN network.

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Synchronisation In Digital Transmission Systems (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明涉及一种面向工业无线与TSN融合的跨网时间同步方法,属于工业互联网领域,包括以下步骤:S1:边界网关的TSN模块以从时钟状态与TSN网络2的TSN交换机进行时钟同步;S2:在边界网关内部,TSN模块作为工业无线模块的主时钟,工业无线模块通过串口与TSN模块进行时钟同步;S3:工业无线网络中的路由设备以从时钟状态和边界网关的工业无线模块通过信标帧同步方式进行时钟同步,同时路由设备又作为节点设备的主时钟,对节点设备进行时钟同步;S4:终端侧转换节点以从时钟状态和工业无线网络的节点设备之间使用免时间戳同步方式进行时钟同步,最后TSN网络1的TSN交换机以从时钟状态和终端侧转换节点进行同步。

Description

一种面向工业无线与TSN融合的跨网时间同步方法 技术领域
本发明属于工业互联网领域,涉及一种面向工业无线与TSN融合的跨网时间同步方法。
背景技术
传统制造工厂其工业控制网络、管理网络等内网可认为是封闭网络,外界无法访问,而随着智能制造的发展与应用推广,需要内部网络与外部网络进行互联互通;同时在旋转类机械装备、老旧工厂升级改造布线不方便和油气管线等生产现场距离过远等工业场景必须采用无线通信技术;更多工业互联用例如机器与机器交互协同、工厂物流、大型制造装备远程维护和AR/VR等则需要无线通信来满足其高效性需求和各方面性能指标。目前,包括有线和无线技术在内的融合网络研究已有一定基础并开始在自动化***、智能城市和家庭区域网络中广泛部署。一方面有线网络由于拓扑稳定性和更小的数据包丢失率而具有更高的可靠性和性能;另一方面在空间范围移动性、能耗、数据吞吐量、及时性、可靠性和其他特征方面网络连接性要求各不相同,因此有必要结合有线和无线网络技术来为工厂网络提供所需的服务质量。
工业控制业务需要由具有确定时延及可靠保障的网络承载,TSN因具有确定时延转发能力且兼容以太网协议而受到广泛关注。然而,有线TSN不能满足海量传感器及AGV等新型设备在智能工厂的广泛部署需求,工业无线网络与TSN协同传输成为工业与自动化领域的研究热点。
时钟同步是指将***中设备的时钟与某一特定的时间信息源之间的偏差限定在一定范围之内的过程,是分布式***中尤为重要的技术,其目的是通过确保任意两个节点之间的有限最大偏移来建立具有预定义精度的全局时间概念。特别是在工业控制网络中,网络中的许多基本操作如对象跟踪、确定性调度以及资源管理等都离不开时钟同步,这意味着网络中的所有节点需要共享一个统一的时间尺度。
当工厂网络将无线网络和有线网络这两种不同通信协议的网络相融合时,如何进行跨网时间同步就成为了一个具有挑战的问题。目前国内外对跨网时间同步的研究还比较少。由于工业异构网络对实时性、低延迟等特征的高要求,直接采用已有的方法满足不了新型工业异构网络的诸多需求。
发明内容
有鉴于此,本发明的目的在于有效地进行跨网时间同步,以异构网络为对象,通过IEEE802.1AS、信标帧同步、免时间戳同步等同步方式同时结合统计学算法,研究异构网络中同步节点之间的时间建模和时间误差补偿模型,提高跨网时间同步精度。
为达到上述目的,本发明提供如下技术方案:
一种面向工业无线与TSN融合的跨网时间同步方法,包括以下步骤:
S1:边界网关的TSN模块以从时钟状态与TSN网络2的TSN交换机进行时钟同步;
S2:在边界网关内部,TSN模块作为工业无线模块的主时钟,工业无线模块通过串口与TSN模块进行时钟同步,从而将时间信息从TSN网络1转发到工业无线网络;
S3:工业无线网络中的路由设备以从时钟状态和边界网关的工业无线模块通过信标帧同步方式进行时钟同步,同时路由设备又作为节点设备的主时钟,对节点设备进行时钟同步,从而确保工业无线网络的时钟同步;
S4:终端侧转换节点以从时钟状态和工业无线网络的节点设备之间使用免时间戳同步方式进行时钟同步,不依赖于专门的同步协议,仅跟随现有网络数据流完成同步参数估计,以此将时间信息从工业无线网络转发到TSN网络1。随后TSN网络1的TSN交换机以从时钟状态和终端侧转换节点完成时钟同步。
进一步,步骤S1中,具体包括以下步骤:
S11:从时钟利用对等延迟测量机制获取时间戳t1、t2、t3、t4;
S12:从时钟求取主从时钟的链路延迟;
S13:从时钟求取处理时间记为time_process;
S14:将链路延迟和处理时间之和记为Δ1。
进一步,步骤S1中,利用延迟测量机制测量主从时钟间的链路延迟,从时钟和主时钟之间通过交换延迟请求报文Delay_rep和延迟响应报文Delay_resp来获取相应的时间戳,Delay_resp(t2)表示主时钟在t3时刻回复给从时钟的延迟响应报文中携带了t2这个时刻;Follow_up_delay_resp(t3)表示主时钟发送给从时钟的跟随延迟响应报文中携带了t3这个时刻;在Follow_up_delay_resp消息发送后,从时钟接收到时间戳t1、t2、t3、t4;从时钟和主时钟之间的链路延迟delay_resp计算如表达式(1)所示:
Figure PCTCN2022091457-appb-000001
处理时间即为本地时间,是主时钟接收时间同步事件消息与主时钟发送下一个后续时间同步事件消息之间的时间间隔的持续时间,处理时间记为time_process,此处即为t3;
将链路延迟和处理时间之和记为Δ1=delay_resp+time_process。
进一步,边界网关包含一个工业无线模块和一个TSN模块,使用SLIP协议对IP数据包进行封装,通过SLIP串口来保证数据传输的透明性,使得数据在两网之间无缝连接。
进一步,步骤S2中,TSN模块时钟输出的同步脉冲信号为1PPS脉冲,工业无线模块通过检测TSN模块时钟输出的同步脉冲信号来调整自己的本地时间;包括以下步骤:
S21:TSN模块将自身的初始时钟T 0通过串口发送给工业无线模块,工业无线模块根据此时间设置自己的本地初始时间local_time0;
S22:启动中断程序开始检测TSN模块输出的1PPS脉冲信号,工业无线模块每检测到1次1PPS信号,记接收到1PPS信号的时间为local_time1;
S23:工业无线模块对自己的本地时钟进行一次校正;
假设工业无线模块在local_time2时刻进行校正,则校正后的时间local_time满足以下表达式:
local_time0=T 0+S/R       (2)
local_time=local_time2+N-local_time1      (3)
其中,S表示串口报文长度,R表示串口波特率,N表示工业无线模块接收到的脉冲个数。
进一步,步骤S3中,将工业无线模块记为A,路由设备记为P 1,节点设备记为P 2
节点A向节点P 1发送时间消息的时刻到节点P 1接收到该消息的时刻为节点A和节点P 1之间交互的一个同步周期,节点A和节点P 1间的这种过程重复N次产生N个同步周期;在每次两节点进行时间信息的单向交互后(即一个同步周期),节点P 1通过时间戳信息计算出时钟频率偏移进行调整;
节点P 1调整后,又向节点P 2发送时间消息,节点P 1向节点P 2发送时间消息的时刻到节点P 2接收到该消息的时刻为节点P 1和节点P 2之间交互的一个同步周期,节点P 1和节点P 2间的这种过程重复N次产生N个同步周期;在每次两节点进行时间信息的单向交互后(即一个同步周期),节点P 2通过时间戳信息计算出时钟频率偏移进行调整。
参考节点A和第一层节点P 1。在节点A和节点P 1的第一个同步周期,节点A将其当前的时间信息
Figure PCTCN2022091457-appb-000002
发送给节点P 1,由P 1记录的接收时间
Figure PCTCN2022091457-appb-000003
表示为:
Figure PCTCN2022091457-appb-000004
式中,
Figure PCTCN2022091457-appb-000005
表示节点P 1在t 0时刻相对于参考节点A的初始时钟偏移;
Figure PCTCN2022091457-appb-000006
表示节点P 1相对于参考节点A的时钟频率偏移;
Figure PCTCN2022091457-appb-000007
表示时间信息收发过程中产生的固定时延,取值范围(1ms~10ms);
Figure PCTCN2022091457-appb-000008
表示时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000009
表示 时间信息收发过程中由于频率偏移而积累的同步误差;
进而,在节点A和节点P 1的第一个同步周期参考节点A和待同步节点P 1之间的收发时间戳的差值
Figure PCTCN2022091457-appb-000010
表示为:
Figure PCTCN2022091457-appb-000011
Figure PCTCN2022091457-appb-000012
为P1第一个周期的校正时刻。随后节点P 1
Figure PCTCN2022091457-appb-000013
时刻利用差值
Figure PCTCN2022091457-appb-000014
来调整自己的本地时钟,P 1调整后的时间用
Figure PCTCN2022091457-appb-000015
表示,且有
Figure PCTCN2022091457-appb-000016
从初始时间t 0到校正时刻
Figure PCTCN2022091457-appb-000017
这段时间内,节点A和节点P 1之间的真实时钟偏差
Figure PCTCN2022091457-appb-000018
为:
Figure PCTCN2022091457-appb-000019
由公式(7)-(5)得节点P 1调整后新的时钟偏移量
Figure PCTCN2022091457-appb-000020
为:
Figure PCTCN2022091457-appb-000021
在节点A和节点P 1的第二个同步周期,节点A将其当前的时间信息
Figure PCTCN2022091457-appb-000022
发送给节点P 1,节点P 1记录的接收时间
Figure PCTCN2022091457-appb-000023
表示为:
Figure PCTCN2022091457-appb-000024
其中,
Figure PCTCN2022091457-appb-000025
为节点A和节点P 1的第二个同步周期时间信息收发过程中产生的随机时延。
进而,在节点A和节点P 1的第二个同步周期,参考节点A和节点P 1的收发时间戳的差值
Figure PCTCN2022091457-appb-000026
表示为:
Figure PCTCN2022091457-appb-000027
设P 1在第一个同步周期调整后的时间为
Figure PCTCN2022091457-appb-000028
该时刻对应的节点A的时间为
Figure PCTCN2022091457-appb-000029
此时由于节点P 1刚刚调整本地时间,节点A和节点P 1的时间很接近,因此设
Figure PCTCN2022091457-appb-000030
公式(10)变为:
Figure PCTCN2022091457-appb-000031
Figure PCTCN2022091457-appb-000032
Figure PCTCN2022091457-appb-000033
代入到公式(11)得:
Figure PCTCN2022091457-appb-000034
Figure PCTCN2022091457-appb-000035
Figure PCTCN2022091457-appb-000036
这段时间内(即节点A和节点P 1的第二个同步周期),节点A和节点P 1之间的真实时钟偏差
Figure PCTCN2022091457-appb-000037
为:
Figure PCTCN2022091457-appb-000038
其中
Figure PCTCN2022091457-appb-000039
为P1在节点A和节点P 1的第二个周期的校正时刻。
由公式(13)-(12)得节点P 1在第二次调整后新的时钟偏移量
Figure PCTCN2022091457-appb-000040
为:
Figure PCTCN2022091457-appb-000041
在节点A和节点P 1的第i个同步周期,得到参考节点A和节点P 1的收发时间戳的差值
Figure PCTCN2022091457-appb-000042
为:
Figure PCTCN2022091457-appb-000043
其中,
Figure PCTCN2022091457-appb-000044
为节点A和节点P 1的第i个周期时间信息收发过程中产生的随机时延。
节点P 1在第i次调整后新的时钟偏移量
Figure PCTCN2022091457-appb-000045
为:
Figure PCTCN2022091457-appb-000046
为了使计算更为简化,公式(15)化为如下表达式:
Figure PCTCN2022091457-appb-000047
Figure PCTCN2022091457-appb-000048
因为随机时延
Figure PCTCN2022091457-appb-000049
Figure PCTCN2022091457-appb-000050
可以看作是均值相同并服从高斯分布的随机变量,易知Q i也是服从高斯分布的随机变量,即Q i~(0,δ 2);使用最大似然估计方法对频率偏移
Figure PCTCN2022091457-appb-000051
进行估计,公式(17)基于观测量
Figure PCTCN2022091457-appb-000052
的对数似然函数表示为:
Figure PCTCN2022091457-appb-000053
Figure PCTCN2022091457-appb-000054
求偏导并令其为0有:
Figure PCTCN2022091457-appb-000055
得其解为:
Figure PCTCN2022091457-appb-000056
P 1仅需要获取式(19)中的参数,即可估算出
Figure PCTCN2022091457-appb-000057
通过选取适当的
Figure PCTCN2022091457-appb-000058
值即可调整
Figure PCTCN2022091457-appb-000059
使节点P 1和节点A的频率逐步接近从而完成P 1和A的时间同步。
进一步,步骤S3中,第一层节点P 1和第二层节点P 2。节点P 1收到来自参考节点A的时间信息后,根据时间戳计算调整自己的本地时钟后又作为节点P 2的参考节点。同样地节点P 2以节点P 1为参考节点,节点P 1向节点P 2发送时间信息;
在节点P 1和节点P 2的第一个同步周期,节点P 1将其当前的时间信息
Figure PCTCN2022091457-appb-000060
发送给节点P 2,节点P 2收到节点P 1发送的时间信息的时间记为
Figure PCTCN2022091457-appb-000061
表示为:
Figure PCTCN2022091457-appb-000062
式中,
Figure PCTCN2022091457-appb-000063
表示节点P 2在t 0时刻相对于参考节点P 1的初始时钟偏移;
Figure PCTCN2022091457-appb-000064
表示节点P 2相对于参考节点P 1的时钟频率偏移;
Figure PCTCN2022091457-appb-000065
表示时间信息收发过程中产生的固定时延,取值范围(1ms~10ms);
Figure PCTCN2022091457-appb-000066
表示时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000067
表示时间信息收发过程中由于频率偏移而积累的同步误差;
在节点P 1和节点P 2的第一个同步周期节点P 1和节点P 2的收发时间戳差值
Figure PCTCN2022091457-appb-000068
为:
Figure PCTCN2022091457-appb-000069
Figure PCTCN2022091457-appb-000070
为P 2节点P 1和节点P 2的第一个同步周期的校正时刻。然后节点P 2
Figure PCTCN2022091457-appb-000071
时刻利用时间戳差值
Figure PCTCN2022091457-appb-000072
调整本地时钟,调整后的本地时间用
Figure PCTCN2022091457-appb-000073
表示,且有
Figure PCTCN2022091457-appb-000074
Figure PCTCN2022091457-appb-000075
从初始时刻到校正时刻这段时间内,节点P 1和节点P 2之间的真实时钟偏差
Figure PCTCN2022091457-appb-000076
为:
Figure PCTCN2022091457-appb-000077
由公式(22)-(21)得节点P 2在节点P 1和节点P 2的第一个同步周期调整后新的时钟偏移量
Figure PCTCN2022091457-appb-000078
为:
Figure PCTCN2022091457-appb-000079
在节点P 1和节点P 2的第二个同步周期,参考节点P 1需要在每个周期内利用节点A和节点P 1在第一个同步周期的收发时间戳的差值
Figure PCTCN2022091457-appb-000080
来校正自己的本地时钟,所以对于节点P 2和P 1,在第二个同步周期其收发时间戳的差值
Figure PCTCN2022091457-appb-000081
表示为:
Figure PCTCN2022091457-appb-000082
其中
Figure PCTCN2022091457-appb-000083
为节点P 1和节点P 2的第二个同步周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000084
为节点P 1和节点P 2的第二个同步周期节点P 2接收到同步消息的接收时间;
Figure PCTCN2022091457-appb-000085
为节点P 1和节点P 2的第二个同步周期时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000086
Figure PCTCN2022091457-appb-000087
代入式(24)得在节点P 1和节点P 2的第二个同步周期节点P 1和节点P 2的收发时间戳差值
Figure PCTCN2022091457-appb-000088
为:
Figure PCTCN2022091457-appb-000089
其中
Figure PCTCN2022091457-appb-000090
表示在节点P 1和节点P 2的第一个同步周期节点P 2收到节点P 1发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000091
表示在节点A和节点P 1的第一个同步周期节点P 1收到节点A发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000092
表示在节点A和节点P 1的第一同步个周期节点A向节点P 1发送时间消息的发送时间。
Figure PCTCN2022091457-appb-000093
为P 2在节点P 1和节点P 2的第二个周期的校正时刻。然后节点P 2
Figure PCTCN2022091457-appb-000094
时刻利用时间戳差值
Figure PCTCN2022091457-appb-000095
调整本地时钟,调整后的时间用
Figure PCTCN2022091457-appb-000096
表示,且有
Figure PCTCN2022091457-appb-000097
Figure PCTCN2022091457-appb-000098
在节点P 1和节点P 2的第一个同步周期节点P 2
Figure PCTCN2022091457-appb-000099
时刻调整后的时间为
Figure PCTCN2022091457-appb-000100
此时间对应的时刻到本次校正时刻
Figure PCTCN2022091457-appb-000101
这段时间内,节点P 2相较于P 1之间真实的时钟偏差
Figure PCTCN2022091457-appb-000102
为:
Figure PCTCN2022091457-appb-000103
由公式(26)-(25)得节点P 2在节点P 1和节点P 2的第二个同步周期调整后的时钟偏移量
Figure PCTCN2022091457-appb-000104
为:
Figure PCTCN2022091457-appb-000105
其中
Figure PCTCN2022091457-appb-000106
为节点P 1和节点P 2的第二个同步周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000107
表示在节点A和节点P 1的第一个同步周期节点A和节点P 1之间的收发时间戳的差值。
在节点P 1和节点P 2的第三个同步周期,参考节点P 1需要在每个周期内利用节点A和节点P 1在第二个同步周期的收发时间戳的差值
Figure PCTCN2022091457-appb-000108
来校正自己的本地时钟,所以对于节点P 2和 P 1,在节点P 1和节点P 2的第三个同步周期其收发时间戳的差值
Figure PCTCN2022091457-appb-000109
表示为:
Figure PCTCN2022091457-appb-000110
Figure PCTCN2022091457-appb-000111
Figure PCTCN2022091457-appb-000112
代入式(28)得在节点P 1和节点P 2的第三个同步周期节点P 1和节点P 2的收发时间戳差值
Figure PCTCN2022091457-appb-000113
为:
Figure PCTCN2022091457-appb-000114
其中
Figure PCTCN2022091457-appb-000115
表示在节点P 1和节点P 2的第二个同步周期节点P 2收到节点P 1发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000116
表示在节点A和节点P 1的第二个同步周期节点P 1收到节点A发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000117
表示在节点A和节点P 1的第二个同步周期节点A向节点P 1发送时间消息的发送时间;
Figure PCTCN2022091457-appb-000118
表示在节点A和节点P 1的第一个同步周期节点P 1收到节点A发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000119
表示在节点A和节点P 1的第一个同步周期节点A向节点P 1发送时间消息的发送时间。
对于节点P 1和节点P 2的第i个同步周期,参考节点P 1需要在每个周期内利用节点A和节点P 1在第i个同步周期的收发时间戳的差值
Figure PCTCN2022091457-appb-000120
来校正自己的本地时钟,所以对于节点P 2和P 1,在第i个同步周期其收发时间戳的差值
Figure PCTCN2022091457-appb-000121
表示为:
Figure PCTCN2022091457-appb-000122
Figure PCTCN2022091457-appb-000123
其中
Figure PCTCN2022091457-appb-000124
表示为节点P 1和节点P 2的第i个同步周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000125
表示在节点P 1和节点P 2的第i-1个同步周期节点P 2收到节点P 1发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000126
表示在节点P 1和节点P 2的第i个同步周期时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000127
表示在节点A和节点P 1的第i-1个同步周期节点P 1收到节点A发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000128
表示在节点A和节点P 1的第i-1个同步周期节点A向节点P 1发送时间消息的发送时间;
Figure PCTCN2022091457-appb-000129
表示为节点P 1和节点P 2的第i-j-1个同步周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000130
表示在节点A和节点P 1的第i-j-1个同步周期节点A向节点P 1发送时间消息的发送时间。
使用最大似然估计方法对频率偏移
Figure PCTCN2022091457-appb-000131
进行估计,公式(30)基于观测量
Figure PCTCN2022091457-appb-000132
的对数似然函数表示为:
Figure PCTCN2022091457-appb-000133
式(31)两边对
Figure PCTCN2022091457-appb-000134
求偏导有:
Figure PCTCN2022091457-appb-000135
令式(32)等于0得其解为:
Figure PCTCN2022091457-appb-000136
P 2仅需要获取式(33)中的参数,即可估算出
Figure PCTCN2022091457-appb-000137
通过选取适当的
Figure PCTCN2022091457-appb-000138
值即可调整
Figure PCTCN2022091457-appb-000139
使节点P 2和节点P 1的频率逐步接近从而完成P 2和P 1的时间同步。
进一步,步骤S4中,以工业无线网络的节点设备为主时钟,记为M,终端侧转换节点作为从时钟,记为S,节点S周期性地向节点M发送普通数据包,通过控制接收方对发送方的响应时间间隔Δ i来隐式地传递同步信息;
在第i个周期发送节点S在
Figure PCTCN2022091457-appb-000140
时刻向接收节点M发送一个普通数据包,接收节点M在
Figure PCTCN2022091457-appb-000141
时刻收到该数据包,并且在等待响应时间间隔Δ i后,在
Figure PCTCN2022091457-appb-000142
时刻向发送节点S返回一个响应报文,发送节点S在
Figure PCTCN2022091457-appb-000143
时刻收到该响应报文,随后记录下收发的本地时刻并存储在本地缓存中,伴随着节点数据包收发过程的重复进行,节点可利用本地缓存的收发时间估计出时钟频率偏移。
进一步,步骤S4中的时间同步次数设为N,有A∈{2,3,4,...N-1},任取素数P、Q、D,其中N=P*Q,当(D*E)mod((P-1)*(Q-1))=1成立,则有C i=(A i EXP D)mod N,其中1<i<N,则C i就是生成的随机数,同时用C i表示第i次同步过程中接收节点需等待的响应时间间隔Δ i。注意,根据RSA算法原理,此处虽然同步N次,但只能得出N-2个同步过程中的响应时间间隔Δ i,但当N足够大的时候,同步N-2次和同步N次可视为是一样的。
设f S、f M分别是节点S和节点M的晶振频率,定义节点S相对于节点M的时钟频率偏移
Figure PCTCN2022091457-appb-000144
第一个同步周期,当节点S在
Figure PCTCN2022091457-appb-000145
时刻向节点M发送一个普通数据包,节点M收到该数据包的时间
Figure PCTCN2022091457-appb-000146
表示为下式:
Figure PCTCN2022091457-appb-000147
其中,
Figure PCTCN2022091457-appb-000148
表示节点S相对于节点M的初始时钟偏移;d 1表示数据包在收发过程中产生的固定时延;X 1表示数据包在收发过程中产生的随机时延;
节点M发送响应报文的时间
Figure PCTCN2022091457-appb-000149
可根据所述RSA算法获得,即
Figure PCTCN2022091457-appb-000150
节点S收到节点M发送的响应报文的时间记为
Figure PCTCN2022091457-appb-000151
则:
Figure PCTCN2022091457-appb-000152
其中,d 2是响应报文在传输中的固定时延,且d 2=d 1-m,这里m的大小与上行的数据包报文长度和下行的响应报文长度差成正比;Y 1是第一个同步周期中下行链路随机时延;
第二个同步周期,节点S在
Figure PCTCN2022091457-appb-000153
时刻向节点M发送一个普通数据包,节点M在
Figure PCTCN2022091457-appb-000154
时刻收到该数据包,根据RSA算法得到第二个同步周期的响应时间间隔Δ 2的值,此时节点M发送响应报文的时间
Figure PCTCN2022091457-appb-000155
则:
Figure PCTCN2022091457-appb-000156
其中,d 2是数据包在第二个同步周期数据包收发过程中产生的固定时延;Y 2第二个同步周期中下行链路随机时延;
第三个同步周期,节点S在
Figure PCTCN2022091457-appb-000157
时刻向节点M发送一个普通数据包,节点M在
Figure PCTCN2022091457-appb-000158
时刻收到该数据包,根据RSA算法得到第三个同步周期的响应时间间隔Δ 3的值,此时节点M发送响应报文的时间
Figure PCTCN2022091457-appb-000159
则:
Figure PCTCN2022091457-appb-000160
其中,Y 3表示第三个同步周期中下行链路随机时延;
第i个同步周期,节点S在
Figure PCTCN2022091457-appb-000161
时刻向节点M发送一个普通数据包,节点M在
Figure PCTCN2022091457-appb-000162
时刻收到该数据包,根据RSA算法可得第i个周期的响应时间间隔Δ i的值,此时节点M发送响应报文的时间
Figure PCTCN2022091457-appb-000163
节点S在
Figure PCTCN2022091457-appb-000164
时刻受到节点M返回的响应报文。可得,
Figure PCTCN2022091457-appb-000165
Figure PCTCN2022091457-appb-000166
Δ i∈(Δ 2,Δ 3,...,Δ N-1)             (40)
Figure PCTCN2022091457-appb-000167
其中X i表示第i个周期数据包在收发过程中产生的随机时延;Y i表示第i个周期下行链路随机时延。
N-1个同步周期后,发送节点S获得一组本地时间戳对
Figure PCTCN2022091457-appb-000168
接收节点M获得一组本地时间戳对
Figure PCTCN2022091457-appb-000169
结合响应时间间隔Δ i,利用最大似然估计方法估计出节点间的频率偏移和固定延迟;
用式(39)-(38)可得:
Figure PCTCN2022091457-appb-000170
因为X i和Y i是均值为0,方差为σ 2的高斯随机分布,令V i=X i+Y i,所以V i是均值为0,方差为2σ 2的高斯随机分布,令
Figure PCTCN2022091457-appb-000171
则式(42)简化为:
Figure PCTCN2022091457-appb-000172
其中,
Figure PCTCN2022091457-appb-000173
则式(43)表示为下列矩阵形式:
Figure PCTCN2022091457-appb-000174
式(44)关于参数(α,d 12)最大似然函数表示为:
Figure PCTCN2022091457-appb-000175
对式(45)两边取对数有:
Figure PCTCN2022091457-appb-000176
对式(46)求关于α的一阶偏导数有:
Figure PCTCN2022091457-appb-000177
对式(46)求关于d 1的一阶偏导数有:
Figure PCTCN2022091457-appb-000178
令式(47)和式(48)为0,联立求解可得ρ (MS)和d 1的最大似然函数为:
Figure PCTCN2022091457-appb-000179
Figure PCTCN2022091457-appb-000180
所以,响应报文在传输中的固定时延d 2=d 1-m。
在终端侧转换节点和工业无线网络的节点设备完成同步后,TSN网络1的TSN交换机以从时钟状态和终端侧转换节点完成同步,其同步方法和步骤S1一致。
本发明的有益效果在于:本方法在无线网络中将时间同步信息从两级网络扩展到三级网络中,每次同步周期后一级节点的时间调整都要考虑上一级节点的校正量,并且在传统的校正量中加入了固定延迟的影响,同时不需要联合估计频率偏移和时钟偏移,只需单独估计出频率偏移就可以实现与网络多级节点的同步,在优化校正量的同时减少了算法复杂度。同时,本发明通过多级转发时间同步信息的方式,使同步信息由融合网络的主时钟经工业无线网络向从时钟逐级进行传递,完成工业无线网络和TSN网络的跨网时间同步。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为本发明所述面向工业无线与TSN融合的跨网网络架构图;
图2为跨网时间同步方法示意图;
图3为TSN网络同步流程图;
图4为对等延迟测量机制图;
图5为边界网关结构示意图;
图6为校正型时间同步消息交互流程图;
图7为基于RSA算法的免时间戳同步协议。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本 发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
其中,附图仅用于示例性说明,表示的仅是示意图,而非实物图,不能理解为对本发明的限制;为了更好地说明本发明的实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
本发明实施例的附图中相同或相似的标号对应相同或相似的部件;在本发明的描述中,需要理解的是,若有术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此附图中描述位置关系的用语仅用于示例性说明,不能理解为对本发明的限制,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
随着工业化与信息化的不断融合发展,在实际的工厂网络中,网络大多是由有线和无线网络融合构成的。在这种异构网络中,网络在进行数据融合、能量管理、确定性调度时要求不同网络间的设备具有较强的时序约束关系。如图1所示,本方案创新性的提出了一种面向工业无线与TSN融合的跨网网络架构,此架构把工业无线网络两端TSN化,由两个TSN网络、一个工业无线网络、一个边界网关以及一个终端侧转换节点组成,其中边界网关包含一个工业无线模块和一个TSN模块。TSN网络1和TSN网络2使用相同的协议,工业无线网络和两个TSN网络使用不同的协议。
TSN网络包括TSN终端和TSN交换机;工业无线网络可以是WIA-PA、无线传感器网络、WirelessHART以及5G等无线通信网络;边界网关和终端侧转换节点用来转换工业无线网络和TSN网络之间的协议,起着沟通两网的桥梁作用。
无线和有线混合网络中的非对称延迟是同步过程不准确的主要原因之一,同步精度受不对称延迟的影响,其影响程度取决于计算出的路径延迟的精度。本方案跨网时间同步方法示意图如图2所示。
整个融合网络以TSN网络2的TSN交换机为主时钟,以TSN网络1的TSN交换机为从时钟,时间信息从TSN网络2的TSN交换机经过边界网关、路由设备、节点设备和终端侧转换节点转发到TSN网络1的TSN交换机,主时钟和从时钟之间跨越了工业无线网络。TSN网络和工业无线网络均使用诸如IEEE 802.1AS、信标帧同步等现有较为成熟的时钟同步协议进行同步,不引入额外的协议开销。
首先,边界网关的TSN模块以从时钟状态与TSN网络2的TSN交换机进行时钟同步; 然后在边界网关内部,TSN模块又作为工业无线模块的主时钟,工业无线模块通过串口与TSN模块进行时钟同步,从而将时间信息从TSN网络2转发到了工业无线网络。
其次,工业无线网络中的路由设备以从时钟状态和边界网关的工业无线模块通过信标帧同步方式进行时钟同步,同时路由设备又作为节点设备的主时钟,从而确保工业无线网络的时钟同步。
最后,终端侧转换节点以从时钟状态和工业无线网络的节点设备之间使用免时间戳同步方式进行时钟同步,不依赖于专门的同步协议,仅跟随现有网络数据流完成同步参数估计,以此将时间信息从工业无线网络转发到TSN网络1。随后TSN网络1的TSN交换机以从时钟状态和终端侧转换节点完成时钟同步。
通过这样一级一级转发时间同步信息的方式,使同步信息由融合网络的主时钟经工业无线网络向从时钟逐级进行传递,完成工业无线网络和TSN网络的跨网时间同步。
TSN网络时间同步过程中以TSN网络2的TSN交换机为主时钟,边界网关内部的TSN模块为从时钟。考虑到有线网络拓扑结构相对稳定,故在进行TSN网络时间同步过程中主要考虑主从时钟的链路延迟和处理时间,其链路延迟根据IEEE 802.1AS计算,其同步流程图如图3所示。
1)链路延迟测量。利用如图4所示的延迟测量机制测量主从时钟间的链路延迟。
图4中从时钟和主时钟之间通过交换延迟请求报文Delay_rep和延迟响应报文Delay_resp来获取相应的时间戳,图中Delay_resp(t2)表示主时钟在t3时刻回复给从时钟的延迟响应报文中携带了t2这个时刻,同理,Follow_up_delay_resp(t3)表示主时钟发送给从时钟的跟随延迟响应报文中携带了t3这个时刻。在跟随延迟响应报文Follow_up_delay_resp消息发送后,从时钟就接收到t1、t2、t3、t4等时间戳,根据IEEE 802.1AS,从时钟和主时钟之间的链路延迟delay_resp计算如表达式(1)所示:
Figure PCTCN2022091457-appb-000181
2)处理时间测量。处理时间即为本地时间(此处为主时钟接收时间同步事件消息与主时钟发送下一个后续时间同步事件消息之间的时间间隔的持续时间)。处理时间记为time_process,此处即为t3。
最后将链路延迟和处理时间之和记为Δ1=delay_resp+time_process。
边界网关结构如图5所示,包含一个工业无线模块和一个TSN模块,为支持WIA-PA等非IP网络的IP化,本方案中使用SLIP协议对IP数据包进行封装,通过SLIP串口来保证数据传输的透明性,使得数据在两网之间无缝连接。
其中TSN模块以从时钟状态同步到TSN网络2的主时钟后,然后又以主时钟状态通过串口通信将时间同步信息转发到工业无线模块,由此将时间信息从TSN网络2转发到了工业无线网络。
TSN模块时钟输出的同步脉冲信号为1PPS脉冲,工业无线模块通过检测TSN模块时钟输出的同步脉冲信号来调整自己的本地时间。
对于两模块之间的时钟同步,首先需要TSN模块将自身的初始时钟T 0通过串口发送给工业无线模块,工业无线模块根据此时间设置自己的本地初始时间local_time0,然后启动中断程序开始检测TSN模块输出的1PPS脉冲信号,每检测到1次1PPS信号,记接收到1PPS信号的时间为local_time1,随后对自己的本地时钟进行一次校正。假设工业无线模块在local_time2时刻进行校正,则校正后的时间local_time满足以下表达式:
local_time0=T 0+S/R             (2)
local_time=local_time2+N-local_time1            (3)
其中,S表示串口报文长度;R表示串口波特率;N表示工业无线模块接收到的脉冲个数。
工业无线网络可运行WIA-PA、ISA100.11a、5G或WirelessHART等协议。在以往的同步校正型算法中,往往只通过发送端和接收端时间戳的差值求取频率偏移从而来调整时钟,并没有考虑固定时延带来的误差影响;或是考虑了固定时延的影响,但只在单跳网络中考虑。为了进一步提升节点时钟的校正精度,本文基于单向消息传播的同步机制,在多跳网络中对频率偏移进行估计和补偿,同时在校正量中加入固定时延。同步过程如下:
在此基础上工业无线模块根据上述方法以从时钟状态已同步到TSN网络2,接下来工业无线模块又以主时钟状态将时间同步信息转发到路由设备和节点设备。
整个工业无线网络首先以边界网关内的工业无线模块作为时钟同步的参考节点,工业无线模块将同步信息传递到路由设备,路由根据此信息调整自己的本地时钟;然后路由设备又作为节点设备的时钟的参考节点并将同步信息转发到节点设备,节点设备以此时间同步信息来调整自己的本地时钟。如图6所示,将工业无线模块记为A,路由设备记为P 1,节点设备记为P 2
节点A向节点P 1发送时间消息的时刻到节点P 1接收到该消息的时刻为节点A和节点P 1之间交互的一个同步周期,节点A和节点P 1间的这种过程重复N次产生N个同步周期;在每次两节点进行时间信息的单向交互后(即一个同步周期),节点P 1通过时间戳信息计算出时钟频率偏移进行调整;
节点P 1调整后,又向节点P 2发送时间消息,节点P 1向节点P 2发送时间消息的时刻到节点P 2 接收到该消息的时刻为节点P 1和节点P 2之间交互的一个同步周期,节点P 1和节点P 2间的这种过程重复N次产生N个同步周期;在每次两节点进行时间信息的单向交互后(即一个同步周期),节点P 2通过时间戳信息计算出时钟频率偏移进行调整。
1、参考节点A和第一层节点P 1
在第一个同步周期,节点A将其当前的时间信息
Figure PCTCN2022091457-appb-000182
发送给节点P 1,由P 1记录的接收时间
Figure PCTCN2022091457-appb-000183
可表示为:
Figure PCTCN2022091457-appb-000184
式中,
Figure PCTCN2022091457-appb-000185
表示节点P 1在t 0时刻相对于参考节点A的初始时钟偏移;
Figure PCTCN2022091457-appb-000186
表示节点P 1相对于参考节点A的时钟频率偏移;
Figure PCTCN2022091457-appb-000187
表示时间信息收发过程中产生的固定时延,取值范围(1ms~10ms);
Figure PCTCN2022091457-appb-000188
表示时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000189
表示时间信息收发过程中由于频率偏移而积累的同步误差。
进而,在第一个同步周期,参考节点A和待同步节点P 1之间的收发时间戳的差值
Figure PCTCN2022091457-appb-000190
可表示为:
Figure PCTCN2022091457-appb-000191
Figure PCTCN2022091457-appb-000192
为P1第一个同步周期的校正时刻。节点P 1
Figure PCTCN2022091457-appb-000193
时刻利用差值
Figure PCTCN2022091457-appb-000194
来调整自己的本地时钟,P 1调整后的时间用
Figure PCTCN2022091457-appb-000195
表示,且有
Figure PCTCN2022091457-appb-000196
从初始时间t 0到校正时刻
Figure PCTCN2022091457-appb-000197
这段时间内,节点A和节点P 1之间的真实时钟偏差
Figure PCTCN2022091457-appb-000198
为:
Figure PCTCN2022091457-appb-000199
由公式(7)-(5)可得节点P 1调整后新的时钟偏移量
Figure PCTCN2022091457-appb-000200
为:
Figure PCTCN2022091457-appb-000201
在第二个同步周期,节点A将其当前的时间信息
Figure PCTCN2022091457-appb-000202
发送给节点P 1。在第二个同步周期,节点P 1记录的接收时间
Figure PCTCN2022091457-appb-000203
可表示为:
Figure PCTCN2022091457-appb-000204
其中,
Figure PCTCN2022091457-appb-000205
为第二个同步周期时间信息收发过程中产生的随机时延。
进而,在第二个同步周期,参考节点A和节点P 1的收发时间戳的差值
Figure PCTCN2022091457-appb-000206
可表示为:
Figure PCTCN2022091457-appb-000207
设P 1在第一个同步周期调整后的时间为
Figure PCTCN2022091457-appb-000208
该时刻对应的节点A的时间为
Figure PCTCN2022091457-appb-000209
此时由于节点P 1刚刚调整本地时间,节点A和节点P 1的时间很接近,因此可设
Figure PCTCN2022091457-appb-000210
所以公式(10)变为:
Figure PCTCN2022091457-appb-000211
Figure PCTCN2022091457-appb-000212
Figure PCTCN2022091457-appb-000213
代入到公式(11)可得:
Figure PCTCN2022091457-appb-000214
Figure PCTCN2022091457-appb-000215
Figure PCTCN2022091457-appb-000216
这段时间内(即第二个同步周期),节点A和节点P 1之间的真实时钟偏差
Figure PCTCN2022091457-appb-000217
为:
Figure PCTCN2022091457-appb-000218
其中
Figure PCTCN2022091457-appb-000219
为P1第二个同步周期的校正时刻。由公式(13)-(12)可得节点P 1在第二次调整后新的时钟偏移量
Figure PCTCN2022091457-appb-000220
为:
Figure PCTCN2022091457-appb-000221
类似地,在第i个同步周期,可以得到参考节点A和节点P 1的收发时间戳的差值
Figure PCTCN2022091457-appb-000222
为:
Figure PCTCN2022091457-appb-000223
其中,
Figure PCTCN2022091457-appb-000224
为第i个同步周期时间信息收发过程中产生的随机时延。
节点P 1在第i次调整后新的时钟偏移量
Figure PCTCN2022091457-appb-000225
为:
Figure PCTCN2022091457-appb-000226
表1给出了节点时钟同步过程中的主要参数,包括参考时间、初始偏移量、校正时间以及校正量等。
表1时间同步中的主要参数
Figure PCTCN2022091457-appb-000227
Figure PCTCN2022091457-appb-000228
为了使计算更为简化,公式(15)可化为如下表达式:
Figure PCTCN2022091457-appb-000229
Figure PCTCN2022091457-appb-000230
因为随机时延
Figure PCTCN2022091457-appb-000231
Figure PCTCN2022091457-appb-000232
可以看作是均值相同并服从高斯分布的随机变量,易知Q i也是服从高斯分布的随机变量,即Q i~(0,δ 2)。使用最大似然估计方法对频率偏移
Figure PCTCN2022091457-appb-000233
进行估计,公式(17)基于观测量
Figure PCTCN2022091457-appb-000234
的对数似然函数表示为:
Figure PCTCN2022091457-appb-000235
Figure PCTCN2022091457-appb-000236
求偏导并令其为0有:
Figure PCTCN2022091457-appb-000237
得其解为:
Figure PCTCN2022091457-appb-000238
P 1仅需要获取式(19)中的参数,即可估算出
Figure PCTCN2022091457-appb-000239
通过选取适当的
Figure PCTCN2022091457-appb-000240
值即可调整
Figure PCTCN2022091457-appb-000241
使节点P 1和节点A的频率逐步接近从而完成P 1和A的时间同步。
2、第一层节点P 1和第二层节点P 2
节点P 1收到来自参考节点A的时间信息后,根据时间戳计算调整自己的本地时钟后又作为节点P 2的参考节点。同样地节点P 2以节点P 1为参考节点,节点P 1向节点P 2发送时间信息。
在节点P 1和节点P 2的第一个同步周期,节点P 1将其当前的时间信息
Figure PCTCN2022091457-appb-000242
发送给节点P 2,节点P 2收到节点P 1发送的时间信息的时间记为
Figure PCTCN2022091457-appb-000243
可表示为:
Figure PCTCN2022091457-appb-000244
式中,
Figure PCTCN2022091457-appb-000245
表示节点P 2在t 0时刻相对于参考节点P 1的初始时钟偏移;
Figure PCTCN2022091457-appb-000246
表示节点P 2相对于参考节点P 1的时钟频率偏移;
Figure PCTCN2022091457-appb-000247
表示时间信息收发过程中产生的固定时延,取值范围(1ms~10ms);
Figure PCTCN2022091457-appb-000248
表示时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000249
为时 间信息收发过程中由于频率偏移而积累的同步误差。
进而可得在节点P 1和节点P 2的第一个同步周期节点P 1和节点P 2的收发时间戳差值
Figure PCTCN2022091457-appb-000250
为:
Figure PCTCN2022091457-appb-000251
Figure PCTCN2022091457-appb-000252
为P 2第一个同步周期的校正时刻。然后节点P 2
Figure PCTCN2022091457-appb-000253
时刻利用时间戳差值
Figure PCTCN2022091457-appb-000254
调整本地时钟,调整后的本地时间用
Figure PCTCN2022091457-appb-000255
表示,且有
Figure PCTCN2022091457-appb-000256
从初始时刻到校正时刻这段时间内,节点P 1和节点P 2之间的真实时钟偏差
Figure PCTCN2022091457-appb-000257
为:
Figure PCTCN2022091457-appb-000258
由公式(22)-(21)得节点P 2在节点P 1和节点P 2的第一个同步周期调整后新的时钟偏移量
Figure PCTCN2022091457-appb-000259
为:
Figure PCTCN2022091457-appb-000260
在第二个同步周期,参考节点P 1需要在每个周期内利用节点A和节点P 1在第一个同步周期的收发时间戳的差值
Figure PCTCN2022091457-appb-000261
来校正自己的本地时钟,所以对于节点P 2和P 1,在节点P 1和节点P 2的第二个同步周期其收发时间戳的差值
Figure PCTCN2022091457-appb-000262
可表示为:
Figure PCTCN2022091457-appb-000263
其中
Figure PCTCN2022091457-appb-000264
为节点P 1和节点P 2的第二个同步周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000265
为第节点P 1和节点P 2的二个同步周期节点P 2接收到同步消息的接收时间;
Figure PCTCN2022091457-appb-000266
为节点P 1和节点P 2的第二个同步周期时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000267
Figure PCTCN2022091457-appb-000268
代入式(24),可得在节点P 1和节点P 2的第二个同步周期节点P 1和节点P 2的收发时间戳差值
Figure PCTCN2022091457-appb-000269
为::
Figure PCTCN2022091457-appb-000270
其中
Figure PCTCN2022091457-appb-000271
表示在节点P 1和节点P 2的第一个同步周期节点P 2收到节点P 1发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000272
表示在节点A和节点P 1的第一个同步周期节点P 1收到节点A发送的时间信 息的接收时间;
Figure PCTCN2022091457-appb-000273
表示在节点A和节点P 1的第一个同步周期节点A向节点P 1发送时间消息的发送时间。
Figure PCTCN2022091457-appb-000274
为P 2第二个同步周期的校正时刻。然后节点P 2
Figure PCTCN2022091457-appb-000275
时刻利用时间戳差值
Figure PCTCN2022091457-appb-000276
调整本地时钟,调整后的时间用
Figure PCTCN2022091457-appb-000277
表示,且有
Figure PCTCN2022091457-appb-000278
在节点P 1和节点P 2的第一个同步周期节点P 2
Figure PCTCN2022091457-appb-000279
时刻调整后的时间为
Figure PCTCN2022091457-appb-000280
此时间对应的时刻到本次校正时刻
Figure PCTCN2022091457-appb-000281
这段时间内,节点P 2相较于节点P 1之间真实的时钟偏差
Figure PCTCN2022091457-appb-000282
为:
Figure PCTCN2022091457-appb-000283
由公式(26)-(25)得节点P 2在节点P 1和节点P 2的第二个同步周期调整后的时钟偏移量
Figure PCTCN2022091457-appb-000284
为:
Figure PCTCN2022091457-appb-000285
其中
Figure PCTCN2022091457-appb-000286
为节点P 1和节点P 2的第二个同步周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000287
表示在节点A和节点P 1的第一个同步周期节点A和节点P 1之间的收发时间戳的差值。
在节点P 1和节点P 2的第三个同步周期,参考节点P 1需要在每个周期内利用节点A和节点P 1在第二个同步周期的收发时间戳的差值
Figure PCTCN2022091457-appb-000288
来校正自己的本地时钟,所以对于节点P 2和P 1,在节点P 1和节点P 2的第三个同步周期其收发时间戳的差值
Figure PCTCN2022091457-appb-000289
可表示为:
Figure PCTCN2022091457-appb-000290
其中
Figure PCTCN2022091457-appb-000291
为节点P 1和节点P 2的第三个同步周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000292
为节点P 1和节点P 2的第三个同步周期节点P 2接收到同步消息的接收时间;
Figure PCTCN2022091457-appb-000293
为节点P 1和节点P 2的第三个同步周期时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000294
Figure PCTCN2022091457-appb-000295
代入式(28)可得在节点P 1和节点P 2的第三个同步周期节点P 1和节点P 2的收发时间戳差值
Figure PCTCN2022091457-appb-000296
为:
Figure PCTCN2022091457-appb-000297
其中
Figure PCTCN2022091457-appb-000298
表示在节点P 1和节点P 2的第二个同步周期节点P 2收到节点P 1发送的时间信息 的接收时间;
Figure PCTCN2022091457-appb-000299
表示在节点A和节点P 1的第二个同步周期节点P 1收到节点A发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000300
表示在节点A和节点P 1的第二个同步周期节点A向节点P 1发送时间消息的发送时间;
Figure PCTCN2022091457-appb-000301
表示在节点A和节点P 1的第一个同步周期节点P 1收到节点A发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000302
表示在节点A和节点P 1的第一个同步周期节点A向节点P 1发送时间消息的发送时间。
类似地,对于节点P 1和节点P 2的第i个同步周期,参考节点P 1需要在每个周期内利用节点A和节点P 1在第i个同步周期的收发时间戳的差值
Figure PCTCN2022091457-appb-000303
来校正自己的本地时钟,所以对于节点P 2和P 1,在节点P 1和节点P 2的第i个同步周期其收发时间戳的差值
Figure PCTCN2022091457-appb-000304
可表示为:
Figure PCTCN2022091457-appb-000305
Figure PCTCN2022091457-appb-000306
其中
Figure PCTCN2022091457-appb-000307
表示为节点P 1和节点P 2的第i个同步周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000308
表示在节点P 1和节点P 2的第i-1个同步周期节点P 2收到节点P 1发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000309
表示在节点P 1和节点P 2的第i个同步周期时间信息收发过程中产生的随机时延;
Figure PCTCN2022091457-appb-000310
表示在节点A和节点P 1的第i-1个同步周期节点P 1收到节点A发送的时间信息的接收时间;
Figure PCTCN2022091457-appb-000311
表示在节点A和节点P 1的第i-1个同步周期节点A向节点P 1发送时间消息的发送时间;
Figure PCTCN2022091457-appb-000312
表示为节点P 1和节点P 2的第i-j-1同步个周期节点P 1向节点P 2发送同步消息的发送时间;
Figure PCTCN2022091457-appb-000313
表示在节点A和节点P 1的第i-j-1个同步周期节点A向节点P 1发送时间消息的发送时间。
使用最大似然估计方法对频率偏移
Figure PCTCN2022091457-appb-000314
进行估计,公式(30)基于观测量
Figure PCTCN2022091457-appb-000315
的对数似然函数表示为:
Figure PCTCN2022091457-appb-000316
式(31)两边对
Figure PCTCN2022091457-appb-000317
求偏导有:
Figure PCTCN2022091457-appb-000318
Figure PCTCN2022091457-appb-000319
令式(32)等于0得其解为:
Figure PCTCN2022091457-appb-000320
P 2仅需要获取式(33)中的参数,即可估算出
Figure PCTCN2022091457-appb-000321
通过选取适当的
Figure PCTCN2022091457-appb-000322
值即可调整
Figure PCTCN2022091457-appb-000323
使节点P 2和节点P 1的频率逐步接近从而完成P 2和P 1的时间同步。
本方法将时间同步信息从两级网络扩展到三级网络中,每次同步周期后一级节点的时间调整都要考虑上一级节点的校正量,并且在传统的校正量中加入了固定延迟的影响,同时不需要联合估计频率偏移和时钟偏移,只需单独估计出频率偏移就可以实现与节点A的同步,在优化校正量的同时减少了算法复杂度。
终端侧转换节点是跨越工业无线网络和TSN网络1的重要节点,起着连接两网的桥梁作用。在以往的跨网时间同步方法中,大多数同步方式是采用通过硬件时钟模块来获取时间戳,这不仅带来了较高的硬件复杂度,还因为不同协议网络使用不同的时间同步机制造成额外的协议开销。
针对上述问题,此处节点设备和终端侧转换节点同步本文使用基于RSA算法的免时间戳同步方法,不依赖专门的同步协议,不受限于不同网络使用不同的时间同步协议,摆脱现有免时间戳同步方法对跟随响应报文的依赖,仅跟随现有网络数据流的发送完成同步参数估计,数据流交互过程中无需携带时间戳信息,从而隐式地获得跨网时间同步所需的同步消息。
如图7所示,以工业无线网络的节点设备为主时钟,记为M,终端侧转换节点作为从时钟,记为S,节点S周期性地向节点M发送普通数据包,通过控制接收方对发送方的响应时间间隔Δ i来隐式地传递同步信息。
在第i个同步周期发送节点S在
Figure PCTCN2022091457-appb-000324
时刻向接收节点M发送一个普通数据包,接收节点M在
Figure PCTCN2022091457-appb-000325
时刻收到该数据包,并且在等待响应时间间隔Δ i后,在
Figure PCTCN2022091457-appb-000326
时刻向发送节点S返回一个响应报文,发送节点S在
Figure PCTCN2022091457-appb-000327
时刻收到该响应报文,随后记录下收发的本地时刻并存储在本地缓存中,伴随着节点数据包收发过程的重复进行,节点可利用本地缓存的收发时间估计出时钟频率偏移。
1、预设响应时间间隔Δ i
在基于RSA算法的免时间戳同步方法中,接收节点M根据RSA算法来确定响应时间间 隔Δ i,RSA算法是一种不对称密码算法,给定一个数,则可以产生这个数范围之内的不同随机数,此处本文只借鉴RSA算法能生成随机数的特点。具体机制如下:
假设此过程同步N次,有A∈{2,3,4,...N-1},任取素数P、Q、D,其中N=P*Q,当(D*E)mod((P-1)*(Q-1))=1成立,则有C i=(A i EXP D)mod N,其中1<i<N,则C i就是生成的随机数,同时用C i表示第i次同步过程中接收节点需等待的响应时间间隔Δ i。注意,根据RSA算法原理,此处虽然同步N次,但只能得出N-2个同步过程中的响应时间间隔Δ i,但当N足够大的时候,同步N-2次和同步N次可视为是一样的。
2、基于RSA算法的免时间戳同步方法建模
设f S、f M分别节点S和节点M的晶振频率,定义节点S相对于节点M的时钟频率偏移
Figure PCTCN2022091457-appb-000328
第一个同步周期,如图7所示,当节点S在
Figure PCTCN2022091457-appb-000329
时刻向节点M发送一个普通数据包,节点M收到该数据包的时间
Figure PCTCN2022091457-appb-000330
可表示为下式:
Figure PCTCN2022091457-appb-000331
其中,
Figure PCTCN2022091457-appb-000332
表示节点S相对于节点M的初始时钟偏移;d 1表示数据包在收发过程中产生的固定时延;X 1表示数据包在收发过程中产生的随机时延。
节点M向节点S发送响应报文的时间
Figure PCTCN2022091457-appb-000333
可根据上述RSA算法获得,即
Figure PCTCN2022091457-appb-000334
节点S收到节点M发送的响应报文的时间记为
Figure PCTCN2022091457-appb-000335
则:
Figure PCTCN2022091457-appb-000336
其中,d 2是响应报文在传输中的固定时延,且d 2=d 1-m,这里m的大小与上行的数据包报文长度和下行的响应报文长度差成正比;Y 1是第一个同步周期中下行链路随机时延。
第二个同步周期,节点S在
Figure PCTCN2022091457-appb-000337
时刻向节点M发送一个普通数据包,节点M在
Figure PCTCN2022091457-appb-000338
时刻收到该数据包,根据RSA算法可得第二个同步周期的响应时间间隔Δ 2的值,此时节点M发送响应报文的时间
Figure PCTCN2022091457-appb-000339
则:
Figure PCTCN2022091457-appb-000340
其中,d 2是数据包在第二个同步同步周期数据包收发过程中产生的固定时延;Y 2第二个同步周期中下行链路随机时延。
第三个同步周期,节点S在
Figure PCTCN2022091457-appb-000341
时刻向节点M发送一个普通数据包,节点M在
Figure PCTCN2022091457-appb-000342
时刻收到该数据包,根据RSA算法可得第三个同步周期的响应时间间隔Δ 3的值,此时节点M发 送响应报文的时间
Figure PCTCN2022091457-appb-000343
则:
Figure PCTCN2022091457-appb-000344
其中,Y 3表示第三个同步周期中下行链路随机时延。
类似地,第i个同步周期,节点S在
Figure PCTCN2022091457-appb-000345
时刻向节点M发送一个普通数据包,节点M在
Figure PCTCN2022091457-appb-000346
时刻收到该数据包,根据RSA算法可得第i个同步周期的响应时间间隔Δ i的值,此时节点M发送响应报文的时间
Figure PCTCN2022091457-appb-000347
节点S在
Figure PCTCN2022091457-appb-000348
时刻受到节点M返回的响应报文。可得,
Figure PCTCN2022091457-appb-000349
Figure PCTCN2022091457-appb-000350
Δ i∈(Δ 2,Δ 3,...,Δ N-1)             (40)
Figure PCTCN2022091457-appb-000351
其中X i表示第i个同步周期数据包在收发过程中产生的随机时延;Y i表示第i个同步周期下行链路随机时延。
N-1个同步周期后,发送节点S可以获得一组本地时间戳对
Figure PCTCN2022091457-appb-000352
接收节点M可以获得一组本地时间戳对
Figure PCTCN2022091457-appb-000353
结合响应时间间隔Δ i,利用最大似然估计方法即可估计出节点间的频率偏移和固定延迟。
用式(39)-(38)可得:
Figure PCTCN2022091457-appb-000354
因为X i和Y i是均值为0,方差为σ 2的高斯随机分布,令V i=X i+Y i,所以V i是均值为0,方差为2σ 2的高斯随机分布。令
Figure PCTCN2022091457-appb-000355
则式(42)可以简化为:
Figure PCTCN2022091457-appb-000356
其中,
Figure PCTCN2022091457-appb-000357
则式(43)可以表示为下列矩阵形式:
Figure PCTCN2022091457-appb-000358
式(44)关于参数(α,d 12)最大似然函数可表示为:
Figure PCTCN2022091457-appb-000359
对式(45)两边取对数有:
Figure PCTCN2022091457-appb-000360
对式(46)求关于α的一阶偏导数有:
Figure PCTCN2022091457-appb-000361
对式(46)求关于d 1的一阶偏导数有:
Figure PCTCN2022091457-appb-000362
令式(47)和式(48)为0,联立求解可得ρ (MS)和d 1的最大似然函数为:
Figure PCTCN2022091457-appb-000363
Figure PCTCN2022091457-appb-000364
所以,响应报文在传输中的固定时延d 2=d 1-m。
在终端侧转换节点和工业无线网络的节点设备完成同步后,TSN网络1的TSN交换机以从时钟状态和终端侧转换节点完成同步,其同步方法和TSN网络时间同步过程一致。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (9)

  1. 一种面向工业无线与TSN融合的跨网时间同步方法,其特征在于:包括以下步骤:
    S1:边界网关的TSN模块以从时钟状态与TSN网络2的TSN交换机进行时钟同步;
    S2:在边界网关内部,TSN模块作为工业无线模块的主时钟,工业无线模块通过串口与TSN模块进行时钟同步;
    S3:工业无线网络中的路由设备以从时钟状态和边界网关的工业无线模块通过信标帧同步方式进行时钟同步,同时路由设备又作为节点设备的主时钟,对节点设备进行时钟同步;
    S4:终端侧转换节点以从时钟状态和工业无线网络的节点设备之间使用免时间戳同步方式进行时钟同步,最后TSN网络1的TSN交换机以从时钟状态和终端侧转换节点进行同步。
  2. 根据权利要求1所述的面向工业无线与TSN融合的跨网时间同步方法,其特征在于:步骤S1中,具体包括以下步骤:
    S11:从时钟利用对等延迟测量机制获取时间戳t1、t2、t3、t4;
    S12:从时钟求取主从时钟的链路延迟;
    S13:从时钟求取处理时间记为time_process;
    S14:将链路延迟和处理时间之和记为Δ1。
  3. 根据权利要求2所述的面向工业无线与TSN融合的跨网时间同步方法,其特征在于:步骤S1中,利用延迟测量机制测量主从时钟间的链路延迟,从时钟和主时钟之间通过交换延迟请求报文Delay_rep和延迟响应报文Delay_resp来获取相应的时间戳,Delay_resp(t2)表示主时钟在t3时刻回复给从时钟的延迟响应报文中携带了t2这个时刻;Follow_up_delay_resp(t3)表示主时钟发送给从时钟的跟随延迟响应报文中携带了t3这个时刻;在Follow_up_delay_resp消息发送后,从时钟接收到时间戳t1、t2、t3、t4;从时钟和主时钟之间的链路延迟delay_resp计算如表达式(1)所示:
    Figure PCTCN2022091457-appb-100001
    处理时间即为本地时间,是主时钟接收时间同步事件消息与主时钟发送下一个后续时间同步事件消息之间的时间间隔的持续时间,处理时间记为time_process,此处即为t3;
    将链路延迟和处理时间之和记为Δ1=delay_resp+time_process。
  4. 根据权利要求1所述的面向工业无线与TSN融合的跨网时间同步方法,其特征在于:边界网关包含一个工业无线模块和一个TSN模块,使用SLIP协议对IP数据包进行封装,通过SLIP串口来保证数据传输的透明性,使得数据在两网之间无缝连接。
  5. 根据权利要求4所述的面向工业无线与TSN融合的跨网时间同步方法,其特征在于:步骤S2中,TSN模块时钟输出的同步脉冲信号为1PPS脉冲,工业无线模块通过检测TSN 模块时钟输出的同步脉冲信号来调整自己的本地时间;包括以下步骤:
    S21:TSN模块将自身的初始时钟T 0通过串口发送给工业无线模块,工业无线模块根据此时间设置自己的本地初始时间local_time0;
    S22:启动中断程序开始检测TSN模块输出的1PPS脉冲信号,工业无线模块每检测到1次1PPS信号,记接收到1PPS信号的时间为local_time1;
    S23:工业无线模块对自己的本地时钟进行一次校正;
    假设工业无线模块在local_time2时刻进行校正,则校正后的时间local_time满足以下表达式:
    local_time0=T 0+S/R      (2)
    local_time=local_time2+N-local_time1      (3)
    其中,S表示串口报文长度,R表示串口波特率,N表示工业无线模块接收到的脉冲个数。
  6. 根据权利要求1所述的面向工业无线与TSN融合的跨网时间同步方法,其特征在于:步骤S3中,将工业无线模块记为A,路由设备记为P 1,节点设备记为P 2
    节点A向节点P 1发送时间消息的时刻到节点P 1接收到该消息的时刻为节点A和节点P 1之间交互的一个同步周期,节点A和节点P 1间的这种过程重复N次产生N个同步周期;在每次两节点进行时间信息的单向交互后,即一个同步周期,节点P 1通过时间戳信息计算出时钟频率偏移进行调整;
    节点P 1调整后,又向节点P 2发送时间消息,节点P 1向节点P 2发送时间消息的时刻到节点P 2接收到该消息的时刻为节点P 1和节点P 2之间交互的一个同步周期,节点P 1和节点P 2间的这种过程重复N次产生N个同步周期;在每次两节点进行时间信息的单向交互后,即一个同步周期,节点P 2通过时间戳信息计算出时钟频率偏移进行调整;
    在节点A和节点P 1的第一个同步周期,节点A将其当前的时间信息
    Figure PCTCN2022091457-appb-100002
    发送给节点P 1,由P 1记录的接收时间
    Figure PCTCN2022091457-appb-100003
    表示为:
    Figure PCTCN2022091457-appb-100004
    式中,
    Figure PCTCN2022091457-appb-100005
    表示节点P 1在t 0时刻相对于参考节点A的初始时钟偏移;
    Figure PCTCN2022091457-appb-100006
    表示节点P 1相对于参考节点A的时钟频率偏移;
    Figure PCTCN2022091457-appb-100007
    表示时间信息收发过程中产生的固定时延;
    Figure PCTCN2022091457-appb-100008
    表示时间信息收发过程中产生的随机时延;
    Figure PCTCN2022091457-appb-100009
    表示时间信息收发过程中由于频率偏移而积累的同步误差;
    进而,在节点A和节点P 1的第一个同步周期参考节点A和待同步节点P 1之间的收发时间 戳的差值
    Figure PCTCN2022091457-appb-100010
    表示为:
    Figure PCTCN2022091457-appb-100011
    Figure PCTCN2022091457-appb-100012
    为P1第一个周期的校正时刻;随后节点P 1
    Figure PCTCN2022091457-appb-100013
    时刻利用差值
    Figure PCTCN2022091457-appb-100014
    来调整自己的本地时钟,P 1调整后的时间用
    Figure PCTCN2022091457-appb-100015
    表示,且有
    Figure PCTCN2022091457-appb-100016
    从初始时间t 0到校正时刻
    Figure PCTCN2022091457-appb-100017
    这段时间内,节点A和节点P 1之间的真实时钟偏差
    Figure PCTCN2022091457-appb-100018
    为:
    Figure PCTCN2022091457-appb-100019
    由公式(7)-(5)得节点P 1调整后新的时钟偏移量
    Figure PCTCN2022091457-appb-100020
    为:
    Figure PCTCN2022091457-appb-100021
    在节点A和节点P 1的第二个同步周期,节点A将其当前的时间信息
    Figure PCTCN2022091457-appb-100022
    发送给节点P 1,节点P 1记录的接收时间
    Figure PCTCN2022091457-appb-100023
    表示为:
    Figure PCTCN2022091457-appb-100024
    其中,
    Figure PCTCN2022091457-appb-100025
    为节点A和节点P 1的第二个同步周期时间信息收发过程中产生的随机时延;
    进而,在节点A和节点P 1的第二个同步周期,参考节点A和节点P 1的收发时间戳的差值
    Figure PCTCN2022091457-appb-100026
    表示为:
    Figure PCTCN2022091457-appb-100027
    设P 1在第一个同步周期调整后的时间为
    Figure PCTCN2022091457-appb-100028
    该时刻对应的节点A的时间为
    Figure PCTCN2022091457-appb-100029
    Figure PCTCN2022091457-appb-100030
    公式(10)变为:
    Figure PCTCN2022091457-appb-100031
    Figure PCTCN2022091457-appb-100032
    Figure PCTCN2022091457-appb-100033
    代入到公式(11)得:
    Figure PCTCN2022091457-appb-100034
    Figure PCTCN2022091457-appb-100035
    Figure PCTCN2022091457-appb-100036
    这段时间内,即节点A和节点P 1的第二个同步周期内,节点A和节点P 1之间的真实时钟偏差
    Figure PCTCN2022091457-appb-100037
    为:
    Figure PCTCN2022091457-appb-100038
    其中
    Figure PCTCN2022091457-appb-100039
    为P1在节点A和节点P 1的第二个周期的校正时刻;
    由公式(13)-(12)得节点P 1在第二次调整后新的时钟偏移量
    Figure PCTCN2022091457-appb-100040
    为:
    Figure PCTCN2022091457-appb-100041
    在节点A和节点P 1的第i个同步周期,得到参考节点A和节点P 1的收发时间戳的差值
    Figure PCTCN2022091457-appb-100042
    为:
    Figure PCTCN2022091457-appb-100043
    其中,
    Figure PCTCN2022091457-appb-100044
    为节点A和节点P 1的第i个周期时间信息收发过程中产生的随机时延;
    节点P 1在第i次调整后新的时钟偏移量
    Figure PCTCN2022091457-appb-100045
    为:
    Figure PCTCN2022091457-appb-100046
    为了使计算更为简化,公式(15)化为如下表达式:
    Figure PCTCN2022091457-appb-100047
    Figure PCTCN2022091457-appb-100048
    Q i是服从高斯分布的随机变量,即Q i~(0,δ 2);使用最大似然估计方法对频率偏移
    Figure PCTCN2022091457-appb-100049
    进行估计,公式(17)基于观测量
    Figure PCTCN2022091457-appb-100050
    的对数似然函数表示为:
    Figure PCTCN2022091457-appb-100051
    Figure PCTCN2022091457-appb-100052
    求偏导并令其为0有:
    Figure PCTCN2022091457-appb-100053
    得其解为:
    Figure PCTCN2022091457-appb-100054
    P 1仅需要获取式(19)中的参数,即可估算出
    Figure PCTCN2022091457-appb-100055
    通过选取适当的
    Figure PCTCN2022091457-appb-100056
    值即可调整
    Figure PCTCN2022091457-appb-100057
    使节点P 1和节点A的频率逐步接近,从而完成P 1和A的时间同步。
  7. 根据权利要求6所述的面向工业无线与TSN融合的跨网时间同步方法,其特征在于:步骤S3中,节点P 1收到来自参考节点A的时间信息后,同样地节点P 2以节点P 1为参考节点,节点P 1向节点P 2发送时间信息;
    在节点P 1和节点P 2的第一个同步周期,节点P 1将其当前的时间信息
    Figure PCTCN2022091457-appb-100058
    发送给节点P 2,节点P 2收到节点P 1发送的时间信息的时间记为
    Figure PCTCN2022091457-appb-100059
    表示为:
    Figure PCTCN2022091457-appb-100060
    式中,
    Figure PCTCN2022091457-appb-100061
    表示节点P 2在t 0时刻相对于参考节点P 1的初始时钟偏移;
    Figure PCTCN2022091457-appb-100062
    表示节点P 2相对于参考节点P 1的时钟频率偏移;
    Figure PCTCN2022091457-appb-100063
    表示时间信息收发过程中产生的固定时延;
    Figure PCTCN2022091457-appb-100064
    表示时间信息收发过程中产生的随机时延;
    Figure PCTCN2022091457-appb-100065
    表示时间信息收发过程中由于频率偏移而积累的同步误差;
    在节点P 1和节点P 2的第一个同步周期节点P 1和节点P 2的收发时间戳差值
    Figure PCTCN2022091457-appb-100066
    为:
    Figure PCTCN2022091457-appb-100067
    Figure PCTCN2022091457-appb-100068
    为P 2节点P 1和节点P 2的第一个同步周期的校正时刻;然后节点P 2
    Figure PCTCN2022091457-appb-100069
    时刻利用时间戳差值
    Figure PCTCN2022091457-appb-100070
    调整本地时钟,调整后的本地时间用
    Figure PCTCN2022091457-appb-100071
    表示,且有
    Figure PCTCN2022091457-appb-100072
    Figure PCTCN2022091457-appb-100073
    从初始时刻到校正时刻这段时间内,节点P 1和节点P 2之间的真实时钟偏差
    Figure PCTCN2022091457-appb-100074
    为:
    Figure PCTCN2022091457-appb-100075
    由公式(22)-(21)得节点P 2在节点P 1和节点P 2的第一个同步周期调整后新的时钟偏移量
    Figure PCTCN2022091457-appb-100076
    为:
    Figure PCTCN2022091457-appb-100077
    在在节点P 1和节点P 2的第二个同步周期,参考节点P 1需要在每个周期内利用节点A和节点P 1在第一个同步周期的收发时间戳的差值
    Figure PCTCN2022091457-appb-100078
    来校正自己的本地时钟,所以对于节点P 2和P 1,其收发时间戳的差值
    Figure PCTCN2022091457-appb-100079
    表示为:
    Figure PCTCN2022091457-appb-100080
    其中
    Figure PCTCN2022091457-appb-100081
    为节点P 1和节点P 2的第二个同步周期节点P 1向节点P 2发送同步消息的发送时间;
    Figure PCTCN2022091457-appb-100082
    为节点P 1和节点P 2的第二个同步周期节点P 2接收到同步消息的接收时间;
    Figure PCTCN2022091457-appb-100083
    为节点P 1和节点P 2的第二个同步周期时间信息收发过程中产生的随机时延;
    Figure PCTCN2022091457-appb-100084
    Figure PCTCN2022091457-appb-100085
    代入上式得在节点P 1和节点P 2的第二个同步周期节点P 1和节点P 2的收发时间戳差值
    Figure PCTCN2022091457-appb-100086
    为:
    Figure PCTCN2022091457-appb-100087
    其中
    Figure PCTCN2022091457-appb-100088
    表示在节点P 1和节点P 2的第一个同步周期节点P 2收到节点P 1发送的时间信息的接收时间;
    Figure PCTCN2022091457-appb-100089
    表示在节点A和节点P 1的第一个同步周期节点P 1收到节点A发送的时间信息的接收时间;
    Figure PCTCN2022091457-appb-100090
    表示在节点A和节点P 1的第一同步个周期节点A向节点P 1发送时间消息的发送时间;
    Figure PCTCN2022091457-appb-100091
    为P 2在节点P 1和节点P 2的第二个周期的校正时刻;然后节点P 2
    Figure PCTCN2022091457-appb-100092
    时刻利用时间戳差值
    Figure PCTCN2022091457-appb-100093
    调整本地时钟,调整后的时间用
    Figure PCTCN2022091457-appb-100094
    表示,且有
    Figure PCTCN2022091457-appb-100095
    Figure PCTCN2022091457-appb-100096
    在节点P 1和节点P 2的第一个同步周期节点P 2
    Figure PCTCN2022091457-appb-100097
    时刻调整后的时间为
    Figure PCTCN2022091457-appb-100098
    此时间对应的时刻到本次校正时刻
    Figure PCTCN2022091457-appb-100099
    这段时间内,节点P 2相较于P 1之间真实的时钟偏差
    Figure PCTCN2022091457-appb-100100
    为:
    Figure PCTCN2022091457-appb-100101
    由公式(26)-(25)得节点P 2在节点P 1和节点P 2的第二个同步周期调整时间后的时钟偏移量
    Figure PCTCN2022091457-appb-100102
    为:
    Figure PCTCN2022091457-appb-100103
    其中
    Figure PCTCN2022091457-appb-100104
    为节点P 1和节点P 2的第二个同步周期节点P 1向节点P 2发送同步消息的发送时间;
    Figure PCTCN2022091457-appb-100105
    表示在节点A和节点P 1的第一个同步周期节点A和节点P 1之间的收发时间戳的差值;
    在节点P 1和节点P 2的第三个同步周期,参考节点P 1需要在每个周期内利用节点A和节点P 1在第二个同步周期的收发时间戳的差值
    Figure PCTCN2022091457-appb-100106
    来校正自己的本地时钟,所以对于节点P 2和P 1,在节点P 1和节点P 2的第三个同步周期其收发时间戳的差值
    Figure PCTCN2022091457-appb-100107
    表示为:
    Figure PCTCN2022091457-appb-100108
    Figure PCTCN2022091457-appb-100109
    Figure PCTCN2022091457-appb-100110
    代入上式得在节点P 1和节点P 2的第三个同步周期节点P 1和节点P 2的收发时间戳差值
    Figure PCTCN2022091457-appb-100111
    为:
    Figure PCTCN2022091457-appb-100112
    其中
    Figure PCTCN2022091457-appb-100113
    表示在节点P 1和节点P 2的第二个同步周期节点P 2收到节点P 1发送的时间信息的接收时间;
    Figure PCTCN2022091457-appb-100114
    表示在节点A和节点P 1的第二个同步周期节点P 1收到节点A发送的时间信息的接收时间;
    Figure PCTCN2022091457-appb-100115
    表示在节点A和节点P 1的第二个同步周期节点A向节点P 1发送时间消息的发送时间;
    Figure PCTCN2022091457-appb-100116
    表示在节点A和节点P 1的第一个同步周期节点P 1收到节点A发送的时间信息的接收时间;
    Figure PCTCN2022091457-appb-100117
    表示在节点A和节点P 1的第一个同步周期节点A向节点P 1发送时间消息的发送时间;
    对于节点P 1和节点P 2的第i个同步周期,参考节点P 1需要在每个周期内利用点A和节点P 1在第i个同步周期的收发时间戳的差值
    Figure PCTCN2022091457-appb-100118
    来校正自己的本地时钟,所以对于节点P 2和P 1,其收发时间戳的差值
    Figure PCTCN2022091457-appb-100119
    表示为:
    Figure PCTCN2022091457-appb-100120
    Figure PCTCN2022091457-appb-100121
    其中
    Figure PCTCN2022091457-appb-100122
    表示为节点P 1和节点P 2的第i个同步周期节点P 1向节点P 2发送同步消息的发送时间;
    Figure PCTCN2022091457-appb-100123
    表示在节点P 1和节点P 2的第i-1个同步周期节点P 2收到节点P 1发送的时间信息的接收时间;
    Figure PCTCN2022091457-appb-100124
    表示在节点P 1和节点P 2的第i个同步周期时间信息收发过程中产生的随机时延;
    Figure PCTCN2022091457-appb-100125
    表示在节点A和节点P 1的第i-1个同步周期节点P 1收到节点A发送的时间信息的接收时间;
    Figure PCTCN2022091457-appb-100126
    表示在节点A和节点P 1的第i-1个同步周期节点A向节点P 1发送时间消息的发送时间;
    Figure PCTCN2022091457-appb-100127
    表示为节点P 1和节点P 2的第i-j-1个同步周期节点P 1向节点P 2发送同步消息的发送时间;
    Figure PCTCN2022091457-appb-100128
    表示在节点A和节点P 1的第i-j-1个同步周期节点A向节点P 1发送时间消息的发送时间;
    使用最大似然估计方法对频率偏移
    Figure PCTCN2022091457-appb-100129
    进行估计,公式(30)基于观测量
    Figure PCTCN2022091457-appb-100130
    的对数似然函数表示为:
    Figure PCTCN2022091457-appb-100131
    Figure PCTCN2022091457-appb-100132
    式(31)两边对
    Figure PCTCN2022091457-appb-100133
    求偏导有:
    Figure PCTCN2022091457-appb-100134
    令式(32)等于0得其解为:
    Figure PCTCN2022091457-appb-100135
    P 2仅需要获取式(33)中的参数,即可估算出
    Figure PCTCN2022091457-appb-100136
    通过选取适当的
    Figure PCTCN2022091457-appb-100137
    值即可调整
    Figure PCTCN2022091457-appb-100138
    使节点P 2和节点P 1的频率逐步接近,从而完成P 2和P 1的时间同步。
  8. 根据权利要求1所述的面向工业无线与TSN融合的跨网时间同步方法,其特征在于:步骤S4中,以工业无线网络的节点设备为主时钟,记为M,终端侧转换节点作为从时钟,记为S,节点S周期性地向节点M发送普通数据包,通过控制接收方对发送方的响应时间间隔Δ i来隐式地传递同步信息;
    在第i个周期发送节点S在
    Figure PCTCN2022091457-appb-100139
    时刻向接收节点M发送一个普通数据包,接收节点M在
    Figure PCTCN2022091457-appb-100140
    时刻收到该数据包,并且在等待响应时间间隔Δ i后,在
    Figure PCTCN2022091457-appb-100141
    时刻向发送节点S返回一个响应报文,发送节点S在
    Figure PCTCN2022091457-appb-100142
    时刻收到该响应报文,随后记录下收发的本地时刻并存储在本地缓存中,伴随着节点数据包收发过程的重复进行,节点可利用本地缓存的收发时间估计出时钟频率偏移。
  9. 根据权利要求8所述的面向工业无线与TSN融合的跨网时间同步方法,其特征在于:步骤S4中的时间同步次数设为N,有A∈{2,3,4,...N-1},任取素数P、Q、D,其中N=P*Q,当(D*E)mod((P-1)*(Q-1))=1成立,则有C i=(A i EXP D)mod N,其中1<i<N,则C i就是生成的随机数,同时用C i表示第i次同步过程中接收节点需等待的响应时间间隔Δ i
    设f S、f M分别是节点S和节点M的晶振频率,定义节点S相对于节点M的时钟频率偏移
    Figure PCTCN2022091457-appb-100143
    第一个同步周期,当节点S在
    Figure PCTCN2022091457-appb-100144
    时刻向节点M发送一个普通数据包,节点M收到该数据包的时间
    Figure PCTCN2022091457-appb-100145
    表示为下式:
    Figure PCTCN2022091457-appb-100146
    其中,
    Figure PCTCN2022091457-appb-100147
    表示节点S相对于节点M的初始时钟偏移;d 1表示数据包在收发过程中产生的固定时延;X 1表示数据包在收发过程中产生的随机时延;
    节点M发送响应报文的时间
    Figure PCTCN2022091457-appb-100148
    可根据所述RSA算法获得,即
    Figure PCTCN2022091457-appb-100149
    节点S收到节点M发送的响应报文的时间记为
    Figure PCTCN2022091457-appb-100150
    则:
    Figure PCTCN2022091457-appb-100151
    其中,d 2是响应报文在传输中的固定时延,且d 2=d 1-m,这里m的大小与上行的数据包报文长度和下行的响应报文长度差成正比;Y 1是第一个同步周期中下行链路随机时延;
    第二个同步周期,节点S在
    Figure PCTCN2022091457-appb-100152
    时刻向节点M发送一个普通数据包,节点M在
    Figure PCTCN2022091457-appb-100153
    时刻收到该数据包,根据RSA算法得到第二个同步周期的响应时间间隔Δ 2的值,此时节点M发送响应报文的时间
    Figure PCTCN2022091457-appb-100154
    则:
    Figure PCTCN2022091457-appb-100155
    其中,d 2是数据包在第二个同步周期数据包收发过程中产生的固定时延;Y 2第二个同步周期中下行链路随机时延;
    第三个同步周期,节点S在
    Figure PCTCN2022091457-appb-100156
    时刻向节点M发送一个普通数据包,节点M在
    Figure PCTCN2022091457-appb-100157
    时刻收到该数据包,根据RSA算法得到第三个同步周期的响应时间间隔Δ 3的值,此时节点M发送响应报文的时间
    Figure PCTCN2022091457-appb-100158
    则:
    Figure PCTCN2022091457-appb-100159
    其中,Y 3表示第三个同步周期中下行链路随机时延;
    第i个同步周期,节点S在
    Figure PCTCN2022091457-appb-100160
    时刻向节点M发送一个普通数据包,节点M在
    Figure PCTCN2022091457-appb-100161
    时刻收到该数据包,根据RSA算法可得第i个周期的响应时间间隔Δ i的值,此时节点M发送响应报文的时间
    Figure PCTCN2022091457-appb-100162
    节点S在
    Figure PCTCN2022091457-appb-100163
    时刻受到节点M返回的响应报文,可得,
    Figure PCTCN2022091457-appb-100164
    Figure PCTCN2022091457-appb-100165
    Δ i∈(Δ 2,Δ 3,...,Δ N-1)       (40)
    Figure PCTCN2022091457-appb-100166
    N-1个同步周期后,发送节点S获得一组本地时间戳对
    Figure PCTCN2022091457-appb-100167
    接收节点M获得一组本地时间戳对
    Figure PCTCN2022091457-appb-100168
    结合响应时间间隔Δ i,利用最大似然估计方法估计出节点间的频率偏移和固定延迟;
    用式(39)-(38)可得:
    Figure PCTCN2022091457-appb-100169
    因为X i和Y i是均值为0,方差为σ 2的高斯随机分布,令V i=X i+Y i,所以V i是均值为0,方差为2σ 2的高斯随机分布,令
    Figure PCTCN2022091457-appb-100170
    则式(42)简化为:
    Figure PCTCN2022091457-appb-100171
    其中,
    Figure PCTCN2022091457-appb-100172
    则式(43)表示为下列矩阵形式:
    Figure PCTCN2022091457-appb-100173
    式(44)关于参数(α,d 12)最大似然函数表示为:
    Figure PCTCN2022091457-appb-100174
    对式(45)两边取对数有:
    Figure PCTCN2022091457-appb-100175
    对式(46)求关于σ的一阶偏导数有:
    Figure PCTCN2022091457-appb-100176
    对式(46)求关于d 1的一阶偏导数有:
    Figure PCTCN2022091457-appb-100177
    令式(47)和式(48)为0,联立求解可得ρ (MS)和d 1的最大似然函数为:
    Figure PCTCN2022091457-appb-100178
    Figure PCTCN2022091457-appb-100179
    所以,响应报文在传输中的固定时延d 2=d 1-m;
    在终端侧转换节点和工业无线网络的节点设备完成同步后,TSN网络1的TSN交换机以从时钟状态和终端侧转换节点完成同步,其同步方法和步骤S1一致。
PCT/CN2022/091457 2021-11-08 2022-05-07 一种面向工业无线与tsn融合的跨网时间同步方法 WO2023077760A1 (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202111323593.1A CN114050884B (zh) 2021-11-08 2021-11-08 一种面向工业无线与tsn融合的跨网时间同步方法
CN202111323593.1 2021-11-08

Publications (1)

Publication Number Publication Date
WO2023077760A1 true WO2023077760A1 (zh) 2023-05-11

Family

ID=80207907

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/091457 WO2023077760A1 (zh) 2021-11-08 2022-05-07 一种面向工业无线与tsn融合的跨网时间同步方法

Country Status (2)

Country Link
CN (1) CN114050884B (zh)
WO (1) WO2023077760A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116456450A (zh) * 2023-06-13 2023-07-18 北京科技大学 一种EtherCAT与5G融合组网的时间同步方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114050884B (zh) * 2021-11-08 2023-05-12 重庆邮电大学 一种面向工业无线与tsn融合的跨网时间同步方法
CN114513274B (zh) * 2022-04-15 2022-07-26 清华大学 面向互联网的多跳高精度时间同步协议实现方法及装置
CN115567139B (zh) * 2022-09-19 2024-04-12 重庆邮电大学 一种面向5g与tsn融合的跨网时间同步方法
CN117478292B (zh) * 2023-12-27 2024-04-02 中国科学院微电子研究所 一种频率同步网辅助时间同步网的时间同步方法及设备
CN117560109B (zh) * 2024-01-11 2024-04-12 鹏城实验室 时间同步方法、设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110943899A (zh) * 2019-12-13 2020-03-31 重庆邮电大学 一种epa工业总线与时间敏感网络适配***及方法
US20200280383A1 (en) * 2019-02-28 2020-09-03 Nxp B.V. Method and system for merging clocks from multiple precision time protocol (ptp) clock domains
CN112637945A (zh) * 2020-12-17 2021-04-09 江苏省未来网络创新研究院 一种工业互联网络时钟同步的方法和***
CN113411153A (zh) * 2021-06-07 2021-09-17 中国人民解放军国防科技大学 一种用于时间敏感网络的集中式控制时间同步方法
CN114050884A (zh) * 2021-11-08 2022-02-15 重庆邮电大学 一种面向工业无线与tsn融合的跨网时间同步方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8775681B2 (en) * 2011-04-27 2014-07-08 GM Global Technology Operations LLC Cross-network synchronization of application S/W execution using flexray global time
CN103580846B (zh) * 2013-08-23 2017-03-01 北京东土科技股份有限公司 一种跨非1588网络传输精密时钟报文的方法及***
CN106452650B (zh) * 2016-12-01 2019-02-22 重庆邮电大学 适用于多跳无线传感器网络的时钟同步频率偏移估计方法
CN107483134B (zh) * 2017-06-15 2019-04-16 西安微电子技术研究所 一种不同网络之间的跨网时钟同步通信装置及方法
CN107888314B (zh) * 2017-11-24 2019-03-08 重庆邮电大学 一种面向工业物联网的跨网络时间同步方法
CN109936519B (zh) * 2019-01-18 2021-06-04 重庆邮电大学 基于定时响应的工业无线传感网时钟频率偏移估计方法
US20210328900A1 (en) * 2020-04-20 2021-10-21 Mellanox Technologies, Ltd. Time-Synchronization Testing in a Network Element
CN111556559B (zh) * 2020-05-09 2021-11-26 重庆邮电大学 基于免时间戳交互与单向消息传播的混合时钟同步方法
CN112904705B (zh) * 2021-01-22 2022-05-17 重庆邮电大学 一种低轨小卫星之间的分级时钟同步方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200280383A1 (en) * 2019-02-28 2020-09-03 Nxp B.V. Method and system for merging clocks from multiple precision time protocol (ptp) clock domains
CN110943899A (zh) * 2019-12-13 2020-03-31 重庆邮电大学 一种epa工业总线与时间敏感网络适配***及方法
CN112637945A (zh) * 2020-12-17 2021-04-09 江苏省未来网络创新研究院 一种工业互联网络时钟同步的方法和***
CN113411153A (zh) * 2021-06-07 2021-09-17 中国人民解放军国防科技大学 一种用于时间敏感网络的集中式控制时间同步方法
CN114050884A (zh) * 2021-11-08 2022-02-15 重庆邮电大学 一种面向工业无线与tsn融合的跨网时间同步方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116456450A (zh) * 2023-06-13 2023-07-18 北京科技大学 一种EtherCAT与5G融合组网的时间同步方法
CN116456450B (zh) * 2023-06-13 2023-10-10 北京科技大学 一种EtherCAT与5G融合组网的时间同步方法

Also Published As

Publication number Publication date
CN114050884A (zh) 2022-02-15
CN114050884B (zh) 2023-05-12

Similar Documents

Publication Publication Date Title
WO2023077760A1 (zh) 一种面向工业无线与tsn融合的跨网时间同步方法
US10056999B2 (en) Time synchronization error compensation method for multi-hop wireless backhaul network based on PTP
CN109565435B (zh) 用于网络同步的***和方法
CN106452650B (zh) 适用于多跳无线传感器网络的时钟同步频率偏移估计方法
CN102244603B (zh) 传输承载时间的报文的方法、设备及***
Gong et al. CESP: A low-power high-accuracy time synchronization protocol
Ren et al. Self-correcting time synchronization using reference broadcast in wireless sensor network
CN103763055A (zh) 一种精确同步时间的方法
CN113055117A (zh) 一种无线分布式网络的时钟同步装置及方法
CN104754722A (zh) 一种面向层次化异构网络的时间同步方法
CN105188126A (zh) 基于平均场的分布式多跳无线网络时钟同步方法
Krummacker et al. Intra-network clock synchronization for wireless networks: From state of the art systems to an improved solution
Pinho et al. A robust gradient clock synchronization algorithm for wireless sensor networks
CN104243079A (zh) 一种实时以太网的微秒级时钟同步方法
CN104836654B (zh) 一种基于Ethernet POWERLINK的时钟同步方法
CN111682919A (zh) 一种基于时间敏感网络的井下传感设备时间同步方法
WO2018098791A1 (zh) 适用于多跳无线传感器网络的时钟同步频率偏移估计方法
CN114584246A (zh) 基于fpga的无线自组网时间同步方法、***及介质
CN114389735A (zh) 一种基于ieee802.1as冗余主时钟的时钟同步方法
KR100999686B1 (ko) 하이브리드 네트워크를 위한 실시간 동기화 방법
Wang et al. A survey of time synchronization algorithms for wireless sensor networks
CN107197514A (zh) 无线传感器网络的可认知全局时钟同步方法及在多跳网络中的应用
US20240243829A1 (en) Cross-network time synchronization method for industrial wireless network and tsn fusion
Gautam et al. A comparative study of time synchronization protocols in wireless sensor networks
TWI421667B (zh) 時鐘同步之方法及應用該方法之網路系統

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 18008268

Country of ref document: US

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22888797

Country of ref document: EP

Kind code of ref document: A1