Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Бендат Дж. Измерение и анализ случайных процессов

.pdf
Скачиваний:
73
Добавлен:
25.10.2023
Размер:
16.48 Mб
Скачать

258 Глава 7

формы с последующим ограничением суммарного сигнала, как* показано на рис. 7.5, в. Ширина получаемого импульса пропор­ циональна значению входного сигнала в соответствующий момент времени. При кодово-импульсной модуляции (КИМ) дискретные отсчеты сигнала преобразуются в последовательность импульсов, как показано на рис. 7.5, г; эти импульсы затем кодируются, обыч­ но в двоичном коде. Сведения о полосе пропускания и отношении сигнал/шум, получаемых при каждом из описанных методов им­ пульсной модуляции, содержатся в табл. 7.1, где указаны также основные достоинства и недостатки этих методов.

На практике число преобразователей сигнала часто превышает число имеющихся каналов для передачи и (или) регистрации сиг-,, нала. Поэтому иногда по одному и тому же каналу для передачи и (или) регистрации приходится передавать сигналы с нескольких преобразователей одновременно. Это достигается тем или иным приемом коммутации или уплотнения сигналов. Для тех типов модуляции, которые требуют снятия отсчетов в дискретные мо­ менты времени, уплотнение сигналов осуществляется достаточно просто, поскольку кодированные сигналы разных преобразова­ телей нетрудно объединить и передать по одному каналу. При

<частотной модуляции уплотнение достигается в результате того, что различным сигналам отводятся неперекрывающиеся полосы частот, представляющие собой участки общей полосы пропуска­ ния канала передачи телеметрии или записи на магнитную ленту. Как показано в табл. 7.1, уплотнение осуществляется наиболее просто при АИМ- и КИМ-сигналах. Однако во всех случаях^ число уплотненных каналов влияет на требования к полосе про,- пускания линии передачи.

7.2.3. Временная ошибка при магнитной регистрации

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

систем.

Чаще всего рассматривается динамическая временная ошибка в одноканальной системе, называемая флаттером. Флаттер можно' определить как отклонение скорости движения магнитной ленты

Общие соображения о сборе и обработке данных

259

, от нормы1*. Он играет особенно важную роль при частотной моду-

'ляцни, когда аналогом амплитуды сигнала является масштаб времени или частота. Интеграл от флаттера и есть динамическая

временная ошибка. Величина флаттера у высококачественных магнитофонов промышленного производства составляет ~0,25% . Как правило, спектр флаттера практически постоянен на всех частотах до некоторой граничной частоты, которая при высокой

.скорости движения ленты равна обычно примерно 1000 Гц. Наиболее неприятная динамическая временная ошибка между

каналами появляется вследствие динамического перекоса, опре­ деляемого как изменение угла, под которым лента проходит маг­ нитную головку. Это явление, очевидно, влияет на результаты взаимного анализа сигналов в многоканальных системах, напри­ мер на взаимные корреляционные и спектральные функции. Динамический перекос можно минимизировать, записывая по возможности предназначенные для взаимного анализа данные через различные головки. Если данные будут анализироваться в цифровой форме, дальнейшего улучшения можно достичь, сме­ шивая в каждом канале синхронизирующий импульс с основным сигналом и затем используя этот импульс для управления анало­ го-цифровым преобразованием2*сигнала соответствующего канала (записанного на соответствующую дорожку).

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

^размеров ленты, что приводит к появлению временной ошибки, У высококачественных серийных магнитофонов с толщиной ленты 0,025 или 0,038 мм статическая временная ошибка состав­ ляет, как правило, от 0,1 до 0,25%.

Статическая временная ошибка в многоканальных системах объясняется скорее допусковыми отклонениями от номинала, чем непостоянством скорости движения ленты. Эти допуски включают в себя угол, под которым лента проходит головку

П В отечественной терминологии английскому термину «flutter» соответ­ ствует термин «коэффициент колебания скорости» и в частном случае для зву­ козаписывающей аппаратуры —«коэффициент детонации». Коэффициент ко­ лебания скорости представляет собой отношение отклонения от средней ско­ рости к значению средней скорости и выражается в процентах. Если скорость 'записи непостоянна, то шкала времени сигналограммы по сравнению с ре­ альным временем, в течение которого происходило изменение регистрируемо­ го сигнала, имеет мгновенную временную ошибку. Более подробно см., на­ пример, Дэвис Г. Л., Применение точной магнитной записи, изд-во «Энергия», 1967 или Атей С., Устройства записи на магнитную ленту, изд-во «Энергия»,' 1969.— Прим. ред.

2) Имеется в виду запуск синхронизирующим импульсом аналого-цифрог вого преобразователя соответствующего канала.— Прим. ред.

17*

260 Глава 7

(статический перекос) и положение головок (допуски в расстоянии между группами головок и между головками в отдельных груп­ пах, а также допуски в азимуте головок). Здесь необходимо отме­ тить, что записывающий и воспроизводящий каналы большинства магнитофонов промышленного производства обслуживаются дву­ мя различными группами головок, т. е. дорожки с четными и не­ четными номерами приходятся иа различные головки. Поэтому допуски от канала к каналу между дорожками при одной и той же головке оказываются меньшими, чем между дорожками при различных головках. Межканальная статическая временная ошиб­ ка становится особенно важной при анализе характеристик сиг­ налов, поступающих с двух или более каналов. Величина меж­ канальной временной ошибки зависит от характеристик исполь­ зуемой аппаратуры, но, как правило, она составляет около 0,025 мм. Поэтому величина временной ошибки зависит от ско­ рости движения магнитной ленты.

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

вработах [12] и [34].

7.3.Подготовка данных

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

чика). На этой стадии необходимо произвести

ряд операций,

цель

которых заключается в

подготовке

к дальнейшему

ана­

лизу

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

электрического напряжения.

Первая из

этих

операций

на

рис. 7.1 названа редактированием.

Под редактированием понимается последовательность опера­ ций, выполняемых до собственно анализа. Цель этих операций состоит в выявлении и исключении аномальных и (или) искажен­ ных сигналов, которые могут возникать при сборе и регистрации данных, например, за счет высокого уровня помех, снижения уровня сигнала, исчезновения сигнала при плохой работе датчи-^_ ка и т. д. Редактирование часто сводится просто к визуальному' анализу данных опытным специалистом. В сравнительно сложных системах сбора и обработки данных для оперативного редактиро­ вания может быть использовано специальное устройство. Чаще

Общие соображения о сборе и обработке данных

261

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

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

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

ленты,

что позволяет обеспечить

непрерывное поступление дан-

. ных на

аналоговый анализатор.

Если же анализ проводится в

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

7.3.1. Преобразование в цифровую форму

Представление сигнала в цифровой форме предполагает вы­ полнение двух независимых операций: а) дискретизации; б) кван­ тования. Дискретизация есть процесс определения точек1), в ко-

j-(t)

Р и с . 7.6. Дискретное представление непрерывноіНреализации.

торых должны быть произведены выборки; квантование есть пре­ образование отсчетоз, сделанных в соответствующие моменты времени, в цифровую форму.

Дискрет изация. При цифровой обработке данных эта опера­ ция обычно производится через равные интервалы времени, как показано на рис. 7.6. Задача заключается в правильном выборе -величины интервала дискретности /г. С одной стороны, при выбор­ ках (отсчетах) в слишком близко расположенных друг к другу точках будут получаться коррелированные данные, число кото­ рых будет избыточно велико, что неизбежно вызовет увеличение как объема, так и стоимости расчетов. С другой стороны, при вы­

1 То есть моментов времени.— Прим, перев.

862 Глава 7

борках в слишком далеко отстоящих друг от друга точках воз-/ можно перепутывание низко- и высокочастотных составляющих исходного процесса. Это явление называется маскировкой (на­ ложением) частот, оно представляет собой источник ошибок, которые отсутствуют при обработке аналогового сигнала; при цифровой обработке, когда сигнал предварительно преобразуется из аналоговой формы в дискретную, проблемы, связанные с мас­ кировкой частот, могу г играть важную роль. Так, например,все, кго видел классические фильмы «вестерны», обращал, несомненно, внимание на то, что при изменении скорости кареты кажется, что меняется и направление вращения ее колес. Это и есть про­ стейший пример ошибок маскировки, вызванных преобразова­ нием непрерывного процесса в дискретный, которое осуществляет­ ся съемочной камерой.

Чтобы пояснить явление маскировки частот, рассмотрим ил­ люстрируемый рис. 7.6 случай преобразования непрерывного сигнала в дискретную форму; интервал времени между последо­ вательными отсчетами1*составляет h с. Следовательно, скоростьдискретизации равна 1 Hi отсчетов в секунду. Для того чтобы дискретная реализация содержала все те же частоты, что и ис­ ходный непрерывный сигнал, на каждый цикл соответствую­ щего колебания должно приходиться по меньшей мере 2 отсчета. Поэтому наиболее высокая частота, которая может быть выделена при дискретизации со скоростью Uh отсчетов в секунду, равна 1/2 h Гц. Содержащиеся в исходном сигнале более высокие часто­ ты будут свернуты в диапазон 0—1/2 /г Гц и смешаются с более низкими частотами этого диапазона (рис. 7.7). Граничная частота.

ь Ч г ( ? - 2 >

называется частотой Найквиста, или частотой свертывания (на­ ложения). Свертывание составляющих исходного процесса в диа­ пазоне частот [О,/с] представляет собой нечто вроде гофрирования частот, как показано на рис. 7.8. В частности, для любой частоты. /, принадлежащей диапазону 0 < / < / с, замаскированными под. частоту f являются высокие частоты

(2fe ± /)

(4/, ± /),..■£ (2л/, ±7)........

(7-3)

Для доказательства этого факта заметим, что при t =

1/2 fe

cos 2k/ / =

cos 2« (2n f ± f) -Jr- =

cos

(7.4).

 

* f c

I C

- ^

l> В отечественной литературе подобный интервал, выбираемый в соот­ ветствии с теоремой Котельникова, называют периодом опросов, или интерва­ лом выборок.— Прим. ред.

Общие соображения о сборе и обработке данных

263

Таким образом, при снятии отсчетов с интервалами 1/2 fc значе­ ния косинуса для частот / и 2nfc ± / одинаковы. Например, если /с = 100 Гд, то составляющие с частотой 30 Гц будут неотличимы

Р и с . 7.7. Пример маскировки частот.

от составляющих с частотами 170, 230, 370, 430 Гц и т. д. Анало­ гичным образом энергия колебаний с такими высокими маскирую­ щими частотами будет неотличима от энергии колебаний с более

fc

2 fc

3fe

4 fc

6 fr

Исходные частоты

Р.'и с. 7.8. Свертывание частот относительно граничной частоты Най­ квиста.

Низкими частотами. Это объясняется тем, что при t = 1/2 fc ха­ рактеризующие энергию величины sina2n// и соі22л// одинаковы для частот / и 2/іД. + /. Следовательно, если граничная частота/, такая, как показано на рис. 7.9, то истинная спектральная плот­ ность (рис. 7.9, а) будет искажена (рис. 7.9, б).

264

Глава 7

Для исключения влияния маскировки на практике исполь-/* зуются два метода. Один из них состоит в выборе настолько малой величины /г, чтобы составляющие с частотами выше /с в исследуе­ мом процессе были физически невозможны. Предположим, что для некоторого процесса желательно получить информацию, от­ носящуюся к частотам ниже 1000 Гц. Это означает, что технически достаточно выбрать величину h ^ 0,50 мс. Предположим, однако, что есть основания полагать, что исследуемый процесс может со­ держать составляющие с частотами до 2000 Гц. Тогда следует при­ нять частоту Найквиста равной 2000 Гц, и, чтобы избежать маски­ ровки, потребуется величина h <1 0,25 мс. В общем можно рекомен-

0(f)

\ А

л

ОГс

Ри с . 7.9. Маскированный за счет свертывания частот энергетический

спектр.

а — истинный спектр; б — маскированный спектр.

довать выбирать частоту fc в 1,5-—2 раза больше максимальной ожидаемой частоты в спектре исследуемого сигнала. При нсполь-* зовании второго метода процесс до снятия отсчетов фильтруется, в результате чего в отфильтрованном процессе составляющие с частотами выше наиболее высокой желаемой граничной частоты содержаться не будут. В этом случае выбор частоты /с, равной мак­ симальной исследуемой частоте, даст точные результаты для ча­ стот ниже /е. Второй метод предпочтительнее первого, так как он позволяет сэкономить машинное время и снизить стоимость рас­ четов.

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

Общие соображения о сборе и обработке данных

265

(большинстве современных аналого-цифровых преобразователях сигнал на выходе дается в двоичном коде, что позволяет непосред­ ственно вводить его в ЦВМ, число уровней квантования также удобно представлять в двоичном коде. Типичные преобразователи

.аналогового сигнала в цифровой дают от 6 до 16 бит информации, что соответствует диапазону от 64 до 65 536 уровней квантования. При идеальном преобразовании ошибка квантования распределена равномерно со стандартным отклонением ~0,29 Ах, где Ах —

.шаг квантования. В этом легко убедиться.

 

отсчета

 

 

 

 

 

 

 

Р и с .

7.10.

Ошибка квантования.

 

Пусть

р(х) — плотность

распределения

ошибки квантования

1

Р М =

1 при —0,50 < х <

0,50,

 

0 при других X.

 

 

 

 

 

 

 

Очевидно,

среднее значение

ошибки равно

нулю,

поскольку

.плотность

р(х) симметрична

относительно

х = 0.

Дисперсия

ошибки

 

 

 

 

 

 

 

 

со

 

 

0,5

 

 

 

С І= J (х—jxj3 р (х) d x= j’ х М х = ~ .

— оо

—0*.5

Стандартное отклонение

ох= Ѵ 1/12 «= 0,29 цены деления.

(7.5)

Это и естьгсреднеквадратичное значение ошибки квантования1*, которое можно рассматривать как среднеквадратичное значение аомехи, наложенной на полезный сигнал. Предположим, напри­ мер, что весь диапазон значений сигнала разбит на 256 уровней.

Вэтом случае отношение пикового значения сигнала к среднеквадра-

ОЭту ошибку принято называть шумом квантования.— Прим. ред.

266

Глава 7

 

тичному значению

помехи составит (256Дх)/(0,29Дх)

1000, или

около 60 дБ.

 

 

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

цесса занимал возможно большую

часть шкалы квантования.

В противном случае разрешающая

способность окажется малой

и ошибка квантования может стать

существенной.'

Практические соображения. Аналого-цифровые преобра­ зователи (АЦП) промышленного производства обычно представ­ ляют собой двоичные либо двоично-десятичные (ДД) системы. Двоичная система преобразует непрерывные изменения электри­ ческого напряжения в двоичный цифровой код, а ДД-система — в цифровой код, который может быть представлен десятью раз­ личными цифрами. Конструкция двоичной системы сравнительна проще, но для обработки данных приходится составлять, как правило, программы в машинном коде. Двоично-десятичная систе­ ма сложнее и относительно менее надежна, но она позволяет про­ изводить непосредственно обработку данных наблюдений по программам, написанным на обычном алгоритмическом языке, например на ФОРТРАНе. Помимо уже рассмотренных ошибок дискретизации и квантования, необходимо указать еще и другие наиболее существенные ошибки, возникающие в системах АЦП.

Кним относятся:

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

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

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

Как правило, скорость преобразования в существующих АЦП составляет от 10ьдо 2,5-ІО5 отсчетов в секунду (с). Однако на практике ограничения, налагаемые цифровыми магнитными устройствами и требуемой разрешающей способностью, приводят

к уменьшению скорости квантования примерно до 5 -104 отсчетов

всекунду. Для определения скорости преобразования в отдель­ ном канале многоканального устройства эту сниженную скорость следует еще разделить ьа число опрашиваемых каналов. Так, n^gi одновременном опросе 10 каналов скорость преобразования в от­ дельном канале составляет, как правило, <~5-103 отсчетов в се­ кунду. Понятно, конечно, что эта скорость дана в реальных единицах времени. Относительная же скорость может быть па

Общие соображения о сборе и обработке данных

267

желанию увеличена или уменьшена за счет уменьшения или уве­ личения скорости подачи данных в систему АЦП.

Необходимо отметить еще одно обстоятельство, касающееся аналого-цифровых преобразователей. При одновременном преоб­ разовании в нескольких каналах наличие интервала времени между опросами в различных каналах1*приводит к появлению межканальной временной ошибки в отсчетах, которая для быстроменяющихся сигналов может оказаться весьма заметной. Ее появ­ ления можно избежать, используя в системах АЦП запоминаю­ щие схемы, в которые данные из всех каналов поступают одновре­ менно и хранятся там до момента опроса соответствующего канала коммутатором.

7.3.2. Предварительная обработка данных

После того как данные наблюдений представлены в цифровой ■форме, необходимо до проведения окончательного их анализа выполнить еще ряд операций. Первая из них — изменение форма­ та. Теоретически эта часть обработки вполне тривиальна, однако ■ее выполнение на практике часто требует значительных усилий. Устройства ввода современных ЦВМ приспособлены для приема данных, записанных только в определенном формате, а именно в том, в каком сама ЦВМ производит вывод данных. Системы преоб­ разования аналогового сигнала в цифровой зачастую не обеспечи­ вают вывода данных в стандартном формате, в результате чего ^оказывается необходимым составление программ в машинном коде.

Следующая операция состоит в преобразовании записанных в дискретном виде данных в реальные физические единицы. Как правило, зависимость между физическими единицами и единица­ ми, в которых представляют данные системы АЦП, линейна. Здесь наиболее распространен, следующий прием: одновременно ■с дискретизацией данных наблюдений производится дискретиза­ ция отсчетного, или калибровочного, сигнала. В этом случае можно непосредственно определить соотношение между отсчетными и физическими единицами измерения.

Наконец, последняя стадия — предварительное оценивание (улучшение) данных. К ней относятся обнаружение и исключе­ ние выбросов, тренда, смещения нуля и других возможных типов ■ошибок. Процедура обнаружения выбросов обычно включает вычисление первой разности. Вычисление скользящих средних позволяет обычно обнаружить смещение нуля. Проблема исклю- ■І^ния тренда рассмотрена ниже в подразд. 9.1.3.

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

Соседние файлы в папке книги из ГПНТБ