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

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

Восстановить пароль Регистрация
Другие темы раздела
C++ Шифр Вернама http://www.cyberforum.ru/cpp-beginners/thread895094.html
Доброе время суток, совсем недавно заинтересовался С++. Скачал исходник - никак не могу понять в чем ошибка, отказывается компилироваться. Прошу вашей помощи с этим вопросом. Возможно что то делаю не так :)
C++ Машина пост Помогите, ни как ни могу вкурить, написать прогу на удаление наимаеньшего массива из двух в машине пост, каретка находится с лева первого массива. http://www.cyberforum.ru/cpp-beginners/thread895092.html
C++ Потоки thread и stream в чем разница
потоки thread и stream в чем разница. правильно? Первое - поток исполнения, он же thread, второе - поток данных - он же stream. правильно? еще интересует вопрос "пул потока" подскажите литературу по thread и ПУЛам.
C++ При работе с вектором изменяются адреса элементов
Здравствуйте. В соответствии с правилами создаю новую тему, хотя вопрос имеет отношение к другой проблеме в этом топике. Ситуация следующая: имеется указатель на нулевой элемент вектора. Именно указатель, а не итератор. При изменении вектора, по всей видимости, меняется адрес этого нулевого элемента. Для ненулевых элементов такое не происходит. В моем случае (см. скриншоты) есть указатель...
C++ Юникод символ в консоли http://www.cyberforum.ru/cpp-beginners/thread895043.html
Здорова! Есть вопрос, от можно ли в консоли как то вывести юникод символы? Допустим есть код символа в юникод U+2660 и как мне его вывести в консоли через программу? Пробовал отак cout <<(char)(+2660)<<endl; выводит не тот символ, а если U добавляешь, то ошибка.
C++ Перевод с радиан в градусы Дано значение α в радианах α > 0 < 2 ∏. Определите значение этого же угла в градусах. В качестве ∏ использовать 3,14. подробнее

Показать сообщение отдельно
FireProoF
17 / 17 / 1
Регистрация: 28.10.2012
Сообщений: 190
08.06.2013, 07:43     QR разложение (есть свои наброски)
Приветствую! Пытаюсь реализовать данный алгоритм по этой схеме.
Вот мой код:
Кликните здесь для просмотра всего текста

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;
}


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