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

функция exp - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 11, средняя оценка - 4.73
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 10:44     функция exp #1
Чем можно заменить функцию exp???у меня степень типа double и вылетает ошибка типо "floatinf point overflow".
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
naymkazp
20 / 20 / 1
Регистрация: 14.11.2010
Сообщений: 52
29.04.2013, 11:00     функция exp #2
Продемонстрируйте код.
Как вариант можно заменить на
C++
1
2
3
4
#include <math.h>
 
double x,y;
y=pow(M_E,x);
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 11:03  [ТС]     функция exp #3
Цитата Сообщение от naymkazp Посмотреть сообщение
Продемонстрируйте код.
Как вариант можно заменить на
C++
1
2
3
4
#include <math.h>
 
double x,y;
y=pow(M_E,x);
да что показывать число типо double является степенью. число е возводим в эту степень. с функцией pow такая же ошибка
naymkazp
20 / 20 / 1
Регистрация: 14.11.2010
Сообщений: 52
29.04.2013, 11:16     функция exp #4
Вот полный код для проверки работоспособности функции с разными числами:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
#include "stdafx.h"
#include <math.h>
#include <iostream>
 
int _tmain(int argc, _TCHAR* argv[])
{
double x,y;
 
for(;;){
    std::cin>>x;
    y=exp(x);
    std::cout<<y<<std::endl;
}
return 0;
}
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 11:36  [ТС]     функция exp #5
Цитата Сообщение от naymkazp Посмотреть сообщение
Вот полный код для проверки работоспособности функции с разными числами:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
#include "stdafx.h"
#include <math.h>
#include <iostream>
 
int _tmain(int argc, _TCHAR* argv[])
{
double x,y;
 
for(;;){
    std::cin>>x;
    y=exp(x);
    std::cout<<y<<std::endl;
}
return 0;
}
такое и у меня работает, а когда степень типа 0.000000011 вылетает
naymkazp
20 / 20 / 1
Регистрация: 14.11.2010
Сообщений: 52
29.04.2013, 11:38     функция exp #6
Цитата Сообщение от Мишутка1 Посмотреть сообщение
такое и у меня работает, а когда степень типа 0.000000011 вылетает
0.000000011 это уже не тип double. попробуйте long double
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 11:40  [ТС]     функция exp #7
Цитата Сообщение от naymkazp Посмотреть сообщение
0.000000011 это уже не тип double. попробуйте long double
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
long double sum = 0;
long double ag2 = 0;
for (int i_n = 0; i_n < 512; i_n++)
{
sum = 0.0;
for (int j_k = 0; j_k < 512; j_k++)
{
 
    if ((i_n == 0)||(j_k == 0)) {
    ag2 = (2 * pi * j_k * i_n)/512.0;
    sum += (Ak[j_k] * 1);
    }
    else
    {
       ag2 = (2 * pi * j_k * i_n)/512.0;
       sum += Ak[j_k] * exp(ag2);//pow(M_E,ag2)
    }
}
//Ak1[i_n] = sum/512;
}
naymkazp
20 / 20 / 1
Регистрация: 14.11.2010
Сообщений: 52
29.04.2013, 11:55     функция exp #8
массив Ak[j_k] чемто заполняеться заранее?
ибо уже на 90 проходе цикла по i_n получаеться сумма в минус бесконечность.
но ошибки нет

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
// 1.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include <math.h>
#include <iostream>
 
int _tmain(int argc, _TCHAR* argv[])
{
long double sum = 0.0;
long double ag2 = 0.;
double pi=3.14;
double Ak[512];
 
 for (int i_n = 0; i_n < 512; i_n++)
 {
     sum = 0.0;
    for (int j_k = 0; j_k < 512; j_k++)
    {
 
    if ((i_n == 0)||(j_k == 0)) {
        ag2 = (2 * pi * j_k * i_n)/512.0;
        sum += (Ak[j_k] * 1.0);
    }
     else
    {
        ag2 = (2.0 * pi * j_k * i_n)/512.0;
        sum += Ak[j_k] * exp(ag2);//pow(M_E,ag2)
        
    }
    }
 //Ak1[i_n] = sum/512;
    std::cout<<i_n<<"  "<<sum<<std::endl;
 }
 
system("pause");
return 0;
}
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 12:00  [ТС]     функция exp #9
Цитата Сообщение от naymkazp Посмотреть сообщение
массив Ak[j_k] чемто заполняеться заранее?
ибо уже на 90 проходе цикла по i_n получаеться сумма в минус бесконечность.
но ошибки нет

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
// 1.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include <math.h>
#include <iostream>
 
int _tmain(int argc, _TCHAR* argv[])
{
long double sum = 0.0;
long double ag2 = 0.;
double pi=3.14;
double Ak[512];
 
 for (int i_n = 0; i_n < 512; i_n++)
 {
     sum = 0.0;
    for (int j_k = 0; j_k < 512; j_k++)
    {
 
    if ((i_n == 0)||(j_k == 0)) {
        ag2 = (2 * pi * j_k * i_n)/512.0;
        sum += (Ak[j_k] * 1.0);
    }
     else
    {
        ag2 = (2.0 * pi * j_k * i_n)/512.0;
        sum += Ak[j_k] * exp(ag2);//pow(M_E,ag2)
        
    }
    }
 //Ak1[i_n] = sum/512;
    std::cout<<i_n<<"  "<<sum<<std::endl;
 }
 
system("pause");
return 0;
}
да он конечно заполнен числами double
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 12:04  [ТС]     функция exp #10
Цитата Сообщение от Мишутка1 Посмотреть сообщение
да он конечно заполнен числами double
вот формула, которая этими циклами реализуется
Миниатюры
функция exp  
naymkazp
20 / 20 / 1
Регистрация: 14.11.2010
Сообщений: 52
29.04.2013, 13:04     функция exp #11
1. что такое i в формуле в степени экспоненты?
2.в самой сумме степень экспоненты может достигать 3215 а это стандартными средствами не посчитать.
3. лучше полное задание выдайте.
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 13:16  [ТС]     функция exp #12
Цитата Сообщение от naymkazp Посмотреть сообщение
в самой сумме степень экспоненты может достигать 3215 а это стандартными средствами не посчитать.
это как понимать?

Добавлено через 4 минуты
Цитата Сообщение от naymkazp Посмотреть сообщение
лучше полное задание выдайте.
это формула обратного ДПФ, там вот нужна экспонента
naymkazp
20 / 20 / 1
Регистрация: 14.11.2010
Сообщений: 52
29.04.2013, 13:17     функция exp #13
Цитата Сообщение от Мишутка1 Посмотреть сообщение
это как понимать?
2*pi*k*n/N=2*3.14*511*512/512=3209.08
это степень экспоненты для 511 элемента суммы для 512 элемента x_n.
а это уже длинная арифметика и высчитывается по своим законам.
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 13:19  [ТС]     функция exp #14
Цитата Сообщение от naymkazp Посмотреть сообщение
2*pi*k*n/N=2*3.14*511*512/512=3209.08
это степень экспоненты для 511 элемента суммы для 512 элемента x_n.
а это уже длинная арифметика и высчитывается по своим законам.
дак вы же сказали, что до 3215 нормально считает?? выше сообщение
naymkazp
20 / 20 / 1
Регистрация: 14.11.2010
Сообщений: 52
29.04.2013, 13:24     функция exp #15
Цитата Сообщение от Мишутка1 Посмотреть сообщение
дак вы же сказали, что до 3215 нормально считает?? выше сообщение
я не сказал что нормально считает)
512 элементов это обязательно условие?
для 241 элементов еще посчитает,выше нет-здесь нужны другие подходы.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
29.04.2013, 13:27     функция exp
Еще ссылки по теме:

C++ вычисление математической функции exp(x)
exp(x) C++
C++ Ряд Тейлора exp(x)

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

Или воспользуйтесь поиском по форуму:
Мишутка1
0 / 0 / 0
Регистрация: 24.04.2013
Сообщений: 47
29.04.2013, 13:27  [ТС]     функция exp #16
Цитата Сообщение от naymkazp Посмотреть сообщение
я не сказал что нормально считает)
512 элементов это обязательно условие?
для 241 элементов еще посчитает,выше нет-здесь нужны другие подходы.
извините я ошибся, что вы так сказали=)) да именно 512 ну по заданию, и меня тоже вот эта буква i в формуле смущает, на wiki там не указано, что она значит. http://ru.wikipedia.org/wiki/%D0%94%...80%D1%8C%D0%B5
Yandex
Объявления
29.04.2013, 13:27     функция exp
Ответ Создать тему
Опции темы

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