RU2138063C1 - Method of single-point range-finding of sources of atmospherics - Google Patents
Method of single-point range-finding of sources of atmospherics Download PDFInfo
- Publication number
- RU2138063C1 RU2138063C1 RU98117369A RU98117369A RU2138063C1 RU 2138063 C1 RU2138063 C1 RU 2138063C1 RU 98117369 A RU98117369 A RU 98117369A RU 98117369 A RU98117369 A RU 98117369A RU 2138063 C1 RU2138063 C1 RU 2138063C1
- Authority
- RU
- Russia
- Prior art keywords
- signal
- range
- atmosphere
- processing
- terrestrial
- Prior art date
Links
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
Изобретение относится к радиотехническим средствам дальнометрии источников электромагнитного излучения, в частности к способам и устройствам пассивной дальнометрии молниевых разрядов облако - земля преимущественно вертикальной поляризации, и может быть использовано в метеорологии и в гражданской авиации для оперативной грозолокации на расстоянии 300 - 2000 км. The invention relates to radio-frequency means of long-range electromagnetic radiation sources, in particular to methods and devices for passive long-range lightning discharges of cloud-to-ground light of predominantly vertical polarization, and can be used in meteorology and in civil aviation for operational lightning location at a distance of 300 - 2000 km.
Известен способ однопунктовой дальнометрии молниевых разрядов, реализованный в [1]. Этот способ основан на использовании зависимости формы атмосферика как от сигнала в источнике, так и от дальности до источника и заключается в том, что принимают атмосферик на вертикальную электрическую антенну, усиливают его, фильтруют в широкой полосе частот и дискретизируют по времени, запоминают полученный сигнал и проверяют его на соответствие каждой из N заранее установленных градаций дальности Ln (n=1 - N) до молниевого разряда, для этого N-кратно параллельно обрабатывают атмосферик, при каждой n-ой обработке образуют невязку Dn * между атмосфериком и сигналом, полученным при прохождении заранее неизвестным излученным сигналом трассы распространения длиной Ln с известной импульсной переходной функцией распространения (функцией преобразования короткого импульсного сигнала при прохождении трассы распространения), определяют сигнал в источнике, минимизирующий величину Dn * с образованием величины Dn, из полученных Dn образуют обратно пропорциональные величины Pn= 1/Dn, запоминают их, сравнивают по величине, находят обработку с максимальным Pn и соответствующую дальность принимают за дальность до молниевого разряда.A known method of single-point ranging lightning discharges, implemented in [1]. This method is based on the use of the dependence of the shape of the atmosphere on the signal in the source, and on the distance to the source and consists in receiving the atmosphere on a vertical electric antenna, amplifying it, filtering it in a wide frequency band and sampling in time, storing the received signal and check it for compliance with each of the N predetermined gradations of the range L n (n = 1 - N) to the lightning discharge, for this the atmosphere is processed N-times in parallel, at each n-th processing they form a residual D n * between the atmosphere and the signal obtained during the passage of a previously unknown radiated signal of the propagation path of length L n with the known pulse transition function of propagation (the function of converting a short pulse signal during the propagation path), determine the signal in the source that minimizes the value of D n * with the formation of the value of D n D n from the received form inversely proportional value P n = 1 / D n, storing them, is compared in magnitude, are processed with maximum P n and the corresponding range taken as the range to the lightning discharge.
Недостатком указанного способа является большая трудоемкость вычислений, связанная с определением в каждом канале обработки одновременно с дальностью до источника также оптимальной формы сигнала в источнике. The disadvantage of this method is the high complexity of the calculations associated with the determination in each processing channel simultaneously with the range to the source also the optimal waveform in the source.
Наиболее близким к заявленному техническим решением, принятым в качестве прототипа, является способ однопунктового определения дальности до грозовых разрядов, реализованный в [2]. Этот способ основан на представлении атмосферика как совокупности земного сигнала и пространственного (состоящего из нескольких ионосферных) сигнала и использовании зависимостей интервалов времени и отношения амплитуд между каждым из ионосферных и земным сигналами атмосферика от дальности до грозового разряда и состоит в том, что 1) принимают атмосферик, содержащий земной и один или несколько ионосферных сигналов, на вертикальную электрическую антенну, 2) усиливают принятый атмосферик, 3) фильтруют его в широкой полосе частот, 4) дискретизируют по времени с шагом Δt на заранее установленном интервале времени T1, 5) запоминают полученный при этом сигнал Zi (i=1 - I1, где I1 = T1/Δt ) и 6) проверяют его на соответствие каждой из N заранее установленных градаций дальности Ln (n=1 - N) до грозового разряда, для этого при каждой n-й обработке 7) компенсируют в атмосферике ожидаемый пространственный сигнал с помощью земного сигнала с использованием ожидаемых при дальности Ln интервалов времени и отношения амплитуд между каждым из ионосферных и земным сигналами, 8) определяют энергию Dn нескомпенсированного остатка пространственного сигнала на заранее установленном интервале времени, 9) образуют выходной сигнал Pn= 1/Dn, обратно пропорциональный полученному сигналу, и 10) запоминают его, при этом 11) полученные при N обработках выходные сигналы сравнивают между собой по величине и определяют обработку с максимальным выходным сигналом, 12) соответствующую дальность принимают за дальность до грозового разряда.Closest to the claimed technical solution, adopted as a prototype, is a one-point method for determining the distance to lightning discharges, implemented in [2]. This method is based on the representation of the atmosphere as a combination of a terrestrial signal and a spatial (consisting of several ionospheric) signal and the use of the dependences of the time intervals and the amplitude ratio between each of the ionospheric and terrestrial signals of the atmosphere from the distance to a lightning discharge and consists in the following: 1) they accept the atmosphere containing terrestrial and one or more ionospheric signals to a vertical electric antenna, 2) amplify the received atmosphere, 3) filter it in a wide frequency band, 4) discrete are timed in increments of Δt on a predetermined time interval T1, 5) remember the signal Z i obtained in this case (i = 1 - I1, where I1 = T1 / Δt) and 6) check it for compliance with each of the N preset range gradations L n (n = 1 - N) before a lightning discharge, for this, at each nth processing 7), the expected spatial signal is compensated in the atmosphere using an earth signal using the time intervals expected at a range of L n and the amplitude ratio between each of the ionospheric and earth signals, 8) determine the energy D n uncompensated of the remaining spatial signal at a predetermined time interval, 9) form the output signal P n = 1 / D n inversely proportional to the received signal, and 10) store it, while 11) the output signals obtained during N processing are compared with each other in magnitude and processing with a maximum output signal is determined, 12) the corresponding range is taken as the distance to a lightning discharge.
Недостатком указанного способа является погрешность дальнометрии, вызванная различием форм ионосферных и земного сигналов в атмосферике из-за различия условий их распространения и приводящая к неполной компенсации пространственного сигнала атмосферика в канале обработки, соответствующем истинной дальности до грозового разряда. The disadvantage of this method is the error of ranging due to the difference in the forms of ionospheric and terrestrial signals in the atmosphere due to the difference in the conditions of their propagation and leading to incomplete compensation of the spatial signal of the atmosphere in the processing channel corresponding to the true range to a lightning discharge.
Целью настоящего изобретения является увеличение дальности действия и повышение точности однопунктовой дальнометрии источников атмосфериков за счет использования для каждой трассы распространения заранее известной импульсной функции перехода от земного и пространственному сигналу и преобразования с помощью указанной функции земного сигнала в атмосферике в сигнал, который используется для компенсации ожидаемого пространственного сигнала. The aim of the present invention is to increase the range and accuracy of single-point ranging of atmospheric sources through the use for each propagation path of the previously known impulse function of the transition from the terrestrial and spatial signal and conversion using the specified function of the terrestrial signal in the atmosphere into a signal that is used to compensate for the expected spatial signal.
Поставленная цель достигается тем, что в известном способе однопунктового определения дальности до молниевого разряда, выключающем прием атмосферика Z(t), содержащего земной и пространственный сигналы, на вертикальную электрическую антенну, усиление принятого атмосферика, его фильтрацию в требуемой широкой полосе частот и дискретизацию по времени с шагом Δt на заранее установленном интервале времени T1, запоминание полученного при этом сигнала Zi(i=1 - I1, где I1 = T1/Δt) и его проверку на соответствие каждой из N заранее установленных градаций дальности Ln (n=1 - N) до молниевого разряда, для этого N-кратную параллельную обработку атмосферика, включающую при каждой n-ой обработке компенсацию в атмосферике ожидаемого пространственного сигнала с помощью земного сигнала этого же атмосферика путем образования сигнала обратной связи Qn,i и вычитания его из сигнала Zi, причем получаемый при этом разностный сигнал Wn,i используется для образования сигнала Qn,i, вычисление энергии Dn сигнала Wn,i на заранее установленном интервале времени i=I0 - I1 (I0 = T0/Δt, где T0 - ожидаемая длительность земного сигнала), образование выходного сигнала обработки Pn=1/Dn, обратно пропорционального полученной энергии, и его запоминание, сравнение полученных значений Pn (n=1 - N) по величине, нахождение обработки с максимальным Pn и определение дальности до источника как дальность, соответствующую обработке, согласно изобретению при каждой n-ой обработке атмосферика компенсирующий сигнал Qn,i (i=1 - I1) образуют из сигнала Wn,i по правилу
где hn,i - заранее известная импульсная функция перехода от земного сигнала Sn,i к пространственному сигналу Un,i при дальности Ln до источника атмосферика
определяемая соотношением
где Gn0(j2πΔfm) и Gnp(j2πΔfm) - заранее известные комплексные коэффициенты передачи трасс распространения соответственно земного и пространственного сигналов, Δf - шаг по частоте, F - максимальная частота спектра анализируемого атмосферика,
Новым в предложенном способе дальнометрии по сравнению с прототипом является использование при каждой n-ой обработке атмосферика для образования компенсирующего сигнала Qn,i заранее известной импульсной функции hn,i перехода от земного к пространственному сигналу.This goal is achieved by the fact that in the known method of single-point determination of the distance to a lightning discharge, turning off the reception of the atmosphere Z (t), containing terrestrial and spatial signals, to a vertical electric antenna, amplification of the received atmosphere, filtering it in the required wide frequency band and time sampling with a step Δt on a predetermined time interval T1, storing the signal Z i obtained in this case (i = 1 - I1, where I1 = T1 / Δt) and checking it for compliance with each of N preset gradations d of the range L n (n = 1 - N) to the lightning discharge, for this an N-fold parallel processing of the atmosphere, including during each nth processing compensation in the atmosphere of the expected spatial signal using the earth signal of the same atmosphere by generating a feedback signal Q n , i and subtracting it from the signal Z i , and the resulting differential signal W n, i used to generate the signal Q n, i , calculating the energy D n of the signal W n, i over a predetermined time interval i = I0 - I1 (I0 = T0 / Δt, where T0 is the expected duration of the Earth sig ala), the output signal formation processing P n = 1 / D n, is inversely proportional to the received energy, and storing, comparing the obtained values P n (n = 1 - N) in size, processing to find a maximum P n and distance to the source as the range corresponding to the processing according to the invention, for each n-th atmospheric treatment, the compensating signal Q n, i (i = 1 - I1) is formed from the signal W n, i according to the rule
where h n, i is the previously known pulse function of the transition from the terrestrial signal S n, i to the spatial signal U n, i at a distance of L n to the atmospheric source
determined by the relation
where G n0 (j2πΔfm) and G np (j2πΔfm) are the previously known complex transmission coefficients of the propagation paths of terrestrial and spatial signals, respectively, Δf is the frequency step, F is the maximum frequency of the spectrum of the analyzed atmosphere,
New in the proposed method of ranging compared to the prototype is the use of an atmosphere at each nth processing to form a compensating signal Q n, i of the previously known impulse function h n, i of the transition from terrestrial to spatial signal.
На фиг. 1 представлена блок-схема устройства однопунктовой дальнометрии, где обозначено: 1 - электрическая антенна, 2 - усилитель, 3 - фильтр, 4 - первый блок задержки, 5 - аналого-цифровой преобразователь, 6 - первый блок памяти, 7 - пороговый блок, 8 - второй блок задержки, 9 - первый одновибратор, 10 - генератор тактовых импульсов, 11 - 13 - N каналов обработки, 14 - блок определения канала обработки с максимальным выходным сигналом. In FIG. 1 is a block diagram of a single-point ranging device, where it is indicated: 1 — an electric antenna, 2 — an amplifier, 3 — a filter, 4 — a first delay block, 5 — an analog-to-digital converter, 6 — a first memory block, 7 — a threshold block, 8 - the second delay unit, 9 - the first one-shot, 10 - the clock generator, 11 - 13 - N processing channels, 14 - block determining the processing channel with the maximum output signal.
На фиг. 2 представлена блок-схема одного канала обработки, где обозначено: 15 - блок вычитания, 16 - квадратор, 17 - блок накопления, 18 - блок вычисления обратной величины, 19 - второй блок памяти, 20 - вычислитель, 21 - блок постоянной памяти, 22 - третий блок памяти, 23 - второй одновибратор. In FIG. 2 is a block diagram of one processing channel, where it is indicated: 15 — subtraction block, 16 — quadrator, 17 — accumulation block, 18 — reciprocal value calculation block, 19 — second memory block, 20 — calculator, 21 — read-only memory block, 22 - the third memory block, 23 - the second one-shot.
На фиг. 3 представлены примеры импульсных функций hn(t) перехода от земного к пространственному сигналу при ночных условиях распространения и дальностях а) L=500 км и б) L=1000 км.In FIG. Figure 3 shows examples of the impulse functions h n (t) of the transition from the terrestrial to a spatial signal under night propagation conditions and ranges a) L = 500 km and b) L = 1000 km.
На фиг. 4 представлены а) нормированная расчетная форма ночного атмосферика Z(t) на расстоянии L= 500 км от излучателя (приведена в [3], рис. 1.4.7.в) и б) зависимость от L обратной невязки P(L) (полученная из Pn=P(Ln) при соединении соседних точек непрерывной линией), вычисленная для этого атмосферика с помощью предлагаемого способа. Как следует из фиг. 4 б), полученная при этом оценка дальности близка к исходной дальности до источника.In FIG. Figure 4 presents a) the normalized calculated shape of the night atmospheric Z (t) at a distance L = 500 km from the radiator (given in [3], Fig. 1.4.7.c) and b) the dependence of the inverse residual P (L) on L (obtained from P n = P (L n when connecting adjacent points by a continuous line) calculated for this atmosphere using the proposed method. As follows from FIG. 4b), the resulting range estimate is close to the original range to the source.
Сущность предлагаемого способа дальнометрии основана на многолучевом распространении атмосферика в волноводном канале земля - ионосфера [3], при котором атмосферик Z(t) описывается как сумма распространяющегося вдоль поверхности земли земного сигнала S(t) и пространственного сигнала U(t), содержащего один или несколько отраженных от ионосферы сигналов, которые распространяются по нескольким путям: отразившийся один раз от ионосферы - первый ионосферный сигнал, отразившийся от ионосферы, затем от земли и еще раз от ионосферы - второй ионосферный сигнал и т.д., и имеет вид
Z(t)=S(t)+U(t). (4)
При этом задержка начала пространственного сигнала относительно начала земного сигнала равна задержке прихода первого ионосферного сигнала относительно земного сигнала
где R= 6370 км - радиус Земли, H - высота нижней отражающей границы ионосферы (днем H=70 км, ночью H=90 км), c=3•105 км/с - скорость распространения электромагнитного сигнала. При этом используется ограниченность по времени существования молниевого разряда, что приводит к ограниченной длительности T0 земного сигнала (практически T0≤120 мкс).The essence of the proposed ranging method is based on the multipath propagation of the atmosphere in the earth – ionosphere waveguide channel [3], in which the atmosphere Z (t) is described as the sum of the earth signal S (t) propagating along the earth’s surface and the spatial signal U (t) containing one or several signals reflected from the ionosphere, which propagate along several paths: reflected once from the ionosphere - the first ionosphere signal reflected from the ionosphere, then from the earth and again from the ionosphere - the second ionosphere ny signal, etc. and has the form
Z (t) = S (t) + U (t). (4)
The delay in the beginning of the spatial signal relative to the beginning of the earth signal is equal to the delay in the arrival of the first ionospheric signal relative to the earth signal
where R = 6370 km is the radius of the Earth, H is the height of the lower reflecting boundary of the ionosphere (during the day H = 70 km, at night H = 90 km), c = 3 • 10 5 km / s is the speed of propagation of the electromagnetic signal. In this case, the limited time of the existence of a lightning discharge is used, which leads to a limited duration T0 of the earth signal (practically T0≤120 μs).
По предлагаемому способу, как и в прототипе, дальность L до молниевого разряда определяется по форме атмосферика Z(t) путем его проверки на соответствие каждой из N заранее установленных градаций дальности Ln (n=1 - N) до источника, для чего атмосферик N-кратно параллельно обрабатывается, причем при каждой n-ой обработке в анализируемом атмосферике осуществляется компенсация ожидаемого пространственного сигнала земным сигналом с помощью ячейки с обратной связью, содержащей блок вычитания и цепь обратной связи, соединяющей выход блока вычитания с его вторым входом. По предлагаемому способу выходной сигнал ячейки с обратной связью Wn(t) преобразуется в цепи обратной связи в компенсирующий сигнал Qn(t) по правилу
с использованием заранее известной импульсной функции hn(t) перехода от земного сигнала к пространственному сигналу (сами эти сигналы заранее неизвестны) при распространении по трассе длиной Ln. При этом hn(t) определяется соотношением
и зависит от длины трассы (hn(t)=0 и Un(t)=0 при t < τn, где τn - задержка прихода первого ионосферного сигнала относительно земного сигнала при дальности Ln до источника (5)).According to the proposed method, as in the prototype, the range L to the lightning discharge is determined by the shape of the atmosphere Z (t) by checking for compliance with each of the N predetermined gradations of the range L n (n = 1 - N) to the source, for which the atmosphere is N it is processed several times in parallel, with each n-th processing in the atmosphere being analyzed, the expected spatial signal is compensated by an earth signal using a feedback cell containing a subtraction unit and a feedback circuit connecting the output of the subtraction unit I'm with his second entrance. According to the proposed method, the output signal of the feedback cell W n (t) is converted in the feedback circuit into a compensating signal Q n (t) according to the rule
using the previously known impulse function h n (t) of the transition from the terrestrial signal to a spatial signal (these signals themselves are not known in advance) during propagation along a path of length L n . Moreover, h n (t) is determined by the relation
and depends on the path length (h n (t) = 0 and U n (t) = 0 for t <τ n , where τ n is the delay in the arrival of the first ionospheric signal relative to the earth signal at a distance of L n to the source (5)).
В результате компенсации из анализируемого атмосферика образуется сигнал Wn(t)
Wn(t)=Z(t)-Qn(t). (8)
Подставляя в правую часть (8) значение Z(t) из (4) и Qn(t) из (6), получаем для Wn(t) рекуррентное соотношение
В канале обработки, настроенном на истинную дальность до источника излучения анализируемого атмосферика, соотношение (9) с учетом (7) выполняется при Wn(t)= Sn(t), откуда Wn(t)=0 при t>T0. Действие цепи обратной связи соответствует тому, что в каждый момент времени, начиная с t = τn (при t < τn сигнала на выходе цепи обратной связи еще нет: Qn(t)=0 и Wn(t)=Sn(t)), из начальной части сигнала Wn(t) с помощью цепи обратной связи образуется сигнал Qn(t), который, поступая на второй вход блока вычитания, компенсирует пространственный сигнал атмосферика U(t), оставляя на выходе блока вычитания только земной сигнал S(t).As a result of compensation, a signal W n (t) is formed from the analyzed atmosphere
W n (t) = Z (t) -Q n (t). (eight)
Substituting the value Z (t) from (4) and Q n (t) from (6) in the right-hand side of (8), we obtain the recurrence relation for W n (t)
In the processing channel tuned to the true distance to the radiation source of the analyzed atmosphere, relation (9), taking into account (7), is fulfilled for W n (t) = S n (t), whence W n (t) = 0 for t> T0. The action of the feedback circuit corresponds to the fact that at each moment of time, starting from t = τ n (for t <τ n, there is no signal at the output of the feedback circuit: Q n (t) = 0 and W n (t) = S n (t)), from the initial part of the signal W n (t), using the feedback circuit, a signal Q n (t) is formed, which, entering the second input of the subtraction unit, compensates the spatial signal of the atmosphere U (t), leaving at the output of the subtraction block Earth signal S (t) only.
В канале, настроенном на дальность, отличную от дальности до анализируемого атмосферика, из сигнала с выхода блока вычитания образуется сигнал обратной связи, отличающийся от пространственного сигнала атмосферика, в результате компенсации пространственного сигнала атмосферика не происходит и на выходе ячейки с обратной связью выходной сигнал не обращается в нуль после окончания земного сигнала. In a channel tuned to a range other than the range to the analyzed atmosphere, a feedback signal is formed from the output of the subtraction unit, which differs from the spatial signal of the atmosphere, as a result of the compensation of the spatial signal of the atmosphere, the output signal does not apply at the output of the feedback cell to zero after the end of the earth signal.
Импульсная функция hn(t) перехода от земного к пространственному сигналу при n-ой обработке вычисляется следующим образом. Обозначая Gn0(jω) и Gпр(jω) - заранее известные коэффициенты передачи трассы соответственно земного и пространственного сигналов при дальности Ln до источника, Sn(jω) и Un(jω) - неизвестные спектры соответственно земного и пространственного сигналов, V(jω) - неизвестный спектр сигнала в источнике, получаем
Исключая в (10) спектр сигнала в источнике V(jω), получаем соотношение между спектрами пространственного и земного сигналов
Un(jω) = Sn(jω)•Gnp(jω)/Gn0(jω). (11)
При этом, осуществляя преобразование Фурье обеих частей соотношения (7), получаем спектр пространственного сигнала также в виде
Vn(jω) = Sn(jω)•Kn(jω), (12)
где Kn(jω) - спектр функции hn(t). Сравнивая соотношения (11) и (12), получаем
K(jω) = Gпр(jω)/Gn0(jω), (13)
в результате функция hn(t), связанная с Kn(jω) преобразованием Фурье, имеет вид
где F - максимальная частота пропускания входного фильтра (F=35 кГц). (Спектры пространственного и земного сигналов можно получить, используя, например, методы расчета, приведенные в [4]).The impulse function h n (t) of the transition from the earth to the spatial signal during the n-th processing is calculated as follows. Denoting G n0 (jω) and G ave (jω) - pre-known transmission path coefficients of terrestrial and space signals at distance L n to the source, S n (jω) and U n (jω) - unknown spectra respectively terrestrial and space signals, V (jω) is the unknown spectrum of the signal in the source, we obtain
Excluding in (10) the spectrum of the signal in the source V (jω), we obtain the relationship between the spectra of spatial and terrestrial signals
U n (jω) = S n (jω) • G np (jω) / G n0 (jω). (eleven)
At the same time, carrying out the Fourier transform of both parts of relation (7), we obtain the spectrum of the spatial signal also in the form
V n (jω) = S n (jω) • K n (jω), (12)
where K n (jω) is the spectrum of the function h n (t). Comparing relations (11) and (12), we obtain
K (jω) = G ave (jω) / G n0 (jω ), ( 13)
as a result, the function h n (t) associated with the K n (jω) Fourier transform has the form
where F is the maximum transmission frequency of the input filter (F = 35 kHz). (Spectra of spatial and terrestrial signals can be obtained using, for example, the calculation methods given in [4]).
При этом в дискретном по времени представлении функция hn(t) равна
где Δf - шаг по частоте (Δf=2 кГц).Moreover, in a time-discrete representation, the function h n (t) is equal to
where Δf is the frequency step (Δf = 2 kHz).
Получаемый при этом скомпенсированный сигнал атмосферика Wn(t) (7) возводится в квадрат и интегрируется на интервале времени T0<t≤T1, не содержащем земного сигнала, с образованием энергии нескомпенсированного остатка Dn
обратно пропорциональная величина которого Pn=1/Dn является показателем соответствия анализируемого атмосферика дальности Ln до источника, а зависимость Pn от дальности L является решающей функцией - за дальность до источника принимается дальность, соответствующая каналу с максимальным выходным сигналом Pn.The resulting compensated atmospheric signal W n (t) (7) is squared and integrated over the time interval T0 <t≤T1 that does not contain the earth signal, with the formation of the energy of the uncompensated residue D n
the inversely proportional value of which P n = 1 / D n is an indicator of the correspondence of the analyzed atmosphere to the range L n to the source, and the dependence of P n on the range L is a decisive function - the range corresponding to the channel with the maximum output signal P n is taken as the distance to the source.
Дальности Ln, соответствующие соседним обработкам (Ln+1>Ln), устанавливаются такими, чтобы их относительное различие было равно ≈5%. При этом в диапазоне дальностей 300 - 2000 км требуемое количество каналов обработки составляет N≈35.The ranges L n corresponding to adjacent treatments (L n + 1 > L n ) are set such that their relative difference is ≈5%. Moreover, in the range of ranges 300 - 2000 km, the required number of processing channels is N≈35.
Для обработки на ЭВМ осуществляется дискретизация выходящих в приведенные соотношения сигналов по времени с заранее установленным шагом Δt=1 мкс, в результате соотношение (6) переходит в (1)
Как видно из приведенных на фиг. 3 а), б) графиков, импульсная переходная функция h(t) состоит из нескольких разнесенных по времени "всплесков", каждый из которых соответствует одному ионосферному сигналу при увеличении дальности до источника сокращаются интервалы времени между "всплесками" (из-за уменьшения разностей длин путей отдельных ионосферных лучей), увеличиваются амплитуды "всплесков" (из-за увеличения как коэффициентов отражения сигналов от иносферы, так и затухания земного сигнала) и изменяются формы "всплесков". Время обработки атмосферика устанавливается равным T1= 400 мкс, при этом во всем диапазоне дальностей обрабатываемой атмосферик содержит земной сигнал и не менее одного ионосферного сигнала (каждая из функций hn(t) (n=1 - N) содержит на интервале 0 - T1 не менее одного "всплеска") для обеспечения зависимости hn(t) от дальности до источника излучения.For processing on a computer, discretization of the signals outputting to the given ratios by time with a predetermined step Δt = 1 μs is performed, as a result, relation (6) goes to (1)
As can be seen from FIG. 3a), b) of the graphs, the pulse transition function h (t) consists of several “bursts” spaced apart in time, each of which corresponds to one ionospheric signal, while increasing the distance to the source, the time intervals between “bursts” are shortened (due to a decrease in the differences the path lengths of individual ionospheric rays), the amplitudes of the “bursts” increase (due to an increase in both the reflection coefficients of signals from the infosphere and the attenuation of the earth signal) and the shapes of the “bursts” change. The processing time of the atmosphere is set to T1 = 400 μs, while in the entire range of ranges of the processed atmosphere it contains a terrestrial signal and at least one ionospheric signal (each of the functions h n (t) (n = 1 - N) contains in the interval 0 - T1 not less than one “burst”) to ensure the dependence of h n (t) on the distance to the radiation source.
Приведенная на фиг. 4 а) форма ночного атмосферика Z(t) (соответствующего дальности L=500 км до источника) содержит существенно перекрывающиеся по времени земной и первый ионосферный сигналы, а также второй ионосферный сигнал. Полученный с использованием описанного способа и приведенный на фиг. 4 б) график решающей функции P(L) имеет максимум при дальности L≈500 км. Referring to FIG. 4 a) the shape of the night atmosphere Z (t) (corresponding to the distance L = 500 km to the source) contains the terrestrial and first ionospheric signals that substantially overlap in time, as well as the second ionospheric signal. Obtained using the described method and shown in FIG. 4 b) the graph of the decisive function P (L) has a maximum at a range of L≈500 km.
Учет различия форм земного и каждого из ионосферных сигналов с использованием заранее известных частотных характеристик трасс распространения для компенсации ожидаемых ионосферных сигналов атмосферика земным сигналом (в отличие от прототипа, где формы земного и ионосферных сигналов считаются одинаковыми), увеличивает дальность действия и уменьшает погрешности измерения дальности до источника. Taking into account the difference in the forms of the terrestrial and each of the ionospheric signals using pre-known frequency characteristics of the propagation paths to compensate for the expected atmospheric ionospheric signals with a terrestrial signal (in contrast to the prototype, where the forms of the terrestrial and ionospheric signals are considered the same), increases the range and reduces the range measurement errors to source.
Предлагаемый способ определения дальности до источников атмосфериков включает следующую последовательность операций: а) принимают атмосферик на вертикальную электрическую антенну 1 высотой 3 м, б) усиливают его в усилителе 2, в) фильтруют в полосе частот 2 - 35 кГц с помощью фильтра 3, г) для сохранения начальной части атмосферика, предшествующей моменту превышения порога, задерживают сигнал на 20 мкс в первом блоке задержки 4 и д) осуществляют в блоке 5 аналого-цифровое преобразование сигнала по времени с шагом Δt = 1 мкм, е) полученный при этом сигнал Zi (i=1 - I1, где I1=400) запоминают в первом блоке памяти 6 и ж) обрабатывают в решающем устройстве, состоящем из N (N=35) параллельных каналов обработки 11 - 13, каждый из которых настроен на фиксированную дальность Ln до источника излучения в установленном заранее диапазоне дальностей 300 - 2000 км, з) сигнал с выхода фильтра, кроме того, сравнивают с установленным заранее пороговым уровнем в пороговом блоке 7, на выходе которого вырабатывают короткий импульсный сигнал в момент первого превышения атмосфериком порогового уровня, полученный при этом сигнал для синхронизации работы подают параллельно и) на второй вход АЦП, к) на вторые входы каналов обработки и л) задерживают на 500 мкс - время регистрации атмосферика во втором блоке задержки 8, м) сигналом с его выхода запускают одновибратор 9, на время действия которого н) запускают генератор тактовых импульсов 10, вырабатывающий 400 импульсов, о) выходные сигналы тактового генератора подают на третьи входы каналов обработки, п) по окончании обработки атмосферика выходные сигналы каналов обработки сравнивают по величине в блоке 14 и р) на одном из выходов этого блока, соответствующем каналу обработки с максимальным выходным сигналом, вырабатывают стандартный сигнал, при этом в каждом n-ом (n=1 - N) канале обработки сигнал атмосферика с первого входа последовательно а1) пропускают через блок вычитания 15, б1) возводят в квадрат в квадраторе 16, в1) накапливают в накопителе 17 и г1) образуют обратно пропорциональную величину накопленного сигнала в блоке 18, причем выходной сигнал этого блока является выходным сигналом канала обработки, д1) сигнал с выхода блока вычитания, кроме того, запоминают в блоке памяти 19 и е1) обрабатывают в вычислителе 20 совместно с выходными сигналами блока постоянной памяти 21, содержащего импульсную функцию перехода от земного к пространственному сигналу при дальности Ln до источника, ж1) образованный при этом на выходе вычислителя сигнал обратной связи подают на второй вход блока вычитания для компенсации в атмосферике пространственного сигнала, при этом сигнал с второго входа канала обработки, соответствующий моменту начала атмосферика, подают параллельно на первые входы з1) накопителя для его обнуления, и1) второго блока памяти и к1) вычислителя, а также л1) задерживают в третьем блоке задержки 22 на время ожидаемого запаздывания прихода пространственного сигнала относительно земного сигнала и м1) подают на вход второго одновибратора 23, сигналом с выхода которого н1) разрешают по второму входу работу накопителя, при этом тактовые импульсы с третьего входа канала обработки подают параллельно на вторые входы о1) второго блока памяти и п1) вычислителя, а также р1) на третий вход постоянной памяти соединены с соответствующими входами вычислителя.The proposed method for determining the distance to atmospheric sources includes the following sequence of operations: a) take the atmosphere to a vertical electric antenna 1 with a height of 3 m, b) amplify it in the amplifier 2, c) filter it in the frequency band 2 - 35 kHz using filter 3, d) in order to save the initial part of the atmosphere, which precedes the moment of exceeding the threshold, a signal is delayed by 20 μs in the first delay block 4 and e) analog-to-digital conversion of the signal in time with a step Δt = 1 μm is carried out in block 5; nal Z i (i = 1 - I1, where I1 = 400) stored in the first memory unit 6 and g) was treated in the deciding device consisting of N (N = 35) parallel processing channels 11 - 13, each of which is tuned to a fixed the distance L n to the radiation source in a predetermined range of ranges 300 - 2000 km, h) the signal from the output of the filter, in addition, is compared with a predetermined threshold level in the threshold unit 7, the output of which produces a short pulse signal at the moment the atmospheric threshold is first exceeded level obtained with this signal Al to synchronize the work, they are fed in parallel: i) to the second input of the ADC, k) to the second inputs of the processing channels and l) they are delayed by 500 μs - the time of registration of the atmosphere in the second delay unit is 8, m) a single-shot signal 9 is launched from its output, for the duration of which n) start the clock generator 10, generating 400 pulses, o) the output signals of the clock are fed to the third inputs of the processing channels, p) at the end of the processing of the atmosphere, the output signals of the processing channels are compared in value in block 14 and p) on one of the outputs of this block corresponding to the processing channel with the maximum output signal produce a standard signal, while in each n-th (n = 1 - N) processing channel, the atmospheric signal from the first input is sequentially a1) passed through the subtraction block 15, b1) the square in the quadrator 16, b1) is accumulated in the drive 17 and d1) form an inversely proportional value of the accumulated signal in block 18, the output signal of this block being the output signal of the processing channel, e1) the signal from the output of the subtraction block, in addition, remember in the memory block 19 and e1) are processed in the calculator 20 together with the output signals of the permanent memory block 21 containing the pulse function of the transition from the earth to the spatial signal at a distance L n to the source, g1) the feedback signal generated at the output of the calculator is fed to the second the input of the subtraction unit to compensate for the spatial signal in the atmosphere, while the signal from the second input of the processing channel, corresponding to the moment of the beginning of the atmosphere, is fed in parallel to the first inputs of the storage device 1) n, i1) of the second memory block and k1) of the calculator, as well as l1) are delayed in the third delay block 22 for the expected delay of the arrival of the spatial signal relative to the terrestrial signal and m1) are fed to the input of the second one-shot 23, the signal from the output of which n1) is allowed to the second input of the drive, while the clock pulses from the third input of the processing channel are fed in parallel to the second inputs o1) of the second memory block and p1) of the calculator, as well as p1) to the third input of the permanent memory are connected to the corresponding inputs la.
Способ определения дальности до молниевых разрядов может быть реализован в устройстве, блок-схема которого приведена на фиг. 1, где обозначено: 1 - электрическая антенна, 2 - усилитель, 3 - фильтр, 4 - первый блок задержки, 5 - аналого-цифровой преобразователь, 6 - первый блок памяти, 7 - пороговый блок, 8 - второй блок задержки, 9 - одновибратор, 10 - генератор тактовых импульсов, 11 - 13 - N каналов обработки, 14 - блок определения канала обработки с максимальным выходным сигналом, а блок-схема одного канала обработки приведена на фиг. 2, где обозначено: 15 - блок вычитания, 16 - квадратор, 17 - накопитель, 18 - блок вычисления обратной величины, 19 - второй блок памяти, 20 - вычислитель, 21 - блок постоянной памяти, 22 - третий блок памяти, 23 - второй одновибратор. A method for determining the range to lightning discharges can be implemented in a device, a block diagram of which is shown in FIG. 1, where it is indicated: 1 - electric antenna, 2 - amplifier, 3 - filter, 4 - first delay block, 5 - analog-to-digital converter, 6 - first memory block, 7 - threshold block, 8 - second delay block, 9 - one-shot, 10 - clock generator, 11 - 13 - N processing channels, 14 - block determining the processing channel with the maximum output signal, and a block diagram of one processing channel is shown in FIG. 2, where it is indicated: 15 - subtraction unit, 16 - quadrator, 17 - drive, 18 - inverse value calculation unit, 19 - second memory unit, 20 - calculator, 21 - permanent memory unit, 22 - third memory unit, 23 - second single vibrator.
В качестве блоков 2 - 4, 6 - 10, 15 - 19 и 21 - 23 могут быть использованы стандартные блоки, в том числе на интегральных микросхемах, приведенные в [5], в качестве блока 5 - стандартный аналого-цифровой преобразователь, в качестве блока 14 - блок, описанный в [6]. As
Для проверки работоспособности метода дальнометрии, в том числе и вычислителя 20, была использована ЭВМ РС-386. To test the operability of the ranging method, including the
При реализации предлагаемого способа дальнометрии устанавливается:
- электрическая антенна 1 - вертикальная штыревая, длиной 3 м (действующая высота 1,5 м),
- усилитель 2 - линейный, переменного тока, широкополосный с регулируемым коэффициентом усиления 50 - 500,
- фильтр 3 - полосовой с полосой пропускания 2 - 35 кГц,
- длительность задержки сигнала в блоке 4 τa= 20 мкс,
- АЦП 5 10-разрядный с шагом по времени Δt = 1 мкс,
- блок памяти 6 объемом 512 слов,
- длительность задержки сигнала в блоке 8 τb= 500 мкс,
- длительность выходного сигнала одновибратора 9 τb= 4 мс,
- шаг по времени тактового генератора 10 τg= 10 мкс,
- блок памяти 19 объемом 512 слов,
- блок постоянной памяти 21 объемом 512 слов,
- длительность обрабатываемого атмосферика T1 400 мкс,
- ожидаемая длительность земного сигнала атмосферика T0 120 мкс,
- шаг по частоте при расчете hn,1 Δf = 2 кГц,
- диапазон дальностей способа 300 - 2000 км,
- относительное расхождение градаций дальности соседних каналов обработки от средней дальности этих каналов ≈5%,
- количество каналов обработки N=35,
- относительное уменьшение погрешностей дальнометрии по сравнению с прототипом ≈20%.When implementing the proposed ranging method, the following is established:
- electric antenna 1 - vertical whip, 3 m long (effective height 1.5 m),
- amplifier 2 - linear, alternating current, broadband with adjustable gain 50 - 500,
- filter 3 - band pass with a passband of 2 - 35 kHz,
- the duration of the delay signal in the block 4 τ a = 20 μs,
- ADC 5 10-bit with a time step Δt = 1 μs,
-
- the duration of the delay signal in the block 8 τ b = 500 μs,
- the duration of the output signal of the one-shot 9 τ b = 4 ms,
- time step of the clock generator 10 τ g = 10 μs,
- a
- a block of read-
- the duration of the treated atmosphere T1 400 μs,
- the expected duration of the terrestrial atmospheric signal T0 120 μs,
- frequency step in the calculation of h n, 1 Δf = 2 kHz,
- the range of the range of the method 300 - 2000 km,
- the relative difference between the gradations of the distance of the adjacent processing channels from the average range of these channels ≈5%,
- the number of processing channels N = 35,
- the relative decrease in the errors of ranging compared with the prototype ≈20%.
Технический результат использования предложенного способа по сравнению с прототипом состоит в расширении диапазона дальностей и повышении точности однопунктовой дальнометрии источников атмосфериков, что может быть использовано в метеорологии и гражданской авиации. The technical result of using the proposed method in comparison with the prototype is to expand the range and increase the accuracy of single-point ranging of atmospheric sources, which can be used in meteorology and civil aviation.
Литература
1. Епанечников В. А. Способ однопунктового определения дальности до молниевого разряда, МКИ G 01 S 11/00, заявка на патент N 97117580/09 с положительным решением от 22.05.98 г.Literature
1. Yepanechnikov V. A. A one-point method for determining the distance to a lightning discharge, MKI G 01 S 11/00, patent application N 97117580/09 with a positive decision of 05/22/98.
2. Александров М.С., Гапонов И.М., Епанечников В.А., Казаров Ю.В. Устройство для однопунктового определения дальности до источника грозового разряда. Патент РФ N 1799155 от 08.10.92, МКИ G 01 S 11/00. 2. Alexandrov M.S., Gaponov I.M., Yepanechnikov V.A., Kazarov Yu.V. Device for one-point determination of the distance to the source of a lightning discharge. RF patent N 1799155 from 08.10.92, MKI G 01 S 11/00.
3. Кононов И. И., Петренко И.А., Снегуров В.С. Радиотехнические методы местоопределения грозовых очагов. Л.: Гидрометеоиздат, 1986. 3. Kononov I.I., Petrenko I.A., Snegurov V.S. Radio engineering methods for locating thunderclouds. L .: Gidrometeoizdat, 1986.
4. Макаров Г.И., Новиков В.В., Рыбачек С.Т. Распространение радиоволн в волноводном канале Земля - ионосфера. -М.: Наука, 1993. 4. Makarov G.I., Novikov V.V., Rybachek S.T. Propagation of radio waves in the Earth’s waveguide channel - the ionosphere. -M .: Science, 1993.
5. "Аналоговые цифровые интегральные микросхемы"// под ред. Якубовского С.В. -М.: Радио и связь, 1985. 5. "Analog Digital Integrated Circuits" // Ed. Yakubovsky S.V. -M .: Radio and communications, 1985.
6. Устройство определения канала с максимальным уровнем сигнала. Авторское свидетельство 1451609, СССР, МКИ 4 G 01 R 19/04. 6. A device for determining a channel with a maximum signal level. Copyright certificate 1451609, USSR, MKI 4 G 01
Claims (1)
где hn,i - заранее известная импульсная функция перехода от земного сигнала Sn,i к пространственному сигналу Un,i при дальности Ln до источника атмосферика, определяемая соотношением:
где Gno(j2πΔfm) и Gnp(j2πΔfm) - заранее известные комплексные коэффициенты передачи трасс распространения соответственно земного и пространственного сигналов;
Δf - заранее установленный шаг по частоте;
F - максимальная частота спектра анализируемого атмосферика The method for determining the range to lightning discharges, which consists in receiving the atmosphere on a vertical electric antenna, amplifying it, filtering it in the required wide frequency band, sampling it in time with a step Δt on a predetermined time interval TI, memorizing the signal Z i ( i = 1 - I1, where I 1 = TI / Δt) and check it for compliance with each of the N predefined gradations of the range L n (n = 1 - N) to the source, for which the atmosphere is N-times processed in parallel, and for each nth compensation processing the spatial signal expected at a range L n to the source is generated in the atmosphere using the terrestrial signal of the same atmosphere, for this a feedback signal Q n, i is formed and subtracted from the signal Z i , and the difference signal W n, i obtained as a result of such subtraction used to generate the signal Q n, i , while calculating the energy D n of the signal W n, i on a predetermined time interval I0 - I1 (where I0 = TO / Δt, T0 is the expected duration of the earth signal), form the processing output signal P n = 1 / D n, is inversely proportional to D n, memorize it and comparing the obtained values P n (n = 1 - N) in size, are processed with maximum P n and determine the range to the lightning discharge as the distance corresponding to this processing, characterized in that for each n-th processing of the atmospheric compensating signal Q n, i (i = 1 - I1) form according to the rule
where h n, i is the previously known pulse function of the transition from the terrestrial signal S n, i to the spatial signal U n, i at a distance of L n to the atmospheric source, defined by the ratio:
where G no (j2πΔfm) and G np (j2πΔfm) are the previously known complex transmission coefficients of the propagation paths of terrestrial and spatial signals, respectively;
Δf is a predetermined frequency step;
F - maximum spectrum frequency of the analyzed atmosphere
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU98117369A RU2138063C1 (en) | 1998-09-15 | 1998-09-15 | Method of single-point range-finding of sources of atmospherics |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU98117369A RU2138063C1 (en) | 1998-09-15 | 1998-09-15 | Method of single-point range-finding of sources of atmospherics |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2138063C1 true RU2138063C1 (en) | 1999-09-20 |
Family
ID=20210564
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU98117369A RU2138063C1 (en) | 1998-09-15 | 1998-09-15 | Method of single-point range-finding of sources of atmospherics |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2138063C1 (en) |
-
1998
- 1998-09-15 RU RU98117369A patent/RU2138063C1/en not_active IP Right Cessation
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6646587B2 (en) | Doppler radar apparatus | |
Nicolson et al. | Measurement of the intrinsic properties of materials by time-domain techniques | |
Webster | A pulsed ultrasonic distance measurement system based upon phase digitizing | |
EP0535957B1 (en) | Noise correction in radars for diffuse targets | |
US4590477A (en) | Automatic calibration system for distance measurement receivers | |
US4624142A (en) | Double reference pulsed phase locked loop | |
Wohlmuth et al. | Radio frequency fluctuation spectra during the solar conjunctions of the Ulysses and Galileo spacecraft | |
MaMullen et al. | A simple rising-edge detector for time-of-arrival estimation | |
RU2549207C2 (en) | Device for detecting hydroacoustic noise signals based on quadrature receiver | |
US4391142A (en) | Frequency tracked gated pulse technique for ultrasonic frequency | |
RU2138063C1 (en) | Method of single-point range-finding of sources of atmospherics | |
RU2293997C1 (en) | Method for correlation processing of signals, reflected from fast-moving targets | |
RU2460093C1 (en) | Method of measuring distance using sonar | |
RU2212685C2 (en) | Way of single-point detection of position of thunderstorms | |
RU2269789C1 (en) | Method for determining position of electric relay and communication lines disruption and device for realization of said method | |
RU2124739C1 (en) | Method for single-point measuring of distance to lightning | |
Shvets et al. | Lightning location and estimate of the lower ionosphere effective height using by dispersion properties of tweek atmospherics | |
RU2654215C1 (en) | Method of measuring distance by range finder with frequency modulation | |
Challis et al. | Ultrasonic attenuation measurements at very high SNR: Correlation, information theory and performance | |
RU2099747C1 (en) | Process of single-point determination of distance to lightning discharge | |
RU2431852C2 (en) | Radiometric registration of weak broadband radio signal | |
RU2550363C1 (en) | Method of measuring ice thickness and determining properties of medium under ice and device therefor (versions) | |
RU2626068C2 (en) | Method for calibration of parametric tract and device for its implementation | |
RU2042958C1 (en) | Method for determining range to lightning discharge source | |
RU2558654C1 (en) | Device determining direction and distance to signal source |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20070916 |