RU2263929C1 - Способ исследования нелинейности преобразования сигналов объектом - Google Patents

Способ исследования нелинейности преобразования сигналов объектом Download PDF

Info

Publication number
RU2263929C1
RU2263929C1 RU2004110640/09A RU2004110640A RU2263929C1 RU 2263929 C1 RU2263929 C1 RU 2263929C1 RU 2004110640/09 A RU2004110640/09 A RU 2004110640/09A RU 2004110640 A RU2004110640 A RU 2004110640A RU 2263929 C1 RU2263929 C1 RU 2263929C1
Authority
RU
Russia
Prior art keywords
signal
test signal
signals
linear
additional test
Prior art date
Application number
RU2004110640/09A
Other languages
English (en)
Inventor
Э.В. Семенов (RU)
Э.В. Семенов
Original Assignee
Томский Государственный Университет Систем Управления И Радиоэлектроники (Тусур)
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Томский Государственный Университет Систем Управления И Радиоэлектроники (Тусур) filed Critical Томский Государственный Университет Систем Управления И Радиоэлектроники (Тусур)
Priority to RU2004110640/09A priority Critical patent/RU2263929C1/ru
Application granted granted Critical
Publication of RU2263929C1 publication Critical patent/RU2263929C1/ru

Links

Images

Landscapes

  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

Изобретение относится к областям радиотехники и акустики и может быть использовано для обнаружения нелинейных объектов. Сущность изобретения состоит в том, что на объект воздействуют тестовым сигналом, принимают от объекта сигнал-отклик на тестовый сигнал, после окончания сигнала-отклика на объект воздействуют дополнительным тестовым сигналом, осуществляют задержку сигнала-отклика на тестовый сигнал в линии задержки, время запаздывания сигнала в которой выбирают равным интервалу времени между началом тестового сигнала и началом дополнительного тестового сигнала, затем выходной сигнал линии задержки пропускают через первый линейный фильтр, а сигнал-отклик объекта на дополнительный тестовый сигнал пропускают через второй линейный фильтр и регистрируют разность сигналов на выходе первого и второго линейных фильтров, при этом отношение передаточных характеристик первого и второго линейных фильтров выбирают равным отношению спектров дополнительного тестового сигнала и тестового сигнала, причем при вычислении спектра тестового сигнала за нулевой момент времени принимают начало тестового сигнала, а при вычислении спектра дополнительного тестового сигнала за нулевой момент времени принимают начало дополнительного тестового сигнала. В частном случае разность выходных сигналов первого и второго линейных фильтров пропускают через третий линейный фильтр. Достигаемый технический результат - исследование нелинейности преобразования сигналов объектами, вносящими сложные линейные искажения сигнала, в том числе изменяющиеся во времени. 1 з.п. ф-лы, 6 ил.

Description

