Программа 26

Программа 26.4

100 'ИТЕРАТИВНОЕ ПРОЕКТИРОВАНИЕ РЕКУРСИВНЫХ ФИЛЬТРОВ
110 '
120 'ИНИЦИАЛИЗАЦИЯ
130 N% = 256 'число точек в БПФ
140 NP% = 8 'число полюсов в фильтре
150 DELTA = .00001 'изменение приращения
160 MU = .2 'размер шага итерации
170 DIM REX[255] 'действительная часть сигнала при БПФ
180 DIM IMX[255] 'мнимая часть сигнала при БПФ
190 DIM T[128] 'желаемая частотная характеристика(только амплитуда)
200 DIM A[8] 'коэффициент рекурсии "a"
210 DIM B[8] 'коэффициент рекурсии "b"
220 DIM SA[8] 'наклон для коэффициентов "a"
230 DIM SB[8] 'наклон для коэффициентов "b"
240 '
250 GOSUB XXXX 'предполагаемая подпрограмма загрузки T[ ]
260 '
270 FOR P% = 0 TO NP% 'инициализация коэффициентов системы
280 A[P%] = 0
290 B[P%] = 0
300 NEXT P%
310 A[0] = 1
320 '
330 ' 'ЦИКЛ ИТЕРАЦИЙ
340 FOR ITER% = 1 TO 100 'цикл для желаемого числа итераций
350 GOSUB 2000 'вычисление новых коэффициентов
360 PRINT ITER% ENEW MU 'печатать текущее состояние
370 IF ENEW > EOLD THEN MU = MU/2 'изменить значение MU
380 NEXT ITER%
390 '
400 '
410 FOR P% = 0 TO NP% 'ПЕЧАТАТЬ КОЭФФИЦИЕНТЫ
420 PRINT A[P%] B[P%]
430 NEXT P%
440 '
450 END