Как интегрировать алгоритм в m файле в симулинк?, Использовать итерационный алгоритм заданый кодом в модели Опции

[des00]

Добрый день!

Есть модель выравнивания двух потоков данных, на основе их корреляции в симулинке. В ручном режиме, с помощью фильтров с переменной задержкой и slider gain я могу выровнять потоки в ручную. Но хотелось бы сделать автоматически. Сложность у меня вот в чем. Алгоритм выравнивания поделен на 2 этапа :
1. Поиск средней задержки, на основе поиска максимума корреляционной функции с помощью перебора.
2. Удержания выравнивания с помощью метода градиентного спуска.
Кодом такая конструкция описывается элементарно, два цикла, переход из одного в другой по условию. Но вот как перенести это в симулинк? Делать графический аналог (рисовать КА, делать case и т.д.) нет никакого желания, а использование встроенных матлабовских функций ограниченно тем, что от вызова к вызову переменные функции не сохраняются. Переносить всю модель в m файлы тоже не хотелось бы, все таки симулинк вещь более наглядная %)
Можно было бы использовать ступенчатую функцию и переключаться между двумя алгоритмами, но во первых время 1 ого этапа может варьироваться (от размера шага), а во вторых коряво это как то %)

Может есть более красивый способ?

Спасибо.

[des00]

С помощью persistent получилось сделать то, что хотел. Конечно не совсем красиво, но тем не менее. Но обнаружил занятный баг, который никак не могу понять почему происходит. В коде есть такое место смены состояния

Код

delta_delay = 0.5;
max_delay   = 9;
...
    if delay < max_delay
        if value > max_value
            max_value = value;
            max_value_delay = delay;
        end
        delay = delay + delta_delay;        
    else
        delay = max_value_delay;
        state = state + 1; % go to hold state
    end
..

так вот, когда delta_delay = 0.5 смена состояний происходит корректно при delay = 9

а вот когда delta_delay = 0.1 смена состояний происходит при delay = 9.1.

Почему, не понятно, моделька в аттаче %)

Прикрепленные файлы

 test_fsm_R14.zip ( 8.3 килобайт ) Кол-во скачиваний: 23

 

[_Anatoliy]

завтра гляну

[_Anatoliy]

Да,странная ошибка,пока не нашёл откуда ноги растут...
Если вместо
1).

Код

if (delay < max_delay)

поставить

Код

dx=floor(delay*10+1e-10);
    if (dx < floor(max_delay*10))

или вместо
2).

Код

delay = delay + delta_delay;

поставить

Код

delay = delay + floor(delta_delay*10);

то всё работает.
У целочисленной переменной delay откуда то возникает дробная часть величиной менее 1е-11.Ошибка представления числа вроде должна быть намного меньше.

[des00]

Да,странная ошибка,пока не нашёл откуда ноги растут...
У целочисленной переменной delay откуда то возникает дробная часть величиной менее 1е-11.Ошибка представления числа вроде должна быть намного меньше.

вот и я тоже не понял почему 10 раз по 0.1 не равно 1 це %) в принципе на результат функционирования схемы это не влияет, но ошибка занятная %)

[Gate]

вот и я тоже не понял почему 10 раз по 0.1 не равно 1 це %) в принципе на результат функционирования схемы это не влияет, но ошибка занятная %)

Скорее всего потому что при преобразовании в двоичный вид 0.1 нельзя точно представить как сумму ряда 2^(-n) - я не проверял. А 0.5 можно - я гарантирую это
Непонятно то, что точность представления double примерно 2^-53=10e-16, а у _Anatoliy получилось 10e-11.

[_Anatoliy]

Непонятно то, что точность представления double примерно 2^-53=10e-16, а у _Anatoliy получилось 10e-11.

Если быть точнее то получается 1е-12.
В процессе вычислений ошибка может накапливаться.

[Gate]

Если быть точнее то получается 1е-12.

Матлаб дает другую ошибку:

Код

>> format long eng
>> 0.1+0.1+0.1+0.1+0.1+0.1+0.1+0.1+0.1+0.1-1.0

ans =

   -111.022302462516e-018

>>