Биквадратное звено на ADSP - 21XX

@Slave · Регистрация: 19.01.2010

Студворк — интернет-сервис помощи студентам

Необходимо было сделать бих фильтр 6 порядка (3 биквадратных звена). Написал код, фильтр не "фильтрует" (выдает значения, близкие к 0 без всяких переходных процессов). Не учитывал прерывания, т.к. писал для самого отладчика. Помогите разобраться, что я сделал не так.

Assembler
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
.module/ram/abs=0 filtr;
.var/dm/circ x_1[6];     {Zi}
.init x_1: 0,0,0,0,0,0;
.var/pm/circ c_f1[15];  {coefficients}
.init c_f1: 0x7CD1 ,     
0x7ae7 ,
0x2a0 ,
0x0000 ,
0xfd60 ,
 
0x83ea ,
0x7d80 ,
0x1383 ,
0xf1f3 ,
0x1381 ,
 
0x7351 ,
0x7d9d ,
0x2c63 ,
0xc5e1 ,
0x2c64 ;
 
i0=^x_1;
l0=%x_1;
 
i4=^c_f1;
l4=%c_f1;
 
m0=0; m1=1; m4=1;
 
do filt_1 until forever;
    
    mr=0,my0=pm(i4,m4);                 {a1}                                       
    mr1=io(0);                                          {input with io(0)}        
    mx0=dm(i0,m1);                      {Zi-1}                                       
    mr=mr-mx0*my0(su), mx1=dm(i0,m0),my0=pm(i4,m4); { mr=, Zi-2, a2}            
    mr=mr+mx1*my0(ss), my0=pm(i4,m4);           {mr = , b0}                 
    mr=mr(rnd);                                 { rounding}                
    if mv sat mr;                               {dont overflow}                  
    dm(i0,m0)=mr1,mr=mr1*my0(ss);                   {Zi-2= mr1, Zi-2 * b0}              
    my0=pm(i4,m4);                  {b1}
    mr=mr+mx0*my0(su),my0=pm(i4,m4);        {mr+Zi-1*b1, b2}
    mr=mr-mx1*my0(ss);mx0=dm(i0,m1);        {mr-Zi-2*b2}
    mr=mr(rnd);                 {rounding}
 
if mv sat mr;
    
    mr=0,my0=pm(i4,m4);                 {a1}                                                                                     
    mx0=dm(i0,m1);                      {Zi-3}                                       
    mr=mr-mx0*my0(su), mx1=dm(i0,m0),my0=pm(i4,m4); { mr=, Zi-4, a2}            
    mr=mr+mx1*my0(ss), my0=pm(i4,m4);           {mr = , b0}                 
    mr=mr(rnd);                                 { rounding}                
    if mv sat mr;                               {dont overflow}                  
    dm(i0,m0)=mr1,mr=mr1*my0(ss);                   {Zi-4= mr1, Zi-4 * b0}              
    my0=pm(i4,m4);                  {b1}
    mr=mr-mx0*my0(su),my0=pm(i4,m4);        {mr-Zi-3*b1, b2}
    mr=mr+mx1*my0(ss);mx0=dm(i0,m1);        {mr+Zi-4*b2}
    mr=mr(rnd);                 {rounding}
if mv sat mr;
 
    mr=0,my0=pm(i4,m4);                 {a1}                                                                                       
    mx0=dm(i0,m1);                      {Zi-5}                                       
    mr=mr-mx0*my0(su), mx1=dm(i0,m0),my0=pm(i4,m4); { mr=, Zi-6, a2}            
    mr=mr+mx1*my0(ss), my0=pm(i4,m4);           {mr = , b0}                 
    mr=mr(rnd);                                 { rounding}                
    if mv sat mr;                               {dont overflow}                  
    dm(i0,m0)=mr1,mr=mr1*my0(ss);                   {Zi-6= mr1, Zi-6 * b0}              
    my0=pm(i4,m4);                  {b1}
    mr=mr-mx0*my0(su),my0=pm(i4,m4);        {mr-Zi-5*b1, b2}
    mr=mr+mx1*my0(ss);              {mr+Zi-6*b2}
    mr=mr(rnd);                 {rounding}
    if mv sat mr;               
filt_1: io(1)=mr1;                                         {output on io(1)}
.endmod;

@Mikl___ · 08.12.2013, 17:12

Slave,
а это программа для какого микропроцессора?

@Slave · 08.12.2013, 17:26 **[ТС]**

Сообщение от Mikl___

Slave,
а это программа для какого микропроцессора?

ADSP 2181

Новые блоги и статьи Все статьи Все блоги /
Переходник USB-CAN-GPIO Eddy_Em 20.03.2026 Достаточно давно на работе возникла необходимость в переходнике CAN-USB с гальваноразвязкой, оный и был разработан. Однако, все меня терзала совесть, что аж 48-ногий МК используется так тупо: просто. . .	Оттенки серого Argus19 18.03.2026 Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .	SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .
[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .	делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .

@Mikl___ Ушел с форума 16373 / 7685 / 1080 Регистрация: 11.11.2010 Сообщений: 13,759
	08.12.2013, 17:12
	Slave, а это программа для какого микропроцессора? 0