Материал: 1503

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Таким образом, для финитных (ограниченных во времени) функций использование усеченного ряда Котельникова предполагает:

задание допустимой относительной

погрешности у и

вычисление по (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) на 12 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 12 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