Патент на изобретение №2165676
|
||||||||||||||||||||||||||
(54) СПОСОБ И УСТРОЙСТВО ДЛЯ ОПРЕДЕЛЕНИЯ ЗАДЕРЖКИ ДИСКРЕТИЗАЦИИ В ВЫБОРКАХ ПОЛУЧЕННОГО СИГНАЛА В СИСТЕМЕ СОТОВОЙ СВЯЗИ
(57) Реферат: Сущность изобретения: восстановление тактовой синхронизации с или без цикловой синхронизации в системе сотовой связи с многопостанционным доступом и временным разделением каналов (TDMA) осуществляется способом, включающим операции оценки косвенных переменных, которые включают информацию о задержке дискретизации линейного комплексного вектора, используя критерий максимального правдоподобия и извлечение задержки дискретизации из оцененных косвенных переменных. Для извлечения задержки дискретизации используется оптимальная процедура и три субоптимальных, которые проще в смысле вычислений. Эти процедуры раскрыты в описании изобретения вместе с описанием устройства для его физической реализации. Для цикловой синхронизации задержка дискретизации при восстановлении синхронизации дополняется рядом интервалов выборок, определенных максимальной амплитудой сигнала множества выборок, вычисленных из оцененных косвенных переменных, и задержкой дискретизации. Техническим результатом является обеспечение улучшенной цикловой синхронизации и восстановление синхронизации в системе связи. 3 с. и 16 з.п. ф-лы, 6 ил., 2 табл. Изобретение относится к восстановлению тактовой синхронизации и цикловой синхронизации в системах связи. Изобретение применимо к любой системе связи, имеющей известную последовательность сигналов или слов синхронизации и, главным образом, к описываемой ниже системе сотовой связи с многопостанционным доступом и временным разделением каналов (TDMA), которая совместима с системой EIA/TIA, описываемой в документе IS-54-B: “Стандарт совместимости при работе мобильной станции с базовой станцией в двойном режиме”. Для удобства и краткости изложения, такого рода система называется ниже просто как система IS-54. В такой системе данные передаются во временных отрезках, каждый из которых содержит слово синхронизации из 14 знаков, за которым следует информационная последовательность. Известно, что в системах связи необходимо восстанавливать синхронизацию при уплотнении с временным разделением каналов (TDM) принятого цифрового сигнала данных с тем, чтобы выборки сигнала были получены при оптимальном временном режиме для дальнейшей обработки и извлечения передаваемой информации. Точно так же хорошо известно, что восстановление тактовой синхронизации, цикловая синхронизация и необходимая обработка выборок затруднена при низком отношении “сигнал-шум” (SNR) и то, что сотовые системы связи часто обладают низким отношением “сигнал-шум”. Обозначая интервал между знаками полученного сигнала и, следовательно, период между последовательными выборками, как T и обозначая задержку дискретизации, т.е. период между оптимальным и реальным временем выборки полученного сигнала как ![]() ![]() ![]() ![]() Согласно одному аспекту, предлагается способ определения задержки дискретизации в выборках принятого сигнала, содержащий следующие стадии: оценка с использованием максимального критерия правдоподобия и косвенных переменных линейного комплексного вектора, который аппроксимирует полученные выборки сигнала, причем косвенные переменные, содержат информацию о задержке дискретизации, и определение задержки дискретизации из оцененных косвенных переменных для восстановление тактовой синхронизации принятого сигнала сообщения. Предпочтительно, каждая составляющая линейного комплексного вектора содержит линейную комбинацию предопределенной пары функций, зависящих от задержки дискретизации, и пару функций, предпочтительно содержащих функции cos ![]() ![]() ![]() ![]() ![]() i) путем умножения оцененных косвенных переменных на действительные элементы предопределенных хранящихся в памяти матриц; ii) путем выполнения одиночного нелинейного преобразования оцененных косвенных переменных; iii) расчетом в соответствии с функцией (1/ ![]() iv) последовательностью операций, включающих первоначальный расчет оценки задержки дискретизации в соответствии с предопределенной функцией оцененной косвенной переменной и итеративно один или несколько раз: оценка комплексного коэффициента замирания выборок полученного сигнала в зависимости от оценки задержки дискретизации и оценки косвенных переменных и повторной оценки задержки дискретизации в зависимости от оцененного комплексного коэффициента замирания и оцененных косвенных переменных. Для того, чтобы обеспечить цикловую синхронизацию и восстановление тактовой синхронизации, указанный способ должен также включать следующие стадии: добавку к определенной задержке дискретизации для восстановления тактовой синхронизации полученного сигнала сообщения, задержку дискретизации, включающую целое число i интервалов дискретизации T для цикловой синхронизации принятого сигнала связи и определение числа i с помощью следующих операций: определение, в зависимости от оцененных косвенных переменных и определенной задержки дискретизации для восстановления синхронизации, амплитуд сигнала множества полученных выборок и выбор числа i, соответствующего максимальной амплитуде выборки. Операция определения амплитуд сигналов может включать итеративное определение амплитуд сигнала в зависимости от итеративного определения задержки дискретизации для восстановления тактовой синхронизации. Другим аспектом предлагаемого изобретения является способ цикловой синхронизации и восстановления тактовой синхронизации в выборках принятого сигнала связи, включающий следующие стадии: используя максимальный критерий правдоподобия, определяются косвенные переменные линейного комплексного вектора, который аппроксимирует выборки принятого сигнала, причем косвенные переменные включают информацию о задержке дискретизации для каждого из множества принятых выборок сигналов; определение первой задержки дискретизации из оцененных косвенных переменных, причем первая задержка дискретизации по времени меньше интервала дискретизации; определение, в зависимости от оцененных косвенных переменных и в соответствии с первой задержкой дискретизации, амплитуды сигнала в каждой из множества полученных выборок сигналов; идентификация индекса максимума найденных амплитуд сигналов относительно текущей выборки; сложение первой задержки дискретизации для выборки, идентифицируемой указанным индексом, с произведением интервала дискретизации, умноженного на индекс для получения итоговой задержки цикловой синхронизации и восстановления тактовой синхронизации. Стадии определения первой задержки дискретизации и амплитуд сигнала могут быть выполнены итеративно. Изобретение также предусматривает устройство для определения задержки дискретизации для выборки принятого сигнала связи, содержащее: множество фильтров с конечной импульсной характеристикой, реагирующих на принятые выборки сигнала для множества косвенных переменных линейного комплексного вектора, который аппроксимирует полученные выборки сигнала в соответствии с максимальным критерием правдоподобия, и калькулятор, который реагирует на множество косвенных переменных для вычисления оцененной задержки дискретизации. Множество фильтров с конечной импульсной характеристикой и расчетный узел выполнены в виде, по меньшей мере, одного цифрового процессора обработки сигнала. Краткое описание чертежей Сущность изобретения будет понятна из следующего описания со ссылкой на соответствующие чертежи, на которых: Фиг. 1 представляет собой блок-схему узлов беспроводного цифрового приемника системы связи; Фиг. 2 представляет собой блок-схему устройства восстановления тактовой синхронизации в соответствии с настоящим изобретением; Фиг. 3, 4 и 5 представляют более подробную схему восстановления тактовой синхронизации в соответствии с одним вариантом изобретения; Фиг. 6 представляет устройство для восстановление тактовой и цикловой синхронизации в соответствии с другим вариантом изобретения. Подробное описание изобретения Ниже описываются модели сигналов и наблюдения для системы IS-54, за которым следует описание процессов, которые могут быть использованы в соответствии со способом изобретения. Далее, со ссылкой на прилагаемые чертежи, подробно описывается физическое осуществление восстановления тактовой и цикловой синхронизации. Хотя это описание относится конкретно к системе IS-54, нужно подчеркнуть, что это лишь один пример осуществления изобретения и что изобретение применимо к другим системам связи с известными последовательностями синхронизации. Модели сигналов и наблюдения В системе IS-54 используется дифференциальная квадратурная фазовая манипуляция (DQPSK) со сдвигом фазы ![]() ![]() где k представляет собой положительное целое число, отождествляющее знак sk, a wk является комплексной величиной, представляющей двоичную информацию с действительными и мнимыми частями ![]() ![]() ![]() ![]() ![]() ![]() ![]() где yi, представляет собой комплексную выборку, идентифицированную индексом i, который может быть целым числом от 1 до 2M+1, T представляет собой знак интервала, ![]() ![]() где ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Способ восстановления тактовой синхронизации Модель наблюдения, описываемая уравнением (2), аппроксимирована линейным вектором модели наблюдения, включающей косвенные переменные, которые содержат информацию о задержке дискретизации ![]() При [] T, представляющем сопряженную транспонированную матрицу в квадратных скобках, выразим (2M+1)-размерный вектор наблюдения выборок принятого сигнала через Yn=[y2n+1, y2n+2,…y2n+2M, y2n+2M+1]T (3) где (M+1)-размерный вектор известных знаков слова синхронизации является S=[s0, s1,…sM-1, sM]T (4) а (2M+1)-размерный вектор наблюдения выборок помех выражен как ![]() Обозначая матрицу импульсной характеристики выражениями ![]() мы можем выразить модель уравнения (2) в векторном виде Yn=G ![]() Для того чтобы линеаризовать эту модель каждый компонент g0( ![]() ![]() ![]() ![]() ![]() ![]() gij( ![]() ![]() ![]() ![]() ![]() где a1ij, a2ij, и a3ij – коэффициенты матрицы для матриц A1, A2 и A3 соответственно. После этого можно выполнить следующую аппроксимацию для матрицы G( ![]() G( ![]() ![]() ![]() ![]() ![]() в интервале от -T/2 до T/2 для задержки дискретизации ![]() Можно использовать несколько пар функций, например представленные ниже следующие пары уравнений от (9) до (13), для линеаризации и обеспечить достаточную точность аппроксимации. В паре уравнений (11) g( ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Хотя могут быть использованы любая эта или другая возможная пара функций, первая пара функции (9) обеспечивает хорошую точность и наименьшую вычислительную сложность и поэтому наиболее предпочтительна. Только эта пара функций далее рассматривается более подробно. Используя аппроксимацию в уравнении (8), модель наблюдения уравнения (7) становится Yn= (A1 ![]() ![]() ![]() ![]() или, соответственно, Yn= A1SUn ![]() ![]() Если Фn представляет собой переменную, которая представляет собой трехмерный комплексный вектор, образованный перестановкой трех косвенных переменных Ф1,n, Ф2,n и Ф3,n т.е. Фn= [Ф1,n, Ф2,n, Ф3,n]T, когда Ф1,n= Un ![]() ![]() Yn=BФn+Hn (16) Оценка косвенных переменных Если V ![]() ![]() ![]() ![]() суффикс ‘ указывает на сопряженную транспонированную матрицу. Представляя оцененные значения диакритическим знаком ^ , максимальная оценка правдоподобия ![]() ![]() C = (B ![]() ![]() ![]() ![]() ![]() ![]() ![]() Восстановление задержки дискретизации Для восстановления задержки дискретизации можно использовать несколько процедур, выбор которых зависит от компромисса между точностью оценки и сложностью вычислений. Ниже описывается оптимальная оценка и три субоптимальных, но вычислительно простых оценки. Оптимальная оценка Из уравнения (17) можно получить новое уравнение наблюдения для косвенных переменных ![]() где Гn представляет собой трехмерный вектор сложных случайных Гауссовых переменных с известной корреляционной матрицей V ![]() ![]() ![]() ![]() где F( ![]() ![]() ![]() ![]() ![]() ![]() и усреднена по комплексной переменной Un, чтобы определить следующее уравнение для функции правдоподобия по отношению к ![]() ![]() где ![]() Согласно максимальному критерию правдоподобия оптимальная оценка задержки дискретизации ![]() ![]() ![]() Выполнив следующую аппроксимацию: ![]() где Q1-Q5 являются матрицами расширения, имеющими действительные элементы, которые могут рассчитываться и храниться в справочной таблице в памяти. Затем, используя пару функций (9), можно вычислить оптимальную задержку дискретизации ![]() ![]() где ![]() ![]() Уравнение (25) содержит нелинейную максимизацию и по этой причине ее практическое осуществление может быть связано со сложными вычислениями. Следующие три субоптимальные альтернативные процедуры устраняют этот недостаток и легко выполнимы. Одноступенчатое нелинейное преобразование Оценка трех косвенных переменных может быть выражена в виде ![]() ![]() ![]() ![]() где функция f() зависит от пары функций ![]() ![]() ![]() Использование этого первого решения для задержки дискретизации описано ниже со ссылкой на фиг. 3 и требует 6М+5 комплексных умножений, 6М комплексных сложений, 1-го вещественного деления и 1-го не скалярного преобразования. Рекурсивная процедура Описываемая ниже альтернативная процедура рекурсивной оценки задержки дискретизации основана на максимизации функции правдоподобия по отношению к переменным ![]() ![]() ![]() ![]() при ![]() ![]() и ![]() где V ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Из уравнения (28) следует, что ![]() ![]() Во-вторых, из уравнения (21), функция правдоподобия для Un, предполагая, что задержка дискретизации ![]() ![]() где Vu= (F( ![]() ![]() ![]() ![]() ![]() отсюда следует, что оценка ![]() Затем рекурсивная процедура суммируется, как содержащая следующий следующие четыре последовательных шага для каждой итерации k: Шаг 1: первоначальная оценка задержки дискретизации в соответствии с уравнением (31): ![]() Шаг 2: оценка коэффициента замирания в соответствии с уравнением (32) ![]() Шаг 3: повторная оценка задержки дискретизации в соответствии с уравнением (33): ![]() Шаг 4: прекращение дальнейших вычислений, если достигнуто предопределенное максимальное число итераций, в противном случае – возврат к шагу 2. Осуществление этой рекурсивной процедуры для задержки дискретизации описано ниже со ссылкой на фиг. 5. Упрошенная процедура Простой вариант “вышеупомянутой рекурсивной процедуры заключается в остановке после шага 1, без выполнения операций 2-4. В этом случае задержка дискретизации оценивается с помощью уравнения: ![]() где Ci, представляет собой i-й элемент известной матрицы C. Осуществление этой простой процедуры для определения задержки дискретизации описано ниже со ссылкой на фиг. 4. Цикловая синхронизация Вышеописанные способы восстановления тактовой синхронизации применяются, когда задержка дискретизации ![]() ![]() ![]() ![]() ![]() Цикловая синхронизация основана на интервале наблюдения, соответствующего неопределенности интервалов знака L задержки дискретизации t, так чтобы это было применено к набору векторов наблюдения Y, каждый из который имеет вид уравнения (2), причем этот набор описывается как ![]() При использовании принципов максимального правдоподобия задача цикловой синхронизации и восстановления тактовой синхронизации сводится к определению максимума функции правдоподобия: ![]() которая может быть разделена на задачу цикловой синхронизации для определения максимума средней функции правдоподобия: ![]() и задачу восстановления синхронизации при определении максимума функции правдоподобия: ![]() где ![]() Вместо прямого вычисления функции правдоподобия в уравнении (37) можно использовать косвенный подход, как показано ниже. Для каждого значения индекса i в диапазоне от -(L-1)/2 до (L-1)/2 уравнения (7), (16), (18) и (20) могут быть выражены как: Yn+2i=G( ![]() ![]() Yn+2iBФ+Hn+2i; ![]() из которых может быть получена следующая функция правдоподобия: ![]() при ![]() Вместо усреднения по переменной ![]() ![]() ![]() где ![]() ![]() Применяя рекурсивную процедуру, описанную выше, для получения оценки ![]() ![]() где ![]() Игнорируя зависимость Vu от ![]() ![]() Таким образом, задача цикловой синхронизации сводится к оценке амплитуды сигнала для каждого значения индекса i и выбора такого значения индекса i, который соответствует максимальной амплитуде сигнала. Комбинированнмй способ цикловой синхронизации и восстановления тактовой синхронизации Сочетая этот способ цикловой синхронизации с описанной выше рекурсивной процедурой для восстановления тактовой синхронизации с описываемой ниже рекурсивной процедурой, которая предполагает использование упомянутой выше пары функций (9) для комбинации цикловой синхронизации и восстановления тактовой синхронизации, выполняются следующие шаги: Шаг 1: для каждого значения индекса i в диапазоне от – (L-1)/2 до (L-1)/2, выполняем следующие операции 1A-1C: 1A: Оценка вектора косвенной переменной в соответствии с уравнением (44): ![]() 1B: Первоначально (для счетчика k = 1) оценивают амплитуду сигнала и задержку дискретизации в соответствии с уравнениями (45) и (46): ![]() ![]() 1C: Итеративно (для k = 2 для максимального числа K итераций) оценивают амплитуду сигнала и задержку дискретизации в соответствии с уравнениями (47) и (48): ![]() ![]() Шаг 2: для каждого значения индекса i в диапазоне от – (L-1)/2 до (L-1)/2 определяют индекс i максимума оцененной амплитуды сигнала, определенной в шаге 1, и определяют общую задержку дискретизации из этого индекса и оцененную в шаге 1 задержку дискретизации в соответствии с уравнениями (49) и (50): ![]() ![]() Сложность этой процедуры зависит от числа итераций K и от количества L. Эта процедура значительно упрощается с небольшим снижением точности аналогично тому, как определяется только восстановление синхронизации, путем исключения шага рекурсии 1, так что шаг 1 включает только операции 1A и 1B. Уравнения (44) – (46), (49) и (50), используемые, как описано выше, с K = 1, и с парой функций (9) уравнение (45) становится уравнением: ![]() Осуществление этой комбинированной и упрощенной процедуры цикловой синхронизации и восстановления тактовой синхронизации описано ниже со ссылкой на фиг. 6. Физическое осуществление изобретения На фиг. 1 представлена блок-схема беспроводного цифрового приемника связи, в котором цифровой радиосигнал проходит через радиочастотную (РЧ) цепь 20 приемника на преобразователь 22 с понижением частоты, чтобы получить сигнал, который дискретизируется схемой выборки 24; выборки преобразуются в цифровую форму аналого-цифровым преобразователем 26. Оцифрованные выборки интерполируются интерполятором 28 в соответствии с восстановленной и оцененной задержкой дискретизации ![]() ![]() ![]() На фиг. 2 представлена блок-схема узлов блока 32, осуществляющих восстановление тактовой синхронизации для выполнения оцененной задержки дискретизации ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() На фиг. 4 представлена блок-схема узлов блока 32, осуществляющих восстановление тактовой синхронизации для получения оцененной задержки дискретизации ![]() ![]() ![]() ![]() На фиг. 5 представлена блок-схема устройства блока 32 для восстановления тактовой синхронизации и извлечения оцененной задержки дискретизации ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Результаты компьютерного моделирования различных описанных выше процессов изложены в табл. 1, приведенной ниже на примере отношения “сигнал-шум” порядка 8 дБ. Выбор слова синхронизации из 6 слов синхронизации, приведенных в IS-54, не оказал существенного влияния на эти результаты. При моделировании использовалась пара функций уравнения (9) в аппроксимации уравнения (8). Коэффициенты постоянных аппроксимирующих матриц Ai. (i = 1 до 3) уравнения (8) и постоянных матриц расширения Qi (i = 1 до 5) уравнения (24) были вычислены раздельно. Для каждой из описанных выше процедур в табл. 1, указано стандартное отклонение оцененной задержки дискретизации ( ![]() На фиг. 6 представлена блок-схема осуществления комбинированной операции цикловой синхронизации и функций восстановления синхронизации узлом 32 для получения оцененной задержка дискретизации ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Компьютерное моделирование, аналогичное описанному выше моделированию восстановления тактовой синхронизации, приведено ниже в табл. 2 на примере отношения “сигнал/шум” порядка 8 дБ. В табл. 2 сравниваются результаты для комбинированной цикловой синхронизации и восстановление тактовой синхронизации, как описано выше, используя пару функций (9) в уравнении аппроксимации (8), с результатами для обычных согласованных фильтров с обозначением для сдвига f нулевой несущей частоты и несущей частоты 300 Гц стандартного отклонения оцененной задержки дискретизации ( ![]() Формула изобретения
![]() ![]() ![]() ![]() ![]() ![]() 9. Способ по любому из пп.1 – 5, отличающийся тем, что определение задержки дискретизации из оцененных косвенных переменных включает стадии первоначального расчета оценки задержки дискретизации в соответствии с предопределенной функцией оцененных косвенных переменных, итеративно один или несколько раз: оценку комплексного коэффициента замирания полученных выборок сигнала в зависимости от оценки задержки дискретизации и оцененных косвенных переменных, повторную оценку задержки дискретизации в зависимости от оцененного комплексного коэффициента затухания и оцененных косвенных переменных. 10. Способ по любому из пп.1 – 5, отличающийся тем, что дополнительно осуществляют сложение найденной задержки дискретизации для восстановления тактовой синхронизации полученного сигнала связи с задержкой дискретизации, содержащей целое число i интервалов выборок для цикловой синхронизации полученного сигнала связи, определение числа i шагами, определение амплитуд сигнала множества полученных выборок сигнала в зависимости от оцененных косвенных переменных, выбор числа i, соответствующего максимальной амплитуде выборки. 11. Способ по п.10, отличающийся тем, что определение амплитуд сигнала включает итеративное определение амплитуд сигнала в зависимости от итеративного определения задержки дискретизации для восстановления тактовой синхронизации. 12. Способ цикловой синхронизации и восстановление тактовой синхронизации путем определения задержки дискретизации принятого сигнала связи, характеризующийся тем, что, используя максимальный критерий правдоподобия, производят оценку косвенных переменных линейного комплексного вектора, который аппроксимирует полученные выборки сигнала, причем косвенные переменные включают информацию о задержке дискретизации, и для каждой из множества полученных выборок сигнала: определяют первую задержку дискретизации из оцененных косвенных переменных, причем первая задержка дискретизации короче интервала выборки, определяют в зависимости от оцененных косвенных переменных и соответствующей первой задержки дискретизации амплитуды сигнала каждого из множества полученных выборок сигнала, идентифицируют индекс максимума определенных амплитуд сигнала относительно текущей выборки, добавляют первую задержку дискретизации для выборки, идентифицированной индексом произведения при умножении интервала выборки на указанный индекс для получения комбинированной задержки дискретизации для цикловой синхронизации и восстановления тактовой синхронизации. 13. Способ по п.12, отличающийся тем, что определение первой задержки дискретизации и амплитуды сигнала выполняют итеративно. 14. Устройство дискретизации полученного сигнала сообщения в системе сотовой связи, характеризующееся тем, что содержит множество фильтров с конечной импульсной характеристикой, реагирующих на полученные выборки сигналов для получения множества косвенных переменных линейного комплексного вектора, который аппроксимирует полученные выборки сигнала по максимальному критерию правдоподобия; калькулятор, реагирующий на множество косвенных переменных для вычисления оцененной задержки дискретизации. 15. Устройство по п.14, отличающееся тем, что калькулятор реализует умножитель для умножения множества косвенных переменных на действительные элементы предопределенных матриц из блока памяти. 16. Устройство по п.14, отличающееся тем, что калькулятор реализует функции для вычисления оцененной задержки дискретизации, как функции вещественных ![]() вещественных ![]() где ![]() ![]() ![]() 17. Устройство по п.14, отличающееся тем, что калькулятор реализует функции расчета оцененной задержки дискретизации как функции atan вещественной ![]() ![]() ![]() ![]() 18. Устройство по любому из пп.14 – 17, отличающееся тем, что калькулятор также реализует функции для расчета оцененного комплексного коэффициента замирания выборок полученного сигнала в зависимости от множества косвенных переменных и оцененной задержки дискретизации и для рекурсивного вычисления оцененной задержки дискретизации в зависимости от множества косвенных переменных и оцененного комплексного коэффициента замирания. 19. Устройство по любому из пп.14 – 18, отличающееся тем, что множество фильтров с конечной импульсной характеристикой и калькулятор выполнены в виде, по меньшей мере, одного цифрового процессора сигнала. Приоритет по пунктам: 23.08.95. по пп.1 – 9, 14 – 19; 23.10.95 по пп.10 – 13. РИСУНКИ
MM4A Досрочное прекращение действия патента Российской Федерации на изобретение из-за неуплаты в установленный срок пошлины за поддержание патента в силе
Дата прекращения действия патента: 27.06.2004
Извещение опубликовано: 20.04.2006 БИ: 11/2006
|
||||||||||||||||||||||||||