СПОСОБ ИЗМЕРЕНИЯ КОЭФФИЦИЕНТА ИНТЕРМОДУЛЯЦИИ СИЛЬНО ЗАШУМЛЕННОГО СИГНАЛА
METHOD FOR MEASURING THE INTERMODULATION COEFFICIENT OF A HIGHLY NOISY SIGNAL

Изобретение относится к измерительной технике и системам обработки информации и может быть использовано для измерения малых и сверхмалых значений коэффициента интермодуляции, определяющего нелинейность амплитудной характеристики сейсмоприемников для сильно зашумленного сигнала.
Известен способ измерения нелинейности амплитудной характеристики сейсмоприемников определением его коэффициента нелинейных искажений. Обеспечивается путем подачи на калибровочный вход сейсмоприемника синусоидального сигнала заданной частоты малой нелинейности и вычислением отношения суммы амплитуд второй и последующих гармоник к амплитуде первой гармоники (Патент РФ №2657116, МПК G01V 1/18, 23.05.2017).
Недостатками способа являются:
- большая погрешность измерения коэффициента нелинейности при его малых и сверхмалых значениях, обусловленная нелинейностью входного калибровочного сигнала;
- большая погрешность измерения коэффициента нелинейности при его малых и сверхмалых значениях, обусловленная влиянием на выходной сигнал микросейсмического фона и аппаратурного шума на частотах второй и последующей гармоник;
- используемый метод вычисления частоты и амплитуды оцифрованного синусоидального сигнала чувствителен к зашумленности сигнала.
Известны способы измерения коэффициента интермодуляции радиоэлектронных и акустических устройств (ГОСТ Р 51472, ГОСТ Р 53576). Недостатками способов является ограниченный диапазон значений измеряемого коэффициента интермодуляции (не менее 6% и не менее 1% соответственно).
Прототипом изобретения является способ измерения коэффициента интермодуляции второго порядка (А.Е. Поляков, Л.В, Стрыгин. Методика измерения IP2 и IP3 двухтонового сигнала. Труды МФТИ. - 2012, - Том 4, издание 2). Способ заключается в подаче на вход нелинейного элемента двухтонового (бигармонического) сигнала на близких по значению частотах и определении отношения энергии гармоники на суммарной частоте к энергии сигнала на основной частоте.
Недостатком способа является отсутствие нормированных показателей точности метода и условий их применения, в том числе отсутствует оценка влияния на точность метода измерений относительного уровня шума в области основной и суммарной гармоник, дрейфа или колебаний частоты и нулевого уровня сигнала, погрешностей применяемых в методе средств измерений.
Техническим результатом является возможность проведения измерений малых и сверхмалых коэффициентов интермодуляции с приемлемой точностью в задаче измерений нелинейности амплитудной характеристики по сильно зашумленным сигналам.
и
1+2
и
1+2
1
2
Технический результат достигается тем, что при определении малых и сверхмалых значений коэффициента интермодуляции К преобразовательного тракта исследуемых приборов, который определяет нелинейность их амплитудной характеристики, измерения проводятся генераторным методом на частоте суммарной гармоники F, коэффициент интермодуляции К определяют как отношение значений оценки амплитуды на частоте суммарной гармоники F к сумме значений оценки значений амплитуд на частотах F и F
(F1+F2)
(F1)
(F2)
1+2
1
2
где A, A, A - оценки амплитуд синусоидальных сигналов на частоте суммарной гармоники F, частотах F и F, определяемые по формуле
1
2
где значения A и А находят решением системы уравнений
11
1
12
2
1
а A + а А=b
21
1
22
2
2
a A + а A=b,
11
2
где а=∑sin(2πk Fc/Fd);
12
a=∑sin(2πk Fc/Fd) cos(2πk Fc/Fd);
1
k
b=∑X sin(2πk Fc/Fd);
21
12
a=a;
22
11
a=N-a;
2
k
b=∑X cos(2πk Fc/Fd);
1+2
1
2
Fc - одна из частот F, F, F;
Fd - частота дискретизации данных регистрации;
N - количество мгновенных значений в данных регистрации;
Xk - мгновенные значения напряжения из набора из N дискретных значений сигнала, при k = 1, …, N;
учитывают влияние шумов, содержащихся в сигнале с выхода канала регистрации исследуемого прибора на результат измерения Ки;
1
2
1+2
оценку уровня шума и определение относительных погрешностей амплитуд синусоидального сигнала на частотах F, F и частоте суммарной гармоники F проводят в три этапа:
1
2
3
4
1+2
- этап 1 - определяют границы участков спектра, для чего проводят оценку уровня шума в области каждого спектрального пика в области спектра справа и слева от пика, лежащей в границах от 0,25 до 0,75 расстояния между соседними пиками или от 0,02 до 0,05 значения частоты гармоники, с учетом ширины пика, далее определяют параметры границ участков для каждого из спектральных пиков на каждой из частот F, F, F, F, F по формуле
ush
ush
где K - коэффициент уширения пика в результате применения окна Блэкмана, K=3;
затем определяют нижнюю (low_) и верхнюю (high_) границы левой (_left) и правой (_right) областей каждого из спектральных пиков по формулам
3
1
где F - частота первой гармоники для частоты F, Гц;
4
2
где F - частота первой гармоники для частоты F, Гц;
- этап 2 - производят оценку среднего уровня шума для частоты Fc с использованием среднего значения спектральной плотности мощности записи бигармонического сигнала на участках спектра, полученных на этапе 1 для данной частоты, по формулам
1
2
где N, N - количество значений спектральной плотности мощности в левой и правой области соответственно;
noise
1
2
3
4
1+2
A - оценка амплитуды приведенного к выходу исследуемого прибора синусоидального сигнала, формируемого шумом на каждой из частот F, F, F, F, F, которая могла быть получена при ее вычислении методом наименьших квадратов, В (далее - амплитуда шума);
1
2
2
,  - наборы значений спектральной плотности мощности в левой и правой области при r = 1, …, N; при t = 1, …, N; соответственно, В/Гц, вычисляют по формуле
k1
N
где X(2πjf) - вектор, отсч., полученный из набора {Xk} при k = 1, …, N дискретным преобразованием Фурье ДПФ{Xk} по всей выборке по формуле
N3 - количество частот в спектре;
Noise
i
i
- этап 3 - непосредственно вычисляют относительные погрешности амплитуд синусоидального сигнала, вносимые шумом, δ(ΔA(F)) на каждой частоте F по формуле
i
i
где A(F) - амплитуда приведенного к выходу исследуемого прибора синусоидального сигнала на частоте F, В, вычисленная по формуле
i
1
2
1+2
F, - одна из частот F, F, F
1
2
- множители при r = N + N - 1 для нормального закона распределения, Р=0,95
U - квантиль для нормального закона распределения, Р=0,95
Затем вычисляют статистическую погрешность счета
где  - множители при N значений зарегистрированных событий для нормального закона распределения, Р=0,95
и определяют относительную погрешность смещения оценки амплитуды δAFc, %, на заданной частоте Fc от относительной неопределенности измерения частоты синусоиды δFc, %, и количества периодов n⋅р заданной частоты Fc, Гц, в анализируемой выборке
1
2
1+2
и
затем по результатам оценки влияния уровней шума и относительных погрешностей амплитуд синусоидального сигнала на частотах F, F и на частоте суммарной гармоники F, статистической погрешности счета и относительной погрешности смещения оценки амплитуды δAFc определяют либо значение коэффициента интермодуляции K,
Kи
и вычисляют границы расширенной относительной неопределенности δ, %, измерения коэффициента интермодуляции
Kи
либо при δ больше установленного предельно допустимого значения,
вычисляют обусловленное шумом верхнее предельное значение коэффициента интермодуляции  по формуле
1
2
3
где A, А, А - амплитуды синусоид, вычисленные методом наименьших квадратов для частот основных и суммарной гармоник соответственно.