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

QR разложение (есть свои наброски) - C++

Восстановить пароль Регистрация
 
FireProoF
17 / 17 / 1
Регистрация: 28.10.2012
Сообщений: 190
08.06.2013, 07:43     QR разложение (есть свои наброски) #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
// QR.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
#include <iostream>
#include "math.h"
using namespace std;
 
void QR(double**, int, double**, double**);
 
int _tmain(int argc, _TCHAR* argv[])
{
    int row=3,col=3;
    double** A = new double*[row];
    for(int i=0;i<row;i++)  A[i]=new double[col];
 
    for(int i=0;i<row;i++)  for(int j=0;j<col;j++)  cin>>A[i][j];
 
    for(int i=0;i<row;i++)
    {
        cout<<endl;
        for(int j=0;j<col;j++)  cout<<A[i][j]<<" ";
    }
    cout<<endl;
 
    double** R = new double*[row];
    for(int i=0;i<row;i++)  R[i]=new double[col];
    R=A;
 
    double** Q = new double*[row];
    for(int i=0;i<row;i++)  Q[i]=new double[col];
    for(int i=0;i<row;i++)  for(int j=0;j<col;j++)  Q[i][j]=1.0;
 
    QR(A,row,Q,R);
 
    for(int i=0;i<row;i++)
    {
        cout<<endl;
        for(int j=0;j<col;j++)  cout<<R[i][j]<<" ";
    }
    cout<<endl;
 
    for(int i=0;i<row;i++)
    {
        cout<<endl;
        for(int j=0;j<col;j++)  cout<<Q[i][j]<<" ";
    }
    cout<<endl;
 
    cout<<endl;
    system("pause");
    return 0;
}
 
void QR(double** A,int N, double** Q, double** R)
{
    double** E = new double*[N];//единичная матрица
    for(int i=0;i<N;i++)    E[i]=new double[N];
 
    for(int i=0;i<N;i++)    
    {
        for(int j=0;j<N;j++)
        {
            if(i==j)    E[i][j]=1;
            E[i][j]=0;
        }
    }
 
    for(int i=0;i<N;i++)    
    {
        cout<<endl;
        for(int j=0;j<N;j++)
        {
            cout<<E[i][j]<<"  ";
        }
    }
 
    double** T = new double*[N];
    for(int i=0;i<N;i++)    T[i]=new double[N];
    
    for(int k=0;k<(N-1);k++)
    {
        for (int i=(k+1);i<N;i++)
        {
            T=E;
            if(R[i][k]!=0)
            {
                T[k][k]=R[k][k]/sqrt(R[k][k] * R[k][k] + A[i][k] * A[i][k]);
                T[i][i]=T[k][k];
                T[k][i]=R[i][k]/sqrt(R[k][k] * R[k][k] + A[i][k] * A[i][k]);
                T[i][k]= -T[k][i];
            }
            
            for (int i=0; i<N; i++)//R=T*R
            {
                for (int j=0; j<N; j++)
                {
                    R[i][j]=0;
                    for (int r=0; r<N; r++)
                    {
                        R[i][j]+=(T[i][r]*R[r][j]);
                    }
                }
            }
 
            for (int i=0; i<N; i++)//Q=T*Q
            {
                for (int j=0; j<N; j++)
                {
                    Q[i][j]=0;
                    for (int r=0; r<N; r++)
                    {
                        Q[i][j]+=(T[i][r]*Q[r][j]);
                    }
                }
            }
 
        }
    }
    return;
}


Укажите на ошибки или покажите готовое решение)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
08.06.2013, 07:43     QR разложение (есть свои наброски)
Посмотрите здесь:

C++ Контрольная работа(( проверь свои знания):D
C++ Разложение в ряд Тейлора(код есть, но ошибочный)
C++ Массив меняет свои значения
Свои пункты в меню панели задач C++
C++ Замена зарезервированных слов и символов в с++ на свои
Класс, знающий все свои экземпляры C++
C++ Си - создать свои функции new, delete
C++ Вставить между одинаковыми символами вставить " * " (звёздочку) (есть наброски)

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

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