Программа 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