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

Разложение матрицы А - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Наименьшее количество квадратов в прямоугольнике http://www.cyberforum.ru/cpp-beginners/thread694046.html
Дан прямоугольник размерностью NxM. Нужно заполнить его наименьшим количеством квадратов. Я погуглил на данную тему, но нигде не нашёл внятного решения. Собственно проблема не в самой реализации, а в понятии алгоритма. В чем он должен заключаться? Я вот подумал, например перебор суммы квадратов чисел, которые дают в сумме площадь прямоугольника и при этом само число не превышает длину наименьшей...
C++ Вычисления выражений Здравствуйте, пожалуйста, помогите составить программы. В первом скрине, вычисление выражения. Во втором, вычисление с применением условного оператора. http://www.cyberforum.ru/cpp-beginners/thread694045.html
C++ Двумерные массивы с использованием функций
вот сама задача: Составить функцию преобразования матрицы A в новую матрицу путем вычеркивания строки с номером L. Составить также функцию печати исходной матрицы, а рядом - полученной. Использовать эти функции для матрицы B и строки с номером 2. вот что придумал: #include "stdio.h" #include "conio.h" void preob(int*b,int *c,int m,int n,int k); void main() { int b; int i,j,m,n,l;
C++ Как просчитать шаги от домов до магазина
Готовлюсь к олимпиаде (уже будет 2 декабря) и на тренировочных курсах попалась такая задача : Не могу понять, как просчитать шаги от домов до магазина. Может надо вычитать индексы массива? Кто знает, пожалуйста, приведите код. Задача C. Дома и магазины На Новом проспекте построили подряд 10 зданий. Каждое здание может быть либо жилым домом, либо магазином, либо офисным зданием. Но...
C++ Сложная задача про полный квадрат http://www.cyberforum.ru/cpp-beginners/thread694037.html
Есть 9 разных чисел a,b,c,d,e,f,g,h,i. Сумма любых двух разных чисел из данного ряда является полным квадратом. При этом a+b+c=d+e+f=g+h+i и сумма (a+b+c) должна быть минимальная. Ничего не могу придумать. Помогите пожалуйста. Добавлено через 1 час 51 минуту (извините не правильно условие понял) ИЗМЕНЕННЫЙ ВАРИАНТ: найти произведение целых чисел: a,b,c,d,e,f,g,h,i; таких что...
C++ Как разобраться с ошибками в коде Помогите пожалуйста разобраться с ошибками в коде вот код: #include <iostream.h> #include <string.h> int main () { struct strc{ подробнее

Показать сообщение отдельно
ejk
0 / 0 / 0
Регистрация: 14.05.2012
Сообщений: 80
10.11.2012, 20:15     Разложение матрицы А
Всем доброй почти ночи!
Работаю над алгоритмом разложения матрицы А на матрицы Q и R
Пользуюсь данным алгоритмом

И вот такой у меня вышел код:

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
#include <stdio.h>
#include <math.h>
#define N 3
void main () {
    int st,i,j,k,l;
    float A[N][N];
    float Q[N][N];
    float R[N][N];
    float T[N][N];
    float X[N][N];
    float B[N];
    float d=0,h=0;
 
 
///Задание матриц///
 
 
    FILE *fp;
    fp = fopen("m.txt","r");
        for (int i=0;i<N;i++){
            for(int j=0;j<N;j++){
                fscanf(fp, "%f", &A[i][j]);
                fscanf(fp, "%f", &R[i][j]);
                }
        }
    printf("\n");
 
    FILE *f;
    f = fopen("s.txt","r");
    for(i=0;i<N;i++){
        fscanf(f, "%f", &B[i]);
        
    }
    printf("\n");
    
    for(i=0;i<N;i++){
        for(j=0;j<N;j++){
            if(i==j){
            T[i][j]=1;
            Q[i][j]=1;
            }
            else {
            T[i][j]=0;
            Q[i][j]=0;
            }
        }
    }
    
    
    //Разложение
 
 
    for (i=0;i<N;i++){
        for(j=i+1;j<N;j++){
            if(R[i][j]!=0){
                T[i][i]=R[i][i]/sqrt(R[i][i]*R[i][i]+A[j][i]*A[j][i]);
                T[j][j]=T[i][i];
                T[i][j]=R[j][i]/sqrt(R[i][i]*R[i][i]+A[j][i]*A[j][i]);
                T[j][i]=-Q[i][j];
            }
        
        d=0;
        
            for(l=0;l<N;l++){
                d+=T[i][l]*Q[l][j];
                Q[i][j]=d;
                        }
    
            for(k=0;k<N;k++){
                d+=T[i][k]*R[k][j];
                R[i][j]=d;
                        }
    
        }
    }
 
    
 
    for(int k=0;k<N;k++){
        for(int u=0;u<N;u++){
            printf("%f ",Q[k][u]);
                }
            printf("\n");
    }
    
    printf("\n");
        
 
}
Что-то идет не так, ибо матрица Q выводится неправильной.
Можете подсказать, где искать ошибку?

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