МИ 2090-90. Рекомендация. Государственная система обеспечения единства измерений. Определение динамических характеристик линейных аналоговых средств измерений с сосредоточенными параметрами. Общие положения

РЕКОМЕНДАЦИЯ

ГОСУДАРСТВЕННАЯ СИСТЕМА ОБЕСПЕЧЕНИЯ ЕДИНСТВА ИЗМЕРЕНИЙ

ОПРЕДЕЛЕНИЕ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ЛИНЕЙНЫХ АНАЛОГОВЫХ СРЕДСТВ ИЗМЕРЕНИЙ С СОСРЕДОТОЧЕННЫМИ ПАРАМЕТРАМИ

ОБЩИЕ ПОЛОЖЕНИЯ

МИ 2090-90

Группа Т.88

Дата введения 1 августа 1991 года

ИНФОРМАЦИОННЫЕ ДАННЫЕ

1. РАЗРАБОТАНА И ВНЕСЕНА Государственным комитетом СССР по управлению качеством продукции и стандартам

ИСПОЛНИТЕЛИ

В.А. Грановский, канд. техн. наук (руководитель темы), Ю.С. Этингер, С.Л. Васильева

2. УТВЕРЖДЕНА ПРИКАЗОМ НПО «Система»

3. ЗАРЕГИСТРИРОВАНА ВНИИМС 04.09.90

4. ВЗАМЕН РД 50-404-83

Настоящая рекомендация распространяется на методы определения динамических характеристик линейных аналоговых средств измерений с сосредоточенными параметрами (далее — исследуемых средств измерений) и устанавливает основные правила выбора и реализации методов. Методы определения динамических характеристик соответствуют одному из следующих общих методов по ГОСТ 8.061:

прямое измерение исследуемым средством измерений испытательного сигнала, воспроизводимого образцовым генератором;

косвенные измерения испытательного сигнала, воспроизводимого генератором.

Настоящий документ рекомендуется применять при контроле выполнения требований ГОСТ 8.256. Термины, используемые в настоящей рекомендации, и пояснения к ним по МИ 1951.

1. ОСНОВНЫЕ ПОЛОЖЕНИЯ

1.1. Для определения динамической характеристики исследуемого средства измерений необходимо (см. рис. 1 а, б):

выбрать испытательный сигнал;

выбрать генератор испытательного сигнала и средство измерений, предназначенное для определения выходного сигнала исследуемого средства измерений;

подать на вход исследуемого средства измерений испытательный сигнал;

определить выходной сигнал исследуемого средства измерений, представляющий собой отклик на испытательный сигнал;

обработать данные по известному алгоритму с целью нахождения оценки требуемой динамической характеристики и погрешности ее определения.

Примечания:

1. При отсутствии необходимого генератора допускается использовать в совокупности источник сигнала, не являющийся мерой, и средство измерений, предназначенное для определения испытательного сигнала. В указанном случае, кроме перечисленных выше операций, необходимо определить испытательный сигнал.

2. Если исследуемое средство измерений — управляемая мера, то генератором испытательного сигнала служит источник управляющих воздействий (задающих кодов).

Рис. 1. Структурные схемы определения динамических характеристик средств измерений: ГИС — образцовый генератор испытательного сигнала (образцовая динамическая мера); СИ — исследуемое средство измерений; УОД — устройство обработки данных, включающее средство измерений, предназначенное для определения выходного сигнала исследуемого средства измерений; ГС — генератор сигнала; ОСИ — образцовое (быстродействующее) средство измерений

1.2 В качестве данных для нахождения динамической характеристики и погрешности ее определения следует использовать:

совокупность оценок испытательного и выходного сигналов исследуемого средства измерений и погрешностей указанных оценок;

сведения об условиях эксперимента по определению динамической характеристики исследуемого средства измерений;

предварительные сведения об исследуемом средстве измерений.

1.2.1 Оценка испытательного сигнала должна быть представлена либо:

номинальным аналитическим выражением и погрешностью его воспроизведения,

номинальными значениями в известные моменты времени и пределами допускаемых отклонений от указанных значений,

пределами допускаемых значений в известные моменты времени (аналитическими выражениями предельных значений),

значениями в известные моменты времени и погрешностями определения указанных значений.

1.2.2. Оценка выходного сигнала исследуемого средства измерений должна быть представлена либо значениями в известные моменты времени и погрешностями измерений указанных значений, либо аппроксимирующим аналитическим выражением и погрешностью аппроксимации.

1.2.3 Предварительные сведения об исследуемом средстве измерений могут включать:

нормы на статический коэффициент преобразования;

структуру модели динамических свойств исследуемого средства измерений (структуру дифференциального уравнения или передаточной функции);

нормы на динамическую характеристику;

данные о возможных входных сигналах;

любые другие данные о средстве измерений и условиях его применения.

2. ВЫБОР МЕТОДОВ ОПРЕДЕЛЕНИЯ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ИССЛЕДУЕМЫХ СРЕДСТВ ИЗМЕРЕНИЙ

2.1. Выбор метода определения динамической характеристики должен быть обусловлен следующими качественными и количественными факторами.

2.1.1 Качественные факторы:

цель исследования;

требуемая динамическая характеристика;

требуемая форма представления динамической характеристики;

предварительная информация об исследуемом средстве измерений;

имеющиеся в наличии средства измерений (генераторы и источники испытательных сигналов, средства определения и регистрации испытательного сигнала и отклика исследуемого средства измерений);

имеющиеся в наличии средства вычислительной техники и их математическое обеспечение.

2.1.1.1 Цель исследования может состоять в следующем:

контроле нормированной динамической характеристики;

определении динамической характеристики типа средств измерений;

определении динамической характеристики экземпляра средств измерений данного типа.

2.1.1.2. Требуемая динамическая характеристика может быть полной или частной.

2.1.1.3. Испытательные сигналы могут быть характеристическими и нехарактеристическими; соответственно, методы определения динамических характеристик могут быть прямыми и косвенными.

2.1.1.4. В зависимости от объема предварительной информации косвенные методы определения полных динамических характеристик могут быть структурно-параметрическими и параметрическими.

2.1.2 Количественные факторы:

требуемая точность определения динамической характеристики;

требуемый диапазон частот, в котором должно быть исследовано средство измерений.

2.2. Выбор метода определения динамических характеристик при контроле обусловлен следующими требованиями.

2.2.1. Контроль следует производить преимущественно по частным динамическим характеристикам.

2.2.2. Контроль полных динамических характеристик допускается производить путем сравнения отклика исследуемого средства измерений на известный испытательный сигнал с допускаемыми откликами. Допускаемые отклики определяют (рассчитывают) с использованием нормированных динамических характеристик.

2.2.3. В качестве испытательных следует использовать преимущественно характеристические сигналы.

2.3. Выбор метода определения динамических характеристик типа средств измерений обусловлен следующими требованиями.

2.3.1. Типовую частную динамическую характеристику следует определять преимущественно с использованием характеристических сигналов.

В случае выполнения указанного требования частную динамическую характеристику следует находить по усредненному отклику группы исследуемых средств измерений.

2.3.2. Типовую полную динамическую характеристику следует определять по динамическим характеристикам экземпляров группы исследуемых средств измерений.

2.4. Выбор метода определения динамической характеристики экземпляра средств измерений данного типа обусловлен следующими требованиями.

2.4.1. Динамическую характеристику экземпляра средства измерений следует определять преимущественно путем уточнения известной типовой характеристики указанного средства измерений.

2.4.2. При необходимости представить полную динамическую характеристику исследуемого экземпляра средства измерений в аналитической форме следует использовать преимущественно параметрические методы.

2.5. При отсутствии требований к форме представления полной динамической характеристики исследуемого средства измерений следует определять указанную характеристику преимущественно прямым методом.

2.6. При необходимости представить полную динамическую характеристику исследуемого средства измерений в аналитической форме следует определять указанную характеристику косвенным методом с использованием преимущественно характеристических испытательных сигналов.

Соотношения для пересчета динамических характеристик из одной в другие для распространенных моделей динамических свойств средств измерений приведены в приложении 1.

2.7. В том случае, когда при выборе метода определения динамической характеристики исследуемого средства измерений требования, обусловленные указанными в п. 2.1.1 качественными факторами, противоречат друг другу, следует считать преимущественными те факторы, которые связаны с наличием технических средств (преимущественно средств измерений).

2.8. Метод определения динамической характеристики должен удовлетворять каждому из требований, обусловленных указанными в п. 2.1.2 количественными факторами.

2.9. В зависимости от требуемой точности определения динамической характеристики метод ее определения должен выбираться таким, чтобы погрешность определения динамической характеристики не превышала допускаемой погрешности.

2.9.1. При определении частной динамической характеристики V характеристика погрешности \Delta_{\epsilon,V} ее определения должна удовлетворять условию

\Delta_{\epsilon,V} \leq \Delta_{V\text{пр}}, (1)

где \Delta_{V\text{пр}} — предел допускаемой погрешности определения V.

2.9.2. При определении полной динамической характеристики в известной аналитической форме выбранный метод должен обеспечивать выполнение условий

\Delta_{v_{i}} \leq \Delta_{v_{i}\text{пр}} , i = 1, ..., m, (2)

где \Delta_{v_{i}} характеристика погрешности определения коэффициента vi, входящего в аналитическое выражение динамической характеристики;

\Delta_{v_{i}\text{пр}} — предел допускаемой погрешности определения коэффициента vi; m — число определяемых коэффициентов.

2.9.3. При определении полной динамической характеристики в табличной или графической форме характеристика погрешности определения \Delta_{\epsilon} должна удовлетворять условию

\Delta_{\epsilon} \leq \Delta_{\text{пр}}, (3)

где \Delta_{\text{пр}} — предел допускаемой погрешности определения динамической характеристики.

2.9.4. Погрешности методов следует вычислять в форме, аналогичной пределу допускаемой погрешности.

2.10. В зависимости от требуемого диапазона частот, в котором должно быть исследовано средство измерений, метод определения динамической характеристики следует выбирать в соответствии с требованиями разд. 3 настоящей рекомендации.

2.11. В зависимости от требуемых метрологических характеристик средств измерений, используемых для определения динамической характеристики исследуемого средства измерений, метод определения динамической характеристики следует выбирать в соответствии с требованиями разд. 4 настоящей рекомендации.

3. ВЫБОР ИСПЫТАТЕЛЬНЫХ СИГНАЛОВ

3.1. Испытательный сигнал для определения динамической характеристики исследуемого средства измерений следует выбирать в зависимости от определяемой характеристики, требуемого диапазона частот, в котором должно быть исследовано средство измерений, и требуемой точности определения динамической характеристики в указанном диапазоне, либо требуемой точности измерений, которые должны проводиться с помощью исследуемого средства измерений.

3.2. В зависимости от определяемой динамической характеристики испытательный сигнал должен выбираться преимущественно следующим:

для определения переходной характеристики и (или) ее параметров (функционалов) — ступенчатый испытательный сигнал, описываемый функцией, достаточно близкой по форме к единичной ступенчатой функции 1(t);

для определения импульсной характеристики и (или) ее параметров (функционалов) — импульсный испытательный сигнал, описываемый функцией, достаточно близкой по форме к дельта-функции \delta_{0}(t) или двоичный псевдослучайный испытательный сигнал;

для определения частотных характеристик и (или) их параметров (функционалов) — сигналы, достаточно близкие к гармоническим.

В случае, если указанные сигналы не могут быть воспроизведены (определены) с требуемой точностью, допускается выбор в качестве испытательных других сигналов, удовлетворяющих требованиям пп. 3.3 - 3.7.

Примечание. Испытательные сигналы аппроксимируют упрощенными функциями (модели сигналов), примеры которых приведены в приложении 3. Точность воспроизведения (определения) испытательного сигнала обусловлена отклонением модели сигнала по форме от характеристического сигнала (далее — неидеальностью испытательного сигнала) и отклонением испытательного сигнала от его модели.

3.3. В зависимости от требуемого диапазона частот, в котором должно быть исследовано средство измерений, амплитудный спектр X(f) испытательного сигнала должен удовлетворять условию

f01 >= fв, (4)

где f01 — частота, представляющая собой наименьший не равный нулю корень уравнения X(f) = \Delta_{\text{в}} (где f — частота); \Delta_{\text{в}} — верхняя граница погрешности определения амплитудного спектра; fв — верхняя граничная частота требуемого диапазона.

Аналитические выражения f01 и \Delta_{\text{в}} для распространенных испытательных сигналов приведены в приложении 2.

3.4. В зависимости от требуемого диапазона частот, в котором должно быть исследовано средство измерений, максимальная скорость x_{m}\prime изменений во времени (производная) импульсного испытательного сигнала xg(t) с длительностью \tau_{\text{и}} должна удовлетворять условию

x_{m}\prime \geq \Delta_{\text{в}}f_{\text{в}}/[\tau_{\text{и}}(sinu)/u], (5)

где u = \pi f_{\text{в}}\tau_{\text{и}} .

Условие (5) заведомо удовлетворяется, если

x_{m} \geq \pi f_{\text{в}}\Delta_{\text{в}}/(sinu/u), (6)

где xm — максимальное значение сигнала xg(t).

3.5. При определении динамической характеристики средства измерений прямым методом испытательный сигнал должен быть выбран таким, чтобы предварительная оценка погрешности определения динамической характеристики, складывающаяся из составляющих, обусловленных неидеальностью испытательных сигналов и ожидаемой неточностью измерений и вычислений, не превышала установленного предела:

для импульсной характеристики g(t)

\Delta_{g\text{ч}} \leq \Delta_{g\text{пр}}, (7)

где \Delta_{g\text{ч}} — числовая характеристика погрешности \Delta_{g}(t) определения импульсной характеристики; \Delta_{g\text{пр}} — предел допускаемой погрешности определения импульсной характеристики;

\Delta_{g}(t) = (y_{g}(t)/y_{\text{инт}}) - g(t) ; y_{g}(t) = K \int_{0}^{t} g(t - \tau)x_{g}(\tau)d\tau ;

y_{\text{инт}} = \int_{0}^{\infty} y_{g}(t)dt = K \int_{0}^{\tau_{\text{и}}} x_{g}(t)dt ;

K — статический коэффициент преобразования исследуемого средства измерений;

для переходной характеристики h(t)

\Delta_{h\text{ч}} \leq \Delta_{h\text{пр}}, (8)

где \Delta_{h\text{ч}} — числовая характеристика погрешности \Delta_{h}(t) определения переходной характеристики; \Delta_{h\text{пр}} — предел допускаемой погрешности определения переходной характеристики; \Delta_{h}(t) = (y_{h}(t)/y_{\text{уст}}) - h(t) ; y_{h}(t) = K \int_{0}^{t} g(t - \tau)x_{h}(\tau)d\tau ; yуст — установившееся значение выходного сигнала, yуст = Kxуст; xуст — установившееся значение ступенчатого испытательного сигнала;

для амплитудно-частотной характеристики A(f)

\Delta_{A\text{ч}} \leq \Delta_{A\text{пр}}, (9)

