Diagnostic method development when weapons characteristics measuring based on spectral analysis for signals phase shift determination Vasyl Lytvyn1,†, Victoria Vysotska1,†, Sergey Tyshko2,∗,†, Oleksandr Lavrut3,†, Tetiana Lavrut3,† and Mariia Nazarkevych1,† 1 Lviv Polytechnic National University, S. Bandera 12, 79013Lviv, Ukraine 2 State Scientific Research Institute of Armament and Military Equipment Testing and Certification, 18001 Cherkasy, Ukraine 3 Hetman Petro Sahaidachnyi National Army Academy, Heroes of Maidan 32, 79026Lviv, Ukraine Abstract The analysis of known methods that have found wide application in measuring technology and are designed to control the technical and operational characteristics associated with the measurement of phase shift during the development, manufacture and operation of weapons and military equipment, has been carried out. Based on this analysis, it was determined that measuring systems designed to determine the phase shift of two harmonic signals have two information transmission channels. A measurement task was set to determine the phase shift of two harmonic signals, using the spectral analysis of the signal obtained by summing the harmonic signals after performing their two-semiperiod transformation. The assumptions necessary list for the analytical ratios synthesis, which establishes the relationship between the phase spectra and the amplitudes (power) of the signal obtained by summing harmonic signals after carrying out their two-semiperiod transformation and phase shift of two harmonic signals, are determined. Analytical relations are proposed that establish the relationship between the above- mentioned characteristics. It is shown that the values of the spectrum of phases and amplitudes calculated using the proposed expressions and ratios for calculating the Fourier series coefficients differ by no more than 0.1%. The application of the proposed approach in the artificial intelligence system to diagnose and determine the state of modern weapons and military equipment will allow us to reduce the requirements for measuring equipment without reducing the accuracy of measurements. Keywords decision-making system, complex system, phase shift, harmonic signal, spectral analysis1 MoDaST-2024: 6th International Workshop on Modern Data Science Technologies, May, 31 - June, 1, 2024, Lviv- Shatsk, Ukraine ∗ Corresponding author. † These authors contributed equally. vasyl.v.lytvyn@lpnu.ua (V. Lytvyn); Victoria.A.Vysotska@lpnu.ua (V. Vysotska); sergeytyshko57@gmail.com (S. Tyshko); alexandrlavrut@gmail.com (O. Lavrut); Lavrut_t_v@i.ua (T. Lavrut); mariia.a.nazarkevych@lpnu.ua (M. Nazarkevych) 0000-0002-9676-0180 (V. Lytvyn); 0000-0001-6417-3689 (V. Vysotska); 0000-0003-3838-2027 (S. Tyshko); 0000-0002-4909-6723 (O. Lavrut); 0000-0002-1552-9930 (T. Lavrut); 0000-0002-6528-9867 (M. Nazarkevych) © 2024 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR ceur-ws.org Workshop ISSN 1613-0073 Proceedings 1. Introduction Today, the Defense Forces of Ukraine actively use various types of weapons and military equipment. An important component of solving the task of effective defence of the state is ensuring the serviceable state of the anti-aircraft defence system, which is part of the units and units of the defence forces, as well as its further modernization, as well as the development of the latest types of weapons. The assessment of the technical condition of the weapons and military equipment (WME) sample, or the decision on admission to the supply of the latest weapons to the Armed Forces, is made based on tests, by analyzing its technical characteristics for compliance with regulatory and technical documents (technical task, technical conditions, operating instructions, repair documentation sets ). Defects are an important technological operation during the capital repair works of WME. As you know, the task of this technological operation is to determine the possibility of using components from the composition of the sample subject to restoration. For mechanical components, various parameters are measured, including geometric dimensions, and physical and chemical properties of the material. To reduce the time and increase the accuracy of measurements, predicting the state of the WME, it is necessary to use intelligent diagnostic systems [16]. The use of artificial intelligence in such systems for processing measurement results, and forecasting (simulation) of the future state of WME is an integral component of the development of this field. The measurement of the above-mentioned characteristics of small arms, artillery and missile-artillery weapons, and wheeled and tracked military equipment is based on non-destructive control methods. Non-destructive control methods include radiography, ultrasonic flaw detection, magnetic resonance research methods, and others. Measuring systems that implement the specified measurement methods widely use phaseometry methods [1, 2, 17, 34]. Also, phase measurement methods are widely used in radar and radio navigation, aviation and space engineering, geodesy, mechanical engineering, communication and many other fields [7, 21-27]. The phase-measuring transformation of various physical processes into the phase shift of harmonic signals ensures high metrological characteristics. Therefore, phaseometry, as a method of transformation and measurement, has long gone beyond the traditional use in radio engineering, navigation and communication and is successfully used in experimental physics, radio physics, experimental medicine, modern fields of science and technology when carrying out precision measurements [3-5, 13-15, 17-20, 37-39]. Based on the above, in [2, 5], a list of WME parameters, which should be converted into a phase shift when measured during sample manufacturing and testing, is defined. These parameters can include: electrical and magnetic conductivity, permeability, geometric dimensions, movement parameters, microdisplacement, capacity, element inductance, liquid level, liquid or gas consumption, angle of rotation, speed of rotation, electric voltage and current, temperature, distance, angle, delay time. Thus, conducting scientific research to find new principles for determining the phase shift, which will allow to reduce the cost of performing work on the control of the characteristics of the latest and modernized weapons at the stages of development and production, is relevant. 2. Related works For harmonic signals, such concepts as phase, initial phase, phase shift and delay time are used in measuring technology. The most complete classification of methods for measuring phase shifts of harmonic signals is given in works [2, 5, 6, 13, 14, 15]. According to the principle of measurement, the methods of phaseometry are divided into compensatory methods and methods of converting the phase shift into other values - voltage, time interval, and geometric parameters of oscillographic images of the investigated signals. These methods differ from each other in technical implementation, complexity and accuracy. Compensation methods [2, 5, 6, 13, 14] are based on the process of balancing (compensating) the phase shift Δφ∈[0,2π) between the measured harmonic signals, that is, reducing the phase shift to zero by adjusting the phase of one of the signals using an adjustable phase shifter (measures of phase shift). This method ensures the achievement of high measurement accuracy, close to the accuracy of a phase shifter. Measurement methods based on the transformation of the phase shift into other signals [2, 5, 6, 13, 14, 15] allow determining the value of the phase shift of the signals after their transformation into other intermediate values that are convenient to use for measurement. These intermediate values include voltage, current, displacement of the electron beam of the oscilloscope, and time values. The disadvantages of known methods include [2, 5, 6]:  a significant impact on the accuracy of measuring the phase shift of the error component, which is caused by the phase asymmetry of the signal transmission channels;  the presence of two channels for conducting analogue-to-digital conversion of input signals, which leads to the need for mutual synchronization of the frequency of clock generators for each of the channels;  significant impact on the accuracy of external and internal noise measurement;  non-linear nature of the scale. At present, one of the effective and widespread ways of reducing the impact of external and internal noise on the quality of solving problems of analyzing and processing signals of various natures is the use of spectral analysis [35-36]. As is known, a periodic signal of any form can be decomposed into harmonic signals whose frequencies are multiples of the frequency (period) of the analyzed signal. A similar research method is called spectral analysis, the mathematical basis of which is the Fourier series [8, 31-33]. Fourier series of arbitrary periodic signals can contain an infinitely large number of terms. One of the advantages of the Fourier transform is that when limiting the Fourier series to any finite number of its terms, it provides the best mean square error of approximation to the original function (for a given number of terms). The convenience of using the frequency representation of signals lies in the fact that harmonic functions are eigenfunctions of operations of transfer, integration, differentiation and other linear operations invariant in coordinates. They pass through linear systems without changing the shape and frequency of the harmonics, only the initial phase and amplitude of the oscillations change. In the general case, when expanding into a Fourier series of a periodic Т Т signal with a period 𝑇𝑠 , it is possible to use an interval [− с⁄2 , с⁄2]. 2𝜋 If we denote the angular frequency through 𝜔с then since 𝜔с = Т , then for the function с Т Т 𝑓(𝑡) on the interval [− с⁄2 , с⁄2], the Fourier series has the form: 𝑛 𝑎0 (1) 𝑓(𝑡) = + ∑{𝑎𝑘 𝑐𝑜𝑠 𝜔𝑐 𝑘𝑡 + 𝑏𝑘 𝑠𝑖𝑛 𝜔𝑐 𝑘𝑡}. 2 𝑘=1 Expressions for determining the coefficients of the Fourier series 𝑎𝑘 , 𝑏𝑘 have the form: 𝑇𝑐 2 (2) 𝑎𝑘 = 𝑇 ∫ 2𝑇𝑐 𝑓(𝑡) 𝑐𝑜𝑠 𝜔𝑐 𝑘𝑡𝑑𝑡 , (𝑘 = 0,1,2, . . . , 𝑛) 𝑐 − 2 𝑇𝑐 2 (3) 𝑏𝑘 = 𝑇 ∫ 2𝑇𝑐 𝑓(𝑡) 𝑠𝑖𝑛 𝜔𝑐 𝑘𝑡𝑑𝑡 (𝑘 = 1,2, . . . , 𝑛). 𝑐 − 2 Series (1) represents the decomposition of the periodic signal f(t) into the sum of real elementary harmonic functions (cosine and sine) with weighting coefficients, the geometric sum of whose values (that is, the values of 𝑎𝑘 and 𝑏𝑘 are nothing but the real amplitudes of the corresponding harmonic oscillations with frequencies 𝑘𝜔с. The set of amplitude values of these harmonics forms a one-sided physically real (only for positive frequencies 𝑘𝜔с signal spectrum. The goal of the work. To propose the scientific and technical basis of the alternative principle of determining the phase shift, the mathematical basis of which is possible to consider analytical relations that establish the relationship between the phase shift of two harmonic signals and the spectral characteristics of the signal obtained when summing them after carrying out a two-semiperiod transformation, which will allow to significantly reduce the error component, due to the phase asymmetry of the signal transmission channels and the influence of external and internal noise during monitoring of the characteristics (parameters) of the WME. 3. Methods and models At present, the determination of the phase shift is of greatest interest to phasometry. Phase shift [9, 28-30] refers to the modulus of the difference between the initial phases of two harmonic signals of the same frequency. As a rule, the measurement of phase shifts of signals is based on the model of a harmonic signal, which is specified without changes in its parameters over an infinite time interval. This model is ideal, and in practice, a model with a finite time window is used, that is, measurements are carried out on a finite time interval. Then, based on the above, we will formulate the problem of determining the phase shift using the spectral analysis of the signal obtained as a result of the summation of two harmonic signals, after carrying out their two-half-period transformation. Let's consider the initial data necessary to solve the specified measurement problem. Let there be two harmonic signals 𝑢1 (𝑡) and 𝑢2 (𝑡), which have a phase shift relative to the other equal to 𝛥𝜑, in the interval from 0 to 2𝜋. Based on the fact that phase shift measurements refer to relative measurements, the mathematical record of changes in signals 𝑢1 (𝑡)and 𝑢2 (𝑡)can be represented in the form: 𝑢1 (𝑡) = 𝑈𝑚1 𝑐𝑜𝑠( 2𝜋𝑓𝑡), (4) 𝑢2 (𝑡) = 𝑈𝑚2 𝑐𝑜𝑠( 2𝜋𝑓𝑡 + 𝛥𝜑), where 𝑈𝑚1 is signal amplitude 𝑢1 (𝑡); 𝑈𝑚2 is signal amplitude 𝑢2 (𝑡); 𝑓 = 1⁄Т is signal frequency; Т is signal period. The module is extracted from signals 𝑢1 (𝑡) and 𝑢2 (𝑡), as a result of which we get: / (5) 𝑢1 (𝑡) = |𝑢1 (𝑡)| = |𝑈𝑚1 𝑐𝑜𝑠( 2𝜋𝑓𝑡)|, / 𝑢2 (𝑡) = |𝑢2 (𝑡)| = |𝑈𝑚2 𝑐𝑜𝑠( 2𝜋𝑓𝑡 + 𝛥𝜑)|, / / After summing the signals 𝑢1 (𝑡) and 𝑢2 (𝑡), we get: / 𝑢Σ (6) 𝑈1𝑚𝑖𝑛 − 𝑈2𝑚𝑖𝑛 / 𝑈1𝑚𝑖𝑛 + 𝑈2𝑚𝑖𝑛 2𝜋𝑓 𝑈1𝑚𝑖𝑛 + 𝑡 + (𝑈1𝑚𝑎𝑥 − ) sin ((2𝑓 − ) 𝑡) , 𝑡1 ≤ 𝑡 < 𝑡2 𝑡1.2 2 𝛥𝜑 =⟨ 𝑈2𝑚𝑖𝑛 − 𝑈1𝑚𝑖𝑛 / 𝑈2𝑚𝑖𝑛 + 𝑈1𝑚𝑖𝑛 2𝜋𝑓 𝑈2𝑚𝑖𝑛 + 𝑡 + (𝑈2𝑚𝑎𝑥 − ) sin (( ) 𝑡) , 𝑡2 ≤ 𝑡 < 𝑡1 𝑡2.3 2 𝛥𝜑 Т where: 𝑈1𝑚𝑖𝑛 and 𝑈2𝑚𝑖𝑛 are breakpoints of the function on the interval from 0 tо ; 2 𝛥𝜑 1 𝛥𝜑 𝑈1𝑚𝑎𝑥 = (𝑈𝑚1 + 𝑈𝑚2 ) 𝑐𝑜𝑠 2 is local maximum on the time interval 𝑡1.2 = 2𝑓 − 2𝜋𝑓 ; ′ 𝛥𝜑 𝛥𝜑 ′ 𝑈2𝑚𝑎𝑥 = (𝑈𝑚1 + 𝑈𝑚2 ) 𝑠𝑖𝑛 2 is local maximum on the time interval 𝑡2.1 = 2𝜋𝑓 . Figure 1: Structural diagram of the transformation of signals into a signal Synthesis of time diagrams that describe the order of formation of the signal / 𝑢∑ (𝑡) using relations (4) and (5) will be carried out in the Matlab environment using the Simulink tool application. To do this, we synthesize the scheme shown in Fig. 1. In this scheme, the synthesis of signal u1 (t) is carried out by the harmonic signal generator "Sine Wave", and the signal u2 (t) by the generator "Sine Wave 1". From the "Sine Wave" output, the signal enters the "Abs" block and the oscilloscope marked "Scope", which is designed to form the time diagram of the 𝑢1 (𝑡) signal. According to the "Sine Wave 1" output, the signal enters the "Abs 1" block and the oscilloscope marked "Scope 2", which is designed to form / the timing diagram of the u2 (t) signal. At the output of the "Abs" block, a signal 𝑢1 (𝑡) is formed, which is fed to the first input of the summation device "Add" and "Scope 1". "Scope / 1" forms a time diagram of the signal 𝑢1 (𝑡). Similarly, at the output of the "Abs 1" block, a / signal 𝑢2 (𝑡) is formed, which is fed to the second input of the summation device "Add" and / "Scope 3", which forms a time diagram of the signal 𝑢2 (𝑡). At the output of the "Add" block, / we will receive the signal 𝑢∑ (𝑡). The signal 𝑢∑ (𝑡) is fed to "Scope 4", which will form its signal time diagram 4. Experiments, results and discussion / / / Fig. 2 shows the time diagrams of the signals 𝑢1 (𝑡)/𝑢2 (𝑡), 𝑢1 (𝑡)/𝑢2 (𝑡) and 𝑢𝛴 (𝑡) at the value 𝛥𝜑 = 0, as well as the "Sine Wave" generator settings tab ” and “Sine Wave1”. Figure 2: Time diagrams of signals 𝑢1 (𝑡) and 𝑢2 (𝑡),, 𝑢1/ (𝑡) and 𝑢2/ (𝑡), and 𝑢𝛴/ (𝑡) at the value of 𝛥𝜑 = 0 According to the settings tabs, generators "Sine Wave" generates a sinusoidal signal with an amplitude equal to 1, and "Sine Wave1" with an amplitude equal to 0.8, other settings of the tabs are as follows: the constant component of the signals is 0, the frequency of the signal is 1 rad/s, the initial phase of the signals is 0 rads, the quantization time is minimal for this device, i.e. 0 on the tab. The timing diagram shown on "Scope" and "Scope 2" fully confirms / the settings on the tab. The time diagrams for the 𝑢∑ (𝑡) signal are shown on the "Scope 4" / tab. The analysis of this time diagram shows that the signal 𝑢∑ (𝑡) has a period equal to π s, the maximum value is equal to 1.8, and the minimum value is equal to 0. The obtained results correspond to the summation of two harmonic signals, after carrying out their two- and-a-half period transformation. / / / Fig. 3 shows the time diagrams of signals 𝑢1 (𝑡) and 𝑢2 (𝑡), 𝑢1 (𝑡) and 𝑢2 (𝑡) and 𝑢𝛴 (𝑡), at the value of 𝛥𝜑 = 0,9 rad, as well as the generator settings tab "Sine Wave" and "Sine Wave1". According to the configuration tabs, the "Sine Wave" generator generates a sinusoidal signal with an amplitude equal to 1, the constant component of the signal is 0, the frequency of the signal is 1 rad/s, the initial phase of the signal is 0 rad, the quantization time is minimal for this device, i.e. 0 on the tab. Figure 3: Time diagrams of signals 𝑢1 (𝑡) and 𝑢2 (𝑡),, 𝑢1/ (𝑡) and 𝑢2/ (𝑡), and 𝑢𝛴/ (𝑡) at the value of 𝛥𝜑 = 0.9 rad According to the settings tabs, the "Sine Wave1" generator generates a sinusoidal signal with an amplitude equal to 0.8, the constant component of the signal is 0, the frequency of the signal is 1 rad/s, the initial phase of the signal is 0.9 rad, the quantization time is minimal / / for this device, i.e. 0 at tab. The time diagrams for signals 𝑢1 (𝑡) and 𝑢2 (𝑡) are shown in the "Scope 1" and "Scope 3" tabs. The analysis of the data of the time diagrams shows that they / / have a period equal to π s, the maximum values of 𝑢1 (𝑡) are equal to 1, and 𝑢2 (𝑡) is 0.8, also / / / the 𝑢2 (𝑡) signal precedes the 𝑢1 (𝑡) about 0.9. The time diagrams for the 𝑢∑ (𝑡) signal is / shown on the "Scope 4" tab. Analysis of this time diagram shows that the signal 𝑢∑ (𝑡) has a period equal to π s. In the time interval from 0 to (π-0.9) s, the maximum value is approximately 1.65, and the minimum value is approximately 0.78. In the time interval from (𝜋 − 0,9) to π, the maximum value is approximately 0.8, and the minimum value is 0.6. / / / Fig. 4 shows the time diagrams of the signals 𝑢1 (𝑡)/𝑢2 (𝑡), 𝑢1 (𝑡)/𝑢2 (𝑡), and 𝑢𝛴 (𝑡) at the value of 𝛥𝜑 = 1.2 rad, as well as the tab for setting the generator "Sine Wave”/“Sine Wave1”. Figure 4: Time diagrams of signals 𝑢1 (𝑡) and 𝑢2 (𝑡), 𝑢1/ (𝑡) and 𝑢2/ (𝑡), and 𝑢𝛴/ (𝑡) at Δφ=1.2 rad According to the configuration tabs, the "Sine Wave" generator generates a sinusoidal signal with an amplitude equal to 1, the constant component of the signal is 0, the frequency of the signal is 1 rad/s, the initial phase of the signal is 0 rad, the quantization time is minimal for this device, i.e. 0 on the tab. According to the settings tabs, the "Sine Wave1" generator generates a sinusoidal signal with an amplitude equal to 0.8, the constant component of the signal is 0, the frequency of the signal is 1 rad/s, the initial phase of the signal is 1.2 rad, the quantization time is minimal for this device, i.e. 0 on the tab. The time / / diagrams for signals 𝑢1 (𝑡) and 𝑢2 (𝑡) are shown in the "Scope 1" and "Scope 3" tabs. The analysis of the data of the time diagrams shows that they have a period equal to 𝜋 s, the / / / maximum values of 𝑢1 (𝑡) are equal to 1, and 𝑢2 (𝑡) is 0.8, also the 𝑢2 (𝑡)signal precedes the / 𝑢1 (𝑡) about 1.2 rad. The time diagrams for the 𝑢∑ (𝑡) signal are shown on the "Scope 4" / tab. Analysis of this time diagram shows that the signal 𝑢∑ (𝑡) has a period equal to π s. In the time interval from 0 to (𝜋 − 1.2) s, the maximum value is approximately 1.5, and the minimum value is 0.95. In the time interval from (𝜋 − 1.2) to π, the maximum value is approximately 1.05, and the minimum value is 0.75. / From Fig. 4, it can be seen, and it is also shown in [10, 11] that 𝑢∑ (𝑡) depending on the value of Δφ of the input signals 𝑢1 (𝑡) and 𝑢2 (𝑡): 1 Т  is a periodic signal with a period of T   == ; 2𝑓 2  has two local maxima 𝑈1𝑚𝑎𝑥 and 𝑈2𝑚𝑎𝑥 , which correspond to the moments of 𝑡1𝑚𝑎𝑥 , 𝑡2𝑚𝑎𝑥 and two breaks 𝑈1𝑚𝑖𝑛 and 𝑈2𝑚𝑖𝑛 , which correspond to the moments of time 𝑡1 , 𝑡2 , respectively;  a change in the value of the phase shift angle 𝛥𝜑 leads to a change in the values of 𝑈1𝑚𝑎𝑥 , 𝑈2𝑚𝑎𝑥 , 𝑈1𝑚𝑖𝑛 end 𝑈2𝑚𝑖𝑛 and time parameters 𝑡1𝑚𝑎𝑥 , 𝑡2𝑚𝑎𝑥 end 𝑡2 . / A change in the values of the above characteristics of the signal 𝑢∑ (𝑡) leads to a change in its shape, which in turn will lead to a change in the values of the coefficients 𝑎𝑘 and 𝑏𝑘 of the Fourier series depending on the change in the magnitude of the phase shift 𝛥𝜑. When conducting spectral analysis, such concepts as amplitude spectrum, power spectrum and phase spectrum are used. When conducting spectral analysis, such concepts as amplitude spectrum, power spectrum and phase spectrum are used. According to [12], the spectrum of amplitudes is understood as the set of absolute values of the coefficients 𝐶𝑘 (𝑘 = 0,1,2, . . . , 𝑛), which are determined from the known values of 𝑎𝑘 and 𝑏𝑘 by the ratio: (7) √𝑎𝑘2 + 𝑏𝑘2 |𝐶𝑘 | = . 2 The value |𝐶𝑘 | indicates the value of the amplitude of the k-th harmonic signal when expanded into a Fourier series. The power spectrum is understood as a set of arguments of the values |𝐶𝑘 |2. The spectrum of phases means the set of arguments ∠𝐶𝑘 (𝑘 = 0,1,2, . . . , 𝑛), which are determined by the known values of 𝑎𝑘 and 𝑏𝑘 by the ratio: 𝑏𝑘 (8) ∠𝐶𝑘 = 𝑎𝑟𝑐𝑡𝑔 . 𝑎𝑘 The value of ∠𝐶𝑘 indicates the value of the initial phase of the k-th harmonic signal when expanded into a Fourier series. Then, based on the above, we will formulate the task of determining the phase shift using Fourier series decomposition, as a synthesis of analytical relations that describe the relationship between the change in the spectrum of amplitudes / (power) and the spectrum of the phases of the signal 𝑢∑ (𝑡) depending on the change in 𝛥𝜑. Synthesis of relations that determine the relationship between the value of 𝛥𝜑 and the characteristics of the spectrum of amplitudes (power) and the spectrum of phases of the / signal 𝑢∑ (𝑡) using the following assumptions. As can be seen from (4), the signal / / / 𝑢∑ (𝑡) is formed by summing two signals 𝑢1 (𝑡) and 𝑢2 (𝑡). These signals are formed by signals 𝑢1 (𝑡) and 𝑢2 (𝑡), while the amplitude of 𝑈𝑚1 differs from the amplitude of 𝑈𝑚2 by v times, and have a phase shift 𝛥𝜑. Then it can be stated that the signal 𝑢1 (𝑡) precedes the signal 𝑢2 (𝑡)by a certain time interval τ, i.e. 𝑢2 (𝑡) = 𝑣𝑢1 (𝑡 − 𝜏). According to the known values of 𝛥𝜑 and T, the time interval τ is determined by the ratio: 𝛥𝜑 (9) 𝜏= Т. 2𝜋 As shown in [12], the shift of the signal in the time domain by some interval τ leads to a change in the phase spectrum, but the spectrum of the signal amplitudes remains unchanged. Taking into account expressions (7) and (8) provided that 𝑈𝑚1.𝑘 and 𝑈𝑚2.𝑘 are / / the amplitude values of the 𝑘 -th harmonics of signals 𝑢1 (𝑡) and 𝑢2 (𝑡), respectively, and 𝜑1.к / / and 𝜑2.к values of the initial phases of the 𝑘 -th harmonics of the signals 𝑢1 (𝑡) and 𝑢2 (𝑡), / respectively, the Fourier series for the signal 𝑢∑ (𝑡) will have appearance: / u (t )  𝑈𝑚1.0 + ∑𝑛 𝑘=1[𝑈𝑚1.𝑘 𝑠𝑖𝑛( 2𝜔𝑘𝑡 + 𝜑1.𝑘 )] + 𝑈𝑚2.0 + (10) 𝑛 + ∑[𝑈𝑚2.𝑘 𝑠𝑖𝑛( 2𝜔𝑘𝑡 + 𝜑2.𝑘 )] = 𝑘=1 = (𝑈𝑚1.0 + 𝑈𝑚2.0 ) + [𝑈𝑚1.1 𝑠𝑖𝑛( 2𝜔𝑡 + 𝜑1.1 ) + 𝑈𝑚2.1 𝑠𝑖𝑛( 2𝜔𝑡 + 𝜑2.1 )] + +[𝑈𝑚1.2 𝑠𝑖𝑛( 4𝜔𝑡 + 𝜑1.2 ) + 𝑈𝑚2.2 𝑠𝑖𝑛( 4𝜔𝑡 + 𝜑2.2 )] + ⋯ + +[𝑈𝑚1.𝑛 𝑠𝑖𝑛( 2𝑛𝜔𝑡 + 𝜑1.𝑛 ) + 𝑈𝑚2.𝑛 𝑠𝑖𝑛( 2𝑛𝜔𝑡 + 𝜑2.𝑛 )] / To synthesize analytical relations for determining the values of the amplitude 𝑈𝑚 ∑ 1 and / / the initial phase 𝜑𝑚 ∑ 1 of the first harmonic of the signal 𝑢∑ (𝑡), we will construct the vector / diagram of the first harmonics of the signals 𝑢1 (𝑡) and 𝑢2 (𝑡)in the presence of a phase shift equal to 𝛥𝜑 of the signals 𝑢1 (𝑡) and 𝑢2 (𝑡), provided that the initial phase of the signal 𝑢1 (𝑡) is zero, taking into account the following remarks and assumptions:  it is known that the frequency of the first harmonic of the signal is equal to the frequency of the signal; / /  it is shown in [10, 11] and it can be seen in Figs. 2...4 that 𝑢1 (𝑡), 𝑢2 (𝑡) and / 𝑢∑ (𝑡)have a period that is two times shorter than 𝑢1 (𝑡) and 𝑢2 (𝑡); / /  it can be asserted that 𝑢2 (𝑡) = 𝑣𝑢1 (𝑡 − 𝜏) and 𝑢2 (𝑡) = 𝑣𝑢1 (𝑡 − 𝜏), that is, the time / / interval of the shift τ between the signals 𝑢1 (𝑡) and 𝑢2 (𝑡), and the shift time interval τ between signals 𝑢1 (𝑡) and 𝑢2 (𝑡) is the same;  based on relation (9) and the above remarks, it can be seen that the value of the / / phase shift of the first harmonics of signals 𝑢1 (𝑡) and 𝑢2 (𝑡)will be equal to 2𝛥𝜑, if there is a phase shift equal to Δφ of signals 𝑢1 (𝑡) and 𝑢2 (𝑡). The above vector diagram is shown in Fig. 5. From Fig. 5, it can be seen that for the calculation of the values of the characteristics of the vector of the first harmonic of the signal / / / 𝑢∑ (𝑡) according to the vectors of the first harmonics of the signals 𝑢1 (𝑡) and 𝑢2 (𝑡), the parallelogram method was used. Figure 5: Vector diagram for determining the value of the initial phase φm/ 1 of the first harmonic of the signal u / (t ) / Thus, the value of the amplitude 𝑈𝑚 ∑ 1 of the first harmonic of the signal / 𝑢∑ (𝑡) according to the known values of the amplitudes of the first harmonics 𝑈𝑚1.1 and / / 𝑈𝑚2.1 = 𝑈 з 𝑚2.1 of the signals 𝑢1 (𝑡) and 𝑢2 (𝑡) in the presence of a phase shift equal to 𝛥𝜑 of signals 𝑢1 (𝑡) and 𝑢2 (𝑡) is determined by the ratio: / (11) 𝑈𝑚 ∑ 1 = √𝑈 2 𝑚1.1 + 𝑈 2 𝑚2.1 + 2𝑈𝑚1.1 𝑈𝑚2.1 𝑐𝑜𝑠 2 𝛥𝜑 / Then the value of the amplitude 𝑈𝑚 ∑ 𝑘 of the k-th harmonic of the signal / 𝑢∑ (𝑡) according to the known values of the amplitudes 𝑈𝑚1.𝑘 and 𝑈𝑚2.𝑘 of the k-th / / harmonics of the signals 𝑢1 (𝑡) and 𝑢2 (𝑡) at a phase shift value equal to 𝛥𝜑 of signals 𝑢1 (𝑡) and 𝑢2 (𝑡) is calculated using the ratio: / (12) 𝑈𝑚 ∑ 𝑘 = √𝑈 2 𝑚1.𝑘 + 𝑈 2 𝑚2.𝑘 + 2𝑈𝑚1.𝑘 𝑈𝑚2.𝑘 𝑐𝑜𝑠 2 𝑘𝛥𝜑. (𝑘 = 0,1,2, . . . , 𝑛) The synthesis of analytical ratios for calculating the values of 𝑈𝑚1.𝑘 will be carried out in / the following order. Let's determine the spectrum of amplitudes for the signal 𝑢1 (𝑡) under the condition 𝑈𝑚1 = 1. The relations for calculating the coefficients of the Fourier series 𝑎1.𝑘 / и 𝑏1.𝑘 of the signal 𝑢1 (𝑡) are as follows: 𝑇′ (13) 2 𝑎1.𝑘 = 𝑇 ′ ∫ 2𝑇′|𝑐𝑜𝑠( 𝜔𝑡)| 𝑐𝑜𝑠( 2𝜔𝑘𝑡)𝑑𝑡. (𝑘 = 0,1,2, . . . , 𝑛) − 2 𝑇′ (14) 2 2 𝑏1.𝑘 = ∫ ′ |𝑐𝑜𝑠( 𝜔𝑡)| 𝑠𝑖𝑛( 2𝜔𝑘𝑡)𝑑𝑡, (𝑘 = 1,2, . . . , 𝑛) 𝑇 ′ −𝑇 2 Based on the known values of 𝑎1.𝑘 and 𝑏1.𝑘 , the value of 𝑈𝑚1(1).𝑘 is calculated by the ratio: 2 2 (15) √𝑎1.𝑘 + 𝑏1.𝑘 𝑈𝑚1(1).𝑘 = . 2 The values of 𝑈𝑚1(1).𝑘 were calculated using the universal mathematical package MathCAD, the results are presented in Fig. 6. Using the property of the linearity of the / Fourier series, the relationship for determining the value of the amplitude 𝑈𝑚 ∑ 𝑘 k-th / harmonic of the signal 𝑢∑ (𝑡) based on the known values of the amplitudes 𝑈𝑚1 and 𝑈𝑚2 at the value of the phase shift of the initial signals 𝑢1 (𝑡) and 𝑢2 (𝑡), which is equal to 𝛥𝜑, is calculated using the relationship: / (16) 𝑈𝑚 ∑ 𝑘 = 𝑈𝑚1(1).𝑘 √𝑈 2 𝑚1. + 𝑈 2 𝑚2. + 2𝑈𝑚1. 𝑈𝑚2. 𝑐𝑜𝑠 2 𝑘𝛥𝜑, (𝑘 = 0,1,2, . . . , 𝑛) / To synthesize analytical relations to determine the values of the initial phase 𝜑𝑚 ∑ 1 of the / first harmonic of the signal 𝑢∑ (𝑡), we will use the theorem of sines for the triangle formed / by the sides 𝑈𝑚1.1 , 𝑈 р 𝑚2.1 and 𝑈𝑚 ∑ 1 on the vector diagram of the first harmonics of the / / signals 𝑢1 (𝑡) and 𝑢2 (𝑡) in the presence of a phase shift equal to 𝛥𝜑 of the signals 𝑢1 (𝑡) and 𝑢2 (𝑡), we write 𝑈𝑚1.1 𝑈 𝑝 𝑚1.1 = 𝑠𝑖𝑛(𝛼) 𝑠𝑖𝑛 (𝜑/ ) 𝑚∑1 / Solving this equation relative to 𝜑𝑚 ∑ 1 based on the fact that 𝑈𝑚1.1 = 𝑈𝑚1(1).1 𝑈𝑚1 , 𝑝 𝑈𝑚2.1 = 𝑈𝑚2.1 = 𝑈𝑚1(1).1 𝑈𝑚2 , / 𝛼 = 2𝛥𝜑 − 𝜑𝑚 ∑ 1 . / We obtain that the value of the initial phase 𝜑𝑚 ∑ 1 of the first harmonic of the signal / 𝑢∑ (𝑡)depending on the value of Δφ is determined by the following relation, provided that the initial phase of the signal u_1 (t) is zero / 𝑈𝑚1 (17) 𝜑𝑚 ∑ 1 = 𝑎𝑟𝑐𝑐𝑡𝑔 [𝑐𝑜𝑠 𝑒 𝑐(2𝛥𝜑) ⋅ ( + 𝑐𝑜𝑠( 2𝛥𝜑))] 𝑈𝑚2 𝑈𝑚1(1).𝑘 = Figure 6: The value of 𝑈𝑚1(1).𝑘 for 𝑈𝑚1 = 1В / Taking into account (9), for the kth harmonic, the initial phase 𝜑𝑚 ∑ 𝑘 on the interval [− 𝜋⁄2 , 𝜋⁄2] of the signal 𝑢∑ (𝑡)will be equal to: / / 𝑈 (18) 𝜑𝑚 ∑ 𝑘 = 𝑎𝑟𝑐𝑐𝑡𝑔 [𝑐𝑜𝑠 𝑒 𝑐(2𝑘𝛥𝜑) ⋅ (𝑈𝑚1 + 𝑐𝑜𝑠( 2𝑘𝛥𝜑))], (𝑘 = 0,1,2, . . . , 𝑛) 𝑚2 Let us consider the possibility of using the obtained relations for the synthesis of the / spectrum of amplitudes and phases of the signal 𝑢∑ (𝑡). To do this, let's compare the / spectrum of amplitudes and phases of the signal 𝑢∑ (𝑡) obtained using relations (17) and (18) and the following expressions: 𝑇′ (19) 2 𝑎𝑘 = 𝑇 ′ ∫ 2𝑇′[|𝑈𝑚1 𝑐𝑜𝑠( 𝜔𝑡)| + |𝑈𝑚2 𝑐𝑜𝑠( 𝜔𝑡 + 𝛥𝜑)|] 𝑐𝑜𝑠 2 𝜔𝑘𝑡𝑑𝑡, (𝑘 = 0,1,2, . . . , 𝑛) − 2 𝑇′ (20) 2 2 𝑏𝑘 = ∫ ′ [|𝑈𝑚1 𝑐𝑜𝑠( 𝜔𝑡)| + |𝑈𝑚2 𝑐𝑜𝑠( 𝜔𝑡 + 𝛥𝜑)|] 𝑠𝑖𝑛 2 𝜔𝑘𝑡𝑑𝑡, (𝑘 = 1,2, . . . , 𝑛) 𝑇 ′ −𝑇 2 2 2 (21) √𝑎1.𝑘 + 𝑏1.𝑘 𝑏𝑘 𝑈𝑚.𝑘 = , ∠𝐶𝑘 = 𝑎𝑟𝑐𝑡𝑔 . 2 𝑎𝑘 As an indicator of approximation of the spectrum of amplitudes, we use the difference / between the value of 𝑼𝒎 ∑ 𝒌 obtained using relation (16) and the value of 𝑼𝒎.𝒌 calculated using relations (19-21). The calculations were performed in the MathCAD mathematical package, and the comparison results for some values of 𝛥𝜑 , 𝑈𝑚1 and 𝑈𝑚2 are presented in Fig. 7a. These results show that the relative deviation of the amplitude spectrum obtained using (17) differs from the values of the amplitude spectrum obtained using (19 ... 21) less than 0.1%. As an indicator of approximation of the phase spectrum, we use the difference / between the value of 𝜑𝑚 ∑ 𝑘 obtained using relation (18) and the value of ∠𝑪𝒌in the interval [-π⁄2,π⁄2] calculated using relation (19, 20, 21). The results of the comparison for some values of Δφ are presented in Fig. 7b. These results show that the relative deviation of the phase spectrum obtained using (19) differs from the values of the phase spectrum obtained using (19, 20, 21) by less than 0.1%. 5. Conclusions An analysis of the known principles of measuring phase shifts was carried out. The analysis showed that a significant contribution to the final measurement error of phase shifts is made by the component caused by the phase asymmetry of signal transmission channels and the influence of external and internal noise. As an alternative approach to determining the phase shift, it is proposed to use a signal obtained as a result of summing harmonic signals after carrying out their two-semiperiod transformation followed by its spectral analysis. Analytical relations are proposed that establish the relationship between the phase shift and the characteristics of the spectrum of amplitudes and phases of the considered signal. The adequacy of the proposed analytical ratios was checked. As a result of the verification, it was established that the relative discrepancy between the characteristics of the spectrum obtained using the proposed ratios and using the ratios for calculating the coefficients of the Fourier series does not exceed 0.1%. a b Figure 7: The difference between the value of a) the quantities U m/  k / U m.k and b) φк / Ck These analytical relations can be considered as a mathematical basis for the synthesis of methods for determining the phase shift of two harmonic signals using the spectral analysis of the signal obtained as a result of their addition after carrying out their two-half-cycle transformation. Further research should be directed to the development of intelligent diagnostic systems, which will be built based on algorithms of artificial intelligence, which, in turn, will reduce the time and increase the accuracy of measurements and forecasting the state of modern WME. References [1] O.M. Kryukov, R.S. Melnikov, V.A. Muzychuk, The method of diagnosing the technical condition of bores and ammunition based on the identification of the characteristics of the ballistic elements of the shot, Collection of scientific works of the National Academy 2(32) (2018) 5-11. doi:10.33405/2409-7470/2018/2/32/155166. [2] Y. V. Kuts, L. M. Shcherbak, Statistical phasometry. Ternopil State Technical University, 2009. [3] O.O. Lavrut, T.V. Lavrut, K.O. Klymovych, Yu.M. Zdorenko, The latest technologies and means of communication in the Armed Forces of Ukraine: the path of transformation and development prospects, Science and technology of the Air Force of the Armed Forces of Ukraine 1(34) (2019) 91-101. doi:10.30748/nitps.2019.34.13. [4] O.O. Lavrut, K.O. Klymovych, M.L. Tarasyuk, O.L. Antonyuk, Status and prospects of the use of modern technologies and means of radio communication in the Armed Forces of Ukraine Weapon, Systems and military equipment 1(49) (2017) 42-49. [5] L Bohdal, L. Kukiełka, S. Legutko, R. Patyk, A. M. Radchenko, Modeling and Experimental Research of Shear-Slitting of AA6111-T4 Aluminum Alloy Sheet, Materials 13(14) (2020) 3175. doi:10.3390/ma13143175. [6] F. Bonavolontà, M. D’Apuzzo, A. Liccardo, G. Mieleb, Harmonic and interharmonic measurements through a compressed sampling approach. Measurement 77 (2016) 1- 15. doi:10.1016/j.measurement. 2015.08.022. [7] M.I. Skolnik, Radar Handbook. Third Edition. The McGraw-Hill Companies, 2008. [8] G. Sun, L. Wu, Z. Kuang, Z. Ma, J. Liu, Practical tracking control of linear motor via fractional-order sliding mode, Automatica 94 (2018) 221-235. doi:10.1016/j.automatica.2018.02.011. [9] Y. Wang, C. Wang, Y. Tao, Fast Frequency Acquisition and Phase Locking of Nonplanar Ring Oscillators, Appl. Sci. 7 (2017) 10-32. doi:10.3390/app7101032. [10] S. Kihong, On the Selection of Sensor Locations for the Fictitious FRF based Fault Detection Method, International Journal of Emerging Trends in Engineering Research 7(7) (2019) 569-575. doi:10.30534/ijeter/2019/277112019. [11] V.V. Pabyrivskyi, N.V. Pabyrivska, P.Y. Pukach, The study of mathematical models of the linear theory of elasticity by presenting the fundamental solution in harmonic potentials, Mathematical Modeling and Computing 7(2) (2020) 259–268. [12] Yu. Kuts, A. Protasov, Y. Lycenco, O. Dugin, O. Bliznuk, V. Uchanin, Using Multidifferential Transducer for Pulsed Eddy Current Object Inspection, in Proceedings of IEEE First Ukraine Conference on Electrical and computer engineering (Ukrcon), May 29 –June 2, 2017. Kyiv, Ukraine, pp.826-829. doi:1109/UKRCON.2017.8100361 [13] S. A. Tyshko, V. G. Smolyar, O. E. Zabula, Analysis of the possibility of using two-half- period conversion for measuring the phase shift of harmonic signals with equal amplitude, Collection of scientific works of the Kharkiv University of the Air Force 2(30) (2013) 42-44. [14] S. V. Gubin, S. O. Tyshko, O. E. Zabula, Yu. M. Chernychenko, Oscillographic method of phase shift measurement based on two-semi-periodic transformation, Radioelectronic and computer systems 4(2019) (2019) 47 – 54. doi:10.32620/reks.2019.4.05. [15] S. Kihong, On the Selection of Sensor Locations for the Fictitious FRF based Fault Detection Method, International Journal of Emerging Trends in Engineering Research 7(7) (2019) 569-575. doi:10.30534/ijeter/2019/277112019. [16] O. Klymovych, V. Hrabchak, O. Lavrut, T. Lavrut, V. Lytvyn, V. Vysotska, The Diagnostics Methods for Modern Communication Tools in the Armed Forces of Ukraine Based on Neural Network Approach, CEUR Workshop Proceedings 2631 (2020) 198-208. [17] S. Tyshko, O. Lavrut, V. Vysotska, O. Markiv, O. Zabula, Y. Chernichenko, T. Lavrut, Compensatory Method for Measuring Phase Shift Using Signals Bisemiperiodic Conversion in Diagnostic Intelligence Systems, CEUR Workshop Proceedings 3312 (2022) 144-154. [18] M. Nazarkevych, B. Yavourivskiy, I. Klyuynyk, Editing raster images and digital rating with software, in Proceedings of IEEE the Experience of Designing and Application of CAD Systems in Microelectronics, 2015, February, pp. 439-441. [19] M. Nazarkevych, Y. Kynash, R. Oliarnyk, I. Klyujnyk, H. Nazarkevych, Application perfected wave tracing algorithm, in Proceedings of IEEE First Ukraine Conference on Electrical and Computer Engineering (UKRCON), 2017, May, pp. 1011-1014. [20] M. Medykovskyy, P. Lipinski, O. Troyan, M. Nazarkevych, Methods of protection document formed from latent element located by fractals, in Proceedings of Xth International Scientific and Technical Conference on Computer Sciences and Information Technologies (CSIT), 2015, September, pp. 70-72. [21] V. Motyka, M. Nasalska, Y. Stepaniak, V. Vysotska, M. Bublyk: Radar Target Recognition Based on Machine Learning, CEUR Workshop Proceedings 3373 (2023) 117-128. [22] L. Mochurad, R. Bliakhar, N. Reverenda, Identification and Tracking of Unmanned Aerial Vehicles Based on Radar Data, CEUR Workshop Proceedings 3426 (2023) 171-181. [23] A. Vasyliuk, T. Basyuk, V. Lytvyn. Specialized interactive methods for using data on radar application models, CEUR Workshop Proceedings 2631 (2020) 1–11. [24] I.N. Garkusha, V.V. Hnatushenko, V. V. Vasyliev, Research of influence of atmosphere and humidity on the data of radar imaging by Sentinel-1, in Proceedings of IEEE 37th International Conference on Electronics and Nanotechnology (ELNANO), Kiev, 2017, pp. 405-408. doi: 10.1109/ELNANO.2017.7939787. [25] V. Hnatushenko, I. Garkusha, V. Vasyliev, Creating soil moisture maps based on radar satellite imagery, Proceedings of SPIE 10426, Active and Passive Microwave Remote Sensing for Environmental Monitoring, 104260J, 3 October 2017. doi:10.1117/12.2278040. [26] O. Kavats, V. Hnatushenko, Y. Kibukevych, Y. Kavats, Flood Monitoring Using Multi- temporal Synthetic Aperture Radar Images, Advances in Intelligent Systems and Computing 1080 (2020) Springer, Cham. doi:10.1007%2F978-3-030-33695-0_5 [27] M. Makaruk, A. Nazarov, I. Shubin, N. Shanidze, Knowledge Representation Method for Object Recognition in Nonlinear Radar Systems, CEUR Workshop Proceedings 2870 (2021) 948-958. [28] V. Lytvyn, V. Vysotska, I. Peleshchak, I. Rishnyak, R. Peleshchak, Time Dependence of the Output Signal Morphology for Nonlinear Oscillator Neuron Based on Van der Pol Model, International Journal of Intelligent Systems and Applications 10 (2018) 8-17. doi:10.5815/ijisa.2018.04.02. [29] A. Semenov, et al., Signal Statistic and Informational Parameters of Deterministic Chaos Transistor Oscillators for Infocommunication Systems, in Proceedings of International Scientific-Practical Conference on Problems of Infocommunications Science and Technology, PIC S and T, 2019, pp. 730–734. [30] A.O. Semenov, A.Y. Savytskyi, O.V. Bisikalo, P.I. Kulakov, Mathematical modeling of the two-stage chaotic colpitis oscillator, in Proceedings of 14th International Conference on Advanced Trends in Radioelectronics, Telecommunications and Computer Engineering, TCSET, 2018, April, pp. 835–839. [31] O. Soprun, M. Bublyk, Y. Matseliukh, V. Andrunyk, L. Chyrun, I. Dyyak, A. Yakovlev, M. Emmerich, O. Osolinsky, A. Sachenko, Forecasting Temperatures of a Synchronous Motor with Permanent Magnets Using Machine Learning, CEUR workshop proceedings 2631 (2020) 95–120. [32] M. Zagirnyak, O. Chorna, O. Bisikalo, O. Chornyi, A model of the assessment of an induction motor condition and operation life, based on the measurement of the external magnetic field, in Proceedings of IEEE 3rd International Conference on Intelligent Energy and Power Systems, IEPS, 2018, January, pp. 316–321. [33] A. Podorozhniak, N. Liubchenko, M. Kvochka, I. Suarez, Usage of intelligent methods for multispectral data processing in the field of environmental monitoring, Advanced Information Systems 5(3) (2021) 97–102. doi:10.20998/2522-9052.2021.3.13. [34] R. R. Imanov, A. A. Bayramov, Development of field signal centers based on the modern telecommunication technologies. Advanced Information Systems 4(1) (2020) 136– 139. doi:10.20998/2522-9052.2020.1.21. [35] T. Toosi, M. Sirola, J. Laukkanen, M. van Heeswijk, J. Karhunen, Method for detecting aging related failures of process sensors via noise signal measurement, International Journal of Computing 18(2) (2019) 135-146. [36] A.B. Lozynskyy, I.M. Romanyshyn, B.P. Rusyn, Intensity Estimation of Noise-Like Signal in Presence of Uncorrelated Pulse Interferences, Radioelectronics and Communications Systems 62(5) (2019) 214–222. [37] K. Smelyakov, M. Volk, I.Ruban, M. Derenskyi, A. Chupryna, Short-Range Navigation Radio System Simulator, CEUR Workshop Proceedings 3403 (2023) 539-554. [38] R.M. Peleshchak, O.V. Kuzyk, O.O. Dan'kiv, Non-linear model of impurity diffusion in nanoporous materials upon ultrasonic treatment, Condensed Matter Physics 17(2) (2014) 23601. [39] R.M. Peleshchak, O.V. Kuzyk, O.O. Dan’kiv, Formation of periodic structures under the influence of an acoustic wave in semiconductors with a two-component defect subsystem, Ukrainian Journal of Physics 61(8) (2016) 741-746.