Научная статья на тему 'Базовый алгоритм синтеза вычислительных структур реального времени для теплоэнергетических измерений'

Базовый алгоритм синтеза вычислительных структур реального времени для теплоэнергетических измерений Текст научной статьи по специальности «Компьютерные и информационные науки»

CC BY
66
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СИНТЕЗ ВЫЧИСЛИТЕЛЬНЫХ СТРУКТУР / ТЕПЛОЭНЕРГЕТИЧЕСКИЕ ИЗМЕРЕНИЯ SYNTHESIS COMPUTING STRUCTURES / HEAT POWER MEASUREMENTS

Аннотация научной статьи по компьютерным и информационным наукам, автор научной работы — Акунович Станислав Иванович, Жиляк Надежда Александровна

В статье идет речь о внедрении унифицированного базового алгоритма синтеза вычислительных структур реального времени в прибор безбумажный самописец-регистратор РМ-2201 электронный самописец, который предназначен для измерения по 6, 9, 12 или 15-ти гальванически изолированным каналам силы и напряжения постоянного тока, а также неэлектрических величин, преобразованных в электрические сигналы постоянного тока, и сопротивления постоянному току, для обработки результатов измерений, их визуализации, сбора, хранения и сигнализации аварийных ситуаций. Одно из основных и наукоемких звеньев рассматриваемого прибора это разработка программно-аппаратного комплекса прибора. Программная часть регистратора РМ-2201 включает в себя один из алгоритмов, позволяющих реализовать наиболее оптимальное распределение машинных ресурсов. В свою очередь программное обеспечение, созданное на основе предложенного алгоритма, обеспечивает обработку и преобразование термо ЭДС в температуру, что ускоряет процесс производства и освобождает другие не менее важные машинные ресурсы прибора.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по компьютерным и информационным наукам , автор научной работы — Акунович Станислав Иванович, Жиляк Надежда Александровна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

BASIC ALGORITHM OF SYNTHESIS COMPUTING STRUCTURES OF REAL TIME FOR HEAT POWER MEASUREMENTS

Discussed are some aspects of implementation of a unified basic algorithm of synthesis of real-time computation structures in the RM -2201 electronic recorder designed for DC current and voltage measurements via 6, 9, 12 or 15 decoupled channels as well as for measurements of non-electrical values converted to electric signals DC current and resistance. Measurement results are further processed, displayed, collected and stored; emergency situations are announced. Development of hardware and software was one of the main tasks in designing the instrument. Basically, software for the RM -2201 recorder comprises one of the algorithms providing for computer resources share. A novel algorithm and software utilizing this algorithm provide for thermo EMF-to temperature processing and conversion. All that resulted in saving labor and instrument computing resources.

Текст научной работы на тему «Базовый алгоритм синтеза вычислительных структур реального времени для теплоэнергетических измерений»

УДК 519.6

С. И. Акунович, канд. техн. наук; Н. А. Жиляк, аспирант

БАЗОВЫЙ АЛГОРИТМ СИНТЕЗА ВЫЧИСЛИТЕЛЬНЫХ СТРУКТУР

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

Discussed are some aspects of implementation of a unified basic algorithm of synthesis of real-time computation structures in the PM-2201 electronic recorder designed for DC current and voltage measurements via 6, 9, 12 or 15 decoupled channels as well as for measurements of non-electric values converted to electric signals of DC current and resistance. Measurement results are further processed, displayed, collected and stored; emergency situations are announced.

Development of hardware and software was one of the main tasks in designing the instrument. Basically, software for the PM-2201 recorder comprises one of the algorithms providing for computer resources share. A novel algorithm and the software utilizing this algorithm provide for thermo EMF-to-temperature processing and conversion. All that resulted in saving labor and instrument computing resources.

