NTT-Gauss.HTM
ЧП Ферма     ▲ Выше     Примеры ТЧП Гаусса

Теоретико-числовое преобразование Гаусса


 









μ = p+jq
N = p2+q2
ρ( Zμ[j] ) = ZN
Она глядит на вас так нежно,
она лепечет так небрежно,
она так тонко весела,
её глаза так полны чувством,
вечор она с таким искусством
из-под накрытого стола
свою мне ножку подала!

Зачем я ею очарован?
Зачем расстаться должен с ней?
Когда б я не был избалован
цыганской жизнию моей,
когда б не смутное влеченье
чего-то жаждущей души,
я б здесь остался — наслажденье
вкушать в неведомой тиши;
забыл бы всех желаний трепет,
мечтою б целый мир назвал...
И всё бы слушал этот лепет,
всё б эти ножки целовал.
Александр Пушкин
сентябрь 1833.


Кольцо ЦКЧ Гаусса
ρ-изоморфизм Гаусса
Пример ПСНВ+ и ПСАНВ± по модулю ЦКЧ Гаусса μ=(5+j2)
Пример ПСНВ+ и ПСАНВ± по модулю ЦКЧ Гаусса μ=(2+j5)
Сравнение ПСНВ+ и ПСАНВ± по модулю μ=(5+j2) и μ=(2+j5)
Примитивные КОРНИ и модулярное ТЧП Гаусса
Выбор параметров кольца Гаусса
Примеры ТЧП Гаусса
Критика ТЧП Гаусса и ограничения его применения в ЦОС
Это не конец
  Волшебник Ферма
  Ферма опять побеждает
ДП Фурье квантованных сигналов в кольце ЦКЧ Гаусса



 
Всё это было бы так просто,
Когда бы не было так сложно.
                            Поговорка XX века

Кольцо ЦКЧ Гаусса.

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

Простейшее после   Z   кольцо, для которого была построена теория чисел, аналогичная обычной,   есть гауссово кольцо   Целых Комплексных Чисел (ЦКЧ)   Zμ[j] .   Гаусс в 1824 г. показал, что в этом кольце бесконечное число простых чисел и что всякое заданное число этого кольца представляется одним и только одним способом в виде произведения конечного числа этих простых и ещё одной из четырех единиц   (+1,+j,-1,-j)   кольца ЦКЧ   Zμ[j] .

Комплексное число будем называть простым комплексным числом, если оно не может быть представлено в виде произведения двух комплексных чисел, отличных от 1.   Комплексное число   γ=r1+jr2   будет кратно комплексному μ=p+jq   (т.е. μ | γ),   если частное   γ:μ   является ЦКЧ, а это будет тогда, когда

( r1·p + r2·q ) 0 mod N  
( r2·p r1·q ) 0 mod N ,
или иначе (что тоже самое)
r1·p + r2·q N = 0         
r2·p r1·q N = 0 ,       

где   N = p2+q2 = Norma(μ) , т.е. число элементов ПСНВ+ по модулю   μ   равно норме   N   комплексного модуля   μ .

Из этих уравнений можно найти зависимости между   r1 и r2

r1 = ПСАНВ± +r2 · q/p N  
r2 = ПСАНВ±r1 · q/p N .


 
Многие вещи нам не понятны не потому,
что наши понятия слабы;
но потому, что сии вещи
не входят в круг наших понятий.
                               Козьма Прутков
                        "Мысли и афоризмы",66.

ρ-изоморфизм Гаусса.

Теорема Гаусса. Всякое ЦКЧ   γ = a+jb   из кольца

Zμ[j] ,   μ=p+jq
при взаимно простых   p и q , т.е.   НбОД(p,q)=1 , сравнимо с одним и только с одним ВЫЧЕТом   h   из ряда   0,1,...,N-1 , т.е.

h = a+ρ·b N ,   N=Norma(μ)=p2+q2 ,   μ=p+jq ,

где коэффициент изоморфизма   ρ   определяется так

ρ = -p/q N = q/p N = (q-p)/(p+q) N .

Изоморфизм Гаусса ρ позволяет заменить выполнение операций над комплексными ВЫЧЕТами выполнением тех же операций над соответствующими им вещественными ВЫЧЕТами по целому модулю, равному норме комплексного модуля, т.е.

ρ ( Zμ[j] ) = ZN .
Обратный переход делается при помощи
ρ1( ZN ) = Zμ[j] ,
где
ρ1(h) = a+jb μ ,   μ = p+jq ,
т.е.
ρ1(h) = a+jb = (p·r1q·r2)/N + j·(p·r2+q·r1)/N ,
r1 = ПСАНВ± h·p N ,   r2 = ПСАНВ±h·q N ,

т.е.   r1 и r2   принадлежат ПСАНВ± по модулю N , при этом   h,   ρ,   p и q одновременно принадлежат либо   ПСНВ+ mod N   либо   ПСАНВ± mod N .

Выше полученные зависимости между   r1 и r2   примут вид

r1 = ПСАНВ± +r2 · ρ N  
r2 = ПСАНВ±r1 · ρ N .



Пример ПСНВ+ и ПСАНВ± по модулю ЦКЧ Гаусса μ=(5+j2)

Пример ПСНВ+ и ПСАНВ± по модулю ЦКЧ Гаусса μ=(2+j5)

Сравнение ПСНВ+ и ПСАНВ± по модулю μ=(5+j2) и μ=(2+j5)

 
Опять скажу — смотри в корень

                Козьма Прутков
          "Мысли и афоризмы",5.

Примитивные КОРНИ и модулярное ТЧП Гаусса.

Теорема. Если   ε = a+jb — ЦКЧ взаимно простое с ЦКЧ   μ , то

εN-1 μ = 1 ,
иными словами ЦКЧ   ε = a+jb   является КОРНЕМ из 1 порядка   (N-1)   по модулю ЦКЧ   μ .

Теорема. Если   ε = a+jb — ЦКЧ взаимно простое с ЦКЧ   μ , а   Tнаименьшее из чисел, при которых

εT μ = 1   (или, что тоже самое,   ε T1 mod μ) ,
то   T|(N-1)   и   T   делит всякий другой показатель K степени КОРНЯ из 1   ε , при котором
εK μ = 1 .

При наименьшем T таком, что T=φ(N), где φ — функция Эйлера, КОРЕНЬ ε будет примитивным (primitive) по модулю ЦКЧ   μ .

Теорема. Если   ε — примитивный КОРЕНЬ по модулю ЦКЧ μ и Norma(μ)=N, то члены ряда (из N1 значений)

1, ε1, ε2 ,..., εN2
будут попарно несравнимы между собой, т.е. этот ряд представляет полную приведенную (по модулю μ ) систему ВЫЧЕТов, если к нему присоединить нулевой элемент.

Пусть   μ=p+jq — ЦКЧ, где НбОД (p,q)=1, а   ε=a+jb — ЦКЧ взаимно простое с ЦКЧ   μ . Тогда ЦКЧ   ε   будет принадлежать мультипликативной группе (приведенной системе ВЫЧЕТов по модулю μ) кольца  Zμ[j] и иметь некоторый период   T   такой, что   T|φ(N) , где   N — норма ЦКЧ   μ .

Построим систему функций на группе   GT   целых чисел t:

χwGauss(t) = εwt ,   где   w,t=0,1,...,T-1 .     
Т.к.
χwGauss(t1+t2) = χwGauss(t1) · χwGauss(t2) ,
χwGauss(0) = 1
и, если в кольце   Zμ[j]   имеет смысл деление на число T (т.е. существует число 1/T), то введенные функции образуют множество характеров группы GT , а значит и ортогональный базис в пространстве   L(GT,Zμ[j]) . Построенный таким образом базис будем называть базисом Гаусса, а соответствующее ему преобразование — ТЧП Гаусса. Естественно, ограничимся такими КОРНЯми
ε T1 mod μ ,
для которых
ρ(ε) = root =2k, k=0,1, ... .

Тогда характеры группы GT в кольце Zμ[j] , преобразованные ρ-изоморфизмом, превращаются в (вещественный целочисленный) базис ЧП Ферма
ρ( χwGauss(t) ) = ρ( εwt ) = ( ρ(ε) )wt = rootwt = χwRader(t) .
Отсюда обратно
χwGauss(t) = εwt = ρ1( ρ(εwt) ) = ρ1( ( ρ(ε) )wt ) =
= ρ1( ( root )wt ) = ρ1( χwRader(t) ) ,
что устанавливает взаимно-однозначное соответствие между ТЧП Гаусса в пространстве   L(GT,Zμ[j])   и   ЧП Ферма   в пространстве   L(GT,ZN) .

NB. Значит гармонический анализ комплексных сигналов   γ[t]   в кольце ЦКЧ Гаусса можно перенести в кольцо целых чисел по модулю N, например, вычисляя спектр Рейдера   SRader[w]   сигнала   ρ(γ[t])   и применяя обратный изоморфизм Гаусса   ρ1(SRader[w])   получим спектр Гаусса   SGauss[w]   комплексного сигнала   γ[t]   (смотри схему) .

Вычисление спектра Гаусса через ЧП Ферма.

ЦКЧ   γ[t] ,

целые   x[t]=ρ( γ[t] )



                    T-1
SγGauss[w] = γ[t] εwt ( mod μ ) ,
                    t=0
             w = 0,1,...,T-1
                   T-1
SxRader[w] = x[t] rootwt (mod N)
                   t=0
            w = 0,1,...,T-1




SγGauss[w] = ρ1( SxRader[w] )


Тоже самое и для свёртки двух сигналов из ЦКЧ.

Вычисление свёртки двух комплексных векторов через ЧП Ферма.
ЦКЧ   γ[t] = ξ[t] * η[t]

целые   x[t]=ρ(ξ[t]), h[t]=ρ(η[t])



γ[t] = ρ1( ρ(γ[t]) )


SxRader[w] = ЧП Ферма ( x[t] )
ShRader[w] = ЧП Ферма ( h[t] )
ρ(γ[t]) = оЧП-Ф ( SxRader[w]·ShRader[w] )

NB. Это позволяет по новому взглянуть на машинную реализацию ДП Фурье комплексных сигналов, а именно, рассматривая их как сигналы в кольце ЦКЧ Гаусса и вычисляя их спектры в кольце целых чисел через ЧП Ферма (или ТЧП Рейдера, в которых умножения заменяются на сдвиги и сложения) с возвратом обратно в кольцо ЦКЧ Гаусса.


 
Нет пределов совершенству.
                          Поговорка

Выбор параметров кольца Гаусса.

Теорема. Если ε принадлежит кольцу ЦКЧ Zμ[j] , и Norma(ε) — простое целое число, то ε — простое число Гаусса.

Теорема. Норма простого числа Гаусса является либо простым числом, либо квадратом простого числа. Соответственно различают простые числа Гаусса I-го и II-го порядка.

Простые натуральные (целые положительные) числа, равные нормам простых ЦКЧ Гаусса, называются разложимыми. Простые натуральные числа, остающиеся простыми в кольце Zμ[j] , называются неразложимыми.

Теорема. Каждое простое число p=4s+1 является разложимым в кольце

Zμ[j] .

Теорема. ТЧП Гаусса по модулю μпримитивного числа Гаусса (т.е. НбОД (p,q)=1), существует тогда и только тогда, когда норма N — нечетное целое и имеет максимальный объём преобразования, выражающийся числом вида 4s.

NB. Наиболее широко используются кольца   Zμ[j]   такие, у которых ПСАНВ± по модулю   μ   геометрически ближе всего к квадрату со сторонами, параллельными осям координат. Это будет в том случае, когда   p и q   сильно отличаются друг от друга по абсолютной величине.

Например, при   p=1 и q=2b/2 ,   где   b=2n, n=1,2,...
норма ЦКЧ   μ = p+jq

N = p2 + q2 = 12 + (2b/2)2 = 1 + 2b = Fn ,
т.е.   N=Fn — число Ферма,   НбОД (p,q)=1 .



Примеры ТЧП Гаусса

 
  Бди!
             Козьма Прутков
      "Мысли и афоризмы",42.

Критика ТЧП Гаусса и ограничения его применения в ЦОС.

Естественно, что модульные операции накладывают ограничения на условия существования связи между спектром Гаусса   SGauss[w]   и Рейдера   SRader[w] .

Найдём условия, при которых использование   Zμ[j]-арифметики   (т.е. по модулю ЦКЧ Гаусса   μ=p+jq ) приводит к тем же результатам, что и вычисление свёртки в   Z[j]-арифметике   (обычных целых комплексных чисел). Очевидно, это возможно в том случае, когда результаты, вычисленные в кольце Гаусса   Zμ[j] , лежат внутри ПСАНВ± по модулю ЦКЧ Гаусса   μ .   Если отвести под мнимую и действительную части ЦКЧ Гаусса одинаковые "динамические" диапазоны, то значения комплексной свёртки   γ[t]   должны попасть внутрь квадрата P1P2P3P4 с центром в начале координат, у которого, например, вершина   P1   имеет координаты

P1(re) = P1(im)   =   ( N/(p+q) )/2     ( Fn-11 )/2 .
Это приводит к необходимости ограничить рабочие области квантованных комплексных сигналов и переходных функций (импульсных характеристик фильтров), т.е. входных данных комплексной свёртки, величиной
( N/(p+q) )/2     ( Fn-11 )/2 .

Поэтому для свёртки   γ[k]   комплексных сигналов   η[t] и ξ[t]   имеем




T-1
ξ[t] · η[ k-t T]
t=0



  ≤  
T-1
ξ[t] ║ · ║ η[ k-t T] ║
t=0
( N/(p+q) )/2 ,
где ║f(t)║ — норма функции   f(t) .   Если норму   f(t)   определить как
f(t) ║ = max |f(t)|        
и предположить, что
max |η(re)[t]| = max |η(im)[t]| = max |ξ(re)[t]| = max |ξ(im)[t]| = (ЦКЧ)A η ,
то наибольшее значение (амплитуда)   (ЦКЧ)A η   входных комплексных данных для свёртки, определится так
2T · ((ЦКЧ)A η)2   ≤   N/((p+q)) /2 ,                            
T · (ЦКЧ)A η    ≤   N/((p+q)) /2     ( Fn-21 )/2 .
Тоже самое и для результата комплексной свёртки, если предположить, что
max |γ(re)[k]| = max |γ(im)[k]| = (ЦКЧ)A γ ,
тогда
(ЦКЧ)A γ   ≤   ( N/(p+q) )/2     ( Fn-11 )/2 .

Пусть, например,   N=Fn=F5=232+1 ,   т.е.   p=1 , q=216 и T=28 . Тогда

T · (ЦКЧ)A η   ≤   ( (232+1) / (1+216) )1/2 /2  
( 232 / 216 )1/2 /2 =
= ( 216 )1/2 /2 = 28/2 ,
(ЦКЧ)A γ   ≤   ( (232+1) / (1+216) )/2  
( 216 )/2 = 216/2 .

Или иначе (точнее) выражение для динамического диапазона комплексных сигналов   η[t] и ξ[t]   при вычислении свёртки    γ[k]   в кольце ЦКЧ Гаусса можно записать так

max |γ[k]|   ≤   ( N/(p+q) )/2     ( Fn-11 )/2 .    
2T · max |η[t]| · max |ξ[t]|   ≤   ( N/(p+q) )/2     ( Fn-11 )/2 ,   
T · max |ξ[t]|   ≤   ( √ N/(p+q) )/2 ( Fn-21 )/2 .    

Точно такие же формулы для оценки динамического диапазона комплексного сигнала из ЦКЧ Гаусса будут и при вычислении ТЧП в кольце Гаусса

max |ЦКЧγ[k]|   ≤   ( N/(p+q) )/2     ( Fn-11 )/2 .   
2T · max |ЦКЧξ[t]|   ≤   ( N/(p+q) )/2     ( Fn-11 )/2 .   

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

T · max |ЦКЧξ(re)[t]|   ≤   ( N/(p+q) )/2     ( Fn-11 )/2 .   


 
Нет вещи столь великой, которую не
превзошла бы величиною ещё большая.
Нет вещи столь малой, в которую не
вместилась бы ещё меньшая.
                           Козьма Прутков
                     "Мысли и афоризмы",4.

Это не конец.


Волшебник Ферма.

NB. Однако, любознательный читатель мог заметить, что в формуле обратного преобразования-изоморфизма Гаусса   ρ1   присутствует деление на норму комплексного модуля   N=Norma(μ) . Казалось бы, это деление, как ложка дёгтя, портит всю бочку мёда с ТЧП Гаусса и ЧП Ферма в том смысле, что при их вычислении можно обойтись без умножений и делений. Но при выборе комплексного модуля из компонент числа Ферма   Fn-1 , т.е. в виде

μ = 1 + j·22n-1
так, что его норма становится равной следующему числу Ферма
N = Norma(μ) = 12 + (22n-1)2 = 1 + 22n = Fn ,

деление на   N=Fn   непостижимым и чудесным образом исчезает из формулы для   ρ1. При этом изоморфизм Гаусса равен числу Ферма   Fn-1 , уменьшенному на 1, т.е.
ρ = q = 22n-1 = Fn-1 1 .

О, здесь (и не только здесь!), чувствуется рука Великого гасконца и мага науки о числах Пьера де Ферма. В этом можно убедиться, заглянув в страну "Зазеркалье", куда повёл белый кролик в жилетке с часами любопытную девочку Алису. "И так, за мной мой читатель..." к следующим чудесам и диковинкам.


Ферма опять побеждает.

Выше были рассмотрены ограничения, накладываемые на "динамический" диапазон входных сигналов, на результат свёртки и на модуль ТЧП Гаусса μ=p+jq ,   от выбора которого многое зависит. Даже при одинаковой норме   N   выбор конкретных   p и q   очень важен. Было показано, что лучше, если   p и q сильно отличаются друг от друга, т.е. p >> q   или   p << q .

Поэтому выбор чисел Ферма   Fn = 1+22n = 1+2b   в качестве нормы   N   просто очевиден.   В этом случае   p=1   и   q=2b/2 .   Значит, если

N = p2 + q2 = 12 + (2b/2)2 = 1+2b = 1+22n = Fn ,
то для того, чтобы результат комплексной свёртки лежал внутри ПСАНВ± по модулю ЦКЧ Гаусса   μ=p+jq   достаточно, чтобы его значения по абсолютной величине (т.е. амплитуда) не превосходили
(ЦКЧ)A γ   ≤   ( N/(p+q) )/2 = ( (1+2b)/(1+2b/2) )/2  
( 2b/2b/2 )/2 = ( 2b/2 )/2 .

Но это же эквивалентно тому требованию, что амплитуда результата свёртки (уже вещественных и целочисленных) векторов   y[t] = ρ( γ[t] )   должна лежать внутри ПСАНВ± по модулю числа Ферма   Fn-1

(ЦЧ)A y   ≤   2b/2 = 22n-1     1+22n-1 = Fn-1 .

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

T · (ЦКЧ)A η   ≤   (1/2) · √ N/((p+q))       ( 2b/2 )1/2/2 = ( 2b/4 )/2 .

Но это же равносильно тому требованию, что амплитуда входных данных h[t] = ρ( η[t] )   и   x[t] = ρ( ξ[t] )   для (уже вещественной и целочисленной) свёртки должна лежать внутри ПСАНВ± по модулю числа Ферма   Fn-2

T · 2·(ЦЧ)A x,h   ≤   2b/4 = 22n-2     1+22n-2 = Fn-2 .

NB. В итоге получается интересное соотношение между величиной модуля преобразования ( при   N=Fn ), диапазоном результата операции (ЧП Ферма или ЦС) и диапазоном входных данных (это при условии, что входные данные поступают из кольца ЦКЧ Гаусса)

( √T · 2·(ЦКЧ)A η ) : (2·(ЦКЧ)A γ ) : N                 
  √ N/(p+q) : N/(p+q) : N  
  2b/4 : 2b/2 : 2b  
( √T · 2·(ЦЧ)A x,h ) : (2·(ЦЧ)A y ) : N    
        Fn-2 : Fn-1 : Fn    
    22n-2 : 22n-1 : 22n

Возможно, что именно из-за этого соотношения исчезает деление на   N=Fn   из формулы для   ρ1 .



ДП Фурье квантованных сигналов в кольце ЦКЧ Гаусса

 

 ...Моя проза: пойми, что пишу для заработка —  ч т е н и е
в с л у х , т.е. усиленно-членораздельного и пояснительного.
Стихи — для себя, прозу — для всех (рифма — "успех"). Моя
вежливость не позволяет мне стоять и читать моим "последним
верным" явно не понятные вещи — за их же деньги. Т.е. часть
моей тщательности (то, что ты называешь анализом) — вызвана
моей сердечностью. Я — отчитываюсь. А Бунин ещё называет мою
прозу "прекрасной прозой, но безумно-трудной", когда она —
для годовалых детей.
                                         Марина Цветаева,
                из письма к Борису Пастернаку, 20-е годы ХХ

ЧП Ферма     ▲ В начало текущей     Примеры ТЧП Гаусса






Последнее обновление 16.09.2013

© 2005 г., Александр Тимофеев, г.Харьков, Украина, eMail: atimopheyev@yahoo.com