Изобретение относится к областям радиотехники и акустики и может быть использовано для обнаружения нелинейных объектов и исследования нелинейности преобразования ими сигналов.
Известны способы исследования нелинейности преобразования сигналов объектом, заключающиеся в том, что исследуемый объект подвергают воздействию синусоидального сигнала с последующей оценкой в сигнале-отклике от объекта спектральных составляющих на частотах гармоник или одновременному воздействию двух синусоидальных сигналов, имеющих разные частоты, с последующей оценкой сигнала на комбинационных частотах [1]. Их недостатком является то, что с их помощью нельзя исследовать нелинейность преобразования объектом импульсных или шумоподобных сигналов (нелинейность объектов проявляется по-разному при воздействии на них разных сигналов).
Известен способ исследования нелинейности преобразования сигналов объектом, заключающийся в том, что исследуемый объект подвергают воздействию тестового узкополосного шума, и в сигнале-отклике от объекта проводят анализ возникающих в результате нелинейного преобразования спектральных составляющих вне полосы частот тестового шума [2]. Недостатком данного способа является то, что тестовый шум занимает только часть анализируемой полосы частот (использование широкополосных, занимающих всю анализируемую полосу частот, тестовых сигналов невозможно).
Известен также способ исследования нелинейности преобразования сигналов объектом, заключающийся в том, что на исследуемый объект воздействуют тестовым сигналом, образованным изменяющимся во времени каким-либо физическим фактором, и в исследуемом отклике определяют величину физического фактора на частоте или частотах одного или нескольких локальных нулей спектральной плотности мощности исходного тестового сигнала [3]. Недостатком этого способа является то, что регистрации доступна только небольшая часть продуктов нелинейного преобразования в окрестности одного или нескольких нулей спектральной плотности мощности исходного тестового сигнала.
Наиболее близким по технической сущности к заявляемому способу исследования нелинейности преобразования сигналов объектом является способ, заключающийся в воздействии на объект тестовым сигналом и регистрации разности тестового сигнала и сигнала-отклика от объекта [4]. Недостатком способа-прототипа является то, что при использовании этого способа требуется компенсация линейных искажений тестового сигнала объектом, что реально только для стабильных во времени линейных искажений с простой частотной зависимостью.
Технический результат, на достижение которого направлено предлагаемое решение, - создание способа, позволяющего исследовать нелинейность преобразования сигналов объектами, вносящими сложные линейные искажения сигнала, в том числе изменяющиеся во времени.
Это достигается тем, что в известном способе исследования нелинейности преобразования сигналов объектом, включающем воздействие на него тестовым сигналом и прием от объекта сигнала-отклика на тестовый сигнал, после окончания сигнала-отклика на объект воздействуют дополнительным тестовым сигналом, осуществляют задержку сигнала-отклика на тестовый сигнал в линии задержки, время запаздывания сигнала в которой выбирают равным интервалу времени между началом тестового сигнала и началом дополнительного тестового сигнала, затем выходной сигнал линии задержки пропускают через первый линейный фильтр, а сигнал-отклик объекта на дополнительный тестовый сигнал пропускают через второй линейный фильтр и регистрируют разность сигналов на выходе первого и второго линейных фильтров, при этом отношение передаточных характеристик первого и второго линейных фильтров выбирают равным отношению спектров дополнительного тестового сигнала и тестового сигнала, причем при вычислении спектра тестового сигнала за нулевой момент времени принимают начало тестового сигнала, а при вычислении спектра дополнительного тестового сигнала за нулевой момент времени принимают начало дополнительного тестового сигнала. В частном случае разность выходных сигналов первого и второго линейных фильтров пропускают через третий линейный фильтр.
Использование двух тестовых сигналов позволяет отличить их нелинейное преобразование от линейного, в том числе и в случае, когда линейные искажения тестового сигнала объектом имеют сложную частотную зависимость и нестабильны во времени. Ниже приводится обоснование этого утверждения.
Нелинейность преобразования сигнала в частотной области проявляется в виде отклонения от уравнения линии:
u(ω)=К(ω)х(ω),
где u(ω) и х(ω) - спектры Фурье сигнала-отклика от объекта и тестового сигнала соответственно, К(ω) - постоянный на каждой конкретной частоте комплексный коэффициент (передаточная характеристика объекта), ω - круговая частота.
Воздействуем на объект двумя сигналами - тестовым x1(t) и дополнительным тестовым x2(t), где t - время. Сигналы-отклики объекта u1(t) и u2(t) (на тестовый сигнал и на дополнительный тестовый сигнал соответственно) пропустим через линейные фильтры с передаточными характеристиками К1(ω) и К2(ω) соответственно и найдем разность ε(ω) выходных сигналов фильтров:
ε(ω)=K1(ω)u1(ω)-K2(ω)u2(ω), (1)
где u1(ω) и u2(ω) - спектры Фурье сигналов-откликов на тестовый сигнал и на дополнительный тестовый сигнал соответственно.
Если
K1(ω)/K2(ω)=x2(ω)/x1(ω), (2)
то ε(ω) можно использовать в качестве характеристики нелинейности преобразования сигналов объектом, поскольку при линейном преобразовании сигналов объектом ε(ω)=0. Если преобразование сигналов линейно, то
ε(ω)=K1(ω)u1(ω)-К2(ω)u2(ω)=K1(ω)K(ω)x1(ω)-K2(ω)K(ω)х2(ω);
ε(ω)+K2(ω)K(ω)x2(ω)=K1(ω)K(ω)x1(ω);
ε(ω)/[K(ω)K2(ω)x1(ω)]+x2(ω)/x1(ω)=K1(ω)/K2(ω).
При выполнении (2)
ε(ω)/[K(ω)K2(ω)x1(ω)]=0. (3)
Если исключить физически нереальный случай K(ω)K2(ω)x1(ω)=∝, то из (3) следует ε(ω)=0.
Ограничений на передаточную характеристику объекта, кроме ее постоянства в течение времени измерения, не накладывается, поэтому предлагаемый способ позволяет отличить нелинейное преобразование тестового сигнала от линейного и исследовать нелинейность преобразования сигнала при наличии сложных линейных искажений сигнала, причем допустимо, чтобы к моменту следующего измерения линейные искажения изменились. На тестовые сигналы не устанавливается ограничений, поэтому тестовые сигналы могут быть широкополосными.
Если почленно взять обратное преобразование Фурье от (1), то получим временнýю зависимость сигнала ε(t):
ε(t)=h1(t)*u1(t)-h2(t)*u2(t). (4)
Здесь h1(t) и h2(t) - импульсные характеристики первого и второго линейных фильтров соответственно, * - символ свертки, h1(t)*u1(t) и h2(t)*u2(t) - сигналы на выходе первого и второго линейного фильтров соответственно при подаче на их входы соответственно сигналов-откликов объекта на тестовый сигнал и на дополнительный тестовый сигнал.
Сигналы-отклики объекта на тестовый и дополнительный тестовый сигналы необходимо регистрировать раздельно, поэтому воздействие дополнительным тестовым сигналом начинают после окончания сигнала-отклика на тестовый сигнал, а для того, чтобы зарегистрировать разность пропущенных через фильтры откликов в соответствии с (4), сигнал-отклик на тестовый сигнал подают на вход первого линейного фильтра через линию задержки, время запаздывания в которой выбирают равным интервалу времени между началом тестового сигнала и началом дополнительного тестового сигнала.
В частном случае сигнал ε(t) пропускают через третий линейный фильтр. С помощью этого фильтра преобразуют форму сигнала ε(t) (например, концентрируют сигнал во времени), для того чтобы он стал более наглядным или лучше заметным на фоне шумов. Также этот фильтр используют для подавления появляющихся в результате нелинейного преобразования спектральных составляющих, лежащих выше исследуемого диапазона частот, наличие которых нежелательно при регистрации сигнала ε(t) с помощью аналого-цифрового преобразователя.
На фиг.1 представлена амплитудная характеристика, на фиг.2 - тестовые сигналы, на примере которых иллюстрируется предложенный способ, на фиг.3 тестовые сигналы представлены последовательно во времени так, как они воздействуют на объект, на фиг.4 представлены отклики объекта на тестовые сигналы фиг.3, на фиг.5 представлены задержанный в линии задержки отклик объекта на тестовый сигнал и отклик объекта на дополнительный тестовый сигнал (входные сигналы первого и второго линейных фильтров соответственно) и на фиг.6 проиллюстрировано получение сигнала ε(t). Сплошной линией на фиг.2 и 3 представлен тестовый сигнал, пунктирной линией - дополнительный тестовый сигнал. Сплошной линией на фиг.4 обозначен отклик объекта на тестовый сигнал, пунктирной линией - отклик объекта на дополнительный тестовый сигнал. На фиг.5 сплошной линией обозначен задержанный в линии задержки отклик объекта на тестовый сигнал, пунктирной линией - отклик объекта на дополнительный тестовый сигнал. На фиг.6 штриховой линией обозначен сигнал на выходе первого линейного фильтра, пунктирной линией - сигнал на выходе второго линейного фильтра и сплошной линией - разность сигналов на выходах первого и второго линейных фильтров.
Ниже приведен пример осуществления изобретения.
В качестве тестового сигнала выбираем импульс Гаусса x1(t)=exp[-((t-t0)/τ1)2]. При τ1=30 мкс, t0=70 мкс получим сигнал, показанный на фиг.2 сплошной линией. В качестве дополнительного тестового сигнала выбираем первую производную импульса Гаусса x2(t)=-[2τ(t-t0)/τ12]exp[-((t-t0)/τ1)2]. При τ=17 мкс получим сигнал, представленный на фиг.2 пунктирной линией. Соотношение спектров дополнительного тестового и тестового сигналов x2(ω)/x1(ω)=jωτ. Для выполнения (2) выбираем передаточную характеристику первого линейного фильтра K1(ω)=jωτ и передаточную характеристику второго линейного фильтра К2(ω)=1.
Тестовые сигналы генерируем в электрическом виде генератором и последовательно во времени так, как представлено на фиг.3. Интервал времени между началом тестового сигнала и началом дополнительного тестового сигнала выбираем равным 0,5 мс. Затем подаем тестовый и дополнительный тестовый сигналы на преобразователь, который преобразует их в физический фактор, воздействующий на объект.
Природа физического фактора может быть различна. Если исследуется нелинейность преобразования сигналов проходными электро- и радиотехническими устройствами (четырехполюсниками), то в качестве физического фактора используется ток или напряжение. Если решается задача дистанционного обнаружения объекта с нелинейными свойствами, то в качестве физического фактора выступает электрическое или магнитное поле. Если исследуются нелинейности преобразования звуковых волн в жидкостях или газах, то в качестве физического фактора используются акустические колебания. Сигналы-отклики объекта преобразуются в электрические сигналы соответствующим преобразователем.
Действие способа рассмотрим на примере объекта с амплитудной характеристикой u(х):
Figure 00000002
представленной на фиг.1, где х и u - тестовый сигнал и отклик объекта на него соответственно.
Принимаем сигналы-отклики от объекта. Преобразованные в электрический вид сигналы-отклики объекта на тестовый сигнал и дополнительный тестовый сигнал (представлены на фиг.4) подаем соответственно на вход линии задержки и вход второго линейного фильтра. Время запаздывания сигнала в линии задержки выбираем равным интервалу времени между началом тестового сигнала и началом дополнительного тестового сигнала, то есть 0,5 мс. На фиг.5 представлены сигнал-отклик на дополнительный тестовый сигнал и задержанный в линии задержки сигнал-отклик на тестовый сигнал. Выходной сигнал линии задержки подаем на вход первого линейного фильтра.
Поскольку передаточная функция второго линейного фильтра выбрана тождественно равной единице, выходной сигнал этого фильтра будет равен входному, то есть u2(t). Дополнительный тестовый сигнал укладывается в линейную часть амплитудной характеристики, поэтому форма сигнала-отклика объекта u2(t) (представленного на фиг.4 пунктирной линией) повторяет форму и амплитуду x2(t). Отсюда следует, что выходной сигнал второго линейного фильтра (фиг.6, пунктирная линия) будет совпадать по форме и амплитуде с x2(t). Выходной же сигнал первого линейного фильтра h1(t)*u1(t) совпадает с x2(t) только в том случае, если его входной сигнал u1(t) совпадает с x1(t):
h1(t)*u1(t)=h1(t)*x1(t)=F-1[K1(ω)x1(ω)]=F-1[(x2(ω)/x1(ω))x1(ω)]=x2(t),
где F-1 - обратное преобразование Фурье.
Однако поскольку амплитудная характеристика нелинейна, u1(t) отличается от x1(t), и выходной сигнал первого линейного фильтра (фиг.6, штриховая линия) отличается по форме от x2(t), то между выходными сигналами первого и второго линейных фильтров появляется разностный сигнал ε(t) (фиг.6, сплошная линия).
Источники информации
1. Шкритек П.Справочное руководство по звуковой схемотехнике: Пер. с нем. - М.: Мир, 1991. С.32.
2. Журавлев В. М. Метод измерения нелинейных искажений с помощью полос шума: Канд. дисс. Л.: ЛИКИ, 1967.
3. Семенов Э.В. Синтез сверхширокополосных тестовых сигналов для обнаружения нелинейного преобразования сигнала путем наблюдения за нулями его спектральной плотности мощности // Сборник докладов Всероссийской научной конференции "Сверхширокополосные сигналы в радиолокации, связи и акустике". Муром, 1-3 июля 2003 г. - Муром: Изд.-полиграфический центр МИ ВлГУ, 2003. - С.99-103.
4. Акулиничев И. Селекция сигнала искажений // Радио, 1983, №10. С.42-44 - прототип.

Claims (2)

1. Способ исследования нелинейности преобразования сигналов объектом, включающий воздействие на объект тестовым сигналом x1(t), где t - время и прием от объекта сигнала-отклика на тестовый сигнал, отличающийся тем, что после окончания сигнала-отклика на объект воздействуют дополнительным тестовым сигналом x2(t), осуществляют задержку сигнала-отклика на тестовый сигнал в линии задержки, время запаздывания сигнала в которой выбирают равным интервалу времени между началом тестового x1(t) сигнала и началом дополнительного тестового сигнала x2(t), затем выходной сигнал линии задержки U1(t) пропускают через первый линейный фильтр с передаточной характеристикой К1(ω), а сигнал-отклик U2(t) объекта на дополнительный тестовый сигнал x2(t) пропускают через второй линейный фильтр с передаточной характеристикой К2(ω) и регистрируют разность сигналов на выходе первого и второго линейных фильтров ε(ω)=K1(ω)·U1(ω)-K2(ω)·U2(ω), где U1(ω) и U2(ω) - спектры Фурье сигналов U1(t) и U2(t), при этом если К1(ω)/К2(ω)=х2(ω)/х1(ω), где х2(ω)/х1(ω) - отношение спектров дополнительного тестового сигнала x2(t) и тестового сигнала x1(t), то ε(ω) используют в качестве характеристики нелинейности преобразования сигналов объектом, причем при вычислении спектра тестового сигнала x1(t) за нулевой момент времени принимают начало тестового сигнала x1(t), а при вычислении спектра дополнительного тестового сигнала x2(t) за нулевой момент принимают начало дополнительного тестового сигнала.
2. Способ по п.1, отличающийся тем, что разность выходных сигналов первого и второго линейных фильтров пропускают через третий линейный фильтр.
RU2004110640/09A 2004-04-07 2004-04-07 Способ исследования нелинейности преобразования сигналов объектом RU2263929C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2004110640/09A RU2263929C1 (ru) 2004-04-07 2004-04-07 Способ исследования нелинейности преобразования сигналов объектом

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2004110640/09A RU2263929C1 (ru) 2004-04-07 2004-04-07 Способ исследования нелинейности преобразования сигналов объектом

Publications (1)

Publication Number Publication Date
RU2263929C1 true RU2263929C1 (ru) 2005-11-10

Family

ID=35865498

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2004110640/09A RU2263929C1 (ru) 2004-04-07 2004-04-07 Способ исследования нелинейности преобразования сигналов объектом

Country Status (1)

Country Link
RU (1) RU2263929C1 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2498341C1 (ru) * 2012-05-23 2013-11-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Нижегородский государственный технический университет им. Р.Е. Алексеева" (НГТУ) Способ обнаружения объектов, содержащих нелинейные элементы

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
АКУЛИНИЧЕВ И. Селекция сигнала искажений. Радио, 1983, №10, с.42-44. *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2498341C1 (ru) * 2012-05-23 2013-11-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Нижегородский государственный технический университет им. Р.Е. Алексеева" (НГТУ) Способ обнаружения объектов, содержащих нелинейные элементы

Similar Documents

Publication Publication Date Title
Blockley et al. Mixer-based, vector-corrected, vector signal/network analyzer offering 300kHz-20GHz bandwidth and traceable phase response
CN101600144B (zh) 采用连续对数扫频信号获得电声产品多个参数的方法及系统
EP0035391A2 (en) Multi frequency eddy current test apparatus and method
US9823281B2 (en) Time frame measuring method with calibration in the frequency range
CN109975771B (zh) 基于信号三阶相位差分的宽带数字信道化方法
RU2263929C1 (ru) Способ исследования нелинейности преобразования сигналов объектом
Morrison et al. Performance of oversampled polyphase filterbank inversion via Fourier transform
Burrascano et al. A swept-sine pulse compression procedure for an effective measurement of intermodulation distortion
JP3672847B2 (ja) レーダ装置及びコヒーレント積分方法
Liou et al. Digital wideband phased array calibration and beamforming using time reversal technique
RU2405163C1 (ru) Способ частотно-временного корреляционного анализа цифровых сигналов
US20150264682A1 (en) Method and apparatus for processing a signal
Gawasawa et al. Vector measurement of nonlinear signals generated in RF SAW/BAW devices
JP3575252B2 (ja) 電磁超音波計測方法および装置
RU2423780C2 (ru) Способ частотно-избирательной фильтрации сигнала, используя пространственно-временной алгоритм проведения измерений энергетической величины волнового пакета на многоотводной линии задержки, обладающей дисперсией для длин волн рабочего диапазона, и суммированием результатов выполненных измерений
JP4334215B2 (ja) 電磁超音波計測方法および装置
Morise et al. Warped‐TSP: An acoustic measurement signal robust to background noise and harmonic distortion
EP3889627A1 (en) Magnetism detection device and magnetism detection method
RU2291463C2 (ru) Способ аналого-дискретной обработки радиолокационных импульсных сигналов
Gunaratne Beamforming of temporally broadband bandpass plane waves using 2D FIR trapezoidal filters
Suzuki et al. Three Dimensional Displays for Demonstrating Transient Characteristics of Loudspeakers
RU2256934C1 (ru) Способ спектрально-временной трансформации сигналов
RU2580410C1 (ru) Устройство для измерения больших токов
Dadic Wiener filtering for anechoic transfer function measurement in acoustics
Shapir Suggestion for a new formula to calculate group-delay from frequency domain measurements

Legal Events

Date Code Title Description
QB4A Licence on use of patent

Effective date: 20070529

QB4A Licence on use of patent

Effective date: 20090728

MM4A The patent is invalid due to non-payment of fees

Effective date: 20150408