Введение. Республика Беларусь обладает развитой структурой машиностроительных предприятий, предприятий стройиндустрии, теплоэнергетического и нефтехимического профилей, предприятий по переработке и производству сельскохозяйственной продукции и фармацевтической индустрии. Для этих отраслей основной формой эксплуатационного контроля является измерение, обработка, визуализация и активный допусковый контроль длительных реализаций технологических процессов. В настоящее время на большинстве предприятий Республики Беларусь и странах СНГ применяются бумажные электромеханические самописцы, находящиеся в эксплуатации с начала 80-х гг. К существенным недостаткам бумажных регистраторов можно отнести большую трудоемкость проведения ретроспективного анализа информации с бумажной ленты, отсутствие возможности оперативного накопления данных в долговременном архиве на машинных носителях информации; сложность проведения оперативного синхронного анализа данных по нескольким каналам с привязкой ко времени и дате съема информации.

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

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

Безбумажный самописец-регистратор РМ-2201 - электронный самописец, осуществляющий сбор, обработку и архивирова-

ние результатов измерения, разработан в ОАО «МНИПИ» с использованием новых технологий и представляет собой современный прибор, не уступающий по своим техническим характеристикам и функциональным возможностям приборам аналогичного класса таких известных фирм, как YOKOGAWA, JUMO и SIEMENS.

Регистратор измерительный многоканальный РМ-2201 предназначен для измерения по 6, 9, 12 или 15-ти гальванически изолированным каналам силы и напряжения постоянного тока, а также неэлектрических величин, преобразованных в электрические сигналы постоянного тока, и сопротивления постоянному току, для обработки результатов измерений, их визуализации, сбора, хранения и сигнализации аварийных ситуаций [1].

Отображение измерений производится по группам каналов. Количество каналов в группе - от 1 до 6. Скорость обновления отображения в каждой группе выбирается пользователем из списка: 1, 2, 3, 5, 10, 15, 30 с, 1, 2, 3, 6 мин.

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

Настройка измерительных (аналоговых) каналов включает в себя:

- набор идентификатора канала (на виртуальной клавиатуре);

- назначение вида сигнала и диапазона;

- назначение уставок;

- выбор шкалы для графиков и гистограмм;

- выбор цвета отображения канала;

- выбор ширины линии для графика.

Для каждого канала (аналогового или математического) можно назначить до четырех

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

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

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

Если значение гистерезиса равно нулю, событие «Вход в зону нормы (из максимума)» регистрируется, когда измеренное значение станет меньше, чем значение уставки.

Если значение гистерезиса больше нуля, то событие «Вход в зону нормы (из максимума)» регистрируется только тогда, когда измеренное значение станет меньше, чем значение уставки минус абсолютное значение гистерезиса.

2. Сигнал ниже уставки. Аналогично пункту 1 назначается значение уставки, время задержки Т и гистерезис.

Если время задержки равно нулю, событие «Превышение уставки по минимуму» регистрируется, когда измеренное значение станет равным или ниже значения уставки.

Если время задержки больше нуля, то событие «Превышение уставки по минимуму» регистрируется через время Т после момента, когда измеренное значение станет равным или ниже значения уставки, и если измерения остаются ниже значения уставки в течение этого периода времени.

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

Если значение гистерезиса равно нулю, событие «Вход в зону нормы (из минимума)» регистрируется, когда измеренное значение станет выше, чем значение уставки.

Если значение гистерезиса больше нуля, то регистрируется событие «Вход в зону нормы (из минимума)» только тогда, когда измеренное значение станет больше, чем значение уставки плюс абсолютное значение гистерезиса.

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

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

