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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 30, средняя оценка - 4.63
greendzol
0 / 0 / 0
Регистрация: 21.11.2010
Сообщений: 10
#1

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

23.12.2011, 23:13. Просмотров 4415. Ответов 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
#include <tchar.h>
#include <iostream>
#include <cmath>
#include <iomanip>
 
using namespace std;
int _tmain(int argc, _TCHAR* argv[])
{
    double *a=NULL, *b=NULL, **sum=NULL;
    const int N=5;
 
    double  x[N]={  -2, -1, 013},
        y[N]=   {   -1, -4, -5, -4, 4};
    int K,i,j,k,m;
    double z,c;
    cout<<"porjadok 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();
 
 
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
23.12.2011, 23:13
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод наименьших квадратов (C++):

Метод наименьших квадратов - C++
Можете около сложных строчек написать что там происходит(комментарии) #include &quot;stdafx.h&quot; #include &lt;tchar.h&gt; #include &lt;iostream&gt; ...

Метод наименьших квадратов - C++
Ребят,вот в чем беда. нужно для курсового решить задачу. нужно апроксимировать данные линейной зависимостью y=ax+b. В ручную все уже...

МНК (Метод наименьших квадратов) - C++
Помогите с реализацией аппроксимации МНК. Есть экспериментальные данные, хочу аппроксимировать их МНК. Не хочу использовать мат пакеты,...

Метод наименьших квадратов! (с Pascal на C++) - C++
Писал прогу на Паскале, так как не владею должными знаниями по С++ приходится просить помочь вас!!! Помогите перевести прогу написанную на...

Метод наименьших квадратов/лаба 3 - аппроксимация/ С++ - C++
Добрый день! Помогите исправить неисправность программы: #include &quot;stdafx.h&quot; #include &lt;iostream&gt; #include &lt;stdio.h&gt; ...

Метод наименьших квадратов, решение системы Гауссом - C++
подскажите, пожалуйста, где ошибка. матрицу Грамма вычисляю правильно, а вот когда после вычисления прямого хода гаусса, пытаюсь вывести...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Sou
0 / 0 / 0
Регистрация: 18.12.2010
Сообщений: 5
08.01.2012, 23:12 #2
Присоединяюсь к просьбе.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
08.01.2012, 23:12
Привет! Вот еще темы с ответами:

Линейная аппроксимация: метод наименьших квадратов (МНК) - C++
#include &quot;stdafx.h&quot; #include &lt;stdio.h&gt; #include &lt;stdlib.h&gt; #include &lt;math.h&gt; // Задание начального набора значений double **...

Линейный Метод Наименьших Квадратов, ошибка в расчетах - C++
Разработка программы аппроксимации табличной функции степенным рядом( линейный метод наименьших квадратов). использование глобальных...

Портировать программу с графиком из Pascal в С++ (метод наименьших квадратов) - C++
имеется код программы , которая реализует метод наименьших квадратов и строит графики. Может ли кто-то переписать код под С++ или ...

Метод Рауз-Болла построения магических квадратов N-ого порядка - C++
Здравствуйте! Занялся проблемой построения всех возможных магических квадратов N-ого порядка. Пока работаю над N = 3 и N = 4. Для N = 3 -...


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

Или воспользуйтесь поиском по форуму:
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru