Таким образом, для финитных (ограниченных во времени) функций использование усеченного ряда Котельникова предполагает:
задание допустимой относительной |
погрешности у и |
|
вычисление по (1.17) значения fmax; |
|
|
по |
формуле t 1/2fmax определение |
искомого интервала |
взятия |
отсчетов. |
|
Вернулись к той же проблеме: заранее надо знать спектральную плотность или корреляционную функцию исследуемого процесса.
Однако, |
если |
частота fmax |
ориентировочно известна, интеграл |
||
|
|
|
|
|
|
S( f )df |
0, |
тогда |
погрешностью у |
можно пренебречь. Этот |
|
fmax |
|
|
|
|
|
подход и получил распространение на практике. |
|||||
Можно привести много примеров, когда ширина спектра, |
|||||
регистрируемого с |
целью |
обработки |
процесса, изменяется во |
||
времени. До |
|
|
|
|
|
99 % мощности ураганного ветра переносится на частотах, расположенных до микрометеорологического максимума. При штиле для переноса такой же мощности необходимо расширить спектр частот на порядок. Поэтому естественным представляется метод дискретизации функции с переменной частотой взятия отсчетов, определяемый динамическими свойствами СП на текущем отрезке времени. Такой метод получил название метод адаптивной дискретизации.
Наибольшее распространение на практике получила следующая процедура. В точке ti берется отсчет функции x(ti ) (рис. 1.10), и начинается вычисление разности зарегистрированного значения отсчета и текущей функции: x x(t) x(ti ). При превышении разностью установленного порога 0, т.е. x 0, берется новый отсчет. Полученная функция имеет вид ступенчатой кривой (рис. 1.10, а).
Рис. 1.10. Пояснение идеи адаптивной дискретизации:
а – без использования производной; б – с использованием производной
В другом варианте (рис. 1.10, б) после взятия очередного отсчета x(ti ) формируется сигнал, x*(t) x(t) x (ti ) t , где x (ti )– производная функции x(t) в момент времени ti . Очередной отсчет производится в точке, определяемой выполнением равенства
x(t) x(t) x(ti) x (ti ) t 0.
Очевидно, при таком подходе к процессу дискретизации СП исследователь в меньшей степени зависит от априорных данных. Необходимость хранения времени взятия отсчетов, т.е. ti , – недостаток адаптивного метода.
Задачи
1. Источник информации генерирует колокольные импульсы вида
e 2t2 , которые через фильтр с частотой среза fmax подлежат записи в память ПЭВМ. Определить интервал дискретизации сигнала t, если допустимая относительная погрешность последующего его восстановления не должна превышать 10%, а =20 с-1.
Решение
В соответствии с (1.4) спектральная плотность колокольного импульса
|
|
2t2 |
e j tdt |
|
|
|
exp 0,5 2 |
2 . |
S f |
e |
|
||||||
|
|
|
||||||
|
|
|
|
|
|
|||
Принимая во внимание известное соотношение
|
1 |
|
|
|
|
|
|
exp 0,5 2 |
2 df , 2 f |
|
|
|
||
|
2 2 |
|||
|
0 |
|
||
и умножая числитель и знаменатель в выражении (1.17) на 1

2 2 , получим
|
|
|
|
|
|
|
|
2 1 Φ 2 fmax |
, |
|
|
S f df |
|
|
|
S f df |
|
||
|
f |
|
|
|
|||||
|
max |
|
0 |
|
|
|
|
||
|
|
|
|
|
|
|
|
||
где Φ z – функция Крампа.
Учитывая, что 0,1, известным способом определяем аргумент
функции Крампа: 2 fmax |
=2,4. При =20 и |
fmax=24/ получим |
t 1
2 fmax 6,56 10 2 с.
2. Найти относительную погрешность представления случайного синхронного двоичного сигнала рядом Котельникова при произвольной граничной частоте. Определить погрешность представления , если граничная частота выбирается равной fЭ .
Решение
Известно, что спектр случайного синхронного телеграфного сигнала описывается формулой
S f h2t |
sin2 fT |
, |
|
fT 2 |
|||
|
|
где Т– длительность посылки; h – ее амплитуда. Подставляя это выражение в формулу (1.17), получим
|
|
2 |
|
sin |
2 |
fT |
|
|
|
||
|
|
|
|
0,5 sin 2 fT |
, |
||||||
|
|
|
|
|
|
||||||
|
|
|
fT |
|
|
|
|
||||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
где sin x sin y/ y dy– интегральный синус. |
|
|
|||||||||
0 |
|
|
|
|
|
|
|
|
|
|
будем иметь = |
При f fЭ , |
учитывая, |
|
что fЭ T 0,5, |
||||||||
2 2 0,5 sin / 0,23. Если f fЭ , то |
2(0,5 sin 2 )/ |
||||||||||
/ 0,096. |
|
|
|
|
|
|
|
|
|
функцией K |
|
Случайный |
процесс |
|
с |
корреляционной |
|||||||
K 0 exp 
дискретизирован с шагом t. Найти погрешность представления такого процесса рядом Котельникова в зависимости от параметров и t.
Решение
Энергетический спектр случайного процесса
S f 2 K 0 / 2 2, |
2 f . |
При ограничении полосы такого спектра частотой fmax
погрешность 1 2
arctg 2 fmax
.
Принимая fmax 0,5 t , найдем 1 2
arctg
t .
1.5. Модели нестационарных случайных процессов
Приступая к обработке результатов наблюдений (измерений), исследователь должен принять решение о модели изучаемого процесса и прежде всего о его стационарности или нестационарности. Классический подход при решении этого вопроса строится по следующей схеме.
При наличии m реализаций (1, 2, …, k, …, m) и разбиении интервала задания СП на n отрезков (i=1, 2, …, n) можно рассчитать
оценки математического ожидания m*(t ), дисперсии |
D* 2 |
(t )* |
|
x |
i |
x |
i |
корреляционной функции Kx*(ti ,tj ). |
Обозначим |
любую |
из |
рассчитанных оценок индексом *(ti ) и построим для каждого сечения ti доверительные интервалы:
Iβm |
m*x(ti ) tβ |
σ |
, |
IβD D*(ti ) tβ σD . |
|
|
m |
|
(1.18) |
|
|
|
|
Что касается доверительного интервала для корреляционной функции, приводить его не будем из-за громоздкости выражения для
дисперсии 2* .
k
Ввыражении (1.18) индекс tβ вычисляется через β –
доверительную вероятность того, что найденная оценка попадет в
доверительный интервал tβ arcΦ |
* 1 |
β |
, где arcΦ |
* |
(x)– функция, |
|||||||||
|
|
|
|
|
|
|
|
|||||||
|
2 |
|
|
|||||||||||
обратная |
интегралу |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
||||
вероятностей |
|
Φ* (x) |
|
1 |
|
|
|
x |
t2 |
|
|
|||
|
|
|
|
|
|
e 2 dt . |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
||||||
Входящие |
в |
формулы |
|
|
|
доверительных |
|
интервалов |
||||||
среднеквадратические |
отклонения |
m* |
|
и |
D* определяются для |
|||||||||
нормального распределения величин по формулам:
m 
D /m ; D 
2D /m 1.
Охарактере связи t f ( ) можно судить по данным табл. 1.2.
Таблица 1.2
Значения используемых на практике доверительных вероятностей и связанных с ними индексов tβ
β |
t |
β |
t |
β |
t |
|
|
|
|
|
|
0,80 |
1,282 |
0,87 |
1,513 |
0,94 |
1,880 |
|
|
|
|
|
|
0,81 |
1,310 |
0,88 |
1,554 |
0,95 |
1,960 |
|
|
|
|
|
|
0,82 |
1,340 |
0,89 |
1,597 |
0,96 |
|