Прибор регистрирует следующие события, которые отображаются в таблице обзора сигнализаций (журнал событий), записываются в файл (если установлен соответствующий ре-

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

Многофункциональность и универсальность прибора достигается, в частности, за счет использования компьютера с тактовой частотой 500-800 МГц (на этих тактовых частотах в настоящее время не требуется принудительного охлаждения процессора, что очень важно для промышленного использования). Также применяются модули с универсальными измерительными входами, позволяющими производить измерение постоянного напряжения и тока, сопротивления постоянному току, что даст возможность подключать почти все используемые на практике первичные датчики, преобразующие неэлектрические сигналы в электрические, в частности термоЭДС в температуру от термоэлектрических преобразователей (ТП).

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

Программная часть регистратора РМ-2201 включает в себя один из алгоритмов, позволяющих реализовать наиболее оптимальное распределение машинных ресурсов. В свою очередь программное обеспечение, созданное на основе предложенного алгоритма, обеспечило обработку и преобразование термоЭДС в температуру, что ускорило процесс производства и освободило другие не менее важные машинные ресурсы прибора [2].

При отладке программного обеспечения была обеспечена работа аналоговых каналов, которые содержат настройки на виды измерений, например на измерение термоЭДС типа B, S, Я, Т, L, N E и т. д. Зависимости выдаваемого ими электрического сигнала от температуры описываются полиномом высокой степени, т. е. номинальной статической характеристикой (НСХ), аппроксимируемой полиномами 814-й степени.

Алгоритм обработки сигналов от ТП создан на основе теории синтеза вычислительных систем реального времени. Теория синтеза вычислительных систем реального времени позволила синтезировать вычислительную систему обработки сигналов от ТП в условиях жестких требований к временным параметрам на основе полиномов, заданных аналитическим выражением во временной области и алгоритмом их реализации, обеспечивающих высокие метрологические показатели регистратора РМ-2201 за счет распараллеливания и конвейеризации реализации его путей. Исходя из предложенной теории, вершинами исходного графа алгоритма представляются полиномы, описывающие зависимость термо-ЭДС от температуры в градусах Цельсия, которые требуются для измерения ЭДС при определенной постановке задачи. При этом с помощью

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

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

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

Y = F(X, Z, C), (1)

где Y, X, Z - множество соответственно выходных, входных и промежуточных переменных; C - множество констант.

Алгоритм решения задачи предполагает выполнение следующих процедур.

1. Формирование графа вычислительного алгоритма:

F(X, ..., С) ^ G (V, Е), (2)

где множество вершин V соответствуют множеству операций отдельных фрагментов алгоритма, множество Е-ребер - информационным связям. Граф G может быть описан матрицей инцидентности, которая имеет I вершин, J дуг и элементы aij.

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

Полным путем Lv графа алгоритма будем называть путь, связывающий вершину V о xn, хп е X с вершиной V:

V о Ут, Ут е Y,

(3)

где V - дуги графа вычислительного алгоритма; Ут - ряд входных переменных.

3. Назначение уровней временной иерархии вершинам графа вычислительного алгоритма.

Уровнем у временной иерархии подмножества вершин графа V = V) называется их приоритет, соответствующий требуемой скорости обработки данных каждой из вершин данного подмножества и определяемый следующим образом:

Vv^ о V : VI о А (Ф;) = А/(у)|у = 1, Г, Д^(у) = тт(Д^(ф с Ф)),

Г о ттД^ф^С, (4)

где Ф = XuYuZ о V; Ф(/) = иф,- о Д(у) - подмножества переменных модели (1) мощностью К.

Подмножество вершин Vv графа алгоритма V = (V, Е), для которого найдено соответствие

V(Y) о Д(у), (5)

называется множеством вершин у-го уровня временной иерархии.

Интервал времени Д^(у), для которого установлено соответствие

А/(у) о V(y) : G = G(V, Е), V(y) с V, (6)

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

4. Назначение функциональных устройств вершинам графа базовой структуры («нагруженный» граф), формирование векторов назначения вида

*=шх. (7)

1 ' /=1 /=1

векторов размером I, элементы которых г/и) = /; йц = 1,-/' = 1, ..., J, I = 1, ..., I, и = 1, ..., w, Я(-п)Ф Ф m)Vn Ф т; п, т = 1, ..., w называются множеством векторов назначения.

Если вершина vi принадлежит множеству вершин V(Y) о Д(у), то данной вершине может быть назначено функциональное устройство в виде термопары, для которого выполняется условие

Т' = (у) о АКТ), (8)

где Т - время обработки /-м устройством операции, отождествленной с вершиной vi.

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

5. Формирование вектора реализации (элемент ТI равен времени выполнения операции, отождествленной с /-й вершиной функционального устройства (в нашем случае - термопарой), назначенного в соответствии с вектором назначения).

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

7. Формирование усеченных путей Т-х уровней временной иерархии. J-м усеченным путем у-го уровня временной иерархии называется путь у-го уровня временной иерархии, для множества вершин V которого справедливо

/-1

V (ЛУ ) = VL (у) \ ^ (к,у). (9)

к=1

8. Определение конвейеризируемых путей и ступеней конвейера.

В ступень конвейер а должны входить ФУ, для которых суммарное время выполнения операций меньше или равно его циклу.

9. Определение множества свертываемых вершин (для неконвейеризированных путей).

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

11. Проверка реализуемости вычислительной структуры реального времени (выполнение условия реализуемости пути графа алгоритма в РМВ).

12. Формирование вектора временной развертки.

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

14. Анализ реализуемости вычислительной структуры (наличие отрицательных элементов в данном векторе указывает на нереализуемость ВС на выбранном векторе назначения).

15. Разработка графа вычислительной структуры.

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

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

Заключение. В результате реализации алгоритма, основанного на теории синтеза вычислительных систем реального времени, будут сформированы:

- граф вычислительной структуры, транспонируемый в структурную схему проектируемой системы, и спецификация его вершин;

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

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

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

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

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

Литература

1. Носенко, В. А. Специализированные информационно-измерительные системы / В. А. Носенко, А. В. Дзюба, А. И. Денис // Доклады Национальной академии наук Беларуси. Минск. -2005. - Т. 49, № 4.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2. Никифорова, З. С. Роль измерительно-информационных систем в повышении технического уровня и качества выпускаемой продукции / З. С. Никифорова, С. А. Попоудина, А. В. Дзюба. - Минск: БелНИИТИ, 1983. - 13 с. -(Обзорная информация).

3. Жиляк, Н. А. Базовый алгор итм синтеза вычислительных структур реального времени / Н. А. Жиляк, А. С. Кобайло // Труды БГТУ. Сер. VI, Физ.-мат. науки и информ. - 2007. -Вып. XV. - С. 147-150.

4. Шпаковский, Г. И. Архитектура параллельных ЭВМ / Г. И. Шпаковский. - Минск: Университетское, 1989. - С. 24-36.

УДК 681.325.3

Д. В. Шиман, ассистент

КОМПЬЮТЕРНАЯ ИМИТАЦИОННАЯ МОДЕЛЬ ПЕРЕДАЧИ ИНФОРМАЦИИ НА ОСНОВЕ ПЕРЕМЕЖИТЕЛЯ И ДВУМЕРНОГО ЛИНЕЙНОГО ИТЕРАТИВНОГО КОДА С ОБЪЕДИНЕННЫМИ ДВОЙНЫМИ ДИАГОНАЛЬНЫМИ ПРОВЕРКАМИ

The article describes the computer simulation model of the information transmissions on the basis of interleaving and two-dimensional linear iterative code with the united dual diagonal checks. This complex information transformation makes it possible to correct repeated errors with smaller expenditures of time and simplicity of the algorithm correction. In model is realized the opportunity of management by throughput of the channel, likelihood characteristics of the errors arising at data transmission, and also definition of amount of not corrected erroneous bit after decoding.

Введение. Как известно, при передаче и хранении информации все больше факторов влияют на ее качество, внося в последовательности данных различного рода ошибки. Для поиска и исправления таких ошибок используются разнообразные методы. Наиболее эффективными из таких методов являются те, которые базируются на применении помехоустойчивых корректирующих кодов, переме-жении символов исходной последовательности и степени группирования [1-3].

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

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

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

Компьютерная имитационная модель передачи информации реализована в среде Delphi 7. За основу исходных данных взят массив, заполненный двоичной информацией случайным образом. Кодер построен на основе ранее разработанного усеченного двумерного линейного итеративного кода с объединенными двойными диагональными проверками [4]. Кодирование информации осуществляется путем добавления проверочных бит к информационной последовательности. Формирование проверочных бит матрицы размерностью k1 строк на k2 столбцов по вертикали, горизонтали и двум диагоналям производится путем суммирования по модулю два всех информационных символов каждого из направления. По своим корректирующим способностям данный код позволяет исправить только одиночные и двойные ошибки. Декодер построен на основе синдромного декодирования, являющегося одним из наиболее распространенных методов коррекции ошибок блочными кодами. Данный декодер считается декодером максимального правдоподобия при декодировании жестких решений.

