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

Метод наименьших квадратов - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 19, средняя оценка - 4.68
corpdelta
0 / 0 / 0
Регистрация: 09.03.2014
Сообщений: 15
29.03.2014, 12:17     Метод наименьших квадратов #1
Можете около сложных строчек написать что там происходит(комментарии)
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
113
114
115
116
117
118
119
120
121
122
123
124
125
#include "stdafx.h"
#include <tchar.h>
#include <iostream>
#include <cmath>
#include <iomanip>
 
using namespace std;
int _tmain()
{
    setlocale(LC_ALL, "rus");
    double *a = NULL, *b = NULL, **sum = NULL;
    const int N = 8;
 
    double  x[N] = { -1, 0, 1, 3, 5, 7, 9, 11 },
            y[N] = { 0, 2, 4, 6, 8, 10, 12, 14 };
    int K, i, j, k, m;
    double z, c;
    cout << "Порядок k= ";
    cin >> K;
    cin.get();
    b = new double[K + 1];
    a = new double[K + 1];
    sum = new double *[K + 1];
    for (int i = 0; i<K + 1; i++)
        sum[i] = new double[K + 1];
 
 
 
    for (i = 0; i<K + 1; i++)
    {
        for (j = 0; j<K + 1; j++)
        {
            sum[i][j] = 0;
            for (k = 0; k<N; k++)
            {
                sum[i][j] += pow(x[k], i + j);
            }
        }
    }
    for (i = 0; i<K + 1; i++)
    {
        b[i] = 0;
        for (k = 0; k<N; k++)
        {
            b[i] += pow(x[k], i) * y[k];
        }
    }
    for (int i = 0; i<K + 1; i++)
    {
        for (int j = 0; j<K + 1; j++)
            cout << sum[i][j] << "\t";
        cout << endl;
    }
 
    for (int i = 0; i<K + 1; i++)
        cout << b[i] << "\t";
    cout << endl;
    for (int i = 0; i<K; i++)
    {
        m = i;
        for (int j = i + 1; j<K + 1; j++)
        {
            if (fabs(sum[m][i])<fabs(sum[j][i]))
                m = j;
        }
        for (int k = i; k<K + 1; k++)
        {
            z = sum[m][k];
            sum[m][k] = sum[i][k];
            sum[i][k] = z;
        }
        z = b[m];  b[m] = b[i];  b[i] = z;
        for (int i = 0; i<K + 1; i++)
        {
            for (int j = 0; j<K + 1; j++)
            {
                cout << sum[i][j] << setw(15);
            }
            cout << "b" << i << "= ";
            cout << b[i] << endl;
        }
        cin.get();
 
        for (int j = i + 1; j<K + 1; j++)
        {
            c = -sum[j][i] / sum[i][i];
            cout<<"\n\n ! "<<c<<endl;
            for (int k = i; k<K + 1; k++)
            {
                sum[j][k] = sum[j][k] + c*sum[i][k];
            }
            b[j] = b[j] + c*b[i];
        }
    }
 
    for (int i = 0; i<K + 1; i++)
    {
        for (int j = 0; j<K + 1; j++)
        {
            cout << sum[i][j] << setw(15);
        }
        cout << "b" << i << "=";
        cout << b[i] << endl;
    }
    cin.get();
 
    a[K] = b[K] / sum[K][K];
 
    for (int i = K + 1 - 2; i >= 0; i--)
    {
        for (int k = i + 1; k<K + 1; k++)
        {
            b[i] = b[i] - a[k] * sum[i][k];
        }
        a[i] = b[i] / sum[i][i];
    }
 
 
    for (int i = 0; i<K + 1; i++)
        cout << "X" << i << "=" << a[i] << endl;
    cin.get();
 
 
    system("Pause");
}
Лучшие ответы (1)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
29.03.2014, 12:17     Метод наименьших квадратов
Посмотрите здесь:

Произведение 2х наименьших из 3х заданых C++
метод наименьших квадратов, решение системы Гауссом C++
метод наименьших квадратов C++
C++ портировать программу с графиком из Pascal в С++(метод наименьших квадратов)
C++ C++ Линейный Метод Наименьших Квадратов. ошибка в расчетах
МНК (Метод наименьших квадратов) C++
C++ Метод наименьших квадратов! (с Pascal на C++)
Метод наименьших квадратов C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
kiborgdelto
70 / 72 / 27
Регистрация: 23.03.2011
Сообщений: 141
29.03.2014, 14:13     Метод наименьших квадратов #2
Сообщение было отмечено автором темы, экспертом или модератором как ответ
10 меняем кодировку символов для нормального вывода сообщений на русском языке
11 создаём указатели на массивы (2 на одномерные и 1 на двумерный)
20 ожидаем нажатия любой клавиши
21-22 создаём 2 динамических одномерных массива
23-25 создаём динамический двумерный массив
77 выводим значение элемента массива и ограничиваем поле вывода данных(т.е. не более 15 символов)

остальные строки это сам метод наименьших квадратов, его думаю знаете
Yandex
Объявления
29.03.2014, 14:13     Метод наименьших квадратов
Ответ Создать тему
Опции темы

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