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

Дискретное преобразование Фурье - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Класс триугольников http://www.cyberforum.ru/cpp-beginners/thread859085.html
Уважемые юзеры) Прошу помощи при написании проги: Нужно создать класс триугольников. Атрибуты задать на свое усмотрение. Зарание благодарен.
C++ Нужно перепрыгнуть строку Проблема такая, есть скрипт он должен открывать клиент игры, а за место этого он открывает браузер. #ifdef _KOR_WEBLAUNCHER_ #ifdef _GMTOOL_ if( lpCmdLine ) { FILE* fpLog = NULL; fpLog = fopen( "./ClientStartLog.txt", "a+" ); if( fpLog ) { http://www.cyberforum.ru/cpp-beginners/thread858998.html
Найти среднее арифметическое элементов для каждой строки заданного двухмерного массива C++
Найти среднее арифметическое элементов для каждой строки заданного двухмерного массива. Значения элементов массива ввести с клавиатуры
C++ Написать программу формировании и вывода массива размером 3х6
Напишите программу формировании и вывода массива размером 3х6, каждые три последовательных элемента в строке которого составлены по следующему правилу: 1-й - равен квадрату суммы текущих индексов; 2-й - равен случайному числу; 3-й - равен полусумме двух предыдущих элементов
C++ выбор ответа http://www.cyberforum.ru/cpp-beginners/thread858988.html
В двухмерном массиве размером n x m хранятся результаты опроса n человек на m вопросов. Ответами на вопросы являются слова «Да», «Нет», «не знаю». Написать программу, указывающую для каждого вопроса проценты ответов «Да», «Нет», «не знаю».
C++ формирование массива по правилу программа формирования и вывода массива размером 6 на 3, каждые три последовательных элемента в строке которого составлены по следующему правилу: 1-й - равен удвоенному числу, введенному с клавиатуры; 2-й - равен квадрату суммы индексов предыдущего элемента; 3-й - равен остатку от деления индекса текущего элемента на 2. подробнее

Показать сообщение отдельно
Женя Курило
2 / 2 / 0
Регистрация: 09.12.2010
Сообщений: 10

Дискретное преобразование Фурье - C++

08.05.2013, 12:01. Просмотров 928. Ответов 0
Метки (Все метки)

Добрый день! Общем задание (2 картинки). Я знаю, не сложно.
Вот код, который я написал.
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
#include <iostream>
#include <cmath>
#include <complex>
#include <iomanip>
 
using namespace std;
const double PI = 3.14159265;
void out_func_value (double f_arr[][80],const int Ngurnal,const double T,const double delta);
int _tmain(int argc, _TCHAR* argv[])
{
    
    const int Ngurnal=16;
    const double T=Ngurnal*1e-6;
    const double delta=T/10;
    //const double PI = 3.14159265;
    // f(x)=Ngurnal*sin(((2*PI)/T)*k*delta)
    double func[2][80];
    out_func_value(func,Ngurnal,T,delta);
    
  
    
    typedef complex <double> C;
    
    C vec[80];
    C i(0,-1);
    C temp;
    for (int k=0;k<80;k++){
        temp=0;
        for (int n=0;n<80;n++){
            temp+=(func[1][n])*(cos(k*n*(2*PI)/80)-sin(k*n*(2*PI/80))*i);
        }
        vec[k]=temp;
    }
    
    
    cout<<endl;
    for (int j=0;j<80;j++)
        cout<<vec[j]<<endl;
    cout<<endl;
    for (int j=0;j<80;j++)
        cout<<abs(vec[j])<<"\t"<<arg(vec[j])<<endl;
    
    cin.get();
    cin.get();
    return 0;
}
void out_func_value (double f_arr[][80],const int Ngurnal,const double T,const double delta)// функция для сохранения значений функции в массив и печати этих значений на экран
{
    cout.setf(ios::left); // устанавливаем выравнивание по левому краю
    for (int i=0;i<80;i++){
        cout<<setw(15)<<(f_arr[0][i]=(((2*PI)/T)*i*delta)); // манипулятор setw  используется для задания ширины вывода
        cout<<setw(20)<<(f_arr[1][i]=(Ngurnal*cos(((2*PI)/T)*i*delta)))<<endl;
    }
        
}
Но когда я проверял свои расчеты в Matlab, я обнаружил, что они неверны. Я предполагаю идет переполнение типа. С long double тоже самое. Как можно решить эту проблему? Заранее спасибо!
Миниатюры
Дискретное преобразование Фурье   Дискретное преобразование Фурье  
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 00:19. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru