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

Система нелинейных уравнений - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Как вставить перевод строки в конец? http://www.cyberforum.ru/cpp-beginners/thread903041.html
имеется строка char str="123456" как в нее грамотно вставить символ перевода строки? пробовал так: char str="123456"; str='\n'; str='\0'; не получается...
C++ Вводы-выводы из файла Приветствую.Изучаю программирование,в общем мне нужно сдать курсовую.В курсовой задачи.Я написал сами коды программ,но к ним по заданию нужны тесты(такие тесты,чтобы заваливали программу).Но проблема заключается в том,что я не могу никак понять,как сделать так чтобы значения брались из файла(например инпут_0.тхт)и выводились в (оутпут_0.тхт). Вот код одной из самых простых моих задач,прошу... http://www.cyberforum.ru/cpp-beginners/thread903039.html
C++ стеки. Переделать под функцию, то что в мейне
Помогите переделать в функцию, то что в мейне. Пытался уже, но ошибки выдаёт. ( //стек.cpp #include "stdafx.h" #include "iostream" #include "stack.h" using namespace std; #define max 20;
C++ Выполнить перестановку в заданном массиве
Выполнить перестановку в данном массиве Е22 так чтобы сначала были записаны подряд 9 первых элементов массива имеют четные индексы й После них разместить элементы массива имеющие нечетные индексы с номерами 1 = 1,3 ..... 27
C++ Проблемой при создании массива через переменную http://www.cyberforum.ru/cpp-beginners/thread902982.html
На C++ столкнулся с проблемой при создании массива через переменную. Таким образом: int list; i,j - переменные типа int на CodeBlocks все прошло без проблем (компилятор MinGW). на Microsoft Visual C++ получил на это 3 ошибки: 1> error C2057: требуется константное выражение 2> error C2466: невозможно выделить память для массива постоянного нулевого размера 3> error C2133: k: неизвестный...
C++ По номеру дня недели определить количество занятий Написать программу, которая по номеру дня недели (натуральному числу от 1 до 7) выдает в качестве результата количество занятий в Вашей группе в этот день(без разницы сколько занятий, разумное от 3-5 занятий сойдет) Выполнить задание с использованием пречисления. подробнее

Показать сообщение отдельно
vlad_light
4 / 4 / 0
Регистрация: 24.09.2012
Сообщений: 178
18.06.2013, 01:56  [ТС]     Система нелинейных уравнений
Проверьте код, пожалуйста. Спасибо!
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
#include <vector>
#include <functional>
#include <iostream>
 
class SystemOfNonLinearEquations 
{
public:
    SystemOfNonLinearEquations ();
    SystemOfNonLinearEquations (const SystemOfNonLinearEquations& sonle);
    SystemOfNonLinearEquations 
        (std::vector< std::function <double (const std::vector<double>)> >& functions);
 
    ~SystemOfNonLinearEquations ();
 
    void solve (const std::vector<double>& initial_approximation, 
                const std::vector<double>& left_boundary, 
                const std::vector<double>& right_boundary,
                const double error);
 
    void addFunction 
        (const std::function< std::function <double (std::vector<double>)> >& function);
    std::vector<double> getSolution () const;
 
private:
    std::vector< std::function <double (std::vector<double>)> > functions_;
    std::vector<double> solution_;
 
    double resudual_ (const std::vector<double> argument) const;
    double distance_ (const std::vector<double>& x, const std::vector<double>& y);
};
Добавлено через 2 часа 19 минут
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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
#include "SonLE.h"
 
SystemOfNonLinearEquations::SystemOfNonLinearEquations 
    (std::vector< std::function <double (const std::vector<double>)> >& functions)
{
    functions_ = functions;
}
 
void SystemOfNonLinearEquations::solve 
    (const std::vector<double>& initial_approximation, 
     const std::vector<double>& left_boundary, 
     const std::vector<double>& right_boundary,
     const double error)
{
    const double golden_const1 = (3 - std::sqrt (5.0)) / 2;
    const double golden_const2 = 1 - golden_const1;
    const double inner_error = error / 4;
    const int nVariables = solution_.size ();
 
    std::vector<double> a, b, x1, x2, previous_solution;
    double fx1, fx2;
        
    do
    {
        previous_solution = solution_;
 
        for (int i = 0 ; i < nVariables; ++i)
        {
            a = b = x1 = x2 = solution_;
 
            a.at (i) = left_boundary.at (i);
            b.at (i) = right_boundary.at (i);
 
            x1.at (i) = a.at (i) + (b.at (i) - a.at (i)) * golden_const1;
            x2.at (i) = a.at (i) + (b.at (i) - a.at (i)) * golden_const2;
 
            fx1 = resudual_ (x1);
            fx2 = resudual_ (x2);
 
            while (std::abs (b.at (i) - a.at (i)) > inner_error)
            {
                if (fx1 < fx2)
                {
                    b.at (i) = x2.at (i);
                    x2.at (i) = x1.at (i);
                    x1.at (i) = a.at (i) + (b.at (i) - a.at (i)) * golden_const1;
                    fx2 = fx1;
                    fx1 = resudual_ (x1);
                }
                else
                {
                    a.at (i) = x1.at (i);
                    x1.at (i) = x2.at (i);
                    x2.at (i) = a.at (i) + (b.at (i) - a.at (i)) * golden_const2;
                    fx1 = fx2;
                    fx2 = resudual_ (x2);
                }
            }
 
            solution_.at (i) = (a.at (i) + b.at (i)) / 2;
        }
    }
 
    while (distance_ (previous_solution, solution_) > error);
}
 
std::vector<double> SystemOfNonLinearEquations::getSolution () const
{
    if (resudual_ (solution_) > 1)
        std::cerr << "getSolution: No solution found";
 
    return solution_;
}
 
double SystemOfNonLinearEquations::resudual_ (const std::vector<double> argument) const
{
    if (functions_.empty ())
    {
        std::cerr << "residual_: No functions initialized";
        return 0;
    }
 
    double result = 0;
    const int nFunctions = functions_.size ();
 
    for (int i = 0; i < nFunctions; ++i)
        result += std::abs (functions_.at (i) (argument));
 
    return result;
}
 
double SystemOfNonLinearEquations::distance_ 
    (const std::vector<double>& x, 
     const std::vector<double>& y)
{
    if (x.size () != y.size ())
    {
        std::cerr << "distance_: Wrong input size";
        return 0;
    }
 
    double distance = 0;
    double temp;
 
    for (auto iterX = x.cbegin (), iterY = y.cbegin (); iterX != x.cend (); ++iterX, ++iterY)
    {
        temp = *iterX - *iterY;
        distance += temp * temp;
    }
 
    return std::sqrt (distance);
}
Добавлено через 18 секунд
 
Текущее время: 09:51. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru