Форум программистов, компьютерный форум CyberForum.ru

Быстрая свёртка м-последовательности - C++

Восстановить пароль Регистрация
 
АлександрКом
11 / 11 / 1
Регистрация: 21.10.2012
Сообщений: 327
19.11.2013, 13:24     Быстрая свёртка м-последовательности #1
Привет =) я пытаюсь сделать быструю свёртку m-последовательности.
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
#include <xdc/std.h>
#include <iostream>
#include <stdio.h>
#include <math.h>
#include "IntMath/Complex.h"
#include "IntMath/Fft.h"
#include "IntMath/Copy.h"
#include "IntMath/Peak.h"
#include "IntMath/Util.h"
#include "Gnss/Const.h"
#include "Mail.h"
#include "LocateBoard.h"
#include "Search.h"
 
 
static Complex32 result[2048];
static Complex32 resultmulti[2048];
static Complex32 obrat[2048];
static int modul[2048];
 
namespace {
 
const UInt prio = 8;
} // namespace
 
 
//------------------------------------------------------
 
Search::Search(Bios::Mailbox& dataBox):
        Handler(dataBox, prio),
        Fourier(2048)
{
}
 
 
Void Search::Process(Bios::Mail* mail) {
 
    CorrMail* dataMail = static_cast<CorrMail*>(mail);
    Complex32* address = dataMail->address;
    count = dataMail->count;
 
    for (int w = 0; w < 2048; w++) {
        resultmulti[w] = 0;
        result[w] = 0;
        obrat[w] = 0;
        modul[w] = 0;
    };
 
    Fourier.DoIt(address, result);
 
 
    for (int s = 0; s < 2048; s++) {
        resultmulti[s] = IntMath::Abs(result[s]);
    };
 
    Fourier.Undo(resultmulti, obrat);
 
    for (int k = 0; k < 2048; k++) {
        modul[k] = IntMath::Abs(obrat[k]);
    };
 
}
Значит, что я делаю: указатель address указывает на 1-ый элемент моей м-последовательности. Её я составил следующим образом - изначально взял 9 регистров и записал в них 1. Затем в цикле записывал значение регистра №7 (это число идёт в мою м-последовательность), складывал по модулю 2 числа из регистров 5 и 9, сдвигал вправо значения всех регистров, а в 1-ый записывал результат складывания 5 и 9 регистров. И так 511 раз, т.е. получил последовательность из 511 чисел. Её я продублировал 4 раза и в конце добавил четыре нуля - получилось 2048 элементов. В этой строчке
C++
1
Fourier.DoIt(address, result);
делаю БПФ, потом
C++
1
2
3
for (int s = 0; s < 2048; s++) {
        resultmulti[s] = IntMath::Abs(result[s]);
    };
умножаю результат на комплексно-сопряжённое число,
C++
1
Fourier.Undo(resultmulti, obrat);
- делаю обратное ПФ, и беру модуль результата -
C++
1
2
3
for (int k = 0; k < 2048; k++) {
        modul[k] = IntMath::Abs(obrat[k]);
    };
.
Когда строю график modul(он на картинке) - там вместо 3-х пиков получаются симметричные между собой и 2 почти одинаковых рядом. В чём загвоздка?
Миниатюры
Быстрая свёртка м-последовательности  
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
19.11.2013, 13:24     Быстрая свёртка м-последовательности
Посмотрите здесь:

Быстрая сортировка C++
Быстрая сортировка C++
Быстрая сортировка C++
Быстрая сортировка C++
C++ Быстрая сортировка
Быстрая сортировка C++
C++ Быстрая сортировка
C++ Быстрая сортировка

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
АлександрКом
11 / 11 / 1
Регистрация: 21.10.2012
Сообщений: 327
19.11.2013, 14:55  [ТС]     Быстрая свёртка м-последовательности #2
Возможно, нужно переместить тему в раздел "Цифровая обработка сигналов"?
Yandex
Объявления
19.11.2013, 14:55     Быстрая свёртка м-последовательности
Ответ Создать тему
Опции темы

Текущее время: 22:30. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru