A flow probe high-precision time detection system, method, device, medium
By employing a dual detection method combining a frequency counter and a B-code signal source, the error problem in flow probe time detection is resolved, achieving a balance between high precision and stability, and ensuring the accuracy and reliability of flow monitoring results.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2022-10-19
- Publication Date
- 2026-04-07
AI Technical Summary
Existing flow probes have errors in high-precision time detection, making it impossible to achieve a balance between stability and accuracy, resulting in inaccurate flow monitoring results.
A dual detection method using a frequency counter and a B-code signal source is employed. By monitoring the error value and error percentage between the flow probe and the B-code signal source in two dimensions, the reliability of high-precision time detection is ensured.
It achieves unified detection of high precision and stability of flow probe time, ensuring the accuracy and reliability of flow monitoring results.
Smart Images

Figure CN115643183B_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the field of network operation and traffic control, and in particular to a high-precision time detection system, method, device and medium for a traffic probe. BACKGROUND
[0002] With the development of business diversification, business integration, network openness and terminal intelligence, the network coverage field has expanded dramatically, new applications have emerged in large numbers, network traffic has surged, and network transmission bandwidth has expanded from hundreds / kilograms to megabytes. However, the current monitoring capability of network traffic is limited to below kilobytes, and the traffic monitoring of large bandwidth transmission links is still in the blank stage. Moreover, the existing technology has not yet reached the fine-grained resources of diversified applications and important businesses. In addition, the traffic monitoring and detection devices currently deployed in the network can only target a single node and local information, and the information obtained from the network is relatively isolated. The correlation analysis and comprehensive presentation of multi-point information along the link are still lacking, and the full-range and full-domain monitoring of traffic cannot be achieved. In terms of rapid positioning of network faults and comprehensive evaluation of network operation status, it is still impossible to provide comprehensive and real-time data basis. At the same time, there is no unified specification suitable for the application characteristics of the business network in the field of traffic monitoring. Various network monitoring methods are independently constructed and have not formed a combined force. Therefore, the research on large bandwidth full-domain traffic sensing and monitoring technology and the formulation of network traffic monitoring specifications are imminent.
[0003] Internationally, the IP performance Metrics Working Group (IPPMWG) of IETF has proposed principles and overall frameworks for defining performance metrics, and has defined some indicators for evaluating the quality, performance and reliability of IP network data transmission services, such as connectivity, one-way packet loss, one-way delay, loopback delay, and other indicators are in the process of standardization, but the specific implementation methods and applications are not within the definition range of IPPMWG.
[0004] ITU-T Study Group 13 (SG13) has also proposed Y.1540 (original L380) recommendations, which define four parameters for measuring IP packet transmission performance on IP networks: speed, accuracy, reliability, and availability. In addition, the y.1541 (original L381) recommendation specifies IP performance and availability indicators and allocation, and classifies IP services into six categories according to QoS.
[0005] Meanwhile, other international organizations have proposed several testing infrastructures. For example, Surveyor, a network testing infrastructure based on the IPPMWG standard proposed by Advanced Network & Services in collaboration with other organizations, can measure the path performance of the Internet between participating organizations; the project also proposed methods and tools for analyzing performance data. MMI, an NSF-initiated and DARPA-funded project, proposed a probe-based distributed, scalable, and dynamic network testing infrastructure. Other projects such as Ripe, AMP, and PingER are also related to network testing.
[0006] In the business network, in order to ensure that all business systems and various business devices maintain a unified time standard, time synchronization devices are widely deployed. These devices obtain the correct real-time clock (absolute time) from multiple clock sources such as Beidou navigation satellites, atomic clocks, and time stations. After encoding, they provide time synchronization signals to the business devices. This encoding adopts the internationally recognized time code standard IRIG-B code encoding format (English name: Inter-Range Instrumentation Group-B). Such time synchronization devices are referred to as "B code signal sources", and the time synchronization signals they emit are referred to as "B code signals".
[0007] A traffic probe is a dedicated device for monitoring network performance metrics of a service network. It captures all data packets flowing through the device from the mirror port of a network switch or router, and obtains performance metrics such as traffic, latency, and packet loss by parsing the IP protocol. The traffic probe simultaneously obtains a time synchronization signal from a B-code signal source device, generating a real-time clock internally to timestamp each captured data packet. Most of the network performance metrics obtained by the traffic probe are time-related; for example, traffic represents the total number of bits transmitted per second on the line. If the time accuracy of the traffic probe is insufficient or inaccurate, the internal time length of one second will be unstable, fluctuating between fast and slow. This will cause data packets from the previous second to be included in the next second, resulting in an underestimation of traffic in the previous second and an overestimation in the next. Therefore, the traffic probe must use an accurate, high-precision real-time time as the basis for performance testing.
[0008] How to detect the high-precision time inside the flow probe is an important technical problem. The detection should include two aspects: (1) Time stability, that is, the flow probe must maintain a stable second value for a considerable period of time, and the interval of each second must be the same. For example, the second value must remain unchanged for 24 consecutive hours, which requires detecting the error between the flow probe and the B code signal source. (2) Time accuracy, that is, the timestamps given by the flow probe are correct, and the calculation results must also be correct and reliable. Summary of the Invention
[0009] The technical problem to be solved by this invention is to determine the accuracy and precision level of time in a high-precision data acquisition card. High-precision time is the basis for flow probe to monitor flow. Whether it is timestamp or flow monitoring index calculation, it is inseparable from a high-precision clock. Therefore, this invention proposes a high-precision time detection system, method, device and medium for flow probe.
[0010] According to one aspect of the present invention, a high-precision time detection system for a flow probe is provided, comprising: a flow probe, a first frequency counter, a second frequency counter, and a B-code signal source, wherein:
[0011] The B-code signal output terminal of the B-code signal source is connected to the B-code signal input terminal of the flow probe;
[0012] The 1PPS signal output terminal of the B code signal source and the 1PPS signal output terminal of the flow probe are respectively connected to channel A and channel B of the first frequency counter.
[0013] The 1MHz signal output terminal of the B-code signal source and the 1MHz signal output terminal of the flow probe are respectively connected to channel A and channel B of the second frequency counter.
[0014] The first frequency counter accumulates and counts the two input pulse signals respectively. The count value of the first frequency counter is incremented by one for each rising edge detected. The count is continuously recorded for a first preset duration to obtain the number of pulses for each of the two input pulse signals.
[0015] The second frequency counter accumulates the counts for the two input square wave signals respectively. The count value of the second frequency counter is incremented by one for each rising edge detected, and the count is continuously recorded for a second preset duration to obtain the number of square waves for each of the two input square wave signals.
[0016] According to another aspect of the present invention, a high-precision time detection method for flow probes is proposed, comprising:
[0017] Step S101: Connect the B code time synchronization signal of the B code signal source to the B code signal input port of the flow probe high-precision data acquisition card; connect the 1PPS signal-S output by the B code signal source to channel A of the first frequency counter; connect the 1PPS signal-D output by the flow probe high-precision data acquisition card to channel B of the first frequency counter.
[0018] Step S102: Start the first frequency counter to accumulate and count the input pulse signals received by channel A and channel B, and continuously record for a first preset time to obtain a first count record. The first count record includes a count value t1 obtained by accumulating and counting the input pulse signals received by channel A, and a count value t2 obtained by accumulating and counting the input pulse signals received by channel B.
[0019] Step S103: Calculate the error value and error percentage between the corresponding count values of channel A and channel B based on the first counting record;
[0020] Step S104: Connect the 1MHz signal-S output from the B code signal source to channel A of the second frequency counter, and connect the 1MHz signal-D output from the flow probe high-precision data acquisition card to channel B of the second frequency counter.
[0021] Step S105: Start the second frequency counter to accumulate and count the input square wave signals received by channel A and channel B, and continuously record for a second preset time to obtain a second count record. The second count record includes the count value t3 obtained by accumulating and counting the input square wave signals received by channel A, and the count value t4 obtained by accumulating and counting the input square wave signals received by channel B.
[0022] Step S106: Calculate the error value and error percentage between the corresponding count values of channel A and channel B based on the second counting record;
[0023] Step S107: Determine the time accuracy and time stability of the flow probe based on the error value dt1, error percentage rt1, error value dt2, and error percentage rt2.
[0024] In step S103, the error value is calculated according to the following formula:
[0025] dt1 = t2 - t1;
[0026] The error percentage is calculated according to the following formula:
[0027] rt1 = (t2 - t1) / t1.
[0028] In step S106, the error value is calculated according to the following formula:
[0029] dt2 = t4 - t3;
[0030] The error percentage is calculated according to the following formula:
[0031] rt2 = (t4 - t3) / t3.
[0032] In step S107, the magnitudes of error value dt1 and error percentage rt1 represent the stability of high-precision time; the smaller the value, the higher the stability. The magnitudes of error value dt2 and error percentage rt2 represent the accuracy of time; the smaller the value, the higher the accuracy.
[0033] In one embodiment of the present invention, step S107 further includes:
[0034] The accuracy level of the high-precision data acquisition card for the flow probe is determined based on the error percentage rt2.
[0035] In one embodiment of the present invention, determining the accuracy level of the high-precision data acquisition card for the flow probe based on the error percentage rt2 includes:
[0036] Based on the magnitude of the error percentage rt2, determine whether the high-precision data acquisition card for the flow probe is of the same precision level as the B-code signal source. If the precision levels are the same, select a higher-level B-code signal source and repeat steps S101 to S106 until the precision level of the high-precision data acquisition card for the flow probe is lower than that of the B-code signal source. This will determine the precision level of the high-precision data acquisition card for the flow probe.
[0037] According to another aspect of the present invention, an electronic device is provided, including a memory and a processor, the memory being used to store one or more computer instructions, wherein the one or more computer instructions are executed by the processor to implement the above-described method steps.
[0038] According to another aspect of the present invention, a computer-readable storage medium is provided having computer instructions stored thereon, wherein the computer instructions, when executed by a processor, implement the above-described method steps.
[0039] The technical solutions provided by the embodiments of the present invention may include the following beneficial effects:
[0040] The advantages and innovations of this invention compared to the prior art are as follows:
[0041] 1. High-precision time detection is performed using a frequency counter and a B-code signal source. The detection error accuracy is directly related to the frequency counter and the B-code signal source. The high precision of the frequency counter and the B-code signal source has national quality certification. Therefore, the high-precision time detection method proposed in this invention has high reliability.
[0042] 2. High reliability of high-precision time detection is ensured through dual-dimensional monitoring, namely, monitoring the accuracy of high-precision time measurement and the stability of high-precision time testing.
[0043] 3. The timing accuracy and time precision of the flow probe are tested using a frequency counter. Two test channels are designed for accuracy and stability. By comparing the frequency changes and count values of the two channels, the error range of the two channels is determined. Then, it is determined whether the time precision of the flow probe is the same as that of the clock source, and finally, the time precision level of the flow probe is determined.
[0044] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and do not limit the embodiments of the present invention. Attached Figure Description
[0045] Other features, objects, and advantages of embodiments of the present invention will become more apparent from the following detailed description of non-limiting implementations, taken in conjunction with the accompanying drawings. In the drawings:
[0046] Figure 1 A schematic diagram of a typical service network deployment scenario for a traffic probe according to an embodiment of the present invention is shown.
[0047] Figure 2 A schematic diagram of the internal structure and main interfaces of a flow probe according to an embodiment of the present invention is shown.
[0048] Figure 3 A schematic diagram of a high-precision time detection system for flow probes according to an embodiment of the present invention is shown.
[0049] Figure 4 A signal connection diagram is shown for a high-precision time detection method for flow probes according to an embodiment of the present invention.
[0050] Figure 5 A flowchart of a high-precision time detection method for flow probes according to an embodiment of the present invention is shown. Detailed Implementation
[0051] Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings to enable those skilled in the art to readily implement them. Furthermore, for clarity, portions unrelated to the description of the exemplary embodiments have been omitted from the drawings.
[0052] In embodiments of the present invention, it should be understood that terms such as “comprising” or “having” are intended to indicate the presence of features, numbers, steps, behaviors, components, portions or combinations thereof disclosed in this specification, and are not intended to exclude the possibility that one or more other features, numbers, steps, behaviors, components, portions or combinations thereof are present or added.
[0053] It should also be noted that, unless otherwise specified, the embodiments and features described in the present invention can be combined with each other. The embodiments of the present invention will now be described in detail with reference to the accompanying drawings and examples.
[0054] Figure 1 This diagram illustrates a typical service network deployment scenario for a traffic probe according to an embodiment of the present invention. Figure 1 As shown, the service network deploys: Application System A device, Application System B device, switch, multiple B-code signal source devices, and traffic probes, among which:
[0055] Application system A and application system B are connected to the switch to enable network connectivity and communication, and to transmit network data packets.
[0056] The traffic probe is connected to the switch and one of the multiple B-code signal source devices to capture all network packets from the mirror port of the switch and obtain B-code signals from the B-code signal source devices to monitor network performance indicators.
[0057] The remaining B-code signal source devices are connected to the application system device nearby and send time synchronization signals in IRIG-B code encoding format, i.e. B-code signals, to the application system device to ensure the time consistency of the entire system.
[0058] Figure 2 A schematic diagram of the internal structure and main interfaces of a flow probe according to an embodiment of the present invention is shown, as follows. Figure 2 As shown, the flow probe is actually a server equipped with a high-precision data acquisition card. It includes two input signals: B code signal and network data packet, and two output signals: 1PPS signal and 1MHz signal. The 1PPS signal means that one pulse signal is output every second, and the 1MHz signal means that one million square wave signals are output every second. The specific internal structure of the flow probe will not be described in detail in this invention.
[0059] Figure 3 A schematic diagram of a high-precision time detection system for flow probes according to an embodiment of the present invention is shown, as follows: Figure 3 As shown, the high-precision time detection system for the flow probe includes: a flow probe, a first frequency counter, a second frequency counter, and a B-code signal source, wherein:
[0060] The B-code signal output terminal of the B-code signal source is connected to the B-code signal input terminal of the flow probe;
[0061] The 1PPS signal output terminal of the B code signal source and the 1PPS signal output terminal of the flow probe are respectively connected to channel A and channel B of the first frequency counter.
[0062] The 1MHz signal output terminal of the B-code signal source and the 1MHz signal output terminal of the flow probe are respectively connected to channel A and channel B of the second frequency counter;
[0063] The first frequency counter accumulates and counts the two input pulse signals respectively. The count value of the first frequency counter is incremented by one for each rising edge detected. The count is continuously recorded for a first preset time, such as 1 day or 24 hours, to obtain the number of pulses for each of the two input pulse signals. The average count value of a 1PPS signal is 1 pulse per second. The theoretical value of the pulse count value in 1 day is 86400. By analyzing the error and error distribution of the corresponding pulse count values of the two input pulse signals, the stability of the flow probe time can be obtained.
[0064] The second frequency counter accumulates and counts the two input square wave signals separately. The count value of the second frequency counter is incremented by one for each rising edge detected. On average, a 1MHz signal can obtain 1 million square wave counts per second, which is a 7-digit number. Due to the limitation of the number of bits of the frequency counter (the maximum number of bits of commonly used counters is 10-11), it can only record continuously for several hours to obtain the number of square waves of each of the two input square wave signals. Taking a second preset duration of 3 hours as an example, after 3 hours, the theoretical value of the square wave count value of the second frequency counter is 10,800,000,000, just exceeding 10 digits. By analyzing the error and error distribution of the corresponding square wave count values of the two input square wave signals, the accuracy of the flow probe time can be obtained.
[0065] In the detection scheme of this invention, the longer the detection time, the more errors accumulate, the easier the detection becomes, and the more accurate the detection results. Since the flow probe is a device customized for the service network, its internal design has resulted in a relatively small error. The detection results of this invention's detection scheme are the error between the flow probe and the B-code signal source, and the percentage of error. The resulting accuracy level depends on the accuracy level of the B-code signal source. Without considering cost or other factors, a higher-precision B-code signal source can be used to perform the same detection on the flow probe and high-precision data acquisition card. Alternatively, a higher-precision frequency counter with more bits can be used, resulting in a longer detection time and a more accurate accuracy level.
[0066] Figure 4 This diagram illustrates the signal connection of a high-precision time detection method for flow probes according to an embodiment of the present invention. Figure 5 A flowchart illustrating a high-precision time detection method for a flow probe according to an embodiment of the present invention is shown. The method detects a flow probe, and the detection parameters are time accuracy and time stability. Figure 5 As shown, the high-precision time detection method for the flow probe includes the following steps S101-S107:
[0067] In step S101, the B code time synchronization signal of the B code signal source is connected to the B code signal input port of the flow probe high-precision data acquisition card; the 1PPS signal-S output by the B code signal source is connected to channel A of the first frequency counter; and the 1PPS signal-D output by the flow probe high-precision data acquisition card is connected to channel B of the first frequency counter.
[0068] In step S102, the first frequency counter is started to accumulate and count the input pulse signals received by channel A and channel B, and continuously record for a first preset time, such as 1 day, to obtain the first count record. The first count record includes the count value t1 obtained by accumulating and counting the input pulse signals received by channel A, and the count value t2 obtained by accumulating and counting the input pulse signals received by channel B.
[0069] In step S103, the error value and error percentage between the corresponding count values of channel A and channel B are calculated based on the first counting record;
[0070] The error value can be calculated using the following formula:
[0071] dt1 = t2 - t1;
[0072] The error percentage can be calculated using the following formula:
[0073] rt1 = (t2 - t1) / t1;
[0074] In step S104, the 1MHz signal-S output from the B code signal source is connected to channel A of the second frequency counter, and the 1MHz signal-D output from the flow probe high-precision data acquisition card is connected to channel B of the second frequency counter.
[0075] In step S105, the second frequency counter is started to accumulate and count the input square wave signals received by channel A and channel B, and continuously record for a second preset time period, such as 4 hours, to obtain a second count record. The second count record includes the count value t3 obtained by accumulating and counting the input square wave signals received by channel A, and the count value t4 obtained by accumulating and counting the input square wave signals received by channel B.
[0076] In step S106, the error value and error percentage between the corresponding count values of channel A and channel B are calculated based on the second counting record;
[0077] The error value can be calculated using the following formula:
[0078] dt2 = t4 - t3;
[0079] The error percentage can be calculated using the following formula:
[0080] rt2 = (t4 - t3) / t3;
[0081] In step S107, the time accuracy and time stability of the flow probe are determined based on the error value dt1, the error percentage rt1, the error value dt2, and the error percentage rt2.
[0082] In this process, the magnitudes of the error value dt1 and the error percentage rt1 represent the stability of the high-precision time measurement; the smaller the value, the higher the stability. The magnitudes of the error value dt2 and the error percentage rt2 represent the accuracy of the time measurement; the smaller the value, the higher the accuracy. Simultaneously, the magnitude of the error percentage rt2 determines whether the high-precision data acquisition card for the flow probe is at the same accuracy level as the B-code signal source. If the accuracy levels are the same, a higher-level B-code signal source can be selected to repeat steps S101 to S106 until the accuracy level of the high-precision data acquisition card for the flow probe is lower than that of the B-code signal source. This determines the accuracy level of the high-precision data acquisition card for the flow probe. For example, the accuracy level of the high-precision data acquisition card for the flow probe can be considered the highest accuracy level among those consistent with the B-code signal source. This completes the high-precision time detection of the flow probe.
[0083] This invention also discloses an electronic device, which includes a memory and a processor; wherein,
[0084] The memory is used to store one or more computer instructions, wherein the one or more computer instructions are executed by the processor to implement any of the above method steps.
[0085] This invention also discloses a computer-readable storage medium storing computer instructions, wherein the computer instructions, when executed by a processor, implement any of the above-described method steps.
[0086] The above description is merely a preferred embodiment of the present invention and an explanation of the technical principles employed. Those skilled in the art should understand that the scope of the invention involved in the embodiments of the present invention is not limited to the technical solutions formed by specific combinations of the above-described technical features, but should also cover other technical solutions formed by arbitrary combinations of the above-described technical features or their equivalent features without departing from the inventive concept. For example, technical solutions formed by substituting the above-described features with (but not limited to) technical features with similar functions disclosed in the embodiments of the present invention.
Claims
1. A high-precision time detection method for flow probes, comprising: Step S101: Connect the B code time synchronization signal of the B code signal source to the B code signal input port of the flow probe high-precision data acquisition card; Connect the 1PPS signal-S output from the B code signal source to channel A of the first frequency counter, and connect the 1PPS signal-D output from the flow probe high-precision data acquisition card to channel B of the first frequency counter. Step S102: Start the first frequency counter to accumulate and count the input pulse signals received by channel A and channel B, and continuously record for a first preset time to obtain a first count record. The first count record includes a count value t1 obtained by accumulating and counting the input pulse signals received by channel A, and a count value t2 obtained by accumulating and counting the input pulse signals received by channel B. Step S103: Calculate the error value and error percentage between the corresponding count values of channel A and channel B based on the first counting record; Step S104: Connect the 1MHz signal-S output from the B code signal source to channel A of the second frequency counter, and connect the 1MHz signal-D output from the flow probe high-precision data acquisition card to channel B of the second frequency counter. Step S105: Start the second frequency counter to accumulate and count the input square wave signals received by channel A and channel B, and continuously record for a second preset time to obtain a second count record. The second count record includes the count value t3 obtained by accumulating and counting the input square wave signals received by channel A, and the count value t4 obtained by accumulating and counting the input square wave signals received by channel B. Step S106: Calculate the error value and error percentage between the corresponding count values of channel A and channel B based on the second counting record; Step S107: Determine the time accuracy and time stability of the flow probe based on the error value dt1, error percentage rt1, error value dt2, and error percentage rt2.
2. The method according to claim 1, wherein in step S103, the error value is calculated according to the following formula: dt1 = t2 - t1; The error percentage is calculated according to the following formula: rt1 = (t2 - t1) / t1.
3. The method according to claim 1 or 2, wherein in step S106, the error value is calculated according to the following formula: dt2 = t4 - t3; The error percentage is calculated according to the following formula: rt2 = (t4 - t3) / t3.
4. In the method according to claim 1 or 2, in step S107, the magnitudes of the error value dt1 and the error percentage rt1 represent the stability of the high-precision time, and the smaller the value, the higher the stability; the magnitudes of the error value dt2 and the error percentage rt2 represent the accuracy of the time, and the smaller the value, the higher the accuracy.
5. The method according to claim 4, wherein step S107 further comprises: The accuracy level of the high-precision data acquisition card for the flow probe is determined based on the error percentage rt2.
6. The method according to claim 5, wherein determining the accuracy level of the high-precision data acquisition card for the flow probe based on the error percentage rt2 includes: Based on the magnitude of the error percentage rt2, determine whether the high-precision data acquisition card for the flow probe is of the same precision level as the B-code signal source. If the precision levels are the same, select a higher-level B-code signal source and repeat steps S101 to S106 until the precision level of the high-precision data acquisition card for the flow probe is lower than that of the B-code signal source. This will determine the precision level of the high-precision data acquisition card for the flow probe.
7. A high-precision time detection system for a flow probe used to implement the high-precision time detection method for a flow probe according to any one of claims 1-6, comprising: The components include a flow probe, a first frequency counter, a second frequency counter, and a B-code signal source, among which: The B-code signal output terminal of the B-code signal source is connected to the B-code signal input terminal of the flow probe; The 1PPS signal output terminal of the B code signal source and the 1PPS signal output terminal of the flow probe are respectively connected to channel A and channel B of the first frequency counter. The 1MHz signal output terminal of the B-code signal source and the 1MHz signal output terminal of the flow probe are respectively connected to channel A and channel B of the second frequency counter.
8. The system according to claim 7, wherein, The first frequency counter accumulates and counts the two input pulse signals respectively. The count value of the first frequency counter is incremented by one for each rising edge detected. The count is continuously recorded for a first preset duration to obtain the number of pulses for each of the two input pulse signals. The second frequency counter accumulates the counts for the two input square wave signals respectively. The count value of the second frequency counter is incremented by one for each rising edge detected, and the count is continuously recorded for a second preset duration to obtain the number of square waves for each of the two input square wave signals.
9. An electronic device comprising a memory and a processor; wherein, The memory is used to store one or more computer instructions, wherein the one or more computer instructions are executed by the processor to implement the steps of the method according to any one of claims 1-6.
10. A computer-readable storage medium having stored thereon computer instructions, wherein, When executed by a processor, the computer instructions implement the steps of the method described in any one of claims 1-6.
Citation Information
Patent Citations
Network performance monitoring and business testing system and testing method
CN109617743A
Method, device and equipment for determining pulse frequency and medium
CN114623939A