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

метод QR разложекния - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ В массив В перенести все элементы массива А, стоящие правее максимального элемента, и имеющие нечетный индекс http://www.cyberforum.ru/cpp-beginners/thread744632.html
Ввести массив А. В массив В перенести все элементы массива А, стоящие правее максимального элемента, и имеющие нечетный индекс. Массив В отсортировать по возрастанию, используя метод шейкер-сортировки.
C++ Нужно совместить два кода Вот первый код, из него нужно организовать вывод результатов подсчета в док файл, либо .txt #include "iostream" #include <windows.h> using namespace std; http://www.cyberforum.ru/cpp-beginners/thread744630.html
C++ Подсчитать количество "счастливых" шестизначных автобусных билетов(сумма первых трех цифр равна сумме трех последних цифр)
Подсчитать количество "счастливых" шестизначных автобусных билетов, т.е таких, в номерах которых сумма трех первых цифр равна сумме трех последних. (Воспользоваться тем, что число "счастливых" билетов равно S02+S12+...+S272, где Sn - кол-во чисел от 0 до 999, сумма цифр которых равна n) вот такое условие,не совсем понимаю как правильно делать
Перегрузка операций C++
#include <iostream> using namespace std; class chislo { public: chislo(); chislo(int a,int b); chislo(const chislo&); chislo(char N,int pls);
C++ Сортировка пузырьком http://www.cyberforum.ru/cpp-beginners/thread744613.html
Нужно отсортировать список, но почему-то выходит за пределы списка void operator++() {int vrem,i,j; tec2=((struct Elsp*)malloc(sizeof(struct Elsp))); tec2->sled=NULL; tec=start; //tec2=tec->sled; for (int i = 0; i <n-1; i++) {tec2=tec->sled;
C++ Помогите с написанием блок-схемы к программе подскажите пожалуйста,как она будет выглядеть к этой программе # include <conio.h> # include <stdio.h> void main () { int i; const n = 7; подробнее

Показать сообщение отдельно
zenid
0 / 0 / 0
Регистрация: 26.12.2011
Сообщений: 149
24.12.2012, 20:11     метод QR разложекния
есть ошибка: нужно поставить begin, не понятно где?

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
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
#include <stdio.h>
#include <conio.h>
#include <math.h>
 
#define _N 10   // число уравнений
#define _M 10   // число неизвестных
 
class matrix {
  public:
   matrix(int n, int m);
   ~matrix();
   void print_slau();
   void print_x();
   void qr();
   void seta(int i, int j, float value);
   void setb(int i, float value);
  private:
   float **a;
   float *b;
   float *x;
   int n,m;
};
// ------------------------------------------------------
// Конструктор
// ------------------------------------------------------
matrix::matrix(int N,int M) {
 n = N; m = M;
 a = new float* [n];
 for (int i=0; i<n; i++) {
      a[i] = new float [m];
      for (int j=0; j<m; j++) a[i][j]=0; }
 b = new float [n];
 x = new float [n];
 for (i=0; i<n; i++) { b[i]=0; x[i]=0; }
}
 
// ------------------------------------------------------
// Деструктор
// ------------------------------------------------------
matrix::~matrix() {
 for (int i=0; i<n; i++) delete a[i];
 delete[] a;
 delete b; delete x;
}
 
// ------------------------------------------------------
// Вывод СЛАУ на экран
// ------------------------------------------------------
void matrix::print_slau() {
 for (int i=0; i<n; i++) {
    for (int j=0; j<m; j++) printf("%.2f  ",a[i][j]);
    printf("| %.3f  ",b[i]);
    printf("\n");
 }
}
 
// ------------------------------------------------------
// Печать решения и невязки решения
// ------------------------------------------------------
void matrix::print_x() {
 float max=0,h;
 for (int i=0; i<n; i++) {
      h=0;
      for (int j=0; j<n; j++) h=h+x[j]*a[i][j];
      if (max<fabs(b[i]-h)) max=fabs(b[i]-h);
      printf("x[%i]=%.3f  ",i,x[i]);
      if (i==5) printf("\n");
 }
 printf("\nМакс. невязка: %f\n",max);
}
// ------------------------------------------------------
// Установить значение матрицы a[][]
// ------------------------------------------------------
void matrix::seta(int i, int j, float value) {
  a[i][j] = value;
}
// ------------------------------------------------------
// Установить значение вектора b[]
// ------------------------------------------------------
void matrix::setb(int i, float value) {
  b[i] = value;
}
 
// ------------------------------------------------------
// QR-алгоритм
// ------------------------------------------------------
void matrix::qr () {
  int l,k;
  float c[_N][_M];
  float s[_N][_M];
  float akk,akl,alk,all,bk,bl;
 
 
 
 
  // Прямой ход
  for (k=0; k<n-1; k++) {   // "Большой" шаг (исключение переменных)
   for (l=k+1; l<n; l++) {     // "Малый" шаг
    c[k][l] = a[k][k] / (sqrt( a[k][k]*a[k][k] + a[l][k]*a[l][k] ));
    s[k][l] = a[l][k] / (sqrt( a[k][k]*a[k][k] + a[l][k]*a[l][k] ));
 
       // Умножение матрицы a[][] на T[k][l]
       akk=a[k][k];  alk=a[l][k];  akl=a[k][l];  all=a[l][l];
       a[k][k] =  akk*c[k][l] + alk*s[k][l];
       a[k][l] =  akl*c[k][l] + all*s[k][l];
       a[l][k] = -akk*s[k][l] + alk*c[k][l];
       a[l][l] = -akl*s[k][l] + all*c[k][l];
 
       // Вектор свободных членов умножается на T[k][l]
       bk = b[k]; bl = b[l];
       b[k] =  bk*c[k][l] + bl*s[k][l];
       b[l] = -bk*s[k][l] + bl*c[k][l];
   }
  }
  // Теперь матрица a[][] -- верхняя диагональная.
 
  // Обратный ход
  float h;
  x[n-1]=b[n-1]/a[n-1][n-1];
  for (l=(n-1);l>=1;l--) {
      h=b[l-1];
      for (k=(l+1);k<=n;k++) h=h-x[k-1]*a[l-1][k-1];
      x[l-1]=h/a[l-1][l-1];
  }
 
}
 
 
// ------------------------------------------------------
// Main
// ------------------------------------------------------
void main() {
 clrscr();
 matrix A(_N,_M);
 for (int i=0; i<_N; i++) for (int j=0; j<_M; j++) A.seta(i,j,i+j);
 for (i=0; i<_N; i++) A.setb(i,1);
 printf("Решение СЛАУ методом вращения (QR-разложение)\n");
 printf("=Исходная система========================================\n");
 A.print_slau();
 printf("\n=Решение системы=========================================\n");
 A.qr();
 A.print_x();
 getch();
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 12:42. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru