◄ ЧП Ферма
▲ Выше
Примеры ТЧП Гаусса ►
|
Теоретико-числовое преобразование Гаусса |
μ = p+jq N = p2+q2 ρ( Zμ[j] ) = ZN |
Она глядит на вас так нежно, она лепечет так небрежно, она так тонко весела, её глаза так полны чувством, вечор она с таким искусством из-под накрытого стола свою мне ножку подала! Зачем я ею очарован? Зачем расстаться должен с ней? Когда б я не был избалован цыганской жизнию моей, когда б не смутное влеченье чего-то жаждущей души, я б здесь остался — наслажденье вкушать в неведомой тиши; забыл бы всех желаний трепет, мечтою б целый мир назвал... И всё бы слушал этот лепет, всё б эти ножки целовал.
Александр Пушкин
сентябрь 1833. |
Всё это было бы так просто, Когда бы не было так сложно. Поговорка XX века |
Многие вещи нам не понятны не потому, что наши понятия слабы; но потому, что сии вещи не входят в круг наших понятий. Козьма Прутков "Мысли и афоризмы",66. |
Пример ПСНВ+ и ПСАНВ± по модулю ЦКЧ Гаусса μ=(5+j2) |
Пример ПСНВ+ и ПСАНВ± по модулю ЦКЧ Гаусса μ=(2+j5) |
Сравнение ПСНВ+ и ПСАНВ± по модулю μ=(5+j2) и μ=(2+j5) |
Опять скажу — смотри в корень Козьма Прутков "Мысли и афоризмы",5. |
Примитивные КОРНИ и модулярное ТЧП Гаусса. Теорема. Если ε = a+jb — ЦКЧ взаимно простое с ЦКЧ μ , то
〈
εN-1
〉μ
= 1 ,
иными словами ЦКЧ
ε
=
a+jb
является КОРНЕМ из 1 порядка
(N-1)
по модулю ЦКЧ
μ
.
Теорема. Если ε = a+jb — ЦКЧ взаимно простое с ЦКЧ μ , а T — наименьшее из чисел, при которых
〈
εT
〉μ
= 1
(или, что тоже самое,
ε ≡ T√1 mod μ) ,
то
T|(N-1)
и
T
делит всякий другой показатель K степени КОРНЯ из 1
ε
,
при котором
〈
εK
〉μ
= 1 .
При наименьшем T таком, что T=φ(N), где φ — функция Эйлера, КОРЕНЬ ε будет примитивным (primitive) по модулю ЦКЧ μ . Теорема. Если ε — примитивный КОРЕНЬ по модулю ЦКЧ μ и Norma(μ)=N, то члены ряда (из N−1 значений)
1,
ε1,
ε2
,...,
εN−2
будут попарно несравнимы между собой, т.е. этот ряд представляет полную
приведенную (по модулю
μ
) систему ВЫЧЕТов, если к нему присоединить нулевой элемент.
Пусть μ=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])
.
Построенный таким образом базис будем называть базисом Гаусса, а
соответствующее ему преобразование — ТЧП Гаусса.
Естественно, ограничимся такими КОРНЯми
ε ≡ T√1 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] (смотри схему) .
Тоже самое и для свёртки двух сигналов из ЦКЧ.
NB. Это позволяет по новому взглянуть на машинную реализацию ДП Фурье комплексных сигналов, а именно, рассматривая их как сигналы в кольце ЦКЧ Гаусса и вычисляя их спектры в кольце целых чисел через ЧП Ферма (или ТЧП Рейдера, в которых умножения заменяются на сдвиги и сложения) с возвратом обратно в кольцо ЦКЧ Гаусса. |
Нет пределов совершенству. Поговорка |
Выбор параметров кольца Гаусса. Теорема. Если ε принадлежит кольцу ЦКЧ Zμ[j] , и Norma(ε) — простое целое число, то ε — простое число Гаусса. Теорема. Норма простого числа Гаусса является либо простым числом, либо квадратом простого числа. Соответственно различают простые числа Гаусса I-го и II-го порядка. Простые натуральные (целые положительные) числа, равные нормам простых ЦКЧ Гаусса, называются разложимыми. Простые натуральные числа, остающиеся простыми в кольце Zμ[j] , называются неразложимыми. Теорема. Каждое простое число p=4s+1 является разложимым в кольце
Zμ[j] .
Теорема. ТЧП Гаусса по модулю μ — примитивного числа Гаусса (т.е. НбОД (p,q)=1), существует тогда и только тогда, когда норма N — нечетное целое и имеет максимальный объём преобразования, выражающийся числом вида 4s.
NB. Наиболее широко используются кольца
Zμ[j]
такие, у которых
Например, при
p=1
и
q=2b/2
,
где
b=2n, n=1,2,...
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]-арифметике
(обычных целых комплексных чисел). Очевидно, это возможно в том случае, когда результаты,
вычисленные в кольце Гаусса
P1(re)
=
P1(im)
=
( N/(p+q) )/2
≈
( Fn-1−1 )/2
.
Это приводит к необходимости ограничить рабочие области квантованных комплексных сигналов
и переходных функций (импульсных характеристик фильтров), т.е. входных данных комплексной
свёртки, величиной
( N/(p+q) )/2
≈
Поэтому для свёртки γ[k] комплексных сигналов η[t] и ξ[t] имеем
где
║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-2−1 )/2
.
Тоже самое и для результата комплексной свёртки, если предположить, что
max |γ(re)[k]|
=
max |γ(im)[k]|
=
(ЦКЧ)A γ
,
тогда
(ЦКЧ)A γ
≤
⌋
( N/(p+q) )/2
⌊
≈
( Fn-1−1 )/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]
в кольце ЦКЧ Гаусса можно записать так
Точно такие же формулы для оценки динамического диапазона комплексного сигнала из ЦКЧ Гаусса будут и при вычислении ТЧП в кольце Гаусса
Однако, для комплексных сигналов с нулевой мнимой частью имеем
|
Нет вещи столь великой, которую не превзошла бы величиною ещё большая. Нет вещи столь малой, в которую не вместилась бы ещё меньшая. Козьма Прутков "Мысли и афоризмы",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
сильно отличаются друг от друга, т.е.
Поэтому выбор чисел Ферма 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
.
Но это же эквивалентно тому требованию, что амплитуда результата свёртки
(уже вещественных и целочисленных) векторов
2·(ЦЧ)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
.
Но это же равносильно тому требованию, что амплитуда входных данных
√T
·
2·(ЦЧ)A x,h
≤
2b/4
=
22n-2
≈
1+22n-2
=
Fn-2
.
NB.
В итоге получается интересное соотношение между величиной модуля преобразования
(
при
N=Fn
),
диапазоном результата операции (ЧП Ферма или ЦС) и
диапазоном входных данных
(это при условии, что входные данные поступают из кольца ЦКЧ Гаусса)
Возможно, что именно из-за этого соотношения исчезает деление на N=Fn из формулы для ρ−1 . |
ДП Фурье квантованных сигналов в кольце ЦКЧ Гаусса |
...Моя проза: пойми, что пишу для заработка — ч т е н и е в с л у х , т.е. усиленно-членораздельного и пояснительного. Стихи — для себя, прозу — для всех (рифма — "успех"). Моя вежливость не позволяет мне стоять и читать моим "последним верным" явно не понятные вещи — за их же деньги. Т.е. часть моей тщательности (то, что ты называешь анализом) — вызвана моей сердечностью. Я — отчитываюсь. А Бунин ещё называет мою прозу "прекрасной прозой, но безумно-трудной", когда она — для годовалых детей. Марина Цветаева, из письма к Борису Пастернаку, 20-е годы ХХ |
◄ ЧП Ферма ▲ В начало текущей Примеры ТЧП Гаусса ► |
Последнее обновление 16.09.2013
© 2005 г., Александр Тимофеев, г.Харьков, Украина, eMail: atimopheyev@yahoo.com |