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

Как запрограммировать экспоненту - C++

Восстановить пароль Регистрация
 
 
Рейтинг: Рейтинг темы: голосов - 58, средняя оценка - 4.78
Rooney
0 / 0 / 0
Регистрация: 22.11.2009
Сообщений: 37
01.08.2010, 17:50     Как запрограммировать экспоненту #1
Добрый вечер!
Помогите пожалуйста с вопросом: как запрограммировать функцию exp(x). Это стандартная функция, мне нужно сделать аналог своей функции.
exp(x) = 1 + x + (x^2)/2!+(x^3)/3!....+(x^n)/n!+.....
Нужно использовать несколько циклов?
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Rooney
0 / 0 / 0
Регистрация: 22.11.2009
Сообщений: 37
02.08.2010, 18:28  [ТС]     Как запрограммировать экспоненту #21
kazak,
alex_x_x
Спасибо большое!!! Всё разобрался
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Rooney
0 / 0 / 0
Регистрация: 22.11.2009
Сообщений: 37
03.08.2010, 22:05  [ТС]     Как запрограммировать экспоненту #22
Подредактировал, вот окончательный вариант. Компилятор Dev C++:
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
/*
Proc 40.
Описать функцию Exp1(x,eps) вещественного типа (параметры x, eps - вещественные,
eps > 0), находящую приближённое значение функции exp(x):
    exp(x) = 1 + x + (x^2)/(2!) + (x^3)/(3!) + ... + (x^n)/(n!) + ...
(n! = 1*2*3*...*n). В сумме учитывать все слагаемые, боольшие eps.
С помощью Exp1 найти приближённое значение экспоненты для данного
x при шести данных eps.
*/
#include<windows.h>
#include<iostream>
#include<math.h>
double exp1( const double x, const double precision );
using namespace std;
int main(){
      SetConsoleCP(1251);
      SetConsoleOutputCP(1251);  
      int nStep = 0;
      double x, eps;
      int i = 0;
      cout << "Введите X" << endl;
      cin >> x;
      cout << "Введите  eps"<< endl;
      cin >> eps;
      cout << "Значение функции exp1(x) = " << exp1(x, eps) << endl;
      cout << "Значение стандартной функции: "<< exp(x) << endl;
      system("pause");
      return 0;
}
 
double exp1( const double x, const double precision ){
       double dVal, dTemp;
       int nStep = 1;
       for( dVal = 1, dTemp = 1; dTemp >= precision ; nStep++ ){
            dTemp *= x/nStep;
            dVal += dTemp;
       }
       return dVal;
}
alex_x_x
бжни
 Аватар для alex_x_x
2441 / 1646 / 84
Регистрация: 14.05.2009
Сообщений: 7,163
09.08.2010, 17:45     Как запрограммировать экспоненту #23
беременным женщинам и детям не показывать..

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
#include <iostream>
#include <iterator>
#include <stdexcept>
#include <numeric>
 
class ExpIterator;
 
ExpIterator begin( const double );
ExpIterator end( const double );
 
class ExpIterator : public std::iterator<std::input_iterator_tag, double> {
public:
  friend ExpIterator begin( const double );
  friend ExpIterator end( const double );
  ExpIterator( const ExpIterator& it ){
    x = it.x;
    dPrecision = it.dPrecision;
    dTemp = it.dTemp;
    bEnd = it.bEnd;
    uStep = it.uStep;
  }
  ExpIterator& operator++() {
    if( bEnd ){
      throw std::out_of_range( "bad idea ++ end iterator" );
    }else{
      dTemp *= x/uStep;
      ++uStep;
    }
    return *this;
  }
  ExpIterator operator++(int){ ExpIterator tmp(*this); operator++(); return tmp; }
  bool operator==(const ExpIterator& rhs) {
    if( !rhs.bEnd && !bEnd ){
      return( uStep == rhs.uStep );
    }else{
      return( bEnd ? rhs.dTemp < dPrecision : dTemp < rhs.dPrecision );
    }
  }
  bool operator!=(const ExpIterator& rhs) {return !operator==(rhs);}
  double operator*() {return dTemp;}
private:
  ExpIterator( const double _x, const double _dPrecision = 0, const bool _bEnd = false ) : 
    x(_x), dPrecision(_dPrecision), uStep(1), bEnd(_bEnd)  {
    if( !bEnd ){ dTemp = 1.0; }
  }
  double x, dPrecision, dTemp;
  unsigned uStep;
  bool bEnd;
};
 
ExpIterator begin( const double x ){ return ExpIterator( x ); }
ExpIterator end( const double dPrecision ){ return ExpIterator( 0, dPrecision, true ); }
 
int main(){
  ExpIterator it = begin( 5.0 ), end_it = end( 1e-5 );
  std::cout << std::accumulate( it, end_it, .0 );
  return 0;
}
nikkka
Мат в 32 хода
 Аватар для nikkka
235 / 170 / 8
Регистрация: 10.09.2009
Сообщений: 1,096
09.08.2010, 18:04     Как запрограммировать экспоненту #24
Цитата Сообщение от fasked Посмотреть сообщение
просто какая-то прогрессия
А разве это не разложение в ряд тейлора?
alex_x_x
бжни
 Аватар для alex_x_x
2441 / 1646 / 84
Регистрация: 14.05.2009
Сообщений: 7,163
09.08.2010, 18:17     Как запрограммировать экспоненту #25
Цитата Сообщение от nikkka Посмотреть сообщение
А разве это не разложение в ряд тейлора?
относительно точки 0, или проще говоря в ряд маклорена
nikkka
Мат в 32 хода
 Аватар для nikkka
235 / 170 / 8
Регистрация: 10.09.2009
Сообщений: 1,096
09.08.2010, 18:54     Как запрограммировать экспоненту #26
alex_x_x,
Цитата Сообщение от Rooney Посмотреть сообщение
exp(x) = 1 + x + (x^2)/2!+(x^3)/3!....+(x^n)/n!+.....
разве это по отношению к 0??
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
09.08.2010, 19:05     Как запрограммировать экспоненту
Еще ссылки по теме:

C++ Запрограммировать на языке С
Найти матричную экспоненту, используя сумму ряда до слагаемого 10-4 C++
Функция вычисляющая экспоненту разложением в ряд Тейлора C++

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

Или воспользуйтесь поиском по форуму:
alex_x_x
бжни
 Аватар для alex_x_x
2441 / 1646 / 84
Регистрация: 14.05.2009
Сообщений: 7,163
09.08.2010, 19:05     Как запрограммировать экспоненту #27
http://www.cyberforum.ru/cgi-bin/latex.cgi?f(x) = \sum_{i=0}^{\infty} f^{(i)}(x-x_0) \frac{(x-x_0)^i}{i!}

при x_0 = 0

http://www.cyberforum.ru/cgi-bin/latex.cgi?f(x) = \sum_{i=0}^{\infty} f^{(i)}(x) \frac{x^i}{i!}
Yandex
Объявления
09.08.2010, 19:05     Как запрограммировать экспоненту
Ответ Создать тему
Опции темы

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