Передача дискретной информации и телеграфия. Бажанов В.Л. - 20 стр.

UptoLike

ВУЗ: 

СамГУПС | Самара

Составители: 

Рубрика: 

20
• выполнять процедуру цифровой фильтрации и фильтрации с
восстановлением;
• анализировать квадратичные оценки I
1
, I
2
, I
3
при различных настройках
фильтров;
• просматривать на мониторе одновременно три графика: X(t), Fx(t) и Xв(t).
Для более полного понимания программы «filtr-w.bas» ниже приводится ее
распечатка и комментарии.
Программа «filtr-w.bas»
(листинг)
5 ' Программа фильтрации с восстановлением 1-го порядка.
10 ' Конечно-разностная оценка производной для восстановления.
15 ' Вывод 3-х графиков.
20 SCREEN 0, 0: KEY OFF ' :COLOR 0,7
30 DIM X(3000), Y(3000), Z(3000)
40 CLS
50 w = 1: Tc = 2 * 3.14 / w: Fx = 0
60 INPUT " Введите Tc/D = (D-интервал оценки производной) ", sh
70 INPUT " Введите Tc/j = (j-шаг дискретизации) ", k2
80 INPUT " Введите Q = (степень фильтpации) ", Q
90 D = Tc / sh: j = Tc / k2: rr = D / j
100 T = Q * j ' Постоянная фильтра
110 DIM L(rr + 2)
120 REM ************************************************************
125 PRINT : PRINT "Формирование массивов X,Y,Z"
130 REM ************************************************************
140 Wp = 0: Wf = 0: Wv = 0 ' Для квадратичных оценок
170 FOR i = 0 TO 2 * k2 + 2 ' Цикл
172 XX = SIN(w * j * i) ' Полезный сигнал
174 REM Сигнал помехи.
176 Fp = SIN(100 * w * j * i) ' + SIN(200 * w * j * i) + (RND(1) - .5)
180 REM Наложение
помехи на истинный сигнал .
182 YY = XX + Fp
184 REM Фильтpация сигнала.
186 Fx = Fx + (YY - Fx) / Q
188 FOR e = 1 TO rr
190 L(rr - e + 1) = L(rr - e)
192 NEXT e
200 L(0) = Fx
202 IF i < k2 THEN GOTO 220 ' Фильтрация должна вработаться.
204 XW = L(0) + j * (Q - 1) * ((L(0) - L(rr)) / D) ' Восстановление.
206 REM На сколько восстановленный сигнал отличается от полезного.
208 Wv = Wv + ((XW - XX) ^ 2) * j
210 REM На сколько исходный сигнал отличается от полезного.
212 Wp = Wp + (Fp ^ 2) * j
214 REM На сколько фильтрованный сигнал отличается от полезного.
216 Wf = Wf + ((Fx - XX) ^ 2) * j

Страницы