Способ измерения изменения профиля поля физической величины

Изобретение относится к способам контактного изменения профиля физической величины в различных средах, в частности профиля температуры в море или атмосфере. При осуществлении способа измерения изменения профиля поля физической величины используют распределенные датчики с переменной погонной функцией чувствительности. Погонную функцию чувствительности этих распределенных датчиков выполняют по весовым функциям пространственных многополосных фильтров-датчиков, границы пропускания которых устанавливают в соответствии с покрытием подмножества номеров орт ряда Фурье v-м значением s-го разряда b-ичного кода номера орты от 0 до N-1. Весовые функции многополосных фильтров-датчиков выполняют равными сумме значений нормированных орт в полосах пропускания. Коэффициенты разложения по ортам изменения профиля за заданное время определяют сверткой за это время выходных сигналов многополосных фильтров-датчиков, соответствующих значениям b-ичного n-го разрядного кода номера орты. Операции свертки производят одновременно по древовидной схеме b-ичного дешифратора коэффициентов орт. Изменение профиля поля физической величины θ(x,t,T) за время Т вычисляют по формуле

где х - пространственная координата на профиле от 0 до L;

t - текущее время;

Т - задаваемое извне время интегрирования, за которое измеряется изменение профиля поля физической величины;

ϕ_i(х) - орта ряда Фурье, аппроксимирующего профиль θ(х);

i - номер орты в b-ичном коде;

a _i(t,T) - текущее изменение коэффициента разложения профиля по орте ϕ_i(х) за время Т. Технический результат изобретения заключается в повышении точности и пространственной разрешающей способности измерений. 5 ил.

 

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

Известен способ измерения профиля температуры на конечном интервале в различных средах с использованием нескольких распределенных датчиков с переменной погонной функцией термочувствительности, модулированной по функциям из ортогонального базиса, реализованный в устройстве [АС СССР 808872, Кл. G01F 23/22. Устройство для измерения температуры / В.А. Гайский. - Опубл. 28.02.81. Бюл. №8].

Этим способом измеряются коэффициенты разложения искомого профиля температуры в ряд Фурье по функциям из этого базиса. При этом необходимо столько распределенных модулированных датчиков, сколько членов разложения (орт) используется. Известно, например, [Латхи Б.П. Системы передачи информации. Перевод с англ., под общей редакцией Б.И. Кувшинова. - М.: «Связь», 1971. - 324 с.], что среднеквадратическая погрешность представления профиля разложением в ряд уменьшается с ростом числа орт.

Например, при представлении функции профиля температуры рядами Уолша на базе L с пространственным разрешением L/N требуется N орт и, соответственно, N распределенных вдоль базы датчиков, являющихся узкополосными фильтрами орт.

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

Во многих областях контроля окружающей среды и технологических процессов представляет интерес не статический профиль поля физической величины, а его изменение во времени [Гайский В.А., Егупов Н.Д., Корнюшин Ю.П. Применение функций Уолша в системах автоматизации научных исследований. - Киев: Наукова думка, 1983. - 211 с]. Например, это относится к измерению параметров поверхностных и внутренних волн в море, параметров акустических и электромагнитных волн, температуры в системах пожарной сигнализации. Для измерения изменений профиля физической величины может быть применен более экономичный или более точный способ, чем для измерения статического профиля.

Способ, заложенный в основу работы измерителя [АС СССР 808872, Кл. G01F 23/22. Устройство для измерения температуры / В.А. Гайский. - Опубл. 28.02.81. Бюл. №8], выбран в качестве прототипа, потому что наиболее близко совпадает с заявленным изобретением по признакам и технической сущности.

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

Поставленная задача решается тем, что для измерения изменения коэффициентов орт представления профиля поля физической величины используют nb распределенных датчиков с погонными функциями чувствительности, модулированными по весовым функциям пространственных многополосных фильтров-датчиков, границы окон пропускания которых устанавливают в соответствии с покрытием подмножества номеров орт v-м значением s-го разряда позиционного b-ичного кода номера орт от 0 до N-1, весовые функции многополосных фильтров-датчиков выполняют равными сумме значений нормированных орт в полосах пропускания, коэффициенты разложения по ортам изменения профиля за заданное время определяют сверкой за это время выходных сигналов датчиков, соответствующих значениям b-ичного n-разрядного номера орты, операции свертки производят одновременно по древовидной схеме b-ичного дешифратора коэффициентов орт, изменение профиля поля физической величины θ(x,t,T) за время Т вычисляют по формуле

где х - пространственная координата на профиле от 0 до L;

t - текущее время;

Т - задаваемое извне время интегрирования, за которое измеряется изменение профиля поля физической величины;

ϕ_i(х) - ортонормированная функция ряда Фурье (орта), аппроксимирующего профиль θ(х);

i - номер орты в b-ичном коде;

a _i(t) - текущее значение коэффициента разложения профиля по орте ϕ_i(x);

a _i(t,T) - текущее изменение коэффициента разложения профиля по орте ϕ_i(х) за время Т:

где

y_s,v(t) - выходной сигнал s, v-гo многополосного фильтра-датчика, в значениях входной величины соответствующего v-му значению s-го разряда i-го номера орты в b-ичном коде;

причем 1 - старший разряд.

При использовании для представления изменений профиля поля физической величины N=2ⁿ ортонормированных функций Уолша-Адамара ϕ_i(x)=had_i(x) [Гайский В.А., Егупов Н.Д., Корнюшин Ю.П. Применение функций Уолша в системах автоматизации научных исследований. - Киев: Наукова думка, 1983. - 211 с], номера i которых кодируют в двоичном коде (b=2, v={0,1}) n разрядами, весовые функции 2n многополосных фильтров-датчиков выполняют равными

g_s,v(x)=0, при x≠1, x≠2^n-s, для ,

и изменение профиля поля физической величины θ(x,t,T) за время T определяют по формуле

при этом a_i(t,T) определяется по выражению (2).

В предлагаемом способе используется m=bn датчиков вместо N=bⁿ, т.е. в раз снижается сложность реализации измерительного кабеля или во столько же раз увеличивается пространственная разрешающая способность (для функций Уолша) и существенно повышается точность при сохранении числа датчиков. Следует отметить, что модуляция весовой функции многополосных фильтров-датчиков может быть сложнее или проще модуляции весовой функции орт. Как показано ниже, для многополосных фильтров-датчиков функций Уолша-Адамара весовые функции проще весовых функций орт. Достигаемый выигрыш в количестве орт N, коэффициенты которых измеряются при реально возможном числе датчиков т, двоичном (b=2) и тройничном (b=3) кодировании, показан в таблице на фиг. 1.

Описание изобретения иллюстрируется чертежами.

На фиг. 2 показана структурно-функциональная схема формирования многополосных фильтров по кодовым шкалам и дешифрация коэффициентов орт свертками для двоичного (а) и троичного (б) кодирования номеров орт.

На фиг. 3 приведен пример применения предлагаемого способа с использованием для представления изменения измеряемого профиля поля физической величины функций Уолша-Адамара.

На фиг. 4 приведена структурно-функциональная схема формирования выходных сигналов y_sv(t) распределенных датчиков с весовыми функциями пространственных фильтров Уолша-Адамара g_sv(x) на базе L из восьми участков в трех разрядах двоичного кода номера i орты и выделения двоичным дешифратором коэффициентов орт свертками в разрядах.

На фиг. 5 показана структурная схема устройства для измерения изменения профиля температуры по предлагаемому способу при использовании фильтров функций Уолша-Адамара и реализации модуляции весовых функций фильтров петлевой укладкой термочувствительных проводов вдоль профиля.

Рассмотрим суть предлагаемого способа.

Предполагаем, что измеряемый мгновенный профиль поля физической величины θ(x,t) вдоль пространственной координаты х на базе [0,L] представляется с достаточной точностью рядом Фурье в виде

где ϕ_i(х) - функции (орты) из ортонормированного базиса, обладающие свойством

Коэффициенты a_i(t) при ортах ϕ_i(х) измеряются известным способом распределенными датчиками по формуле

Если за некоторое время Т изменится профиль физической величины на θ(x,t,T), то изменятся и коэффициенты на a_i(t,T), т.е.

Измерение a_i(t,T) известным способом с использованием одного распределенного датчика на каждую орту, реализующего формулу (7), требует N датчиков, что ограничивает возможности повышения точности и разрешающей способности за счет увеличения N из-за возрастания сложности реализации.

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

В этих случаях, когда информативным является не статический мгновенный профиль θ(x,t) физической величины, а его изменение за заданное время T, т.е. θ(x,t,T), возможно уменьшение числа необходимых распределенных датчиков при том же числе N измеряемых изменений коэффициентов a_i(t,T) орт.

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

Упорядочим орты в их естественной последовательности . Кодируем номера i орт в b-ичной позиционной системе счисления n-разрядным словом, в которой каждый разряд может принимать "b" значений v={0,1, …, (b-1)}. Каждому значению v s-го разряда кода номера орты поставим в соответствие подмножество G_sv номеров орт, покрываемых шкалой кодирования так, как показано на фиг. 2. В каждом подмножестве G_sv будет N/b орт, а подмножеств будет nb. Столько же многополосных фильтров-датчиков используется, причем подмножества G_sv определяют границы многополосных фильтров, а весовые функции g_sv(x) датчиков выполняются суммированием значений орт из подмножеств G_sv по формуле

При этом выходные сигналы многополосных фильтров-датчиков y_sv{t) будут изменяться с изменением θ(x,t) во времени

Таким образом, каждому номеру i орты в b-ичном коде соответствует свой набор пространственных многополосных фильтров-датчиков. Если бы это были временные фильтры и временные орты, то их коммутация в последовательной цепочке согласно с кодом номера орты позволяла бы сформировать окно пропускания только для орты с этим номером. Поскольку здесь пространственные многополосные фильтры-датчики включены параллельно, то для выделения приращения коэффициента избранной орты, как общей составляющей в выходных сигналах всех датчиков, необходимо произвести последовательно свертку выходных сигналов датчиков y_sv(t) на заданном интервале времени T, соответствующих v элементам b-ичного кода в s-м разряде номера орты.

Для выделения сигналов изменения коэффициентов всех орт свертка выходных сигналов всех датчиков выполняется по древовидной схеме дешифратора b-ичного кода по выражению

где i=[1,2,3,…n_v] номер орты представлен в 6-ичном коде и индексация разрядов кода номера орты совпадает с индексацией весовых функций g_sv датчиков и их выходных сигналов y_sv(t).

Подстановка выражения (11) в выражение (12) дает

В выражении (13) произведения сумм всех орт дадут выражения вида

В итоге из выражения (13) получим

Окончательно получим

Таким образом, для измерения изменений N коэффициентов орт представления пространственного профиля физической величины потребовалось logbN распределенных многополосных фильтров-датчиков вместо N. Поскольку для функций Уолша N=2ⁿ, то для них целесообразно кодирование номеров орт двоичным кодом (b=2). Как раз здесь имеет место случай, когда известная экономичность физического представления слов при троичном кодировании (b=3) дает минимальное число распределенных датчиков для измерения. Целесообразность кодирования с другими основаниями позиционной системы счисления остается не ясной.

Рассмотрим реализацию способа с использованием представления профиля поля физической величины функций Уолша-Адамара had_i{x) [Гайский В.А., Егупов Н.Д., Корнюшин Ю.П. Применение функций Уолша в системах автоматизации научных исследований. - Киев: Наукова думка, 1983. - 211 с]. На фиг. 3 показана матрица функций Уолша-Адамара для N=8 (а), двоичная шкала кодирования номеров орт (б), формирование подмножеств орт G_sv в соответствии с разрядами кодирования и весовых функций g_sv(х) многополосных фильтров суммой орт G_vs (в) и, наконец, сводка весовых функций многополосных фильтров g_vs(x) Уолша-Адамара (г).

Отметим особенности весовых функций g_sv(x) многополосных функций Уолша-Адамара. Во-первых, они имеют уже не два (как had_i(х)) значения, а три (+,0,-). Во-вторых, поскольку амплитуда всех g_sv(x) равна |N/2|, то общий множитель N12 может быть вынесен и модуляция g_sv(x) может выполняться значениями (+1,0,-1). В третьих, модуляция весовых функций g_sv(х) проще, чем модуляция had_i(x), поскольку реально только две координаты на базе [0,L] имеют не нулевые значения. В-четвертых, значения весовых функций многополосных фильтров Уолша-Адамара определяются выражениями

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

Легко убедиться в работоспособности способа на конкретных числовых примерах, представленных на фиг. 3 и 4. При задании значения Т-интервала времени изменения профиля поля физической величины, необходимо учитывать следующее:

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

- при цифровой регистрации выходных сигналов y_sv(t) датчиков число отсчетов каждого из них во времени для реализации свертки их значений должно быть не менее N;

- при реализации операций сверки интервал интегрирования T, как в реальном времени, так и при апостериорной обработке данных, может быть выбран скользящим от t-Т до t или с синхронизацией нуля от t=0 до t=Т;

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

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

Устройство содержит измерительный кабель (ИК) 1, блок измерительных преобразователей (БИП) 2, вычислитель (В) 3.

Измерительный кабель 1 состоит из распределенных многополосных фильтров-датчиков с весовыми функциями g_sv(x) термочувствительных проводов, погонная модуляция термочувствительности которых выполнена продольной петлевой укладкой проводов так, как это показано на фиг. 5, в соответствии с фиг. 3 и 4.

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

Вычислитель 3 предназначен для восстановления мгновенных значений температуры на входе датчиков и выполнения операций свертки за время Т с вычислением изменения коэффициента a_i(t,T) при had_i(x) и вычисления изменений профиля θ(х,t,Т) по формуле (8).

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

Измерительный кабель 1 размещается по траектории х профиля θ(x,t), распределенные многополосные фильтры-датчики g_sv(x) отслеживают изменение температуры θ(х,t) на профиле соответствующим изменением сопротивлений, которые преобразуются блоком измерительных преобразователей в цифровые последовательности, поступающие далее в вычислитель 3. Вычислитель 3, на основании градуировочных характеристик измерительных каналов, включающих датчик и измерительные преобразователи, а также динамики выходных сигналов, восстанавливает известным способом мгновенное значение средневзвешенной по профилю температуры для каждого из датчиков, что и является выходным сигналом y_sv{t) соответствующего датчика. Далее вычислитель производит свертку выходных сигналов датчиков за заданное время Т, вычисляя сигналы z_i(t,T) по формуле (12), определяет изменения коэффициентов a_i(t,T) орт по формуле (16) и восстанавливает значения изменения профиля температуры θ(х,t,Т) по формуле (8).

Способ измерения изменения профиля поля физической величины с использованием распределенных датчиков с переменной погонной функцией чувствительности, отличающийся тем, что используют распределенные датчики, погонную функцию чувствительности которых выполняют по весовым функциям пространственных многополосных фильтров-датчиков, границы пропускания которых устанавливают в соответствии с покрытием подмножества номеров орт ряда Фурье v-м значением s-го разряда b-ичного кода номера орты от 0 до N-1, весовые функции многополосных фильтров-датчиков выполняют равными сумме значений нормированных орт в полосах пропускания, коэффициенты разложения по ортам изменения профиля за заданное время определяют сверткой за это время выходных сигналов многополосных фильтров-датчиков, соответствующих значениям b-ичного n-го разрядного кода номера орты, операции свертки производят одновременно по древовидной схеме b-ичного дешифратора коэффициентов орт, изменение профиля поля физической величины θ (x,t,T) за время Т вычисляют по формуле

где х - пространственная координата на профиле от 0 до L;

t - текущее время;

Т - задаваемое извне время интегрирования, за которое измеряется изменение профиля поля физической величины;

ϕ_i (x) - орта ряда Фурье, аппроксимирующего профиль θ (х);

i - номер орты в b-ичном коде;

а_i (t) - текущее значение коэффициента разложения профиля по орте ϕ_i (x);

а_i (t, T) - текущее изменение коэффициента разложения профиля по орте ϕ_i (x) за время T:

где

y_s,v (t) - выходной сигнал s, v-го многополосного фильтра-датчика, в значениях входной величины соответствующего v-му значению s-го разряда i-го номера орты в b-ичном коде;

причем 1 - старший разряд,

при этом в случае использования базиса функций Уолша-Адамара ϕ_i (x)=had_i(x), номера i которых в естественном порядке кодируют в двоичном коде (b=2, v={0, 1}) n разрядами, используют 2n многополосных фильтров-датчиков, весовые функции g_s,_v (x) которых выполняют по формулам

при v={0, 1};

при v=1;

при v=0;

g_s,v (x)=0, при x≠1, х≠2_n-s + 1,

и изменение профиля поля физической величины θ (x,t,T) за время Т вычисляют по формуле