caxapa.ru :: А как истинно православно квантовать синус?!

IBAH (10.11.2015 19:34, просмотров: 10296)

А как истинно православно квантовать синус?! сгенерил таблицу signed short sin[1024] к целым переходил округлением на постоянную составляющую проверил - ОК! первую гармонику проверил е-13! ортогональность синус-косинус проверил - ОК! коэф. гармоник по 10 гармоникам 120Дб ну вот какая жопа: s=0; for(k=256;k<(256+512);k++) s+=sin[k]; if(s==0) print("OK!"); else print("ZOPA!"); результат ZOPA! то есть интеграл за пол периода косинуса не равен нулю!!! нутром чую, что надо как то не так квантовать, а как понять не могу...

Ответить

- А у Вас какие значения в sin[0] и sin[1023]? - Tiro(12.11.2015 00:53)
- A*sin(2pi*(n+0.5)/N), где N - кол-во отсчётов на период, должно делиться на 4. n = 0..N-1 - argus98(11.11.2015 09:47)
  - Чтобы не вводить дополнительное условие можно записать так A*sin(2pi*(n+0.5)/(4*N)) - Codavr(11.11.2015 10:51)
  - Отблин не успел! думал я один такой умный... всем спасибо IBAH(11.11.2015 09:53)
- А я что-то полюбил кусочками прямых аппроксимировать, если есть время на одно лишнее умножение... всё-таки точнее, чем ступеньки. - Николай Коровин(10.11.2015 23:06)
- Потому что интеграл очень уж хреновый. Вычмата не было? мимокрокодил(147 знак., 10.11.2015 22:19, )
  - что-то мне этот вычмат сильно подгонку под ответ напоминает... - IBAH(11.11.2015 08:37)
    - Простейший метод трапеций - мимокрокодил(11.11.2015 10:04, )
- и сколько дельта от нуля? - RED_DRAGON(10.11.2015 21:47)
  - Порой от "OK" до "ZOPA" всего один шаг) - Олдфаг(10.11.2015 21:59, )
    - Именно один шаг! если брать сумму от 256 до <=(256+512), все ОК - IBAH(10.11.2015 22:47)
      - А как же иначе? Должна быть симметрия отрицательных и положительных четвертинок косинуса. А при условии < самая первая выборка цикла не была компенсирована. - Олдфаг(10.11.2015 23:22, )
        
        Именно это счастье отсутствует!!! поэтому я и спрашиваю - как квантовать, чтобы стать счасливым, как выяснилось округлять не катит IBAH(11.11.2015 09:43)
- А чему именно равен? Чтоб нам легче гипотезы генерить. - Крок(10.11.2015 21:41)
  - таблица int16, от -32767 до +32767, в остатке s=32767 - IBAH(10.11.2015 22:50)
    - Один лишний отчёт. .... надо подумать. - Крок(10.11.2015 23:33)
    - А вы чего ждали? Возьмите массив покороче, 4 или 8 точек на период, и посмотрите глазами, какие значения вы суммируете. Заодно увидите, почему добавление одного отсчёта даёт сумму 0. - как дети, чесслово(10.11.2015 23:21, )
      - Или сдвинуть все отсчеты на полкванта. - Крок(11.11.2015 08:57)
        
        каюсь, не въехал, почему то показалось, что это не решение, а методика поиска решения - IBAH(11.11.2015 10:31)
        
        чего там сдвигать..., я думал задача обсосана, кто нибудь скажет "аппроксимировать надо МНК" или что нибудь подобное IBAH(11.11.2015 09:20)
- Пардон, это просто из любви к искусству или что-то сделать надо? Но вообще кагбэ напрашивается квантование по четверти периода и дальнейшее размножение, не? - SciFi(10.11.2015 20:15 - 20:20)
  - делать ничего не надо! надо получить целочисленный синус максимально похожий на настоящий - IBAH(10.11.2015 22:54)
    - К.О. сообщает: у Вас противоречие в соседних предложениях. "Ничего не надо" - "надо". - SciFi(10.11.2015 22:55)
      - вот такая я противоречивая личность :), это я к тому что задача чисто математическая - IBAH(11.11.2015 08:41)
        
        Да какая там математика... Больше похоже, что на единицу ошибся в каком-то индексе. - SciFi(11.11.2015 08:50)
        
        см. выше - IBAH(11.11.2015 09:25, ссылка)
- Какой у вас тип у s? Может быть, просто сумма не влезает в short? - Ксения(10.11.2015 19:54)
  - это я для понятливости написал, на самом деле я в столбик складываю :) - IBAH(10.11.2015 22:55)

Микроконтроллеры