где \Delta_{A\text{ч}} — числовая характеристика погрешности \Delta_{A}(f) определения амплитудно-частотной характеристики; \Delta_{A\text{пр}} — предел допускаемой погрешности определения амплитудно-частотной характеристики; \Delta_{A}(f) = (y_{A,m}/(Kx_{A,m})) - A(f) ; y_{A,m} = ma_{t \in}x[_{0}y_{,1/Af}(]t,f) ; x_{A,m} = ma_{t}x_{\in}[_{0}x_{,1A/}/_{f}(]t,f) ; y_{A}(t,f) = K \int_{- \infty}^{\infty} g(t - \tau)x_{A}(\tau,f)d\tau ; xA(t, f) — периодический испытательный сигнал, близкий к гармоническому с периодом повторений, равным 1/f; x_{A}(t,f) = sin2\pi ft + \sum_{n = 2}^{N_{\Sigma}} \alpha_{n}sinn(2\pi ft + \varphi_{n}) ; n — номер гармонической составляющей сигнала; \alpha_{n} — коэффициент гармонической составляющей с номером n >= 2; N_{\Sigma} — общее число гармонических составляющих (включая основную); \varphi_{n} — фазовый сдвиг гармонической составляющей с номером n >= 2 относительно основной;

для фазочастотной характеристики Ф(f)

\Delta_{\Phi\text{ч}} \leq \Delta_{\Phi\text{пр}}, (10)

где \Delta_{\Phi\text{ч}} — числовая характеристика погрешности \Delta_{\Phi}(f) , определения фазочастотной характеристики; \Delta_{\Phi\text{пр}} — предел допускаемой погрешности определения фазочастотной характеристики; \Delta_{\Phi}(f) = \epsilon_{0} - \zeta_{0} - \Phi(f) ; \epsilon_{0} , \zeta_{0} — соответственно, корни уравнений

sin\epsilon + \sum_{n = 2}^{N_{\Sigma}} a_{n}sinn(\epsilon + \varphi_{n}) = 0; (11)

sin(\zeta + \Phi(z)) + \sum_{n = 2}^{N_{\Sigma}} a_{n}(\frac{A(nz)}{A(z)})sin(n(\zeta + \varphi_{n}) + \Phi(nz)) = 0, (12)

которым соответствуют положительные производные по времени сигналов xA(t, f) и yA(t, f); z = 2\pi fT ; T — коэффициент или среднее геометрическое двух коэффициентов передаточной функции исследуемого средства измерений.

Выражения z для наиболее распространенных моделей динамических свойств исследуемых средств измерений приведены в приложении 1.

Примечание. Допускается определять фазочастотную характеристику прямым методом по ненулевым синфазным значениям испытательного сигнала и отклика на него с использованием соотношений, отличных от приведенных.

Составляющие погрешности определения прямыми методами переходной, импульсной, амплитудно-частотной, фазочастотной характеристик, обусловленные неидеальностью испытательных сигналов, в зависимости от параметров указанных сигналов и коэффициентов передаточной функции исследуемого средства измерений для распространенных моделей динамических свойств средств измерений и испытательных сигналов приведены в приложениях 3 - 8.

3.6. В зависимости от требуемой точности определения косвенными методами динамических характеристик исследуемых средств измерений испытательные сигналы должны выбираться на основе соотношений для оценивания погрешностей, входящих в состав конкретных алгоритмов обработки данных. Примеры алгоритмов приведены в приложении 8.

3.7. В зависимости от требуемой точности измерений, которые должны проводиться с помощью исследуемого средства измерений, испытательный сигнал должен быть выбран таким, чтобы удовлетворялись следующие неравенства:

во временнуй области

\Delta_{\text{ч}} \geq K_{\text{з}}\Delta_{n,t}, (13)

где \Delta_{\text{ч}} — числовая характеристика временнуй зависимости \Delta(t) искажений испытательного сигнала x(t); \Delta_{n,t} — предел допускаемой погрешности измерений; \Delta(t) = (y(t)/K) - x(t) ; y(t) — выходной сигнал исследуемого средства измерений; Kз — коэффициент запаса;

в частотной области

\Delta_{f\text{ч}} \geq K_{\text{з}}\Delta_{n,f}, (14)

где \Delta_{f\text{ч}} — числовая характеристика частотной зависимости \Delta_{f}(f) искажений испытательного сигнала x(t); \Delta_{n,f} — предел допускаемой погрешности измерений; \Delta_{f}(f) = (Y(f)/K) - X(f) ; Y(f) — амплитудный спектр выходного сигнала исследуемого средства измерений.

Коэффициент запаса следует установить в диапазоне Kз = 3 — 10, имея в виду, что при исследованиях более точных средств измерений Kз должно быть большим.

При необходимости уточнения коэффициент запаса Kз для конкретного метода определения динамических характеристик средств измерений следует устанавливать путем последовательных приближений в соответствии с приложением 9.

4. ТРЕБОВАНИЯ К СРЕДСТВАМ ИЗМЕРЕНИЙ, ИСПОЛЬЗУЕМЫМ ПРИ ОПРЕДЕЛЕНИИ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СРЕДСТВ ИЗМЕРЕНИЙ

4.1. Генераторы испытательных сигналов должны обеспечивать требуемые согласно разд. 3 характеристики испытательных сигналов.

4.2. Погрешности средств измерений, используемых для воспроизведения или определения испытательного сигнала и для регистрации выходного сигнала исследуемого средства измерений, должны удовлетворять условию

\mu_{\text{ч}} \leq 0,2\Delta_{\text{ч}}, (15)

где \mu_{\text{ч}} — числовая характеристика погрешности \mu(t) ; \mu(t) = \Delta_{x}(t) + \Delta_{y}(t)/K ; \Delta_{x}(t) — погрешность генератора испытательного сигнала (исключая неидеальность испытательного сигнала) или регистрирующего прибора, используемого для определения испытательного сигнала; \Delta_{y}(t) — погрешность регистрирующего прибора, используемого для определения выходного сигнала исследуемого средства измерений.

4.3. Импульсная характеристика gр(t) регистрирующего прибора, используемого для определения испытательного сигнала или выходного сигнала исследуемого средства измерений, должна удовлетворять условию

E_{x_{\text{р}}} + (E_{y_{\text{р}}}/K_{\text{р}}^{2}) - \frac{2}{K_{\text{р}}} \int_{0}^{\tau_{\text{н}}} x_{\text{р}}(t)y(t)dt \leq \Delta_{\text{р}}^{2}\tau_{\text{н}}, (16)

где ; E_{y_{\text{р}}} = \int_{0}^{\tau_{\text{н}}} y_{\text{р}}^{2}(t)dt ; Kр — статический коэффициент преобразования регистрирующего прибора; \tau_{\text{н}} — длительность интервала времени наблюдений выходного сигнала.

При регистрации испытательного сигнала xр(t) = x(t); y_{\text{р}}(t) = K_{\text{р}} \int_{0}^{t} g_{\text{р}}(t - \tau)x(\tau)d\tau — результат регистрации x(t); \Delta_{\text{р}} = \Delta_{x\text{ч}} .

При регистрации выходного сигнала xр(t) = y(t); y_{\text{р}}(t) = K_{\text{р}} \int_{0}^{t} g_{\text{р}}(t - \tau)y(\tau)d\tau — результат регистрации y(t); \Delta_{\text{р}} = \Delta_{y\text{ч}} .

4.3.1. Известные номинальная импульсная характеристика регистрирующего прибора gр.ном(t) и предел \Delta_{g} допускаемого отклонения от указанной характеристики должны удовлетворять условиям:

при выборе прибора

\Delta_{g} \leq \frac{E_{x_{\text{р}}} + (E_{y_{\text{р}.\text{ном}}}/K^{2}) - I_{xy} - \Delta_{\text{р}}^{2}\tau_{\text{н}}}{\int_{0}^{\tau_{\text{н}}} [x_{\text{р}}(t) - (y_{\text{р}.\text{ном}}(t)/K_{\text{р}})]J(t)dt} ; (17)

при проверке после регистрации

\Delta_{g} \leq \frac{E_{x_{\text{р}}} + (E_{y_{\text{р}}}/K_{\text{р}}^{2}) - I_{xy} - \Delta_{\text{р}}^{2}\tau_{\text{н}}}{\int_{0}^{\tau_{\text{н}}} x_{\text{р}}(t)J(t)dt} , (18)

где E_{y_{\text{р}.\text{ном}}} = \int_{0}^{\tau_{\text{н}}} y_{\text{р}.\text{ном}}^{2}(t)dt ; y_{\text{р}.\text{ном}}(t) = K_{\text{р}} \int_{0}^{t} g_{\text{р}.\text{ном}}(t - \tau)x_{\text{р}}(\tau)d\tau ; I_{xy} = \frac{2}{K_{\text{р}}} \int_{0}^{\tau_{\text{н}}} x_{\text{р}}(t)y_{\text{р}.\text{ном}}(t)dt ; J(t) = 2 \int_{0}^{t} x_{\text{р}}(\tau)d\tau .

4.4. Частотные характеристики Aр(f), Фр(f) регистрирующего прибора, используемого для определения испытательного сигнала или выходного сигнала исследуемого средства измерений, должны удовлетворять условию

E_{x_{\text{р}}} + (E_{y_{\text{р}}}/K_{\text{р}}^{2}) - 4\pi \int_{f_{\text{н}}}^{f_{\text{в}}} A_{\text{р}}(f)X_{\text{р}}^{2}(f)cos\Phi_{\text{р}}(f)df \leq \Delta_{\text{р}}^{2}\tau_{\text{н}}, (19)

где fн — нижняя граница требуемого диапазона частот, в котором исследуется средство измерений, Xр(f) — амплитудный спектр сигнала x_{\text{р}}(\tau) .

4.4.1. Известные номинальные частотные характеристики Aр.ном(f), Фр.ном(f) регистрирующего прибора и пределы допускаемых отклонений \Delta_{A} , \Delta_{\Phi} от указанных характеристик должны удовлетворять условиям:

при выборе прибора

E_{x_{\text{р}}} + (E_{y_{\text{р}.\text{ном}}}/K^{2}) - I_{\text{с}} + \Delta_{A}(I_{y} - I) + \Delta_{\text{ф}}I_{\text{с}} \leq \Delta_{\text{р}}^{2}\tau_{\text{н}}; (20)

при проверке после регистрации

E_{x_{\text{р}}} + (E_{y_{\text{р}.\text{ном}}}/K_{\text{р}}^{2}) - I_{\text{с}} - \Delta_{A}I + \Delta_{\text{ф}}I_{s} \leq \Delta_{\text{р}}^{2}\tau_{\text{н}}, (21)

где I_{\text{с}} = \frac{2}{K_{\text{р}}} \int_{f_{\text{н}}}^{f_{\text{в}}} X_{\text{р}}(f)Y_{\text{р}.\text{ном}}(f)cos\Phi_{\text{р}.\text{ном}}(f)df ; Yр.ном(f) = KрAр.ном(f)Xр(f); I_{y} = \frac{2}{K_{\text{р}}} \int_{f_{\text{н}}}^{f_{\text{в}}} X_{\text{р}}(f)Y_{\text{р}.\text{ном}}(f)df ; I = 2 \int_{f_{\text{н}}}^{f_{\text{в}}} X_{\text{р}}^{2}(f)cos\Phi_{\text{р}.\text{ном}}(f)df ; I_{s} = \frac{2}{K_{\text{р}}} \int_{f_{\text{н}}}^{f_{\text{в}}} X_{\text{р}}(f)Y_{\text{р}.\text{ном}}(f)sin\Phi_{\text{р}.\text{ном}}(f)df .

5. АЛГОРИТМЫ ОБРАБОТКИ ДАННЫХ ПРИ ОПРЕДЕЛЕНИИ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК

5.1. Алгоритм обработки данных при определении динамической характеристики обусловлен определяемой характеристикой и испытательным сигналом. Исходными данными служат испытательный и выходной сигналы исследуемого средства измерений и предварительные сведения об исследуемом средстве измерений, представленные в соответствии с требованиями п. 1.2.

5.2. При определении прямым методом переходной и импульсной характеристик исследуемого средства измерений оценки указанных характеристик следует находить с помощью соотношений:

переходной характеристики

\tilde{h}(t) = y_{h}(t)/y_{\text{уст}}; (22)

импульсной характеристики

\tilde{g}(t) = y_{g}(t)/y_{\text{инт}}. (23)

5.3. При определении прямым методом амплитудно-частотной и фазочастотной характеристик исследуемого средства измерений оценки указанных характеристик при частоте f следует находить с помощью соотношений:

амплитудно-частотной характеристики

\tilde{A}(f) = y_{A,m}/(Kx_{A,m}); (24)

фазочастотной характеристики

\tilde{\Phi}(f) = t_{y}(f) - t_{x}(f), (25)

где ty(f), tx(f) — значения времени, представляющие собой наименьшие по модулю корни уравнений yA(ty, f) = 0 и xA(tx, f) = 0, удовлетворяющие условию ty(f) — tx(f) < 1/f.

Примечание. Допускается в случае определения фазочастотной характеристики прямым методом по ненулевым синфазным значениям испытательного сигнала и отклика на него использовать соотношения, отличные от приведенных.

5.4. В том случае, когда либо не задана требуемая точность определения динамической характеристики или требуемая точность измерений с помощью исследуемого средства измерений, либо отсутствует возможность воспроизвести выбранный в соответствии с требованиями разд. 3 испытательный сигнал, кроме указанных в пп. 4.2 - 4.4 операций, следует также оценивать погрешность определения найденной динамической характеристики.

5.4.1. Оценивание погрешностей определения динамических характеристик, найденных прямыми методами, должно производиться в соответствии с требованиями разд. 6 рекомендаций.

5.5. Обработка данных с целью определения динамических характеристик косвенным методом с использованием испытательных сигналов, достаточно близких к характеристическим, должна включать:

регистрацию отклика исследуемого средства измерений на испытательный сигнал и отсчитывание значений отклика;

пересчет отсчитанных значений в значения требуемой динамической характеристики (в аналитическое выражение требуемой характеристики);

оценивание погрешностей рассчитанных значений (погрешностей коэффициентов рассчитанного аналитического выражения).

Соотношения для пересчета для распространенных моделей динамических свойств средств измерений приведены в приложении 1.

5.6. Обработка данных с целью определения динамических характеристик с использованием испытательных сигналов, отличных от указанных в п. 5.5, должна включать:

расчет по совокупности испытательного и выходного сигналов требуемой динамической характеристики в форме таблицы значений или аналитического выражения;

оценивание погрешностей рассчитанных значений или коэффициентов аналитического выражения.

При определении динамической характеристики в аналитической форме по совокупности рассчитанных значений, кроме указанных выше операций, должны быть также произведены аппроксимация указанной совокупности аналитическим выражением и оценивание погрешностей найденных коэффициентов этого выражения.

5.7. Алгоритм обработки данных при определении динамической характеристики средства измерений косвенным методом должен удовлетворять следующим требованиям:

обеспечивать получение возможно более простых аналитических выражений для динамических характеристик;

предусматривать возможность упрощений вычислений и получаемого результата за счет использования предварительных сведений об исследуемом средстве измерений;

включать в себя методику оценивания погрешностей найденной динамической характеристики, учитывающую неидеальность испытательного сигнала, погрешности регистрации отклика (выходного сигнала) средства измерений и погрешности вычислений;

обеспечивать регулярную зависимость найденной характеристики от исходных данных (при уменьшении погрешностей определения испытательного сигнала и отклика на него исследуемого средства измерений погрешность оценки характеристики не должна возрастать).

Алгоритмы обработки данных при определении полных динамических характеристик косвенными методами приведены в приложении 8.

6. ОЦЕНИВАНИЕ ПОГРЕШНОСТЕЙ ОПРЕДЕЛЕНИЯ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК

6.1 При оценивании погрешностей определения динамических характеристик следует учитывать:

погрешность определения испытательного сигнала, в том числе погрешность, обусловленную неидеальностью испытательного сигнала, и погрешность генератора испытательного сигнала или

погрешность источника сигнала и регистрирующего прибора и погрешность отсчитывания значений зарегистрированного сигнала (при аналоговой регистрации);

погрешность, обусловленную взаимодействием генератора (источника) испытательного сигнала и исследуемого средства измерений;

погрешность определения выходного сигнала исследуемого средства измерений, в том числе:

погрешность регистрирующего прибора, предназначенного для регистрации выходного сигнала, и погрешность отсчитывания зарегистрированных значений (при аналоговой регистрации),

погрешность, обусловленную взаимодействием регистрирующего прибора и исследуемого средства измерений,

погрешность, обусловленную помехами, воздействующими на выход исследуемого средства измерений и на вход регистрирующего прибора;

погрешность обработки данных, в том числе:

погрешность, обусловленную неадекватностью принятой модели исследуемому средству измерений,

погрешность, вызванную несовершенством выбранного критерия сравнения экспериментальных и расчетных данных,

погрешность, обусловленную приближенным характером используемых формул,

погрешность округления чисел,

погрешность аппроксимации экспериментальных данных,

погрешности используемых средств обработки.

6.2. Погрешность определения переходной (импульсной) характеристики прямым методом следует оценивать путем суммирования составляющих на основе следующего выражения, представляющего структуру модели указанной погрешности:

\hat{\Delta}(t) = \hat{\Delta}_{\text{с}}(t) + \hat{\Delta}_{\text{сл}}(t) + (\Delta_{\text{отсч}}(t) + \Delta_{y}(t))/\hat{y}, (26)

где при определении h(t): \hat{\Delta}(t) = \Delta_{h}(t) ; \hat{y} = y_{\text{уст}} ; при определении g(t): \hat{\Delta}(t) = \Delta_{g}(t) ; \hat{y} = y_{\text{инт}} ; \hat{\Delta}_{\text{с}}(t) = (K/\hat{y}) \int_{0}^{t} g(t - \tau)\hat{\Delta}_{x\text{с}}(\tau)d\tau ; \hat{\Delta}_{\text{сл}}(t) = (K/\hat{y}) \int_{0}^{t} g(t - \tau)\hat{\Delta}_{x\text{сл}}(\tau)d\tau ; \hat{\Delta}_{x\text{с}}(t) — систематическая составляющая погрешности \hat{\Delta}_{x}(t) ; \hat{\Delta}_{x\text{сл}}(t) — случайная составляющая погрешности \hat{\Delta}_{x}(t) ; \hat{\Delta}_{x}(t) — погрешность определения испытательного сигнала; \Delta_{\text{отсч}}(t) — погрешность отсчитывания значений выходного сигнала (при аналоговой регистрации).

При известных характеристиках составляющих границы погрешности \hat{\Delta}(t) следует находить по ГОСТ 8.207.

Выражения и числовые характеристики систематической составляющей \hat{\Delta}_{\text{с}}(t) погрешности \hat{\Delta}(t) , обусловленной неидеальностью испытательного сигнала, для средств измерений с распространенными моделями динамических свойств и распространенных испытательных сигналов приведены в приложениях 3 - 5.

6.3. Погрешности определения частотных характеристик прямым методом следует оценивать путем суммирования составляющих на основе следующих выражений, отражающих структуру указанных погрешностей:

для амплитудно-частотной характеристики

\Delta_{A}(f) = \frac{1}{K} (K \frac{A(f)\hat{\Delta}_{xf}(f)}{y_{A,m}} + \frac{\Delta_{x_{A,m}}(f)}{x_{A,m}} + \frac{\Delta_{y_{A,m}}(f)}{y_{A,m}}) \frac{y_{A,m}}{x_{A,m}} , (27)

где \Delta_{x,f}(f) — частотная зависимость погрешности воспроизведения (определения) испытательного сигнала; \Delta_{x_{A,m}}(f) — погрешность определения xA, m при частоте f; \Delta_{y_{A,m}}(f) — погрешность определения yA, m при частоте f;

для фазочастотной характеристики

\Delta_{\text{ф}}(f) = \Delta_{\Phi\text{с}}(f) + \Delta_{\Phi\text{сл}}(f) + \Delta_{\Phi,x}(f) + \Delta_{\Phi,y}(f), (28)

где \Delta_{\Phi\text{с}}(f) — систематическая составляющая, обусловленная неидеальностью испытательного сигнала, при частоте f; \Delta_{\Phi\text{сл}}(f) — суммарная случайная составляющая, обусловленная случайными погрешностями определения испытательного и выходного сигналов, при частоте f; \Delta_{\Phi,x}(f) , \Delta_{\Phi,y}(f) — соответственно, погрешности, обусловленные неточностью фиксации нулевых уровней (при определении фазочастотной характеристики по ненулевым синфазным значениям — уровней, соответствующих указанным значениям) испытательного и выходного сигналов, при частоте f.

Числовые характеристики составляющих погрешностей, обусловленных неидеальностью испытательных сигналов, для распространенных моделей динамических свойств средств измерений и распространенных испытательных сигналов приведены в приложениях 6, 7.

6.4. Погрешности определения динамических характеристик косвенными методами следует оценивать с помощью операций, входящих в состав конкретных алгоритмов обработки данных.

Примеры алгоритмов обработки данных при определении полных динамических характеристик косвенными методами приведены в приложении 8.

6.5. Погрешность определения динамической характеристики косвенным методом с использованием испытательных сигналов, достаточно близких к характеристическим, следует оценивать путем суммирования составляющих погрешности, указанных в пп. 6.2, 6.3, и погрешности пересчета значений отклика исследуемого средства измерений в значения требуемой характеристики (в аналитическое выражение требуемой характеристики).

Примеры определения динамических характеристик средств измерений приведены в приложении 10.

Приложение 1 Справочное

СООТНОШЕНИЯ ДЛЯ ПЕРЕСЧЕТА ПОЛНЫХ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СРЕДСТВ ИЗМЕРЕНИЙ

Обозначения динамических характеристик и коэффициентов	Аналитические выражения полных динамических характеристик h(t), g(t), A(z), Ф(z) и формулы для пересчета их коэффициентов через коэффициенты W(p) для
	W(p) = \frac{1}{Tp + 1}	W(p) = \frac{T_{4}^{2}p^{2} + T_{3}p + 1}{(T_{1}p + 1)^{2}(T_{2}p + 1)} T1 > 0; T2, T3, T4 >= 0	W(p) = \frac{T_{3}^{2}p^{2} + T_{2}p + 1}{(T_{1}p + 1)^{3}} T1 > 0; T2, T3 >= 0	W(p) = \frac{T_{5}^{2}p^{2} + T_{4}p + 1}{(T_{1}p + 1)(T_{2}p + 1)(T_{3}p + 1)} T1 > T2 > 0; T2 > T3 >= 0; T4, T5 >= 0	W(p) = \frac{T_{5}^{2}p^{2} + T_{4}p + 1}{(T_{1}^{2}p^{2} + 2\gamma T_{2}p + 1)(T_{3}p + 1)} T1 > 0; 0<\gamma<1 ; T3, T4, T5 >= 0
h(t)	1 - e^{-}^{\frac{t}{T}}	1 - (D_{0} + D_{1}t)e^{-}^{\frac{t}{T_{1}}} - E_{0}e^{-}^{\frac{t}{T_{2}}}	1 - (1 + D_{1}t + D_{2}t^{2})e^{-}^{\frac{t}{T_{1}}}	1 - D_{0}e^{-}^{\frac{t}{T_{1}}} - E_{0}e^{-}^{\frac{t}{T_{2}}} - F_{0}e^{-}^{\frac{t}{T_{3}}}	\omega = \sqrt{1 - \gamma^{2}} /T_{1}
g(t)	A_{0}e^{-}^{\frac{t}{T}}	(A_{0} + A_{1}t)e^{-}^{\frac{t}{T_{1}}} + B_{0}e^{-}^{\frac{t}{T_{2}}}	(A_{0} + A_{1}t + A_{2}t^{2})e^{-}^{\frac{t}{T_{1}}}	A_{0}e^{-}^{\frac{t}{T_{1}}} + B_{0}e^{-}^{\frac{t}{T_{2}}} + \frac{F_{0}}{T_{3}} e^{-}^{\frac{t}{T_{3}}}	\omega = \sqrt{1 - \gamma^{2}} /T_{1}
D0	-	\frac{T_{1}^{2} - T_{4}^{2} + T_{2}(T_{3} - 2T_{1})}{(T_{1} - T_{2})^{2}}	-	\frac{T_{5}^{2} + T_{1}(T_{1} - T_{4})}{(T_{1} - T_{2})(T_{1} - T_{3})}	\frac{1}{T_{1}(T_{1}^{2} + T_{3}(T_{3} - 2\gamma T_{1})) \sqrt{1 - \gamma^{2}}}
D1	-	\frac{T_{4}^{2} + T_{1}(T_{1} - T_{3})}{T_{1}^{2}(T_{1} - T_{2})}	\frac{T_{1}^{2} - T_{3}^{2}}{T_{1}^{3}}	-
D2	-	-	\frac{T_{3}^{2} + T_{1}(T_{1} - T_{2})}{2T_{1}^{4}}	-	(T_{1}^{2} - T_{5}^{2} + T_{3}T_{4} - 2\gamma T_{1}T_{3})T_{1} \sqrt{1 - \gamma^{2}}
E0	-	\frac{T_{4}^{2} + T_{2}(T_{2} - T_{3})}{(T_{1} - T_{2})^{2}}	-	\frac{T_{5}^{2} + T_{2}(T_{2} - T_{4})}{(T_{1} - T_{2})(T_{3} - T_{2})}	\frac{T_{5}^{2} + T_{3}(T_{3} - T_{4})}{T_{1}^{2} + T_{3}(T_{3} - 2\gamma T_{1})}
F0	-	-	-	\frac{T_{5}^{2} + T_{3}(T_{3} - T_{4})}{(T_{1} - T_{3})(T_{2} - T_{3})}	-
A0	1/T	-D1 + D0/T1	-D1 + 1/T1	D0/T1	D0/T1
A1	-	D1/T1	-2D2 + D1/T1	-	\gamma D_{1} + D_{2} \sqrt{1 - \gamma^{2}}
A2	-	-	D2/T1	-	\gamma D_{2} - D_{1} \sqrt{1 - \gamma^{2}}
B0	-	E0/T2, T2 > 0	-	E0/T2	-
A(z)	\sqrt{(B_{1}^{2}z^{2} + (1 + B_{3}^{2}z^{2})^{2})/((1 + B_{2}^{2}z^{2})((1 - z^{2})^{2} + 4\Gamma^{2}z^{2}))}
Ф(z)	arctg(z(C1 + C2z2 — C3z4)/(1 + C4z2 + C5z4))
z	2\pi fT	2\pi fT_{1}	2\pi fT_{1}	2\pi f \sqrt{T_{1}T_{2}}	2\pi fT_{1}
Г	1	1	1	(T_{1} + T_{2})/(2 \sqrt{T_{1}T_{2}})	\gamma
B1	1	T3/T1	T2/T1	T_{4}/ \sqrt{T_{1}T_{2}}	T4/T1
B2	0	T2/T1	1	T_{3}/ \sqrt{T_{1}T_{2}}	T3/T1
B3	0	T4/T1	T3/T1	T_{5}/ \sqrt{T_{1}T_{2}}	T5/T1
C1	B1 - B2 - 2Г
C2	C_{3} - C_{1} - 2\Gamma(1 - B_{3}^{2} + B_{1}B_{2})
C3	B_{2}B_{3}^{2}
C4	B_{1}B_{2} - (B_{3}^{2} + 1) + 2\Gamma(B_{1} - B_{2})
C5	B_{3}^{2} - B_{1}B_{2} + 2\Gamma C_{3}

T1 — T5 — коэффициенты передаточной функции W(p), имеющие смысл постоянных времени;

\gamma — коэффициент затухания колебаний средств измерений, содержащих колебательное звено.

Приложение 2 Справочное

АНАЛИТИЧЕСКИЕ ВЫРАЖЕНИЯ f01 И \Delta_{\text{в}} ДЛЯ РАСПРОСТРАНЕННЫХ ИСПЫТАТЕЛЬНЫХ СИГНАЛОВ

Испытательный сигнал		Уравнение для определения f01	Выражение для \Delta_{\text{в}}
Наименование	Аналитическое выражение
Сигнал с линейным нарастанием до установившегося значения		\frac{x_{\text{уст}}}{2\pi f_{01}} \frac{sinu}{u} = \Delta_{\text{в}} , u = \pi f_{01\tau_{\text{ф}}}	\Delta_{\text{в}} = \delta_{\text{в}}x_{\text{уст}}\tau_{\text{ф}}
Сигнал с экспоненциальным нарастанием до установившегося значения	x_{h}(t) = x_{\text{уст}}(1 - e^{- \beta t}) , t >= 0	f_{01} = \frac{\beta}{2\pi \sqrt{2}} \sqrt{\sqrt{1 + \frac{4x_{\text{уст}}^{2}}{\beta^{2}\Delta_{\text{в}}^{2}}} - 1}	\Delta_{\text{в}} = \frac{x_{\text{уст}}}{\beta} \delta_{\text{в}}
Сигнал с косинусоидальным нарастанием до установившегося значения		\frac{x_{\text{уст}}}{2\pi f_{01}} \frac{cos\pi f_{01}\tau_{\text{ф}}}{1 - (2f_{01}\tau_{\text{ф}})^{2}} = \Delta_{\text{в}}	\Delta_{\text{в}} = \delta_{\text{в}}x_{\text{уст}}\tau_{\text{ф}}
Импульсный сигнал прямоугольной формы		x_{m}\tau_{\text{и}} \frac{sinu}{u} = \Delta_{\text{в}} , u = \pi f_{01}\tau_{\text{и}}	\Delta_{\text{в}} = \delta_{\text{в}}x_{m}\tau_{\text{и}}
Импульсный сигнал симметричной треугольной формы			\Delta_{\text{в}} = \delta_{\text{в}}x_{m}\tau_{\text{и}}/2
Импульсный сигнал синусоидальной формы		\frac{2}{\pi} x_{m}\tau_{\text{и}} \frac{cos\pi f_{01}\tau_{\text{и}}}{1 - (2f_{01}\tau_{\text{и}})^{2}} = \Delta_{\text{в}}	\Delta_{\text{в}} = 2\delta_{\text{в}}x_{m}\tau_{\text{и}}/\pi
Импульсный сигнал косинусоидальной формы			\Delta_{\text{в}} = \delta_{\text{в}}x_{m}\tau_{\text{и}}/2

xуст — установившееся значение испытательного сигнала xh(t) с ненулевым установившимся значением;

xm — максимальное значение импульсного сигнала xg(t);

\tau_{\text{ф}} — длительность фронта нарастания сигнала xh(t);

\tau_{\text{и}} — длительность сигнала xg(t);

\delta_{\text{в}} — безразмерная величина, изменяющаяся в диапазоне от 0 до 1.

Приложение 3 Справочное

ПОГРЕШНОСТИ ОПРЕДЕЛЕНИЯ ПРЯМЫМИ МЕТОДАМИ ПЕРЕХОДНОЙ И ИМПУЛЬСНОЙ ХАРАКТЕРИСТИК ИССЛЕДУЕМОГО СРЕДСТВА ИЗМЕРЕНИЙ, ОБУСЛОВЛЕННЫЕ НЕИДЕАЛЬНОСТЬЮ ИСПЫТАТЕЛЬНЫХ СИГНАЛОВ, В ЗАВИСИМОСТИ ОТ ПАРАМЕТРОВ УКАЗАННЫХ СИГНАЛОВ, ВИДА И КОЭФФИЦИЕНТОВ ПЕРЕДАТОЧНОЙ ФУНКЦИИ СРЕДСТВА ИЗМЕРЕНИЙ (СИГНАЛ С ЛИНЕЙНЫМ НАРАСТАНИЕМ ДО УСТАНОВИВШЕГОСЯ ЗНАЧЕНИЯ ВМЕСТО СИГНАЛА, ОПИСЫВАЕМОГО ФУНКЦИЕЙ 1(t), — ДЛЯ ПЕРЕХОДНОЙ ХАРАКТЕРИСТИКИ; ИМПУЛЬСНЫЙ СИГНАЛ ПРЯМОУГОЛЬНОЙ ФОРМЫ ВМЕСТО СИГНАЛА, ОПИСЫВАЕМОГО ФУНКЦИЕЙ \delta_{0}(t) , — ДЛЯ ИМПУЛЬСНОЙ ХАРАКТЕРИСТИКИ)

W(p)	\Delta_{hc}(t)	\Delta_{gc}(t)
\frac{1}{Tp + 1}	m_{t}ax\|\Delta_{hc}(t)\| \leq \frac{\tau_{\text{ф}}}{2T} ,\text{при} \frac{\tau_{\text{ф}}}{T} \leq 0,3
\frac{T_{2}p + 1}{(T_{1}p + 1)^{2}}
\frac{T_{3}p + 1}{(T_{1}p + 1)(T_{2}p + 1)}
\frac{T_{3}p + 1}{T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1}

ПРИЛОЖЕНИЕ 4 Справочное

ПОГРЕШНОСТИ ОПРЕДЕЛЕНИЯ ПРЯМЫМ МЕТОДОМ ПЕРЕХОДНОЙ ХАРАКТЕРИСТИКИ ИССЛЕДУЕМОГО СРЕДСТВА ИЗМЕРЕНИЙ, ОБУСЛОВЛЕННЫЕ НЕИДЕАЛЬНОСТЬЮ ИСПЫТАТЕЛЬНОГО СИГНАЛА, В ЗАВИСИМОСТИ ОТ ПАРАМЕТРА УКАЗАННОГО СИГНАЛА, ВИДА И КОЭФФИЦИЕНТОВ ПЕРЕДАТОЧНОЙ ФУНКЦИИ СРЕДСТВА ИЗМЕРЕНИЙ (СИГНАЛ С ЭКСПОНЕНЦИАЛЬНЫМ НАРАСТАНИЕМ x_{h}(t) = 1 - e^{- \beta t} ВМЕСТО СИГНАЛА, ОПИСЫВАЕМОГО ФУНКЦИЕЙ 1(t)

Погрешности \Delta_{hc}(t) определения переходной характеристики или их числовые характеристики представляются следующими выражениями в зависимости от передаточных функций W(p):

W_{1}(p) = \frac{1}{Tp + 1} ;

m_{t}ax|\Delta_{hc}(t)| = (\beta T)^{- (1 - (\beta T)^{- 1})^{- 1}};

при \beta T \geq 10 можно пользоваться выражением

max|\Delta_{hc}(t)| \leq (\beta T)^{- 1};

W_{2}(p) = \frac{T_{3}p + 1}{(T_{1}p + 1)^{2}} ;

\Delta_{hc}(\tau) = [\frac{\chi m_{3}}{1 - \chi} + \frac{1 - m_{3}}{(1 - \chi)^{2}} + \frac{\tau(1 - m_{3})}{1 - \chi}]e^{- \tau} + \frac{\chi m_{3} - 1}{(1 - \chi)^{2}} e^{- \chi\tau},

где \chi = \beta T_{1} , m_{3} = T_{3}T_{1}^{- 1} , \tau = tT_{1}^{- 1} .

W_{3}(p) = \frac{T_{3}p + 1}{(T_{1}p + 1)(T_{2}p + 1)} ;

\Delta_{hc}(\tau) = \frac{\chi m_{3} - 1}{(1 - \chi)(1 - \chi m_{2})} e^{- \chi\tau} + \frac{m_{3} - m_{2}}{(1 - m_{2})(1 - \chi m_{2})} e^{- m_{2}\tau} + \frac{1 - m_{3}}{(1 - \chi)(1 - m_{2})} e^{- \tau},

где m_{2} = T_{2}T_{1}^{- 1} .

W_{4}(p) = \frac{T_{3}p + 1}{T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1} ;

Для W2(p), W3(p), W4(p) отсутствуют аналитические выражения для m_{t}ax|\Delta_{hc}(t)| , поэтому даны графические зависимости приведенных максимального модуля погрешности \frac{|\Delta_{hc}|m}{h_{m}} , среднего квадратического значения \frac{\Delta_{hc}|_{\text{с}.\text{кв}}}{h_{m}} и среднего модуля \frac{|\Delta_{hc}|_{\text{с}p}}{h_{m}} от параметра сигнала \chi , полученные численным способом. Здесь {}_{h_{m} = m_{t}ax}{_{h(t)}} . Указанные характеристики рассчитаны на интервале времени с длительностью 5T1, разделенном на 100 частей.

Указанные зависимости числовых характеристик погрешности представлены: для W3(p) — на рис. 1 - 3, для W4(p) — на рис. 4 - 6. Зависимости для W2(p) отдельно не представлены, так как последняя является частным случаем W2(p) при m2 = 1.

Рис 1. Зависимость |\Delta_{hc}|_{m}/h_{m} от \chi для W3(p): 1 — 3 - m2 = 0,06; m3 = 30; 1; 0,5; 4 — m3 = m2 и m2 = 0,05; m3 = ₀; 5, 6, 8 — m2 = 0,3; m3 = 30; 1; ₀; 7 — m2 = 0,03; m3 = 0,5 и m2 = 0,7; m3 = 1; 9, 8, 10 — m2 = 0,7; m3 = 30; 0,5; 0; 11, 12 — m2 = 1; m3 = 30; 0

Рис. 2. Зависимость |\Delta_{hc}|_{cp}/h_{m} от \chi для W3(p): 1 — m2 = 0; 1; m2 = 0,03; 0,3; 1; 2 — m2 = 0,03; m3 = 3; 3 — m3 = 3; m2 = 0,3; 4 — m3 = 3; m2 = 1; 5 — m3 = 30; m1 = 0,03; 6 — m3 = 10; m2 = 1 и m3 = 30; m1 = 0,3; 7 — m1 = 30; m2 = 1

Рис. 3. Зависимость |\Delta_{hc}|_{c.\text{кв}}/h_{m} от \chi для W3(p): 1, 2 — m3 = 1; m2 = 0,03; 0,3; 3 — m2 = 1 и m3 = 0; m1 = 0,03; 0,3; 4 — m2 = 0; m2 = 1; 5, 6, 7 — m3 = 30; m2 = 0,03; 0,3; 1

Рис. 4. Зависимость |\Delta_{hc}|_{m}/h_{m} от \chi для W4(p): 1, 2 — m3 = 1; \gamma = 0,05 ; 0,9; 3, 4 — m3 = 0; \gamma = 0,05 ; 0,9; 5, 6 — m3 = 30; \gamma = 0,05 ; 0,9

Рис. 5. Зависимость |\Delta_{hc}|_{cp}/h_{m} от \chi для W4(p): 1 — m3 = 1; \gamma = 0,05 и m3 = 30; \gamma = 0,7 ; 2 — m3 = 0; 1; \gamma = 0,7 ; 3 — m3 = 0; \gamma = 0,05 ; 4 — m3 = 30; \gamma = 0,05

Рис. 6. Зависимость |\Delta_{hc}|_{c.\text{кв}}/h_{m} от \chi для W4(p): 1, 2 — m3 = 1; \gamma = 0,05 ; 0,03; 3, 4 — m3 = 0; \gamma = 0,05 ; 0,3; 5, 6 — m3 = 30; \gamma = 0,05 ; 0,3

Для получения оценок \frac{|\Delta_{hc}|_{m}}{h_{m}} при m3 <= 0,5 и \gamma = 0,05 ; 0,9 (для W4(p)) можно пользоваться приближенным соотношением

\frac{|\Delta_{hc}|_{m}(\chi,\gamma,m_{3})}{h_{m}} |_{m_{3} = 0,5} \leq \frac{2}{3} \frac{|\Delta_{hc}|_{m}(\chi,\gamma,0)}{h_{m}} + \frac{|\Delta_{hc}|_{m}(\chi,\gamma,1)}{3h_{m}} .

Значения hm приведены в табл. 1 и 2.

Таблица 1

Значения hm в зависимости от параметров W3(p) m2 и m3

m2	m3
	0	0,03	0,05	0,1	0,3	0,5	1	3	5	10	30
0,03	0,993	0,993	0,993	0,994	0,995	0,996	1,000	2,756	4,549	9,029	26,950
0,05	0,993	0,993	0,993	0,994	0,995	0,996	1,000	2,667	4,365	8,633	25,705
0,1	0,993	0,993	0,993	0,993	0,995	0,996	1,000	2,486	4,021	7,835	23,369
0,3	0,990	0,991	0,991	0,991	0,993	0,995	1,000	2,050	3,227	6,302	18,119
0,5	0,987	0,987	0,987	0,988	0,991	0,993	1,000	1,800	2,778	5,263	15,253
0,7	0,978	0,980	0,980	0,981	0,985	0,989	0,999	1,628	2,470	4,627	13,317
0,9	0,967	0,968	0,969	0,970	0,976	0,982	0,996	1,499	2,241	4,157	11,890
1,0	0,960	0,961	0,961	0,963	0,970	0,976	0,993	1,446	2,146	3,963	11,306

Таблица 2

Значения hm в зависимости от параметров W4(p) \gamma и m3

\gamma	m3
	0	0,03	0,05	0,1	0,3	0,5	1	3	5	10	30
0,03	1,910	1,910	1,911	1,915	1,951	2,019	2,298	3,963	5,811	10,538	29,594
0,05	1,854	1,855	1,856	1,859	1,893	1,958	2,226	3,839	5,632	10,218	28,723
0,1	1,729	1,729	1,730	1,733	1,763	1,821	2,064	3,558	5,226	9,495	26,718
0,3	1,372	1,372	1,373	1,374	1,391	1,426	1,592	2,728	4,024	7,345	20,754
0,5	1,163	1,163	1,163	1,164	1,172	1,191	1,298	2,192	3,243	5,944	16,851
0,7	1,046	1,046	1,046	1,046	1,049	1,056	1,112	1,822	2,701	4,967	14,123
0,9	0,988	0,988	0,989	0,990	0,994	0,998	1,011	1,554	2,305	4,252	12,117

Приложение 5 Справочное

ПОГРЕШНОСТИ ОПРЕДЕЛЕНИЯ ПРЯМЫМ МЕТОДОМ ПЕРЕХОДНОЙ ХАРАКТЕРИСТИКИ ИССЛЕДУЕМОГО СРЕДСТВА ИЗМЕРЕНИЙ, ОБУСЛОВЛЕННЫЕ НЕИДЕАЛЬНОСТЬЮ ИСПЫТАТЕЛЬНОГО СИГНАЛА, В ЗАВИСИМОСТИ ОТ ПАРАМЕТРА СИГНАЛА, ВИДА И КОЭФФИЦИЕНТОВ ПЕРЕДАТОЧНОЙ ФУНКЦИИ СРЕДСТВА ИЗМЕРЕНИЙ [СИГНАЛ С КОСИНУСОИДАЛЬНЫМ ФРОНТОМ ВМЕСТО СИГНАЛА, ОПИСЫВАЕМОГО ФУНКЦИЕЙ 1(t)]

В настоящем приложении приведены полученные расчетным путем зависимости следующих числовых характеристик погрешности \Delta_{hc}(t) :

\frac{|\Delta_{hc}|_{m}}{h_{m}} ,| \frac{\Delta_{hc}}{h} |_{m}, \frac{|\Delta_{hc}|_{cp}}{h_{m}} , \frac{|\Delta_{hc}|_{c.\text{кв}}}{h_{m}}

от параметра сигнала \chi_{1} = \tau_{\text{ф}}/\text{Т}_{1} .

Указанные характеристики рассчитаны на интервале времени с длительностью 5T1, разделенном на 100 частей. Зависимости характеристик погрешности представлены на рис. 1 - 4 для

W_{3}(p) = \frac{T_{3}p + 1}{(T_{1}p + 1)(T_{2}p + 1)} ,

а на рис. 5 - 9 для

W_{4}(p) = \frac{T_{3}p + 1}{T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1} .

Значения hm приведены в табл. 1 и 2 приложения 4.

Рис. 1. Зависимость |\Delta_{hc}|_{m}/h_{m} от \chi_{1} для W3(p): 1 — m3 = 1; m2 = 0,1; 2 — m3 = 1; m2 = 0,3 и m2 = 30; m2 = 1; 3 — m3 = 1; m2 = 0,7 и m3 = 0,5; m2 = 0,1; 4 — m3 = 30; m2 = 0,1; 5 — m3 = 30; m2 = 0,3; 6 — m3 = 30; m2 = 0,7; 7 — m3 = 0,5; m2 = 0,1; 8 — m3 = 0; m2 = 0,1 и m3 = 0,1; 9, 10, 11 — m3 = 0; m2 = 0,3; 0,7; 1

Рис. 2. Зависимость |\Delta_{hc}/h_{m}|_{m} от \chi_{1} для W3(p): 1 — m3 = 0; m2 <= 1; 2 — m3 = 0,03; m2 = 1; 3 — m3 = 30; m2 = 1 и m3 = 0,03; m2 = 0,05; 4 — m3 = 30; m2 = 0,05

Рис. 3. Зависимость |\Delta_{hc}|_{cp}/h_{m} от \chi_{1} для W3(p): 1 — 0 <= m3 <= 1; 0 < m2 <= 1; 2 — m3 = 3; m2 = 1; 3 — m3 = 3; m2 = 0,3 и m3 = 5; m2 = 1; 4 — m3 = 30; 0 < m2 <= 1

Рис. 4. Зависимость |\Delta_{hc}|_{c.\text{кв}}/h_{m} от \chi_{1} для W3(p): 1 — m3 = 1; m2 = 0,1; 2 — m3 = 1; m2 = 1 и m3 <= 0,1; m2 = 0,1; 3 — m3 = 0,5; m2 = 0,1; 4 — m3 <= 0,5; m2 = 1; 5 — m3 = 30; m2 = 0,1; 6 — m3 = 30; m2 = 1

Рис. 5. Зависимость |\Delta_{hc}|_{m}/h_{m} от \chi_{1} для W4(p): 1 — 10 — \gamma = 0,3 ; m3 = 30; 10; 5; 3; 1; 0,5; 0,3; 0,1; 0,05; 0

Рис. 6. Зависимость |\Delta_{hc}/h|_{m} от \chi_{1} для W4(p): 1 — 5 - \gamma = 0,03 ; m3 = 0,3; 0,1; 0,05; 0,03; 0

Рис. 7. Зависимость |\Delta_{hc}|_{cp}/h_{m} от \chi_{1} для W4(p): 1 — 6 - \gamma = 0,03 ; m3 = 30; 3; 1; 0,5; 0,1; 0; 7 — 12 — \gamma = 0,3 ; m3 = 30; 3; 1; 0,5; 0,1; 0

Рис. 8. Зависимость |\Delta_{hc}|_{cp}/h_{m} от \chi_{1} для W4(p): 1 — 10 — \gamma = 0,7 ; m3 = 30; 10; 5; 3; 1; 0,5; 0,3; 0,1; 0,05; 0

Рис. 9. Зависимость |\Delta_{hc}|_{c.\text{кв}}/h_{m} от \chi_{1} для W4(p): 1 — 6 - \gamma = 0,03 ; m3 = 30; 3; 1; 0,5; 0,1; 0; 7 — 12 — \gamma = 0,3 ; m3 = 30; 3; 1; 0,5; 0,1; 0

Приложение 6 Справочное

ПОГРЕШНОСТИ ОПРЕДЕЛЕНИЯ ПРЯМЫМ МЕТОДОМ АМПЛИТУДНО-ЧАСТОТНОЙ ХАРАКТЕРИСТИКИ ИССЛЕДУЕМОГО СРЕДСТВА ИЗМЕРЕНИЙ, ОБУСЛОВЛЕННЫЕ НЕИДЕАЛЬНОСТЬЮ ИСПЫТАТЕЛЬНОГО СИГНАЛА В ЗАВИСИМОСТИ ОТ ПАРАМЕТРОВ СИГНАЛА, ВИДА И КОЭФФИЦИЕНТОВ ПЕРЕДАТОЧНОЙ ФУНКЦИИ СРЕДСТВА ИЗМЕРЕНИЙ (СУММА ОСНОВНОЙ И ВТОРОЙ ГАРМОНИЧЕСКИХ СОСТАВЛЯЮЩИХ x_{A}(t) = sin2\pi ft + \alpha_{2}sin2(2\pi ft + \varphi_{2}) , 0 \leq \alpha_{2} \leq 0,2 , ВМЕСТО ГАРМОНИЧЕСКОГО СИГНАЛА)

Погрешность \Delta_{A}(f) по определению представляет собой разность между максимальным значением приведенного ко входу отклика исследуемого средства измерений на неидеальный испытательный сигнал, отнесенный к своему максимальному значению, и значением искомой амплитудно-частотной характеристики указанного средства при частоте f (см. п. 3.5).

Для испытательного сигнала, представляющего собой сумму гармонических функций с частотами f и 2f, соответственно, максимальные значения сигнала и отклика на него имеют вид:

В данном приложении представлены зависимости приведенных средней квадратической и максимальной погрешностей и максимальной относительной погрешности:

\frac{|\Delta_{A}|_{\text{с}.\text{кв}}}{A_{max}} = \delta_{11}(\alpha_{2},\varphi_{2},m_{2},m_{3}); \frac{|\Delta_{A}|_{max}}{A_{max}} = \delta_{21}(\alpha_{2},\varphi_{2},m_{2},m_{3});

| \frac{\Delta_{A}}{A} |_{max} = \delta_{31}(\alpha_{2},\varphi_{2},m_{2},m_{3})\text{для}W_{3}(p) = \frac{T_{3}p + 1}{(T_{1}p + 1)(T_{2}p + 1)} ,

m_{2} = \frac{T_{2}}{T_{1}} ,m_{3} = \frac{T_{3}}{T_{1}} ;

\frac{|\Delta_{A}|_{\text{с}.\text{кв}}}{A_{max}} = \delta_{21}(\alpha_{2},\varphi_{2},\gamma,m_{3}); \frac{|\Delta_{A}|_{max}}{A_{max}} = \delta_{22}(\alpha_{2},\varphi_{2},\gamma,m_{3});

| \frac{\Delta_{A}}{A} |_{max} = \delta_{32}(\alpha_{2},\varphi_{2},\gamma,m_{3})\text{для}W_{4}(p) = \frac{T_{3}p + 1}{T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1} ,

m_{3} = \frac{T_{3}}{T_{1}} .

Значения сдвига \varphi_{2} между гармоническими составляющими испытательного сигнала приняты равными: 0, - \pi/8 , - \pi/4 .

Значения A(z) и Ф(z) для указанных передаточных функций рассчитывались по соотношениям, приведенным в приложении 1. Указанные выше числовые характеристики погрешности определялись в диапазоне изменений безразмерной частоты z от 0,01 до 100. Для коэффициентов передаточных функций были приняты следующие значения: m2 = 0; 0,1; 0,3; m3 = 0; 0,1; 0,3; 1,0 для W3(p); \gamma = 0,1 ; 0,7; 0,9; 1,0; m3 = 0; 0,1; 0,3; 1,0 для W4(p).

Построенные по результатам расчетов графики зависимостей числовых характеристик погрешностей проанализированы и обобщены следующим образом. Близкие кривые, отличающиеся не более чем на 30%, заменены мажорантными; бульшая часть кривых аппроксимирована кусочно-линейными функциями, зависящими от параметра \alpha_{2} и состоящими из двух отрезков прямых.

Ниже приведены приближенные аналитические выражения указанных характеристик и таблицы значений входящих в них коэффициентов, а также графические зависимости.

1.

Значения \alpha^{*} , \chi_{1} , \chi_{2} , \chi_{3} приведены в табл. 1.

Таблица 1

\alpha^{*} , \chi_{j}	m3				m2	\varphi_{2}
	0	0,1	0,3	1,0
\alpha^{*}	0,075	0,075	0,075	0,075	0,1	0
\chi_{1}	0,047	0,055	0,077	0,149
\chi_{2}	0,064	0,072	0,097	0,202
\chi_{3}	-0,0013	-0,0013	-0,0015	-0,004
\alpha^{*}	0,075	0,075	0,075	0,075	0,3
\chi_{1}	0,041	0,041	0,055	0,092
\chi_{2}	0,054	0,058	0,072	0,123
\chi_{3}	-0,001	-0,0013	-0,0013	-0,0023
\alpha^{*}	0,070	0,070	0,075	0,075	0,9
\chi_{1}	0,032	0,032	0,037	0,057
\chi_{2}	0,047	0,078	0,054	0,076
\chi_{3}	-0,0011	-0,0032	-0,0013	-0,0014
\alpha^{*}	0,138	0,138	0,090	0,100	0,1	- \frac{\pi}{8}
\chi_{1}	0,068	0,076	0,100	0,215
\chi_{2}	0,052	0,055	0,076	0,169
\chi_{3}	0,0022	0,0026	0,0022	0,0046
\alpha^{*}	0,125	0,125	0,120	0,116	0,3
\chi_{1}	0,060	0,064	0,076	0,130
\chi_{2}	0,053	0,053	0,059	0,103
\chi_{3}	0,0009	0,0014	0,002	0,0031
\alpha^{*}	0,125	0,125	0,125	0,100	0,9
\chi_{1}	0,054	0,054	0,060	0,081
\chi_{2}	0,046	0,047	0,049	0,066
\chi_{3}	0,001	0,0009	0,0014	0,0015
\alpha^{*}	0,100	0,100	0,088	0,076	0,1	- \frac{\pi}{4}
\chi_{1}	0,060	0,064	0,075	0,184
\chi_{2}	0,046	0,044	0,055	0,128
\chi_{3}	0,0014	0,002	0,0018	0,0043
\alpha^{*}	0,088	0,089	0,075	0,075	0,3
\chi_{1}	0,057	0,060	0,067	0,113
\chi_{2}	0,044	0,045	0,048	0,082
\chi_{3}	0,0011	0,0013	0,0014	0,0023
\alpha^{*}	0,750	0,750	0,750	0,750	0,9
\chi_{1}	0,056	0,056	0,056	0,070
\chi_{2}	0,038	0,039	0,042	0,049
\chi_{3}	0,0014	0,013	0,0011	0,0016

2. \delta_{21\text{н}}(\alpha_{2},\varphi_{2}) \leq \delta_{21}(\alpha_{2};\varphi_{2};m_{2};m_{3}) \leq \delta_{21\text{в}}(\alpha_{2},\varphi_{2}) ; 0 \leq \alpha_{2} \leq 0,2 ; m2 = 0,1; 0,3; 0,9; m3 = 0; 0,1; 0,3; 1,0; \varphi_{2} = 0 ; - \pi/8 ; - \pi/4 .

Зависимости \delta_{21\text{н}}(\alpha_{2}\varphi_{2}) и \delta_{21\text{в}}(\alpha_{2}\varphi_{2}) представлены: для \varphi_{2} = 0 на рис. 1, для \varphi_{2} = - \pi/8 на рис. 2, для \varphi_{2} = - \pi/4 на рис. 3.

Рис. 1. Зависимость |\Delta_{A}|_{max}/A_{m} от \alpha_{2} для W3(p): 1 — \delta_{21\text{в}}(\alpha_{2},0) ; 2 — \delta_{21\text{н}}(\alpha_{2},0)

Рис. 2. Зависимость |\Delta_{A}|_{max}/A_{m} от \alpha_{2} для W3(p): 1 — \delta_{21\text{в}}(\alpha_{2}, - \pi/8) ; 2 — \delta_{21\text{н}}(\alpha_{2}, - \pi/8)

Рис. 3. Зависимость |\Delta_{A}|_{max}/A_{m} от \alpha_{2} для W3(p): 1 — \delta_{21\text{в}}(\alpha_{2}, - \pi/4) ; 2 — \delta_{21\text{н}}(\alpha_{2}, - \pi/4)

3. Зависимости |\Delta_{A}/A|_{max} = \delta_{31}(\alpha_{2},0,m_{2},m_{3}) для m2 = 0,1; 0,3; 0,9 при 0,1 <= m3 <= 1,0 представлены на рис. 4.

Рис. 4. Зависимость |\Delta_{A}/A|_{max} от \alpha_{2} для W3(p): 1 — \delta_{31}(\alpha_{2};0;0,1;0) ; 2 — \delta_{31}(\alpha_{2};0;0,3;0) ; 3 — \delta_{31}(\alpha_{2};0;0,9;0) ; 4 -

4. \delta_{31\text{н}}(\alpha_{2}, - \pi/8) \leq \delta_{31}(\alpha_{2}, - \pi/8,m_{2},m_{3}) \leq \delta_{31\text{в}}(\alpha_{2}, - \pi/8) ,

где

m2 = 0,1; 0,3; 0,9; m3 = 0; 0,1; 0,3; 1,0.

5.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 2.

Таблица 2

\alpha^{*} , \kappa_{j}	m2				m2
	0	0,1	0,3	1,0
\alpha^{*}	0,063	0,075	0,075		0,1
\kappa_{1}	1,23	0,947	0,893
\kappa_{2}	0,982	0,655	0,616
\kappa_{3}	0,016	0,022	0,021
\alpha^{*}	0,063	0,063	0,050		0,3
\kappa_{1}	1,23	1,17	1,02
\kappa_{2}	0,989	0,880	0,717
\kappa_{3}	0,015	0,018	0,015
\alpha^{*}	0,075		0,075	0,063	0,9
\kappa_{1}	1,24		1,07	0,992
\kappa_{2}	0,96		0,752	0,669
\kappa_{3}	0,021		0,024	0,020

6.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 3.

Таблица 3

\alpha^{*} , \kappa_{j}	m2				\gamma
	0	0,1	0,3	1,0
\alpha^{*}	0,100			0,100	0,1
\kappa_{1}	0,034			0,034
\kappa_{2}	0,051			0,048
\kappa_{3}	-0,0017			-0,0014
\alpha^{*}	0,075		0,070	0,075	0,7
\kappa_{1}	0,032		0,036	0,053
\kappa_{2}	0,052		0,057	0,074
\kappa_{3}	-0,0015		-0,0015	-0,0016
\alpha^{*}	0,070		0,063	0,070	0,9
\kappa_{1}	0,031		0,035	0,054
\kappa_{2}	0,048		0,052	0,075
\kappa_{3}	-0,0012		-0,0011	-0,0015
\alpha^{*}	0,075	0,075	0,090	0,100	1,0
\kappa_{1}	0,032	0,032	0,038	0,059
\kappa_{2}	0,046	0,047	0,052	0,073
\kappa_{3}	-0,0011	-0,0011	-0,0013	-0,0014

7.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 4.

Таблица 4

\alpha^{*} , \kappa_{j}	m3				\gamma
	0	0,1	0,3	1,0
1	2	3	4	5	6
\alpha^{*}	0,145			0,150	0,1
\kappa_{1}	0,076			0,072
\kappa_{2}	0,070			0,068
\kappa_{3}	0,0009			0,0006
\alpha^{*}	0,130	0,128	0,115	0,125	0,7
\kappa_{1}	0,063	0,064	0,069	0,082
\kappa_{2}	0,056	0,057	0,060	0,069
\kappa_{3}	0,0009	0,0009	0,001	0,0016
\alpha^{*}	0,113	0,125	0,115	0,105	0,9
\kappa_{1}	0,055	0,056	0,061	0,081
\kappa_{2}	0,049	0,049	0,052	0,067
\kappa_{3}	0,0007	0,0009	0,001	0,0015
\alpha^{*}	0,125	0,125	0,125	0,110	1,0
\kappa_{1}	0,052	0,054	0,058	0,076
\kappa_{2}	0,045	0,045	0,047	0,062
\kappa_{3}	0,0009	0,0011	0,0014	0,0015

8.

9.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 5.

Таблица 5

\alpha^{*} , \kappa_{j}	m3				\gamma
	0	0,1	0,3	1,0
\alpha^{*}	0,075		0,075	0,065	0,7
\kappa_{1}	0,068		0,071	0,081
\kappa_{2}	0,054		0,055	0,059
\kappa_{3}	0,0011		0,0012	0,0014
\alpha^{*}	0,075			0,063	0,9
\kappa_{1}	0,059			0,073
\kappa_{2}	0,041			0,054
\kappa_{3}	0,0014			0,0012
\alpha^{*}	0,075			0,050	1,0
\kappa_{1}	0,056			0,072
\kappa_{2}	0,041			0,051
\kappa_{3}	0,0011			0,0011

10.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 6.

Таблица 6

\kappa_{j}	m3
	0; 0,1; 0,3	1,0
\kappa_{1}	0,956	0,937
\kappa_{2}	0,790	0,739
\kappa_{3}	0,017	0,0204

11.

Зависимости \delta_{22\text{в}}(\alpha_{2},0) и \delta_{22\text{в}}(\alpha_{2},0) представлены на рис. 5.

Рис. 5. Зависимость |\Delta_{A}|_{max}/A_{m} от \alpha_{2} для W4(p): 1 — \delta_{22\text{н}}(\alpha_{2},0) ; 2 — \delta_{22\text{в}}(\alpha_{2},0)

12.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 7.

Таблица 7

\alpha^{*} , \kappa_{j}	m3		\gamma
	0; 0,1; 0,3	1,0
\alpha^{*}	0,075	0,075	0,1
\kappa_{1}	0,747	0,760
\kappa_{2}	0,632	0,668
\kappa_{3}	0,009	0,007
\alpha^{*}	0,125	0,088	0,7
\kappa_{1}	0,500	0,423
\kappa_{2}	0,453	0,373
\kappa_{3}	0,006	0,004
\alpha^{*}	0,125	0,075	0,9
\kappa_{1}	0,416	0,387
\kappa_{2}	0,353	0,327
\kappa_{3}	0,008	0,0045
\alpha^{*}	0,120	0,090	1,0
\kappa_{1}	0,396	0,356
\kappa_{2}	0,330	0,290
\kappa_{3}	0,008	0,006

13.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 8.

Таблица 8

\alpha^{*} , \kappa_{j}	m3		\gamma
	0; 0,1; 0,3	1,0
\alpha^{*}	0,070		0,1
\kappa_{1}	0,910
\kappa_{2}	0,735
\kappa_{3}	0,0012
\alpha^{*}	0,070	0,070	0,7
\kappa_{1}	0,521	0,432
\kappa_{2}	0,408	0,344
\kappa_{3}	0,008	0,006
\alpha^{*}	0,080	0,080	0,9
\kappa_{1}	0,403	0,334
\kappa_{2}	0,394	0,252
\kappa_{3}	0,007	0,007
\alpha^{*}	0,075	0,075	1,0
\kappa_{1}	0,373	0,297
\kappa_{2}	0,290	0,220
\kappa_{3}	0,006	0,006

14.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 9.

Таблица 9

\alpha^{*} , \kappa_{j}	m3
	0; 0,1; 0,3	1,0
\alpha^{*}	0,075	0,06
\kappa_{1}	3,80	4,96
\kappa_{2}	3,14	3,93
\kappa_{3}	0,05	0,06

15. Зависимости \delta_{32}(\alpha_{2};0;\gamma;m_{3}) для \gamma = 0,7 ; 0,9; 1,0 и m3 = 0; 0,1; 0,3; 1,0 представлены на рис. 6.

Рис. 6. Зависимость |\Delta_{A}/A|_{max} от \alpha_{2} для W4(p): 1 — \delta_{32}(\alpha_{2};0;0,7;0) ; 2 — \delta_{32}(\alpha_{2};0;0,9;0) ; 3 — \delta_{32}(\alpha_{2};0;1,0;0) ; 4 -

16.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 10.

Таблица 10

\alpha^{*} , \kappa_{j}	m3
	0; 0,1	0,3	1,0
\alpha^{*}	0,090	0,090	0,113
\kappa_{1}	2,89	3,00	3,65
\kappa_{2}	2,27	2,64	3,03
\kappa_{3}	0,055	0,032	0,07

17. \delta_{32\text{н}}(\alpha_{2}, - \pi/8) \leq \delta_{32}(\alpha_{2}, - \pi/8,\gamma,m_{3}) \leq \delta_{32\text{в}}(\alpha_{2}, - \pi/8) ;

где

18.

Значения \alpha^{*} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 11.

Таблица 11

\alpha^{*} , \kappa_{j}	m3
	0; 0,1	0,3	1,0
\alpha^{*}	0,030	0,070	0,070
\kappa_{1}	1,92	2,00	2,47
\kappa_{2}	1,97	2,02	2,73
\kappa_{3}	-0,0015	-0,0014	-0,018

19. Зависимости \delta_{32}(\alpha_{2}; - \pi/4;\gamma;m_{3}) для \gamma = 0,7 ; 0,9; 1,0 представлены на рис. 7.

Рис. 7. Зависимость |\Delta_{A}/A|_{max} от \alpha_{2} для W4(p): 1 - 2 — \delta_{32}(\alpha_{2}; - \pi/4;1,0;0,3) ; 3 — \delta_{32}(\alpha_{2}; - \pi/4;0,9;0,3) ; 4 — \delta_{32}(\alpha_{2}; - \pi/4;0,7;0,3) ; 5 -

Пересчет приведенных погрешностей в абсолютные для передаточных функций W3(p) не требуется, поскольку для этих функций Amax = 1. Пересчет для передаточных функций W4(p) может быть сделан на основе значений Amax, приведенных в табл. 12.

Таблица 12

\gamma	m3
	0	0,1	0,3	1,0
0,1	5,025	5,050	5,242	7,080
0,7	1,000	1,000	1,002	1,107
0,9	1,000	1,000	1,000	1,000
1,0	1,000	1,000	1,000	1,000

Приложение 7 Справочное

ПОГРЕШНОСТИ ОПРЕДЕЛЕНИЯ ПРЯМЫМ МЕТОДОМ ФАЗОЧАСТОТНОЙ ХАРАКТЕРИСТИКИ ИССЛЕДУЕМОГО СРЕДСТВА ИЗМЕРЕНИЙ, ОБУСЛОВЛЕННЫЕ НЕИДЕАЛЬНОСТЬЮ ИСПЫТАТЕЛЬНОГО СИГНАЛА, В ЗАВИСИМОСТИ ОТ ПАРАМЕТРОВ СИГНАЛА, ВИДА И КОЭФФИЦИЕНТОВ ПЕРЕДАТОЧНОЙ ФУНКЦИИ СРЕДСТВА ИЗМЕРЕНИЙ (СУММА ОСНОВНОЙ И ВТОРОЙ ГАРМОНИЧЕСКИХ СОСТАВЛЯЮЩИХ x_{A}(t) = sin2\pi ft + \alpha_{2}sin2(2\pi ft + \varphi_{2}) , 0<\alpha_{2} \leq 0,2 , ВМЕСТО ГАРМОНИЧЕСКОГО СИГНАЛА)

Погрешность \Delta_{\text{Ф}}(f) в соответствии с п. 3.5 представляет собой разность между корнями \epsilon_{0} и \zeta_{0} уравнений (11) и (12) (которым соответствуют положительные производные по времени испытательного и выходного сигналов исследуемого средства измерений) и искомой фазочастотной характеристикой при частоте f.

Для испытательного сигнала, представляющего собой сумму гармонических функций с частотами f и 2f, соответственно, уравнения (11) и (12) имеют вид: sin\epsilon + \alpha_{2}sin2(\epsilon + \varphi_{2}) = 0 , 0 \leq \epsilon<2\pi ; sin(\zeta + \text{Ф}(z)) + \alpha_{2}(A(2z)/A(z))sin(2\zeta + \varphi_{2}) + \text{Ф}(2z)) = 0 , 0 \leq \zeta<2\pi .

В данном приложении представлены полученные расчетным путем зависимости максимальной погрешности:

|\Delta_{\text{Ф}}|_{max} = \Delta_{1}(\alpha_{2},\varphi_{2},m_{2},m_{3}) для W_{3}(p) = \frac{T_{3}p + 1}{(T_{1}p + 1)(T_{2}p + 1)} ,

и m_{2} = \frac{T_{2}}{T_{1}} , m_{3} = \frac{T_{3}}{T_{1}}

|\Delta_{\text{Ф}}|_{max} = \Delta_{2}(\alpha_{2},\varphi_{2},\gamma,m_{3}) для W_{4}(p) = \frac{T_{3}p + 1}{(T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1)} ,

m_{3} = \frac{T_{3}}{T_{1}} .

Значения сдвига \varphi_{2} между гармоническими составляющими сигнала приняты равными: 0, - \pi/8 , - \pi/4 .

Указанная выше характеристика погрешности определялась в диапазоне изменений безразмерной частоты z от 0,01 до 100:

Частотные характеристики A(z) и Ф(z), соответствующие передаточным функциям W3(p) и W4(p), рассчитывались по соотношениям, приведенным в справочном приложении 2.

Для коэффициентов передаточных функций были приняты значения: m2 = 0,1; 0,3; 0,7; m3 = 0; 0,1; 0,3; 0,7; 1,0 для W3(p) и \gamma = 0,1 ; 0,3; 0,5; 0,7; 1,0; m3 = 0; 0,1; 0,3; 0,7; 1,0 для W4(p).

Бульшая часть построенных по результатам расчетов зависимостей максимальной погрешности определения фазочастотной характеристики была аппроксимирована кусочно-линейными функциями, зависящими от параметра \alpha_{2} и состоящими из двух отрезков прямых.

Ниже приведены приближенные аналитические выражения и таблицы значений входящих в них коэффициентов, а также графики указанных выше зависимостей.

1.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 1.

Таблица 1

\kappa_{j}	m3
	0	0,1	0,3	0,7; 1,0
\kappa_{1}	0,48	0,58	0,56	0,48
\kappa_{2}	0,4275	0,495	0,4775	0,4275
\kappa_{3}	0,005225	0,0085	0,00825	0,00525

2.

3.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 2.

Таблица 2

\kappa_{j}	m3
	0	0,1	0,7; 1,0
\kappa_{1}	0,45	0,5	0,495
\kappa_{2}	0,4	0,485	0,466
\kappa_{3}	0,005	0,0015	0,003

4.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 3.

Таблица 3

\kappa_{j}	m3
	0	0,1; 0,3	0,7; 1,0
\kappa_{1}	0,43	0,517	0,495
\kappa_{2}	0,375	0,483	0,475
\kappa_{3}	0,0055	0,0034	0,002

5.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 4.

Таблица 4

\kappa_{j}	m3
	0	0,1	0,3; 0,7; 1,0
\kappa_{1}	0,965	0,934	0,8945
\kappa_{2}	0,711	0,6904	0,645
\kappa_{3}	0,023	0,022	0,02201

6.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 5.

Таблица 5

\kappa_{j}	m3
	0	0,1	0,3; 0,7; 1,0
\kappa_{1}	0,9445	1,0334	0,9445
\kappa_{2}	0,7545	0,7364	0,7545
\kappa_{3}	0,0171	0,0268	0,0171

7.

8.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 6.

Таблица 6

\kappa_{j}	m3
	0	0,3; 0,7; 1,0
\kappa_{1}	0,9612	1,0
\kappa_{2}	0,67	0,714
\kappa_{3}	0,0235	0,0258

9.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 7.

Таблица 7

\kappa_{j}	m3			m2
	0	0,1	0,3; 0,7; 1,0
\kappa_{1}	1,26	0,93	0,87	0,1
\kappa_{2}	1,09	0,795	0,79
\kappa_{3}	0,017	0,0135	0,008
\kappa_{1}	1,25	1,045	0,96	0,3
\kappa_{2}	1,13	0,895	0,87
\kappa_{3}	0,012	0,015	0,009
\kappa_{1}	1,27	1,17	0,99	0,7
\kappa_{2}	1,10	1,03	0,88
\kappa_{3}	0,017	0,014	0,011

10. Зависимость |\Delta_{\text{Ф}}|_{max}(\alpha_{2}) = \Delta_{2}(\alpha_{2};0;0,1;m_{3}) для 0<\alpha_{2} \leq 0,1 , m3 = 0; 0,1 представлена на рис. 1.

Рис. 1. Зависимость |\Delta_{\text{Ф}}|_{max} от \alpha_{2} для W4(p)

11. Зависимости |\Delta_{\text{Ф}}|_{max}(\alpha_{2}) = \Delta_{2}(\alpha_{2};0;0,3;m_{3}) для 0<\alpha_{2} \leq 0,2 , m3 = 0; ₀,1; 0,3 представлены на рис. 2.

Рис. 2. Зависимость |\Delta_{\text{Ф}}|_{max} от \alpha_{2} для W4(p): 1 — m3 = 0; 0,1; 2 — m2 = 0,3

12. Зависимости |\Delta_{\text{Ф}}|_{max}(\alpha_{2}) = \Delta_{2}(\alpha_{2};0;0,5;m_{3}) для 0 <= m3 <= 1,0 даются следующими соотношениями:

12.1.

12.2.

0,1 <= m3 <= 1,0

13. \Delta_{2}(\alpha_{2};0;0,7;0) = 0,27\alpha_{2} .

14. 0,49\alpha_{2} \leq \Delta_{2}(\alpha_{2};0;0,7;m_{3}) \leq 0,5\alpha_{2} .

0,1 <= m3 <= 1,0

15. \Delta_{2}(\alpha_{2};0;1,0;0) = 0,315\alpha_{2} .

16. \Delta_{2}(\alpha_{2};0;1,0;m_{3}) = 0,5\alpha_{2} .

0,1 <= m3 <= 1,0

17.

18.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 8.

Таблица 8

\kappa_{j}	m3
	0	0,1; 0,3
\kappa_{1}	0,875	1,1125
\kappa_{2}	0,7	0,91
\kappa_{3}	0,014	0,0164

19.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 9.

Таблица 9

\kappa_{j}	m3
	0	0,1; 0,3	0,7; 1,0
\kappa_{1}	0,907	1,114	1,03
\kappa_{2}	0,72	0,917	0,806
\kappa_{3}	0,01404	0,0152	0,0168

20.

Значения \alpha_{\text{п}} , \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 10.

Таблица 10

\alpha_{\text{п}} , \kappa_{j}	m3
	0	0,1	0,3	0,7; 1,0
\alpha_{\text{п}}	0,09	0,1	0,1	0,09
\kappa_{1}	0,917	1,1	1,056	1,0
\kappa_{2}	0,7043	0,845	0,814	0,7455
\kappa_{3}	0,01915	0,0255	0,0242	0,023

21.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 11.

Таблица 11

\kappa_{j}	m3
	0	0,1	0,3	0,7; 1,0
\kappa_{1}	0,947	1,12	1,067	1,02
\kappa_{2}	0,7358	0,888	0,8238	0,764
\kappa_{3}	0,016	0,0174	0,01824	0,0192

22. Зависимость |\Delta_{\text{Ф}}|_{max}(\alpha_{2}) = \Delta_{2}(\alpha_{2}; - \pi/4;0,1;m_{3}) для m3 = 0; 0,1 представлена на рис. 3.

Рис. 3. Зависимость |\Delta_{\text{Ф}}|_{max} от \alpha_{2} для W4(p): 1 — m2 = 0; 2 — m3 = 0,1

23.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 12.

Таблица 12

\kappa_{j}	m3
	0	0,1	0,3
\kappa_{1}	1,234	1,313	1,234
\kappa_{2}	1,112	1,143	1,112
\kappa_{3}	0,0147	0,0204	0,0147

24.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 13.

Таблица 13

\kappa_{j}	m3
	0	0,1	0,3	0,7	1,0
\kappa_{1}	1,24	1,31	1,21	1,04	0,991
\kappa_{2}	1,107	1,1434	1,0712	0,929	0,8855
\kappa_{3}	0,01467	0,01834	0,01528	0,01223	0,01162

25.

Значения \kappa_{1} , \kappa_{2} , \kappa_{3} приведены в табл. 14.

Таблица 14

\kappa_{j}	m3				m2
	0	0,1	0,3	0,7; 1,0
\kappa_{1}	1,250	1,31	1,187	1,015	0,7
\kappa_{2}	1,15	1,14	1,023	0,87
\kappa_{3}	0,01	0,017	0,0164	0,0145
\kappa_{1}	1,245	1,28	1,145	0,98	0,9
\kappa_{2}	1,135	1,12	0,975	0,91
\kappa_{3}	0,011	0,016	0,017	0,007

Приложение 8 Справочное

ПРИМЕРЫ АЛГОРИТМОВ ОБРАБОТКИ ДАННЫХ ПРИ ОПРЕДЕЛЕНИИ ПОЛНЫХ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК КОСВЕННЫМИ МЕТОДАМИ

1. Метод дискретизации уравнения свертки

1.1 Исходные данные.

1.1.1. N* значений испытательного сигнала, взятых с периодом дискретизации t_{1}^{*} = \tau_{\text{н}}/N^{*} , x_{\nu}^{*} = x(t_{\nu}^{*}) , t_{\nu}^{*} = \nu t_{1}^{*} , \nu = 1,2...,N^{*} - 1 ; x_{N^{*}}^{*} = x(t_{N^{*}}^{*})/2 .

1.1.2. Значения выходного сигнала y_{\nu}^{*} = y(t_{\nu}^{*}) , \nu = 1,2,...,N^{*} .

1.1.3. Верхние границы погрешностей определения значений испытательного сигнала d_{\nu}^{*} = \Delta_{\text{в}x_{\nu}}^{*} , \nu = 1,2,...,N^{*} .

1.1.4. Верхние границы погрешностей определения значений выходного сигнала \lambda_{\nu}^{*} = \Delta_{\text{в}y_{\nu}}^{*} , \nu = 1,2,...,N^{*} .

1.2 Последовательность вычислений.

1.2.1. Определяют момент времени t_{\nu}^{~} , для которого выполняется условие (\tau_{\text{н}}/t_{\nu}^{*})(d_{\nu}^{*}/x_{\nu}^{*})<1 , \nu = 1,2,... ;

вычисляют N = \tau_{\text{н}}/t_{\nu}^{*} ;

принимают новый период дискретизации: t_{1} = t_{\nu}^{*} ;

выбирают из указанных в п. 1.1 совокупностей x_{\nu}^{*} , y_{\nu}^{*} , d_{\nu}^{*} , \lambda_{\nu}^{*} ;

значения, соответствующие моментам времени ti = it1: [xi = x(ti), i = 1, ..., N — 1; xN = x(tN)/₂; yi = y(ti), i = 1, ..., N; d_{i} = \Delta_{\text{вх}_{i}} , i = 1, ..., N; \lambda_{i} = \Delta_{\text{вх}_{i}} , i = 1, ..., N.

1.2.2 Составляют матрицу значений испытательного сигнала

и вектор значений выходного сигнала

вычисляют обратную матрицу

где z1 = 1/x1; z_{j} = - \frac{1}{x_{1}} \sum_{i = 1}^{j - 1} (- 1)^{j - i + 1} \frac{x_{i + 1}^{j - i}}{x_{1}^{j - i}} , j = 2, ..., N;

вычисляют вектор оценок значений импульсной характеристики

где gi = g(ti), i = 0, 1, ..., N — 1.

1.2.3 Составляют матрицу абсолютных значений

и вектор погрешностей значений испытательного сигнала

вычисляют вектор

C = |X-1|Dx;

составляют матрицу

где ci — элементы вектора C;

вычисляют обратную матрицу

где w1 = 1/(1 — c1);

w_{j} = \frac{1}{c_{1} - 1} \sum_{i = 1}^{j - 1} (- 1)^{j - i} \frac{c_{j + 1}^{j - i}}{(1 - c_{1})^{j - i}} , j = 2, ..., n;

составляют вектор абсолютных значений

и вектор погрешностей значений выходного сигнала

вычисляют вектор оценок погрешностей определения значений импульсной характеристики RG = R-1(C·[G] + X-1Dy).

2. Многократное интегрирование входного и выходного сигналов

2.1 Исходные данные.

2.1.1 Те же, что в п. 1.1.

2.1.2. Структура дифференциального уравнения исследуемого средства измерений y(t) + \sum_{k = 1}^{m} a_{k}y^{(k)}(t) = \sum_{j = 1}^{l} b_{j}x^{(j)}(t) + Kx(t) , l < m.

2.1.3 Статический коэффициент преобразования средства измерений K.

2.2 Последовательность вычислений.

2.2.1 Вычисляют интегралы от входного и выходного сигналов на промежутках

[0,\tau_{r}] , r = 1, ..., m + l; ;

Для приближенного вычисления интеграла Im-j(s·2m-j) кратности m — j в точке {}^{\vee}t = s \cdot 2^{m - j} безразмерного времени {}^{\vee}t = t/t_{1} можно использовать рекуррентную формулу Симпсона:

где j = m — 1, ..., 1,0;

I₀(s) обозначает подынтегральную функцию в точке {}^{\vee}t = s . Для вычисления интеграла заданной кратности j можно использовать рекуррентную формулу

I_{j}(2s) = I_{j}(2(s - 1)) + \frac{2^{j - 1}}{3} (f_{2(s - 1)} + 4f_{2s - 1} + f_{2s}),

где fs — значение подынтегральной функции в точке {}^{\vee}t = s . Верхнюю границу промежутка интегрирования \tau_{m + i} следует выбирать равной времени наблюдений \tau_{\text{н}} . Верхние границы промежутков \tau_{r} , r = 1, ..., m + l — 1, следует выбирать приближенно равномерно в интервале [0, \tau_{m + 1} ], при этом должны быть выполнены условия \tau_{r - 1}<\tau_{r} . В случае использования формулы Симпсона количество значений x_{v}^{*}(y_{v}^{*}) на каждом промежутке интегрирования должно быть кратно значению 2m+1.

2.2.2 Составляют систему алгебраических уравнений

\psi_{m + 1}(r) + \sum_{k = 1}^{m} a_{k}\psi_{m - k + 1}(r) = K\text{Ф}_{n + 1}(r) + \sum_{j = 1}^{l} b_{j}\text{Ф}_{m - j + 1}(r) , r = 1, ..., m + l, или RиA = F,

где матрица системы

вектор свободных членов

вектор искомых коэффициентов

2.2.3 Вычисляют определитель матрицы Rи.

Если detR_{\text{и}} \neq 0 , то составляют обратную матрицу

где \overline{R_{ij}} — алгебраическое дополнение ij-го элемента матрицы Rи.

Если detRи = ₀, то изменяют один из промежутков интегрирования, вычисляют интегралы Фm+1-j(r), \psi_{m + 1 - k}(r) на новом промежутке и вычисления по пп. 2.2.2; 2.2.3 повторяют.

2.2.4. Составляют матрицу |R_{\text{и}}^{- 1}| абсолютных значений элементов R_{\text{и}}^{- 1} .

2.2.5 Составляют матрицу погрешностей элементов Rи

и вектор погрешностей элементов F

где \tau_{j}^{\vee} = \tau_{j}/t_{1} ; ; .

2.2.6. Вычисляют матрицу критической обусловленности системы алгебраических уравнений C_{\text{и}} = |R_{\text{и}}^{- 1}|D_{R_{\text{и}}} .

Вычисляют фактор критической обусловленности \text{Ф}_{\text{кр}} = \sum_{i = 1}^{m + l} c_{li} , где cli — диагональный элемент матрицы Cи.

2.2.7. Если Фкр < 1, то вычисляют вектор \tilde{A} оценок элементов A: \tilde{A} = R_{\text{и}}^{- 1}F .

Если Фкр >= 1, то увеличивают период дискретизации при разбиении промежутков интегрирования, и вычисления по пп. 2.2.2 - 2.2.7 повторяют.

2.2.8. Вычисляют вектор погрешностей D_{A} = (I - C)^{- 1}(C_{\text{и}}\tilde{A} + |R_{\text{и}}^{- 1}|D_{F}) , где I — единичная матрица.

3. Адаптивный метод определения полных динамических характеристик.

3.1 Исходные данные.

3.1.1. Те же, что в п. 1.1, кроме x_{N*}^{*} , равного x_{N*}^{*} = x(t_{v}^{*}) . Испытательный сигнал может быть задан также аналитическим выражением.

3.1.2 Совокупность типовых испытательных сигналов, представленных на рисунке:

	x0(t) = xуст1(t);
	x_{1}(t) = x_{m}e^{- \beta t}1(t) ;
	x_{2}(t) = x_{m} \frac{t}{\tau_{\text{и}}} [1(t) - 1(t - \tau_{\text{и}})] ;
	x_{3}(t) = x_{q}sin \frac{\pi t}{q\tau_{\text{и}}} [1(t) - 1(t - \tau_{\text{и}})] ;
	x_{4}(t) = x_{q}(1 - cos \frac{\pi}{q\tau_{\text{и}}} t)[1(t) - 1(t - \tau_{\text{и}})] ,

3.1.4 Иерархическая система моделей динамических свойств средств измерений с передаточными функциями Wk(p):

\frac{1}{(T_{1}p + 1)(T_{2}p + 1)} ; \frac{1}{(Tp + 1)^{3}} ; \frac{1}{(Tp + 1)^{2}} ; \frac{1}{Tp + 1} ;

\frac{T_{2}p + 1}{(T_{1}p + 1)^{2}} ; \frac{1}{(T_{1}p + 1)^{2}(T_{2}p + 1)} ; \frac{1}{T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1} ;

\frac{1}{(T_{1}p + 1)(T_{2}p + 1)(T_{3}p + 1)} ; \frac{(T_{2}p + 1)^{2}}{(T_{1}p + 1)^{3}} ; \frac{T_{2}p + 1}{(T_{1}p + 1)^{3}} ;

\frac{T_{3}p + 1}{T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1} ; \frac{T_{3}p + 1}{(T_{1}p + 1)(T_{2}p + 1)} ; \frac{T_{4}p + 1}{(T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1)(T_{3}p + 1)} ;

\frac{T_{3}^{2}p^{2} + T_{2}p + 1}{(T_{1}p + 1)^{3}} ; \frac{(T_{3}p + 1)^{2}}{(T_{1}p + 1)^{2}(T_{2}p + 1)} ; \frac{T_{3}p + 1}{(T_{1}p + 1)^{2}(T_{2}p + 1)} ;

\frac{T_{4}p + 1}{(T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1)(T_{3}p + 1)} ; \frac{T_{4}p + 1}{(T_{1}p + 1)(T_{2}p + 1)(T_{3}p + 1)} ;

\frac{(T_{4}p + 1)^{2}}{(T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1)(T_{3}p + 1)} ; \frac{(T_{4}p + 1)^{2}}{(T_{1}p + 1)(T_{2}p + 1)(T_{3}p + 1)} ;

\frac{T_{5}^{2}p^{2} + T_{4}p + 1}{(T_{1}p + 1)(T_{2}p + 1)(T_{3}p + 1)} ; \frac{T_{4}^{2}p^{2} + T_{3}p + 1}{(T_{1}p + 1)^{2}(T_{2}p + 1)} ;

\frac{T_{5}^{2}p^{2} + T_{4}p + 1}{(T_{1}^{2}p^{2} + 2\gamma T_{1}p + 1)(T_{3}p + 1)} .

3.2. Принцип алгоритма состоит в минимизации квадратических функционалов от разности откликов исследуемого средства измерений и моделей его динамических свойств (выбираемых последовательно из указанной иерархии) на заданный испытательный сигнал и нахождении той модели, которая дает значение функционала, не превышающее заданного.

Комплекс программ, реализующих указанный выше принцип вычислений, включен в отраслевой фонд алгоритмов и программ Госстандарта СССР (ФАП АИУС ГИВЦ Госстандарта СССР). Комплекс включает программы:

для нахождения простейшей по структуре оптимальной передаточной функции,

для оценивания погрешностей найденных коэффициентов передаточной функции,

для пересчета найденной передаточной функции в другие полные динамические характеристики.

Приложение 9 Рекомендуемое

УСТАНОВЛЕНИЕ КОЭФФИЦИЕНТА ЗАПАСА Kз

Выбор испытательного сигнала и установление коэффициента Kз производят путем последовательных приближений:

принимают Kз0 = 1;

проверяют, удовлетворяет ли сигнал x(t) неравенствам (13), (14) п. 3.7 рекомендации; при нарушении неравенства выбирают более широкополосный сигнал;

находят требуемую динамическую характеристику исследуемого средства измерений;

оценивают погрешность определения динамической характеристики средства измерений;

оценивают неточность \Delta_{\Delta}(t) оценки временной зависимости \Delta(t) искажений испытательного сигнала, обусловленную погрешностью определения динамической характеристики исследуемого средства измерений; при условии \Delta_{\Delta\text{ч}}/\Delta_{\text{ч}} \leq 0,3 операции выбора заканчивают;

при условии \Delta_{\Delta\text{ч}}/\Delta_{\text{ч}}>0,3 принимают Kз1 > Kз0 и операции выбора повторяют.

Приложение 10 Справочное

ПРИМЕРЫ ОПРЕДЕЛЕНИЯ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СРЕДСТВ ИЗМЕРЕНИЙ

1. Контроль нормированной динамической характеристики

1.1 Исходные данные

Исследуемое средство измерений — индукционный расходомер с самопишущим измерительным (вторичным) устройством (по ГОСТ 11988).

Требуемая динамическая характеристика — переходная характеристика.

Требуемая форма представления динамической характеристики — графическая.

Предварительная информация — модель динамических свойств расходомера представляет собой динамическое звено 1-го порядка с постоянной времени T, номинальное значение которой Tн = 1 с; допускаемые отклонения переходной характеристики от номинальной по модулю не должны превышать \Delta h = 15\% .

Имеющийся в наличии генератор испытательного сигнала — мера расхода жидкости, обеспечивающая линейное нарастание расхода от 0 до 0,03 м3/ч за время \tau_{\text{ф}} , не превышающее 0,02 с; погрешность воспроизведения сигнала не превышает 2%.

Требуемая точность определения динамической характеристики — не установлена.

Требуемый диапазон частот, в котором должен быть исследован расходомер, — не установлен.

1.2 Последовательность операций

1.2.1 Выбор метода контроля.

По п. 2.2.1 рекомендации контроль осуществляется по постоянной времени T.

Допускаемые относительные отклонения \delta_{\text{т}} параметра T от Tн определяются условием \delta_{h} \leq e\Delta_{h}/1 - e\Delta_{h} . Для \Delta_{h} = 15\% получается \delta_{\text{т}} \leq 69\% , чему соответствуют пределы допускаемых значений T:Tmin = 0,31 с, Tmax = 1,69 с.

1.2.2 Выбор испытательного сигнала.

По п. 2.2.3 рекомендации в качестве испытательного выбирается сигнал, достаточно близкий по форме к единичной ступенчатой функции.

1.2.3 Оценка неидеальности испытательного сигнала.

Погрешность от неидеальности испытательного сигнала (приложение 4, 1-я строка таблицы) вследствие выполнения условия \tau_{\text{ф}max}/T_{min} = 0,06<0,3 оценивается по формуле {}_{t}^{max}|\Delta_{h.\text{с}}(t)| \leq \tau_{\text{ф}}/2T = 0,01 .

Значение {}_{t}^{max}|\Delta_{h,\text{с}}(t)| пренебрежимо мало по сравнению с \Delta_{h} и испытательный сигнал можно считать близким к характеристическому.

1.2.4 Проверка требований п. 4.2 рекомендации.

Искажение испытательного сигнала \Delta(t) = 1/ky(t) - x(t) , где y(t) выходной сигнал индукционного расходомера; x(t) — испытательный сигнал,

m_{t}ax\Delta(t) = (1 - e^{- \tau_{\text{ф}}/T}) \frac{T}{\tau_{\text{ф}}} = 0,93 при t = \tau_{\text{ф}}.

В качестве регистрирующего используется вторичное устройство расходомера. Таким образом,

\frac{max\mu(t)}{x_{max}} \leq \frac{max\Delta_{x}(t)}{x_{max}} + \frac{max\Delta_{y}(t)}{kx_{max}} = 3\%,

\frac{max\mu(t)}{x_{max}} \leq 0,2m_{t}ax\Delta(t) = 19,8\%.

1.2.5 Нахождение контролируемого параметра.

Исследуемый расходомер подключается к мере, и выходной сигнал записывается на ленте самопишущего прибора. За начало отсчета времени принимается момент, смещенный относительно момента включения меры на время запаздывания. Время запаздывания определяется скоростью движения жидкости по трубопроводу от меры к расходомеру и длиной трубопровода.

Оценка \tilde{T} параметра T находится как длительность интервала времени от начала отсчета до момента, когда выходной сигнал достигнет уровня, составляющего 63% от установившегося значения. Оценка \tilde{T} = 1,42c .

Найденная оценка сравнивается с границами диапазона допускаемых значений T. Пределы допускаемых отклонений могут быть установлены с учетом и без учета погрешности контроля. В первом случае найденная оценка удовлетворяет неравенству T_{min}<\tilde{T}<T_{max} , что свидетельствует о соответствии динамической характеристики установленным нормам. Во втором случае диапазон допускаемых значений T должен быть сужен на величину погрешности контроля, которая складывается из погрешности от неидеальности испытательного сигнала и инструментальной погрешности \mu(t) и не превышает 5%. Поэтому суженный диапазон для T имеет границы при \delta_{T}\prime \leq e10\%/1 - e10\% = 37\% ; T_{m}\prime_{in} = 0,63c ; T_{m}\prime_{ax} = 1,37c . Оценка \tilde{T} оказывается вне диапазона допускаемых значений и, следовательно, динамические свойства расходомера не отвечают установленным нормам.

2. Определение динамической характеристики экземпляра средств измерений данного типа

2.1 Исходные данные

Исследуемое средство измерений: осциллографический гальванометр типа М1013-1800 № 001184.

Требуемая динамическая характеристика — переходная характеристика.

Требуемая форма представления динамической характеристики — аналитическая.

Предварительная информация — модель динамических свойств гальванометра представляет собой колебательное звено 2-го порядка со степенью успокоения \gamma в пределах 0,1 — 0,4 и собственной частотой f0 не менее 1800 Гц.

Имеющийся в наличии генератор испытательного сигнала — генератор импульсов точной амплитуды Г5-75, создающий на согласованной нагрузке прямоугольные импульсы напряжения с длительностью 50 нс — 1 с (в режиме внутреннего запуска) и фронтом не более 10 нс. Погрешность установки амплитуды импульса: от дискретности установки амплитуды не более 0,5 мВ на поддиапазоне 1 — 9,9999 В; в динамическом режиме не более (0,01·103 + 2) мВ = 12 мВ; от выбросов на вершине основного импульса — не более 1% от амплитуды импульса, что соответствует 10 мВ; от паразитной модуляции параметров — не более 1 мВ. Результирующая погрешность установки амплитуды импульса с номинальным значением 1 В равна

\delta_{\text{и}} = \frac{2,35}{10^{3}} \cdot 100\% = 2,35\%.

Требуемая точность определения динамической характеристики — погрешность определения \gamma и f0 не должны превышать 30 и 25%, соответственно.

Требуемый диапазон частот, в котором исследуется гальванометр — от 0 до 1800 Гц.

2.2 Последовательность операций

2.2.1 Выбор метода и испытательного сигнала

По п.п. 2.2.3, 2.4.2, 2.6 рекомендации выбирается прямой метод определения переходной характеристики с использованием испытательного сигнала, достаточно близкого к единичной ступенчатой функции. Длительность импульса выбирается максимальной для данного генератора и равной 1 с.

2.2.2 Оценка неидеальности испытательного сигнала

Погрешность от неидеальности испытательного сигнала (сигнал с линейным нарастанием до установившегося значения вместо сигнала, описываемого функцией 1(t) - приложение 3, последняя строка таблицы):

m_{t}ax|\Delta_{h,c}(t)| \leq \gamma \sqrt{1 - \gamma^{2}} \cdot \frac{\tau_{\text{ф}}}{T} l^{\frac{\gamma}{\sqrt{1 - \gamma^{2}}}} \cdot arctg(-^{\frac{\gamma}{\sqrt{1 - \gamma^{2}}}} {}^{)}.

Подставляя T = \frac{\sqrt{1 - \gamma^{2}}}{2\pi f_{0}} , получим

Погрешность \Delta_{h,c}(t) пренебрежимо мала, и можно считать, что испытательный сигнал практически совпадает с характеристическим.

2.2.3. Проверка требований п. 3.3 рекомендации (приложение 2, 1-я строка таблицы). Поскольку погрешность \Delta_{\text{в}} не задана, принимается максимальное значение \delta_{\text{в}} = 1 . При этом f01 ~= 1 МГц. Таким образом, условие f01 > fв выполнено.

2.2.4 Проверка требований п. 4.2 рекомендации

Искажения испытательного сигнала

\Delta(t) = \frac{1}{k} y(t) - x(t),

где y(t) — выходной сигнал;

k — статический коэффициент преобразования;

x(t) — испытательный сигнал исследуемого средства измерений.

\Delta(t) = - (\frac{\gamma}{\sqrt{1 - \gamma^{2}}} sin2\pi f_{0}t + cos2\pi f_{0}t)e^{-}^{\frac{2\pi f_{0}\gamma}{\sqrt{1 - \gamma^{2}}}} {}^{t}.

Максимум функции \Delta(t) получим при t = 1/2f0

В качестве регистрирующего прибора используется светолучевой осциллограф типа Н115, амплитудная погрешность которого не превышает 10%. Учитывая амплитуду входного сигнала, равную 1 В, получим

\frac{max\mu(t)}{1} \leq \frac{max\Delta_{x}(t)}{1} + \frac{max1/k\Delta_{y}(t)}{1} = 12,35\%.

\frac{max\mu(t)}{1} \leq 0,2 \frac{max\Delta(t)}{1} = 14,6\%.

2.2.5 Нахождение параметров \gamma и f0

Отклик осциллографического гальванометра на ступенчатый испытательный сигнал

Ввиду относительно малого успокоения гальванометра параметры находятся по отклику (см. рисунок) простейшим графоаналитическим способом с помощью выражений:

f_{0} = \frac{1}{2(t_{2} - t_{1})} , \gamma = \frac{1}{\sqrt{1 + (\frac{\pi}{H})2}} ,

где H = ln\xi_{1}/\xi_{2} ; \xi_{1} = y_{max} - y_{\text{уст}} — амплитуда первого выброса;

ymax, yуст, ymin — максимальное, установившееся и минимальное (не равное нулю) значение отклика;

t1, t2 моменты времени, соответствующие первому и второму пересечениям кривой отклика установившегося значения.

Значения отклика: yуст = ₂,5 В; ymax = 4,5 В; yуст = 1,55 В; на основании чего \gamma = 0,173 ; f0 = 2857 Гц при t1 = 0,075 мс; t2 = 0,25 мс.

2.2.6 Оценивание погрешностей \gamma и f0

На основании выражений для \gamma и f0 погрешности их определения равны

\Delta f_{0} = 2f_{0}^{3}\Delta_{\text{т}}, \Delta_{\gamma} = \pi^{2} \frac{\gamma^{3}}{H^{3}} (\frac{1}{\xi_{1}} \Delta\xi_{1} + \frac{1}{\xi_{2}} \Delta\xi_{2}),

где \Delta_{\text{т}} — погрешность определения интервала времени [t₁, t₂];

\Delta\xi_{1} — погрешность определения амплитуды первого выброса \xi_{1} ;

\Delta\xi_{2} — погрешность определения амплитуды \xi_{2} .

Поскольку из ряда осциллографических гальванометров, входящих в комплект Н115, был выбран гальванометр типа М014-7000 с рабочей полосой частот 0 — 5000 Гц и собственной частотой fог = 7000 Гц (что намного больше верхнего предела частотного диапазона, в котором исследуется гальванометр), причем перед испытаниями он был настроен на оптимальное \gamma , выбранное в соответствии с [1], пренебрегаем погрешностью, вносимой неравномерностью амплитудно-частотной характеристики гальванометра как пренебрежимо малой по сравнению со статической погрешностью.

Так как за начало отсчета при определении отклика осциллографического гальванометра принято установившееся значение отклика (т.е. измерение амплитуды отклика проводится относительно установившегося значения), то погрешность измерения амплитуды по осциллограмме равна

\Delta\xi_{1} = \Delta\xi_{\text{ампл}} + 2\Delta_{\text{отсч}.\text{у}},

где \Delta\xi_{\text{ампл}} — погрешность, зависящая от амплитудной погрешности гальванометра, \Delta\xi_{\text{ампл}} = 0,1\xi ;

\Delta_{\text{отсч}.\text{у}} — погрешность отсчитывания по оси y.

Поскольку масштаб по оси y My = 0,05 В/мм (максимальное регистрируемое напряжение 6 В при ширине фотоленты 120 мм), а погрешность отсчитывания по оси y при использовании оптических средств достаточно мала, в данном случае она равна 0,125 мм, то \Delta_{\text{отсч}.y} = 0,125 \times 0,05\text{В}/\text{мм } = 6,25\text{мВ} . Представляя значение для \xi_{1} и \xi_{2} получим

\Delta\xi_{1} = 0,1\xi_{1} + 2\Delta_{\text{отсч}.y} = 0,1775\text{В},

\Delta\xi_{2} = 0,1\xi_{2} + 2\Delta_{\text{отсч}.y} = 0,1075\text{В},

При этом \Delta_{j} = 0,06 ; \delta_{j} = 35,2\% .

Так как полученное значение \delta_{j} = 35,2\% превышает допустимое значение погрешности определения и нет возможности использовать осциллограф с меньшим значением амплитудной погрешности, то возникает необходимость уменьшить статическую погрешность прибора. Это можно сделать путем коррекции систематической составляющей погрешности осциллографа: установки нулевой линии, калибровки, учета погрешности нелинейности прибора [2]. В результате статическая погрешность не превышает 5%, а погрешность определения \gamma

\delta_{\gamma} \leq 20,5\%. \Delta_{\gamma} \leq 0,035;

Погрешность \Delta_{\text{т}} определения интервала времени [t₁, t₂] определяется, в основном, погрешностью отсчитывания \Delta t_{\text{отсч}} и погрешностью развертывания \Delta_{\text{р}} , т.е. \Delta_{\text{т}} = \Delta t_{\text{отсч}} + \Delta_{\text{р}} . Погрешность отсчитывания \Delta t_{\text{отсч}} зависит от \Delta t_{1\text{отсч}} , \Delta t_{2\text{отсч}} — погрешностей отсчитывания моментов времени t1, t2, определяющих интервал [t1, t2], т.е. \Delta t_{\text{отсч}} = \Delta t_{1\text{отсч}} + \Delta t_{2\text{отсч}} . Применение оптических увеличительных устройств позволяет уменьшить погрешность отсчитывания до значения 0,125 мм, что при скорости v = 104 мм/с соответствует \Delta t_{1\text{отсч}} = \Delta t_{2\text{отсч}} = 12,5\text{мкс} . Отсюда \Delta t_{\text{отсч}} = 25\text{мкс} .

Определенная экспериментально [3] неравномерность движения фотоленты прибора Н115 при скорости v = 104 мм/с на участке l = 10 мм составляет \delta \leq 6,5\% , что на интервале [t1, t2] обуславливает появление погрешности развертывания \Delta_{\text{р}} = (t_{2} - t_{1}) \cdot \delta = 0,175 \cdot 0,065 = 0,011\text{мс} . При этом \Delta_{\text{т}} = 0,036\text{мс} , \Delta f_{0} = 588\text{Гц} , \delta f_{0} = 20,6\% . Таким образом, точность определения параметров переходной характеристики гальванометра (с учетом коррекции систематической составляющей статической погрешности регистрирующего прибора) удовлетворяет требованиям.

Следует заметить, что экспериментатор не всегда располагает данными о неравномерности движения носителя для конкретного типа прибора и соответствующей скорости.

Зная период повторения меток tмет и период повторений T систематической составляющей \Delta_{\text{рс}} погрешности \Delta_{\text{р}} [4]: T = \pi d/v , где d — диаметр ведущего вала лентопротяжного механизма, v — скорость движения фотоленты, можно получить достаточно близкую оценку «снизу» погрешности \Delta_{\text{т}} в следующих случаях:

1. когда определяемый интервал t₂ - t₁ намного меньше периода T(t2 — t1 << T) и метки времени не позволяют уменьшить систематическую \Delta_{\text{рс}} составляющую погрешности развертывания на интервале [t1, t2];

2. когда tмет << T и метки времени позволяют существенно уменьшить \Delta_{\text{рс}} на исследуемом интервале. В примере время затухания переходного процесса составляет несколько миллисекунд и максимальная для прибора Н115 частота меток fмет = 500 Гц слишком мала для такого процесса, к тому же с целью получения удобной для расшифровки осциллограммы выбрана наибольшая скорость движения фотоленты v = 104 мм/с, поэтому на всей осциллограмме будет не более двух меток, т.е. имеет место первый случай. Очевидно, что при выполнении условия t2 — t1 << T, систематическая составляющая погрешности развертывания на интервале [t1, t2] много меньше максимального значения погрешности \Delta_{\text{р}} . Пренебрегая случайной и систематической составляющей \Delta_{\text{р}} на интервале [t1, t2], оценим \Delta_{\text{т}} «снизу» \tilde{\Delta}_{\text{т}} = \Delta_{\text{тотсч}} = 25\text{мкс} . При этом получим оценку «снизу» для погрешности определения частоты \Delta f_{0} = 409\text{Гц} , \delta_{f0} = 14,3\% .

Рассмотрим второй случай, когда существенное значение для уменьшения систематической составляющей погрешности развертывания имеет соотношение периодов tмет и T; т.е. когда выполняется условие tмет << T. Поскольку наличие меток времени позволяет уменьшить систематическую составляющую погрешности \Delta_{\text{р}} до значения, обусловленного нестабильностью скорости протяжки на интервале tмет, то при условии tмет << T этой нестабильностью скорости на интервале [t1, t2] можно пренебречь и тогда, когда t2 — t1 > tмет (при этом систематическая составляющая погрешности \Delta_{\text{р}} на интервале [t1, t2] равна систематической составляющей погрешности развертывания на интервале tмет, т.е. пренебрежимо мала по сравнению с максимальным значением \Delta_{\text{р}} ), и тогда, когда t2 — t1 < tмет << T. Как и в первом случае, погрешность \Delta_{\text{т}} можно оценить «снизу»: \Delta_{\text{т}} = \Delta t_{\text{отсч}} .

ЛИТЕРАТУРА

1. Бренер М.Д., Шкабардня М.С. О погрешности измерений с помощью светолучевых осциллографов в динамическом режиме. — Метрология, 1980, № 8, с. 22 — 31.

2. Шкабардня М.С. Теория и принципы построения быстродействующих самопишущих приборов. — М.: Энергоатомиздат, 1984. — 214 с.

3. Артамонов В.М. и др. Светолучевые осциллографы. — Л.: Энергоиздат. Ленинградское отделение, 1982. — 104 с.

4. Грановский В.А. Динамические измерения: Основы метрологического обеспечения. — Л.: Энергоатомиздат. Ленинградское отделение, 1984. — 224 с.

5. ГОСТ 8.256-77. ГСИ. Нормирование и определение динамических характеристик аналоговых средств измерений.

6. МИ 1951-88. ГСИ. Динамические измерения. Термины и определения.

7. ГОСТ 8.207-76. ГСИ. Прямые измерения с многократными наблюдениями. Методы обработки результатов наблюдений. Основные положения.

8. ГОСТ 8.061-80. ГСИ. Содержание и построение поверочных схем.

---