Канал

а

б

Рис. 1. Схемы передачи информации: а - на основе кодера и перемежителя; б - на основе кодера

На основе ранее разработанного псевдослучайного метода перемежения построены блоки перемежитель и деперемежитель [5]. Процедура перемежения состоит в перестановке символов кодированной последовательности до ее модуляции и восстановлении исходной последовательности после демодуляции. Перестановка позволяет так разнести рядом стоящие символы, чтобы они оказались разделены группой других символов, передаваемых в том же блоке данных. Данная операция не вносит избыточности, а только изменяет порядок следования символов или бит. В разработанном методе пе-ремежение (деперемежение) основано на использовании матрицы, строки и столбцы которой формируются последовательностью входных символов. Для удобства реализации метода матричные расчеты преобразованы в линейные:

(( с ь. \ У1 ^

adr.

int

n2 -1

1

mod n1

-1

(((к mod(n2 -1) +1)

х (n2 -1) = v,

mod n 2) -

)-О"

где int - целочисленное деление; к - адрес памяти (0 < к < (N -1)); n2 - длина столбцов; n1 -длина строк; d1, d2 - степенные коэффициенты. В приведенном соотношении учитывается тот факт, что считывание символов осуществляется по столбцам. Данная формальная запись означает, что v-й бит входной последовательности будет считан на выходе под к-м номером (adrk), что позволяет определить местоположение бит в выходной последовательности без перемежения строк и столбцов. Данная запись используется и при деперемежении с небольшой модификацией (заменой степенных коэффициентов d1, d2 на e1, e2). Реализация определения местоположения бита выполнено следующей функцией:

function TPsevdo.Pos(const i : integer; const R : TRec) : integer; var

PosA, PosB : Integer;

S : AnsiString; begin

PosA := (Power((i div (R.n)) + 1, R.e, (R.nr + 1)) mod (R.nr + 1)) - 1; // определение строкового местоположения бит после перемежения строк PosB := (Power((i mod (R.n)) + 1, R.d, (R.n + 1)) mod (R.n + 1)); // определение столбчатого местоположения бит после перемежения столбцов

PosB := (PosB - 1) * R.nr;

Result := PosA + PosB; //определение позиции выходного бита end;

Power - функция, реализующая возведение в степень типа xy:

function TPsevdo.Power(const A, B : Integer; const AMod : Integer = 0) : Integer;

var i : integer; begin Result := 1; if B > 0 then begin

for i := 1 to B do begin

Result := Result * A; if AMod <> 0 then Result := Result mod AMod; end; end; end;

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

Перед началом моделирования необходимо определить начальные параметры модели. В окне Входные данные нужно осуществить выбор пропускной способности моделируемого канала передачи данных из наиболее распространенных на практике. Одними из таких каналов являются волоконно-оптические каналы, характеризующиеся малым ослаблением сигнала и высокой степенью защищенности от внешних электромагнитных помех [6]. Далее следует определить объем передаваемой информации в мегабайтах, выбрать длину информационного слова для усеченного двумерного линейного итеративного кода с объединенными диа-тональными проверками в битах. В окне Зашумленность канала необходимо указать вероятность появления ошибки для каждого вида ошибок (одиночных, двойных, тройных и т. д.).

В начале моделирования, исходя из значения длины информационного слова кода, определяется оптимальное количество строк и столбцов проверочной матрицы с точки зрения избыточности. Оптимальной матрицей является квадратная матрица, в которой к1 строк равно к2 столбцов. Данное утверждение доказано в работе [1]. В случае, если невозможно сформировать квадратную проверочную матрицу, то необходимо воспользоваться матрицей, близкой к квадратной. Данное решение реализовано следующим образом:

S := $FFFFFFFF; while i >= 1 do begin

if fInputData.edBlockSize.Value mod i = 0 then if fInputData.edBlockSize.Value * 3 + i + i * fInputData.edBlockSize.Value < St hen // определение блока с минимальным количеством проверочных символов

begin k1 := i;

k2 := fInputData.edBlockSize.Value div i; S := k2 * 3 + k1 + k1 * k2; Data_and_Likdd_k1 := k1; Data_and_Likdd_k2 := k2; end; // определение значений k1 строк и k2 столбцов Dec(i);

fProcessInfo.pbProcess.Progress := fProcessInfo.pbProcess.Progress + 1; end;

Рис. 2. Интерфейс компьютерной имитационной модели

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

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

Процесс моделирования можно контролировать путем нажатия на кнопки Начать моделирование, Пауза/продолжить, Остановить моделирование. В процессе моделирования выводится в окно Моделирование информация о количестве переданной информации, времени, затраченном на ее передачу, количестве переданных проверочных символов. Количество и кратность ошибок, имитированных во время

эксперимента, отображается в таблице окна Процесс моделирования. Также в этой таблице показывается количество неисправленных ошибок и их кратность для каждой схемы моделирования (после likdd - на основе кодера и без перемежителя, после перемеж - на основе кодера и перемежителя). В нижней части окна ведется история выбранных начальных параметров модели, имитированных и неисправленных ошибок, их кратность и позиция. Данная информация записывается в файл.

Рис. 3. Окно Размер блока перемежителя

В ходе проведения экспериментов выяснилось, что не всегда выбранная длина информационного слова кода и параметры процесса перемежения способны корректировать ошибки, возникающие в процессе передачи информации. Так, при средней вероятности появления ошибки, равной 10-5, на 1 Мбайт передаваемой информации с размером блока перемежения 32x32, пропускной способностью 1 кбит/с, длиной информационного слова 100 моделью имитируется в среднем 300 ошибок. Из общего количества ошибок в среднем 90 не исправляется при использовании схемы передачи данных без перемежителя и 20 ошибок - с переме-жителем. Изменяя в модели только размер блока перемежения на 32x320 и пропускную способность канала на 10 кбит/с, получим уменьшение неисправленных ошибок в схеме с использованием перемежителя в среднем до 1 при неизменном объеме передаваемой информации. Другим способом сокращения неисправленных ошибок является уменьшение длины информационного слова усеченного двумерного линейного итеративного кода с объединенными диагональными проверками, но при этом увеличивается объем передаваемой информации. Данный факт объясняется не совсем оптимальным выбором начальных параметров модели.

Заключение. Разработанная компьютерная модель передачи информации показала, что выбранный вариант комплексного преобразования передаваемой информации на основе пе-ремежителя и двумерного линейного итератив-

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

Литература

1. Урбанович, П. П. Избыточность в полупроводниковых интегральных микросхемах памяти / П. П. Урбанович, В. Ф. Алексеев, Е. А. Верниковский. - Минск: Наука и техника, 1995. - 263 с.

2. Мак-Вильямс, Ф. Теория кодов, исправляющих ошибки / Ф. Мак-Вильямс, Н. Слоэн; пер . с англ.; под ред. Л. А. Басалыго. - М.: Связь, 1979. - 746 с.

3. Скляр, Б. Цифровая связь. Теоретические основы и практическое применение / Б. Скляр; пер. с англ. - 2-е изд. - М.: Издательский дом «Вильямс», 2003. - 1104 с.

4. Шиман, Д. В. Метод псевдоравномерно-го перемежения символов турбокода / Д. В. Шиман // Труды БГТУ. Сер VI, Физ.-мат. науки и информ. - 2005. - Вып. XIII. - С. 168-170.

5. Шиман, Д. В. Свойства и параметры линейных итеративных кодов с двойными диагональными проверками / Д. В. Шиман, Д. М. Романенко // Труды БГТУ. Сер. VI, Физ.-мат. науки и информ. - 2007. - Вып. XV. -С.151-154.

6. Гуров, И. П. Основы теории информации и передачи сигналов / И. П. Гуров. - СПб.: BHV, 2000. - 97 с.

i Надоели баннеры? Вы всегда можете отключить рекламу.