Статистическая обработка в любительских дозиметрах

/*
* (PCINT5/RESET/ADC0/dW) PB5 |1 OS 8| VCC
* (PCINT3/XTAL1/CLKI/OC1B/ADC3) PB3 |2 oled 7| PB2 (SCK/USCK/SCL/ADC1/INT0/PCINT2)
* (PCINT4/XTAL2/CLKO/OC1B/ADC2) PB4 |3 MOD GM 6| PB1 (MISO/DO/AIN1/OC0B/OC1A/PCINT1)
* GND |4 oled 5| PB0 (MOSI/DI/SDA/AINO/OC1A/AREF/PCINT0)
*
*/

#include <Tiny4kOLED.h>
#include <font16x32digits.h>
#include <avr/interrupt.h>

#define F_CPU 8000000UL
#define MOD 4 /* Управление ВВ ключем */
#define OS A0 /* Напряжение на ВВ части */
#define GM 1 /* Вход импульсов с GM */
#define MASSIVE_SIZE 10 /* Размер массива интервалов */
#define FAST_INTERVAL 2000
/* СБМ-19 -- 9000 */
#define SLOW_INTERVAL 9000
/* Доверительная вероятность
100 - 1 sigma 68.3%
200 - 2 sigma 95.4%
300 - 3 sigma 99.7%
*/
#define SIGMA 300
#define POWER 400 /* Напряжение ВВ */
#define CALIBR 0.4736 /* Вольт на разряд (0.4985) */

uint8_t width = 128;
uint8_t height = 32;

long /*cmax = 0, cmin = 1000000,*/ numm1, numm = 0;
long summ1, summ = 0;
int perc;
volatile int ms_old = 0, cnt, ind, interval;
volatile int inter[MASSIVE_SIZE];
int active, ms, refr, refr2, ii;

/* Обработчик прерываний от GM */
ISR(PCINT0_vect)
{
if ( ms_old != 0 ) /* Первое значение */
interval = millis() - ms_old;
else
ms_old = millis();
if (interval > 0) /* Прерывание по фронту игнорируем */
{
inter[ind++] = interval;
if ( ind > MASSIVE_SIZE - 1 ) /* Массив интервалов */
ind = 0;
ms_old = millis();
cnt++;
}
}


void setup() {
pinMode(MOD, OUTPUT);
analogReference(INTERNAL);
pinMode(OS, INPUT);
for ( ind = 0; ind < MASSIVE_SIZE; ind++)
inter[ind] = 0;
ind = 0;
pinMode(GM, INPUT_PULLUP);
oled.begin(width, height, sizeof(tiny4koled_init_128x32br), tiny4koled_init_128x32br);
//oled.setFont(FONT16X32DIGITS);
oled.setFont(FONT6X8);
oled.clear();
oled.on();
GIMSK = 0b00100000; /* Включаем прерывания для PCINT, PCIE */
PCMSK = 1 << GM; /* Разрешаем прерывания от PCINT1 */
sei();
}

/*
* Вывод на OLED
* rw -- Строка дисплея, считать от 0
* str -- строка для вывода
*/
void writeOled(int rw, String str)
{
oled.setCursor(0, rw);
oled.print(str);
oled.clearToEOL();
}

void loop() {
ms = millis();
if (ms - refr > SLOW_INTERVAL) // Значения измерения
{
refr = millis();
//if ( cnt > cmax)
// cmax = cnt;
//if (cnt < cmin)
// cmin = cnt;
summ = summ + cnt;
cnt = 0;
numm++;
if (summ > 0)
{
perc = round((float) (SIGMA /sqrt(summ)));
if ( perc > 100 )
perc = 100;
} else {
perc = 100;
}
writeOled(0, (String) round((float) summ / numm) + "uR " + (String) perc + "% " + (String) summ);

}

if (ms - refr2 > FAST_INTERVAL) // Значение поиска и высокого напряжения.
{
refr2 = millis();
summ1 = 0;
for ( ii = 0; ii < MASSIVE_SIZE; ii++)
{
summ1 = summ1 + inter[ii];
}
if (summ1 > 0)
{
writeOled(1, (String) round((float) MASSIVE_SIZE * SLOW_INTERVAL / summ1) + "uR " + (String) round((float) FAST_INTERVAL * MASSIVE_SIZE / summ1) + "cps");
} else {
writeOled(1, "0cps 0uR ");
}
writeOled(2, (String) round((float) analogRead(OS) * CALIBR) + "V " + (String) (active / (FAST_INTERVAL / 1000)) + " pps");
active = 0;
}

/* Упрравление ВВ частью */
if (analogRead(OS) < POWER / CALIBR)
{
active++;
PORTB |= (1 << MOD);
asm("nop");
asm("nop");
asm("nop");
asm("nop");
asm("nop");
PORTB &= ~(1 << MOD);
//digitalWrite(MOD, HIGH);
//digitalWrite(MOD, LOW);
}
}