Hardella IDE

[Владимир Ситников]

Мою формулу можно свести в итоге к такой: T = T0 + ΣdT. Приращение dT можно вычислить в программе CoDeSys и передать при обращении к программе PRU. В итоге все сведется к суммированию заранее известного числа раз. Или я не прав ?

Это будет другая формула, и она будет неким приближением исходной. dT нужно будет корректировать каждые 100мкс.

[Вольд]

Да я не верно написал. Должно быть так: 1/T(t) = 1/T0 + a * t;

1/T0 - это константа, она вычисляется однажды в основной программе и передается в PRU.

[Вольд]

T(t) = T0/(1+T0*a*t); или T(t) = T0/(1+T0*aΣdt);

Если делать вычисления в лоб, то на каждом шаге надо будет суммировать, умножать, и делить. Не фонтан.

Надо попробовать разложить функцию в степенной ряд, возможно полегчает.

Если ограничиться двумя членами, то в окрестности точки 133 мкс (7500 Гц) получается так:

T(t) = T0/(1+T0*a*t) = 200/(1+200*0,000000833*t) = 196 - 0,0325*(t-133)

Вроде не плохо получилось. Только что делать с коэффициентом 0,0325 ? Владимир, что скажите ?

[Владимир Ситников]

Вроде не плохо получилось. Владимир, что скажите ?

Я бы предложил сделать вычисления без PRU, проверить, что получаются нужные цифры и тогда уже переходить к кодированию.

И ещё: время лучше не в мкс измерять, а 5нс интервалах. PRU работает на частоте 200МГц, и 1 такт занимает 5нс. Счётчик тактов доступен из PRUграммы.

[Вольд]

И ещё: время лучше не в мкс измерять, а 5нс интервалах.

А это для чего нужно ?

[Владимир Ситников]

А это для чего нужно ?

Блок PRU_CURRENT_TIME: https://hardella.com/docs/pru/standa...u_current_time

[Вольд]

Блок PRU_CURRENT_TIME: https://hardella.com/docs/pru/standa...u_current_time

Зачем мерить время в 5 нс интервалах ? Что это дает ?

[Вольд]

Владимир, операция целочисленного деления сколько у вас времени занимает ?

[Вольд]

С коэффициентом 0,0325 в степенном ряду что делать ? Как с дробными числами работать ?

Может так избавиться от дробного числа: 196 - 0,0325*(t-133) = 196 - (t-133)/31

Тогда погрешность при целочисленном делении будет большая.

Может так сделать T(t)*31 = 196*31 - (t-133), т.е. в результате вычислений получить T(t)*31, а потом результат поделить на 31.

[IVM]

А оператор CASE в Hardella есть ?