Сейсмодатчик своими руками - Страница 6

BolshoyK · 21.07.2013, 07:05

niXto, Попробовал план...

последние x,y,z в серии из 30 замеров, среднее значение x,y,z, условия - датчик прижат пальцем к столу на время замера. Старался не трястись...

датчик настроен на 8g:

10:54:02.968› -2;-5;-15;-2.67;-5.40;-14.83
10:54:03.031› -1;-5;-14;-2.40;-5.60;-15.03
10:54:03.031› -4;-7;-16;-3.10;-5.67;-15.13
10:54:03.093› -4;-7;-16;-2.47;-5.53;-14.53
10:54:03.156› -2;-4;-13;-2.80;-5.33;-14.40
10:54:03.156› -3;-7;-15;-2.57;-5.40;-14.47
10:54:03.218› -3;-7;-15;-2.50;-5.93;-14.40
10:54:03.281› -2;-5;-14;-2.70;-5.27;-14.97
10:54:03.343› -2;-5;-14;-2.50;-5.43;-14.50
10:54:03.343› -3;-6;-15;-2.73;-5.63;-14.90
10:54:03.406› -2;-5;-14;-2.87;-5.63;-14.80
10:54:03.468› -3;-5;-14;-2.97;-5.13;-14.50
10:54:03.468› -2;-5;-15;-2.50;-5.40;-14.80
10:54:03.531› -4;-7;-16;-2.50;-5.30;-14.67
10:54:03.593› -3;-6;-14;-2.73;-5.30;-14.27
10:54:03.656› -3;-5;-14;-2.77;-5.40;-14.73
10:54:03.718› -1;-4;-16;-2.47;-5.50;-14.57
10:54:03.718› -3;-7;-16;-2.53;-5.50;-14.90
10:54:03.781› -2;-5;-14;-2.63;-5.80;-14.70
10:54:03.843› -1;-4;-13;-2.57;-5.17;-14.60
10:54:03.843› -4;-7;-16;-2.70;-5.17;-14.83
10:54:03.906› -2;-4;-13;-2.87;-5.57;-14.67
10:54:03.968› -1;-4;-13;-2.07;-4.97;-14.27
10:54:04.031› -3;-7;-16;-2.40;-5.47;-14.40

Датчик настроен на 2g:

11:03:29.968› -9;-21;-56;-8.37;-20.27;-57.27
11:03:30.031› -10;-21;-60;-9.00;-20.60;-56.13
11:03:30.093› -11;-21;-58;-9.23;-20.67;-57.70
11:03:30.156› -10;-22;-58;-9.23;-22.03;-56.73
11:03:30.156› -10;-20;-57;-8.77;-20.90;-57.00
11:03:30.218› -7;-20;-57;-8.07;-20.40;-57.10
11:03:30.281› -9;-21;-57;-8.97;-20.47;-56.27
11:03:30.281› -9;-20;-58;-9.20;-20.23;-57.27
11:03:30.343› -7;-19;-55;-8.73;-20.30;-56.43
11:03:30.406› -8;-22;-57;-9.43;-21.47;-57.43
11:03:30.468› -11;-20;-58;-8.73;-21.13;-57.30
11:03:30.468› -8;-20;-56;-8.83;-20.63;-56.80
11:03:30.531› -9;-19;-56;-9.27;-21.03;-56.97
11:03:30.593› -12;-22;-57;-9.10;-20.73;-57.33
11:03:30.593› -9;-20;-57;-8.37;-21.00;-56.47
11:03:30.656› -6;-20;-56;-9.17;-20.60;-56.13
11:03:30.718› -10;-21;-58;-9.43;-20.83;-56.43
11:03:30.781› -7;-19;-56;-8.53;-20.13;-57.37
11:03:30.843› -8;-18;-56;-8.57;-19.43;-56.93
11:03:30.843› -8;-21;-58;-8.73;-20.57;-57.17
11:03:30.906› -11;-21;-59;-8.87;-20.23;-58.43
11:03:30.968› -7;-20;-56;-7.77;-20.07;-56.53
11:03:30.968› -10;-21;-57;-9.27;-21.10;-57.23
11:03:31.031› -8;-22;-56;-8.13;-20.93;-55.80

Код:

void loop()
{
    x = y = z = 0;
    for (i=0; i‹30; i++){
      xVal = mySensor.readAxis('x');
      yVal = mySensor.readAxis('y');
      zVal = mySensor.readAxis('z');
      x = x + xVal;
      y = y + yVal;
      z = z + zVal;
    }
    x = x / 30;
    y = y / 30;
    z = z / 30;

    Serial.print(xVal,DEC);
    Serial.print(";");
    Serial.print(yVal,DEC);
    Serial.print(";");
    Serial.print(zVal,DEC);
    Serial.print(";");

    Serial.print(x, 2);
    Serial.print(";");
    Serial.print(y, 2);
    Serial.print(";");
    Serial.println(z, 2);
}

Как выглядят графики 8g режима:

Как выглядят графики 2g режима:

niXto · 21.07.2013, 08:27

Ну как бы работает, только это простое усреднение, а не свертка... Заведи кольцевой буфер на 30 значений, после окончания измерения записывай новое значение на место самого старого в буфере, высчитывай среднее значение в буфере и клади его в другой буфер на 30 значений. По окончании нового измерения снова клади его на место самого старого в кольцевом, считай среднее в кольцевом и клади в другой буфер. Через 30 измерений считай среднее во втором буфере, очищай его и выводи получившееся число на график. Кольцевой буфер никогда не очищай - в нем постоянно крутятся 30 последних сырых значений

Таким макаром при той же латентности сгладишь пики ещё в квадрат раз

BolshoyK · 21.07.2013, 11:50

Так ?

15:45:06.921› -8;-20;-55;-8.63;-20.43;-56.70
15:45:07.046› -8;-19;-58;-8.64;-20.55;-56.58
15:45:07.109› -9;-21;-56;-8.43;-20.61;-56.94
15:45:07.171› -10;-22;-58;-8.07;-20.67;-56.47
15:45:07.296› -9;-21;-56;-8.48;-20.16;-56.39
15:45:07.359› -12;-24;-58;-8.54;-20.79;-56.18
15:45:07.421› -12;-21;-57;-8.95;-20.85;-56.26
15:45:07.484› -10;-22;-56;-8.82;-20.17;-56.08
15:45:07.609› -6;-22;-55;-8.83;-20.82;-56.47
15:45:07.671› -9;-21;-58;-8.12;-20.73;-56.56
15:45:07.734› -10;-23;-57;-8.51;-20.66;-56.25
15:45:07.859› -10;-21;-58;-8.46;-20.62;-56.43
15:45:07.921› -11;-21;-57;-8.89;-20.17;-56.68
15:45:07.984› -9;-20;-57;-8.69;-19.52;-56.32
15:45:08.046› -8;-22;-57;-8.42;-20.33;-56.10
15:45:08.171› -7;-20;-56;-8.30;-20.80;-56.39
15:45:08.234› -9;-21;-57;-8.74;-20.58;-56.63
15:45:08.296› -8;-22;-58;-9.38;-20.80;-57.01
15:45:08.359› -9;-21;-57;-8.65;-20.76;-56.63
15:45:08.484› -6;-22;-54;-8.43;-20.62;-56.52
15:45:08.546› -8;-20;-58;-8.39;-20.62;-56.56
15:45:08.609› -7;-20;-57;-8.84;-20.88;-56.94
15:45:08.734› -9;-19;-56;-8.83;-20.60;-56.83
15:45:08.796› -8;-20;-58;-8.81;-19.94;-56.76
15:45:08.859› -8;-19;-55;-8.40;-20.82;-57.11
15:45:08.921› -10;-19;-57;-9.02;-20.82;-56.52
15:45:09.046› -10;-23;-55;-9.31;-20.42;-56.53
15:45:09.109› -8;-20;-57;-9.02;-20.45;-56.47
15:45:09.171› -8;-22;-55;-8.19;-20.40;-56.56
15:45:09.296› -7;-20;-57;-8.38;-20.92;-56.51
15:45:09.359› -6;-18;-55;-8.53;-20.35;-56.19
15:45:09.421› -9;-20;-55;-8.45;-20.19;-56.14
15:45:09.484› -7;-19;-56;-8.46;-20.65;-56.12
15:45:09.609› -9;-20;-57;-8.25;-20.42;-55.62
15:45:09.671› -9;-20;-57;-8.06;-20.53;-56.20
15:45:09.734› -9;-19;-56;-8.39;-20.28;-56.45
15:45:09.796› -7;-20;-55;-8.52;-20.33;-56.55
15:45:09.921› -11;-21;-59;-9.03;-21.08;-56.31
15:45:09.984› -9;-21;-54;-8.93;-20.74;-57.07
15:45:10.046› -7;-20;-58;-8.25;-20.77;-56.46

Код:

void loop()
{
    for (i=0; i‹30; i++){
      buf1[i][0] = mySensor.readAxis('x');
      buf1[i][1] = mySensor.readAxis('y');
      buf1[i][2] = mySensor.readAxis('z');
      x = y = z = 0;
      for (j=0; j‹30; j++){
        x = x + buf1[j][0];
        y = y + buf1[j][1];
        z = z + buf1[j][2];
      }
      buf2[i][0] = x / 30;
      buf2[i][1] = y / 30;
      buf2[i][2] = z / 30;
      
    }
    x = y = z = 0;
    for (j=0; j‹30; j++){
      x = x + buf2[j][0];
      y = y + buf2[j][1];
      z = z + buf2[j][2];
    }
    x = x / 30;
    y = y / 30;
    z = z / 30;

    Serial.print(buf1[29][0],DEC);
    Serial.print(";");
    Serial.print(buf1[29][1],DEC);
    Serial.print(";");
    Serial.print(buf1[29][2],DEC);
    Serial.print(";");

    Serial.print(x, 2);
    Serial.print(";");
    Serial.print(y, 2);
    Serial.print(";");
    Serial.println(z, 2);

}

BolshoyK · 21.07.2013, 14:34

С экзелем не устроил результат оффлайновый. Написал прогу, рисует в реалтайме данные с акселя.
Первые 2 - возил по осям. прижав к столу. 3-й - поднимал и опускал. Ну... что сказать, руки у меня не могут соблюсти только вертикальное перемещение.

tempora · 21.07.2013, 14:51

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

Ну... что сказать, руки

С руками трудно сверяться, может попробовать "возить" движком с шатунным механизмом? И размеры перемещений будут известны, и другие значения, да и форма д.б. синусоидальной - даже визуально м.б. оценить соответствие, а уж для "алгеброй поверить" - просто раздолье.

akegor · 21.07.2013, 15:00

BolshoyK, это Вы с целью "для проэталонировать"?
Лучше использовать гравитацию, покрутив по осям. Результат будет гораздо точнее.
А вот "возить" после этого можно будет по линейке и проверять перемещения.

BolshoyK · 21.07.2013, 16:07

akegor, Да не, это просто оценить как отображаются мои действия.
Тут положил архив с прогой на C#, ардуино скетч и библиотеку для акселерометра. Вдруг кто захочет подсобить...http://serg22.sibgtu.ru/radio/GraphAccel01.rar
Прога настроена на COM6.

akegor · 21.07.2013, 16:23

BolshoyK, с Шарпом и Ард... не общаюсь (с Шарпом - пока, да и не очень хочется, пока не подопрет). Потому - практически не помогу. Разве только теорией, как сам это вижу.

BolshoyK · 21.07.2013, 16:28

akegor, ну у меня ардуина тоже валялась, в первый раз подключил и это первая прога на ней у меня. А шарп - это просто что под рукой было. Достаточно просто писать несложные вещи, одноразовые.
С теорией самые сложности. Дальше то что делать ?

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

akegor · 21.07.2013, 16:37

Теперь бы определиться, что в итоге нужно получить, что конкретно интересует?
Стоит ли так сильно заглаживать?
Считаю, что нужно перейти непосредственно к перемещениям в нормальных единицах, в которых можно оценить результат, и посмотреть результат от сырых данных и от сглаженных.
Мне так видится.
ПС. Если от сырых данных получится дрожание +/- мм - стоит ли выглаживать и терять быстрые смещения?