■D
□------------------
Запропоновано метод передискретиза-ції дискретних сигналів, який не потребує поновлення безперервних сигналів по його відлікам з подальшою його дискретизацією на новій частоті, який не потребує ніяких фільтрів за для передискретизації та дискретного перетворення Фур’є, а також великих обчислювальних витрат
Ключові слова: дискретний сигнал, кортеж, образ сигналу, багаторівневе перетворення, матриця
□----------------------------------□
Предложен метод передискретизации дискретных сигналов, не требующий восстановления непрерывного сигнала по его отсчётам с последующей дискретизацией его на новой частоте, не требующий никаких фильтров для передискретизации и дискретного преобразования Фурье, а также больших вычислительных ресурсов
Ключевые слова: дискретный сигнал, кортеж, образ сигнала, многоуровневое преобразование, матрица
□----------------------------------□
The method resampling of discrete signals not requiring restoration of a continuous signal on its readout with subsequent resampling it on new frequency, which is not requiring any filters for resampling and discrete transformation Fourier, and large computing resources is offered
Keyword: discrete signals, tuple, image signal, multilevel transformation, matrix -------------------□ □---------------------
УДК 004.383.3
МНОГОУРОВНЕВЫЙ
МЕТОД
ПЕРЕДИСКРЕТИЗАЦИИ
ДИСКРЕТНЫХ
СИГНАЛОВ
И.Г. Филиппенко
Доктор технических наук, профессор, заслуженный деятель науки и техники Украины Украинская государственная академия железнодорожного транспорта пл. Фейербаха, 7, г. Харьков, Украина, 61050 Контактний тел.: (057) 702-59-60
И.О. Филиппенко
Харьковский национальный университет радиоэлектроники пр. Ленина, 14, г. Харьков, Украина, 61166
1. Введение
Передискретизация (aHra.resampling) в обработке сигналов - это изменение частоты дискредитации дискретного сигнала. Алгоритмы передискретизации широко применяются при обработке звуковых сигналов, радиосигналов и изображений. Теоретически ’’Идеальная передискретизация эквивалентна восстановлению непрерывного сигнала по его отсчётам с последующей дискретизацией его на новой частоте [2]”. Её вычисление производится по выражению, приведенному в [4, р. 3],
sinW^t -1,)
x(t)=Е -------------,
1 WSt-t,)
где x(tj) - значение i-го отсчёта сигнала, t, - момент времени, соответствующий этому отсчёту, rnj = 2nf - циклическая частота дискретизации, x(t) - интерполированное значение сигнала в момент времени t. Это выражение ”...не реализуемо на практике [2] ” из-за того, что для осуществления вычисления ’необходимо обработать бесконечное число его отсчётов (как в прошлом, так и в будущем) [2]”.
Интерполированное значение сигнала x(t) в момент времени t может быть вычислено при помощи интеграла свертки [4, р. 21]:
x(t) =Х x(t,)h(t -1,),
i
где h(t) - импульсная характеристика восстанавливающего фильтра.
”Прямое вычисление новых отсчётов сигнала по вышеприведённым формулам требует значительных вычислительных ресурсов... [2]”. Существуют более простые частные случаи передискретизации, к которым относятся децимация и интерполяция с целыми коэффициентами (изменение частоты дискретизации в целое число раз) и изменение частоты дискретизации в рациональное число раз. В этих случаях используется широкий класс цифровых фильтров (см. [2]).
Существует алгоритм передискретизации при помощи дискретного преобразования Фурье (ДПФ), который включает в себя вычисление ДПФ, вставки в середину спектра необходимого числа нулевых компонент и вычисление обратного ДПФ с нормировкой [2].
Целью работы является разработка метода передискретизации дискретных сигналов не требующего восстановления непрерывного сигнала по его отсчётам с
у5
последующей дискретизацией его на новой частоте, не требующего никаких фильтров, дискретного преобразования Фурье и больших вычислительных ресурсов.
2. Постановка задачи
Задан непрерывный сигнал, подвергшийся дискретизации, и представленный в виде ряда значений N точек
X = (x0,x1,...,xi,...,xN-1) (1)
N-мерным вектором, компоненты которого представляют собой значения i-ых точек отсчета дискретного сигнала для заданной частоты дискретизации fd.Необходимо осуществить передискретизацию вектора (1) методом, не требующим восстановления непрерывного сигнала по его отсчётам с последующей дискретизацией его на новой частоте, не требующим никаких фильтров, дискретного преобразования Фурье и больших вычислительных ресурсов.
3. Основная часть
Процедуре передискретизации дискретного сигнала предшествует прямое многоуровневое преобразование вектора X (1) [1], суть которого сводится к следующему. Вектор X (1) вначале преобразуется к униполярному (если сигнал биполярный) виду
X ^ X+ =
x+ ,x+V
.,xi , —,XN—1
(2)
где X+ I i=[0N-I] = xi + lmin Znachnl ,
а затем к нормализованному виду
—+ X
/ + + + + \ (xo,X1 ,...,Xi,...,Xn-1 ) ,
где
—+ Xi
i=[0,N-1]
Xi /maxZnachn
+
xi
*
1.0
maxZnachn
(3)
если maxZnachn > 1 если maxZnachn < 1.0
Здесь minZnachn и maxZnachn - минимальное значение точки отсчета биполярного и максимальное значение точки отсчета униполярного сигнала, соответственно.
Затем осуществляется преобразование нормированного вектора в матричный вид
X =(Xo, X1,..., Xi,..., Xn-1 ) =
x0;0 X0;1 ... X0;i ... x0;N -1
X1;0 X1;1 ... X1;i ... x1;N -1
xr;0 Xr1 ... Xr;i ... Xr;N-1
xL-1;0 xL-1;1 ... XL-1;i ... XL-1;N-1
(4)
при помощи скал_ярно-векторного преобразования компонент вектора X следующим образом.
Xi
^ X,
r=[0,L-1]
= (X0;i,X1;i,...,Xr;i,...,XL-1;i)T
(5)
где Xr;i I r=[0, L-1] i=[0,N-1]
d, если (xi - q ■ r) > 0; 0, если (xi - q ■ r) < 0.
Здесь: r є (0, 1, ,...,N-1); L - число уровней скалярновекторного преобразователя; d = 1; q - шаг квантования xi сигнала по уровню, который определяется как q = 1.0/L. ^+
Каждый элемент матрицы X может принимать значение нуля или единицы.
Далее осуществляется преобразование матрицы (4) в кортеж P, названный в работе [1] образом униполярного нормированного сигнала (3),
/S0,
S1,
P =
Sr,
(6)
S
L -1
для чего каждой X(r) строке (4) ставится в соответствие кортеж Sr, элементами которого являются кортежи длиной 2 (пары)
X
+
(r)
r=[0,L-1]
^ Sr
/ < ileft r;0, iright r;0 \< ileftr;j,irightr;j
>, < heftr;1,irightr;1 >,.. >,..., < ileft r;nr, iright r;nr
>
, (7)
где каждая пара (<*, *>) представляет собой два номера точек отсчета X(r) строки, в которых произошли изменение нулевого значения xri отсчета на единичное значение, а затем наоборот - изменение единичного значение отсчета на нулевое значение.
Многоуровневый метод передискретизации однополярных нормированных дискретных сигналов
Известно[2], что вид передискретизации определяется частотой передискретизации fresampı и частотой дискретизации fd исходного сигнала, предназначенного для передискретизации. При этом, начальное число отсчетов не передискретизированного сигнала - N, пе-редискретизированного сигнала - Nresampı. Если fresampı > fd, то речь идет о повышении частоты дискретизации, т.е. об интерполяции (англ. interpolation). Если fresampı < fd, то речь идет о понижении частоты дискретизации, т. е. о децимации (англ. decimation) [2].
Суть предлагаемого метода передискретизации дискретных униполярных нормированных сигналов сводится к преобразованию каждой компоненты пары образа не передискретизированного сигналу в компоненты новой пары образа Presampl по следующей формуле:
resampl
ц.ч.ч.1 если д.ч.ч.1 < 0.5 ц.ч.ч.1 +1 если д.ч.ч.1 > 0.5 ц.ч.ч.2 если д.ч.ч .1 < 0.5 ц.ч.ч.2 +1 если д.ч.ч .1 > 0.5
, если N > Nresampl
, (8)
, если N < Nresampl
где:
ц.ч.ч.1 = ц.ч.ч. (icomp /((N- 1)/(Ndecim -1))) , (9)
д.ч.ч.1 = д.ч.ч. (icomp/ ((N - 1)/(Ndecim -1))), (10)
ц.ч.ч.2 = ц.ч.ч. (icomp x((Nmterpoi - 1)/(N-1))), (11)
д.ч.ч.2 = д.ч.ч. (icomp x((Ni„terpoi - 1)/(N-1))), (12)
где iresampl - значение компоненты передискрети-зированного образа; icomp - значение компоненты не
Системы управления
передискретизированного образа; ц.ч.ч. - целая часть числа, заключенного в круглые скобки; д.ч.ч. - дробная часть числа, заключенного в круглые скобки: Ndecim, Ninterpol - число отсчетов передискретизированного дискретного сигнала в случае децимации и в случае интерполяции, соответственно.
В вышеприведенных формулах (9 - 12) выражения
((N - 1) / (Ndecim -1)), (13)
((Ninterpol - 1)/(N -1)) (14)
это коэффициенты передискретизации, в дальнейшем именуемые как ’’коэффициент децимации” и ’’коэффициент интерполяции”, соответственно.
В приведенных выше коэффициентах передискретизации все значения чисел точек отсчетов уменьшены на единицу, поскольку нумерация точек отсчета начинается с нуля (1).
Коэффициенты передискретизации, выраженные через частоты передискретизации, легко определяются из (13) и (14), и имеют следующий вид для коэффициента децимации как:
((f/k - 1)/(fdec,mtk - 1)) , (15)
для коэффициента интерполяции как;
while(in)
//#26 Restoration of a discrete signal on an image &&& from &&& void Resmp::restorDscrtnUnplrSgnlFromImage( char *nameImage,
double* pt_ArrRstrSgnl, unsigned long int stlbArrRstrSgnl )
{
unsigned int left = 0; //index nachala edinichnogo potoka unsigned int right = 0; //index kontsa edinichnogo potoka unsigned int count = 0; // dlia while
unsigned int nmbrElmnt = 0; //kolichestvosimvolov v obraze signala
ifstream in(namelmage); // otkrytie file dlia vvoda
if(!in){
cout << "namelmage otkryt' nevozmozhno\n";
}
in.close(); //zakryt' potok in.open(nameImage);
nmbrElmnt++;
count++;
if(count
{
}
else
{
}
}
1)
in >> left; count++;
in >> right;
//cout << "left = " << left << 'W << "right = " << right << endl; for(unsigned int stlbBkTrnsfSgnl = 0; stlbBkTrnsfSgnl < stlbArrRstrSgnl;
stlbBkTrnsfSgnl++)
{
if( (left <= stlbBkTrnsfSgnl) && (stlbBkTrnsfSgnl <= (right ) ) ) //- 1 *(pt_ArrRstrSgnl + stlbBkTrnsfSgnl) += q;
}
//outFile << "left = " << left << " right = " << right << endl; left = 0; right = 0; count = 0;
cout << "nmbrElmnt = " << nmbrElmnt - 1 << endl;
in.close();
}
((Le^k - 1)/(fdtk -1)), (16)
где tk - длительность дискретного сигнала.
В результате образуется передискретизированный образ униполярного нормированного дискретного
сигнала Presampl.
Восстановление передискретизированного дискретного сигнала по передискретизированному образу с целью оценки степени подобия передискретизированного и_со ответствующего идеального сигналов.
Каждая Xi компонента вектора X (3) восстанавливается по формуле, взятой из [1]
Xi
i=[0,N-1] ^ Xi
q! fr,
(17)
где Xi
элемент восстановленного вектора
X
f
Vre[0,L-1]
1, если 1, если
^left;*, ^right;* t )^Sr fp ■ “ İright;*
IA
ч ^left;*, iright;* t HSr w >>i| r ^ İright;
Читается fVre[0L-1] следующим образом. Значение функции fVre[0L-1] для всех г, принадлежащих [0, L-1], равно единице, если в Sr существует пара (кортеж) (iıeft.„iright.„) , удовлетворяющая условию f < i|r < iright.*, т.е. i-ый отсчет г-го уровня дискретного многоуровневого сигнала находится в пределах [iıeft;*, iright;*]. В противном случае значение функции fr равно нулю.
На листинге № 1 представлена функция на СИ++ класса Resmp, восстанавливающая дискретный сигнал по заданному образу сигнала.
Листинг № 1.
Рассмотрим описанный выше метод передискретизации на примере.
Пример передискретизации дискретных сигналов
Задан непрерывный сигнал xı(t)
Xj(t) = sin(2n- 2f ■ t) + sin(2n- 3f ■ t) + sin(2n- 64f ■ t), (18) где fi - частота, равная 1 Hz,
и непрерывный треугольный периодический сигнал x2(a, b, c, d, u, t)
x2(a, b, c,d,u,t) =
0 для t < 0, t-(a+k(t))
b + k(t)-(a + k(t))
для (a + k(t)) < t < (b + k(t)),
(c + k(t+))-((b)+ k(t)) для (b + k(t))<' <(c + k(t)), , <19>
-(d+tt(tC))+,k(1)) для (c + k(t>)<1 <(d + k(t»,
(u + k(t)) -t (u + k(t))- (d + k(t))
для (d + k(t))< t <(u + k(t)).
где k(t) - ступенчатая возрастающая функция, показанная на рис. 1, определяется как
k(t) = n ■ u|t=n-u+r, (20)
где n є I - целые числа; u,r eR - действительные числа; 0 < г < u .
Рис. 1. Ступенчатая возрастающая функция
Сформирован исходный непрерывный сигнал x(t) как произведение двух сигналов xj(t) (18) и x2(t) (19)
x(t) = x1(t)' x2(t). (21)
Осуществим дискретизацию сигнала x(t) путем замещения t на idT, где i - номер отсчета дискретного сигнала, dT - период квантования непрерывного сигнала x(t) по времени. В результате дискретизации имеем исходный дискретный сигнал
Далее производится преобразование дискретного сигнала, представленного на рис. 2, к униполярному нормированному виду, показанному на рис. 3.
Рис. 3. Исходный униполярный нормированный дискретный сигнал
0
x = x . x Перед передискретизацией исходного униполяр-
ного нормированного дискретного сигнала (УНДС) где функции x1i и x2i принимают следующий вид: создадим его образ согласно формулам (4) - (7).
Фрагменты образа не передискретизи-x1i =(sin(2n.2.i.dT) + sin(2n.3.i.dT) + sin(2n.64.i.dT)).xtmg_j,(22a) рованного исходного УНДС приведены на
листинге № 2.
Передискретизация УНДС осуществляется путем преобразования его образа и восстановления сигнала по преобразованному образу. Восстановленный сигнал и будет передискре-тизированным дискретным сигналом.
x2i (a, b,c,d,u,i • dT) = -
i • dT-(a+k(i • dT))
b + k(i • dT)-(a + k(i • dT)) (c + k(i • dT))-i • dT
(c + k(i • dT)) - (b + k(i • dT)) i • dT-(c + k(i • dT))
- (d + k(i • dT)) - (c + k(i • dT))
(u + k(i • dT))-i • dT
- (u + k(i • dT)) - (d + k(i • dT))
для i • dT < 0,
для (a + k(i • dT)) < i • dT < (b + k(i • dT)), для (b + k(i • dT)) < i • dT < (c + k(i • dT)), (23a) для (c + k(i • dT)) < i • dT < (d + k(i • dT)), для (d + k(i • dT)) < i • dT < (u + k(i • dT)).
где
k(i.dT) = n. uİ!.dT=n.n+r .
В данном примере принято: период квантования по времени исходного дискретного сигнала dT = 0.000725 сек.; частота дискретизации fd = 1380,07 Hz; число точек отсчета N = 2070; длительность сигнала tk = dTN = =1.50075 сек.
График зависимости исходного дискретного сигнала от i показан на рис. 2.
Рис. 2. График зависимости исходного дискретного сигнала от i
Листинг № 2.
Фрагмент не передискретизированного УНДС
P2070 -
1223 1225 2069 1222 1225 2069 1222 1225 2069 1201 1203 1222 1201 1204 1222 1201 1204 1222 1201 1204 1221 1200 1204 1221 1200 1204 1221 1200 1205 1221 1200 1205 1221 1200 1205 1221 1200 1205 1221
1225 2069
1226 2069 1226 2069 1226 2069 1226 2069 1226 2069 1226 2069 1226 2069
1226 1244 1246 2069
1227 1244 1246 2069
Пусть необходимо изменить частоту дискретизации УНДС с fd = 1380,07 Hz на fdecim = 384,76 Hz, т.е. осуществить децимацию сигнала.
Для этого осуществим преобразование образа исходного УНДС с частотой дискретизации fd = 1380,07 Hz, представленного на листинге № 2 в образ переди-скретизированного УНДС с частотой дискретизации fdecim = 384,76 Hz путем вычисления каждой компоненты новой пары кортежа Р577из2070 по формуле (8).
Фрагменты образа передискретизированного исходного УНДС приведены на листинге № 3.
Листинг № 3. Фрагмент передискретизированного образа исходного УНДС
//#62 Creation of the other image out of given image
void Resmp::crtnOtherImgOutOfGvnImg( char *nameGvnImg, char *nameOtherImg, double f_dT_Gvn, double f_dT_Other )
{
int count = 0; long nmbrElmnt = 0; double left = 0; double leftTo = 0; double right = 0; double right512 = 0;
ofstream out(nameOtherImg, ios::out); //sozdanie faeila dlia vyvoda if(!out){cout << "File nameOtherImg otkryt' nevozmozhnoW;} out.clear();
ifstream in(nameGvnImg, ios::in); // otkrytie file dlia vvoda if(!in){
cout << "File nameGvnImg otkryt' nevozmozhnoW;} in.close(); //zakryt' potok in.open(nameGvnImg); //otkryt' potok
0 340 341 576 / 0 340 341 576 / 0 340 341 576 0 334 335 340 341 576 0 334 335 340 341 576 0 334 335 340 341 576 Р577из2070 = ( 0 334 33 5 340 34 1 576
0 334 335 340 341 576 \ 0 334 335 340 341 576
\ 0 334 335 340 341 576 \ 0 334 335 340 341 576 \ 0 334 335 340 341 346 347 576 \0 334 335 340 342 346 347 576
Из сравнения строк образов P2070 и Р577из2070 можно увидеть процедуру преобразования одного образа в другой образ как преобразование пары одного образа в пары другого образа.
Разместить образы, а не их фрагменты, на листингах формата А4 не представляется возможным потому что образ P2070 содержит 5338 пар, образ Р577из2070 содержит 5182 пар.
Процесс восстановления сигнала по вновь полученному образу Р577из2070 осуществляется по формуле (17) и функции, приведенной на листинге № 1. Восстановленный сигнал автоматически становится передискре-тизированным сигнал, показанным на рис. 4.
Рис. 4. Передискретизированный сигнал
На листинге № 4 приведена функция на СИ++ преобразования образа, заданного char *nameDvnImg в образ, заданный *nameOtherImg.
Листинг № 4. Функция преобразования образа
while(in) //
{
nmbrElmnt++;
count++;
if(count == 1)
{
in >> left; count++;
}
{
in >> right;
if(left == 0)
{
leftTo = 0.0;
}
{
leftTo = trnsfrmElmntOtherImgOutOfElmntDvnImg( f_dT_Other, f_dT_Gvn, (double)left ); //#61 sin206
}
if(right == 0)
{
leftTo = 0.0;
}
{
right512 = trnsfrmElmntOtherImgOutOfElmntDvnImg( f_dT_Other, f_dT_Gvn, (double)right ); //#61 sin206
}
//cout << leftTo << ' ' << right512 << ' ' << endl; out << leftTo << ' ' << right512 << ' ' << endl; left = 0; right = 0; count = 0;
}
}
//nmbrParnameOtherlmgnameGvnImg = (double)(nmbrElmnt - 1)/2.0;
cout << "nmbrElmntImage nameGvnImg = " << (double)(nmbrElmnt - 1)/2.0 << endl;
in.close();
out.close();
}
Естественно возникает вопрос, насколько передискретизация сигала идеальная. ’’Идеальная передискретизация эквивалентна восстановлению непрерывного сигнала по его отсчётам с последующей дискретизацией его на новой частоте.[2]”. В рассматриваемом примере нет необходимости в восстановлении непрерывного сигнала, поскольку для чистоты компьютерного эксперимента исходный непрерывный сигнал задан (21). Нам необходимо только передискретизировать его на частоте дискретизации fdecim = 384,76 Hz, что и было сделано. Идеальный сигнал, с которым сравнивается передискретизированный сигнал, показан на рис. 5.
Рис. 5. Идеальный сигнал с частотой дискретизации fd = 384,76 Hz
сигнала
Визуально количественно оценить степень подобия (идеальности) дискретных сигналов, приведенных на рис. 4 и рис. 5 практически невозможно, Они практически идентичны друг другу. Для количественной оценки степени подобия передискретизированно-го (рис. 4) и идеального сигнала (рис. 5) воспользуемся коэффициентом корреляции [3, стр. 52]:
г =
м
N-1
I fkgk
I f
k=0
(24)
где (f,g) - скалярное произведение векторов f и g; ||f||||g|| - нормы векторов f и g; fk и gk - значения точек отсчета векторов f и g.
Функция на СИ++, реализующая вычисление коэффициента корреляции, приведена на листинге 5.
Листинг № 5. Функция на СИ++, реализующая вычисление коэффициента корреляции двух векторов
//#28 sin32 Coefficient correlation
double Resmp::coefCorrl( unsigned size, double* pt_Vctr_1, double* pt_Vctr_2 )
{
double skaliarProizvdn = 0; double summKvdrtRaspoznSgnl = 0; double modRaspoznSgnl = 0; double summKvdrtObrazSgnl = 0; double modObrazSgnl = 0;
for(unsigned int i = 0; i < size; i++)
{
skaliarProizvdn += (*(pt_Vctr_1 + i) * *(pt_Vctr_2 + i)); summKvdrtRaspoznSgnl += pow(*(pt_Vctr_1 + i), 2.0); summKvdrtObrazSgnl += pow(*(pt_Vctr_2 + i), 2.0);
}
//cout << "skaliarProizvdn = ” << skaliarProizvdn << endl;
//cout << "ssummKvdrtRaspoznSgnl = " << summKvdrtRaspoznSgnl << endl; //cout << "summKvdrtObrazSgnl = " << summKvdrtObrazSgnl << endl;
modRaspoznSgnl = sqrt(summKvdrtRaspoznSgnl); modObrazSgnl = sqrt(summKvdrtObrazSgnl); //summKvdrtObrazSgnl
//cout << "modRaspoznSgnl = " << modRaspoznSgnl << endl;
//cout << "modObrazSgnl = " << modObrazSgnl << endl;
//cout << "modRaspoznSgnl*modObrazSgnl = " << modRaspoznSgnl*modObrazSgnl << endl;
}
return skaliarProizvdn / (modRaspoznSgnl*modObrazSgnl);
Коэффициент корреляции г577иэ2070и577 между пере-дискретизированным (рис. 4) и идеальным УНДС (рис. 5) равен 0.911933.
С целью более детальной проверки предложенного метода передискретизации дискретных сиг-
налов были выполнены, описанным выше методом, следующие передискретизации: f2070^577, f2070m6100, f610fa32070. Произведены оценки степени подобия пере-дискретизированных и соответствующих идеальных сигналов. Результаты полученных оценок приведены в табл. 1, где s - символ подобия сигналов, заданных векторами.
Таблица 1
Коэффициенты корреляции идеальных и соответствующих
передискретизированных сигналов
Объекты Коэффициент корреляции (г) сигналов
корреляции униполярных биполярных
Xf = Xf 0.997465 0.911933
Xf = Xf 0.996559 0.947199
f2070 12070из577
Xf = Xf 0.999706 0.980699
16100 f6100 из 2070
Xf = Xf 0.999782 0.988586
f2070 12070из6100
В табл. 1 под объектами корреляции понимаются идеальный дискретный сигнал (например X^ ) и передискретизированный (например Xf ) сигнал.
Идеальный дискретный сигнал, получен из известного непрерывного сигнала путем его дискретизации с частотой дискретизации fd = 384,76 Hz. Передискретизированный сигнал (например Xf ) получен в ре-
зультате передискретизации с частоты дискретизации 1380.07 Hz на частоту 384.76 Hz . При частоте дискретизации 384,76 Hz число отсчетов N = 577. При частоте дискретизации 1380.07 Hz число отсчетов N = 2070. Длина сигнала (tk) по прежнему равна 1.50075 сек.
4. Выводы
- Предложен новый метод передискретизации дискретных сигналов.
- Степень ’’идеализации” передискретизации униполярных дискретных сигналов для примеров, приведенных данной работы, лежит в пределах (0.999782 - 0.996559), для биполярных дискретных сигналов в пределах (0.988586 - 0.911933).
- Предложенный многоуровневый метод передискретизации дискретных сигналов не требует больших вычислительных ресурсов.
Литература
1. Филиппенко И.Г. Прямое и обратное преобразование дискретных сигналов / И. Г. Филиппенко, Филиппенко И. О. // Восточно-Европейский журнал передовых технологий. - 2010. - №4/3 (46). - С. 29-32.
2. Передискретизация. [Электронный ресурс] wikipedia - Режим доступа:http://гu.wikipedia.oгg/wiki/ %D0%9F%D0%B5%D1% 80%D0%B5%D0%B4%D0%B8%D1%81%D0%BA%D1%80%D0%B5%D1%82%D0%B8%D0%B7%D0%B0%D1%86%D0%B8%D 1%8F - 29.08.2010 г. - Загл. с экрана.
3. Сато Ю. Обработка сигналов. Первое знакомство. / Сато Ю. ; пер с яп.; под ред. Есифуми Амемия. - М.: Издательский дом «Додэка XXI”, 2009. - 176 с.
4. Crochiere Ronald E. Multirate digital signal processing / Ronald E. Crochiere, Lawrence R. Rabiner. - Prentice-Hall, 1983. - 411 p.